CodeForces 75C Modified GCD （二分+數論）

題目連結：http://codeforces.com/problemset/problem/75/C

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))

#define pii pair<ll,ll>
#define mk(x,y) make_pair(x,y)

const int  maxn =1e5+5;
const int mod=998244353;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
題目大意：求在給定範圍内的兩個數的最大公約數。

先預處理出其gcd的所有因子，
然後排序，對每個查詢在容器裡二分找位置，取最大的那個，
判斷一下存在性即可。

時間複雜度：O（sqrt(n)+qloglogn）
*/
ll n,a,b,m;
int q;
vector<ll> p;
int main()
{
    scanf("%lld%lld",&n,&m);

    ll gd=gcd(n,m);
    for(int i=1;1LL*i*i<=gd;i++) if(gd%i==0){
        ll x=i,y=gd/i;
        p.push_back(x);
        if(x!=y) p.push_back(y);
    }

    sort(p.begin(),p.end());
    scanf("%d",&q);
    while(q--){
        scanf("%lld%lld",&a,&b);
        int pa=lower_bound(p.begin(),p.end(),a)-p.begin();
        int pb=upper_bound(p.begin(),p.end(),b)-p.begin();pb--;
        ///cout<<pa<<" "<<pb<<endl;
        if(pb>=pa) printf("%lld\n",p[pb]);
        else puts("-1");
    }
    return 0;
}