4199. 公约数 - AcWing题库 

思路：

最大整数x一定是最大公约数的因数，所以先用__gcd(a,b)求出a和b的最大公因数，再用O(log(n))的算法求出最大公因数的因数，放到vector中，并将vector排序。利用STL中的upper_bound(res.begin(),res.end(),num)找到res中大于num的最小位置，如果该位置的元素大于等于l，则成立。

推荐博客：关于lower_bound( )和upper_bound( )的常见用法_lowerbound和upperbound-CSDN博客

for(auto i : v)遍历容器元素_for auto 遍历-CSDN博客

【C++】__gcd(x,y)函数_∑ y∈c gcd(x,y)-CSDN博客

C++代码：

#include<bits/stdc++.h>
using namespace std;int a,b,q,l,r;
vector<int> res;
int main(){cin>>a>>b>>q;int temp=__gcd(a,b);for(int i=1;i*i<=temp;i++){if(temp%i==0){res.push_back(i);res.push_back(temp/i);}}sort(res.begin(),res.end());/*for(auto i:res){cout<<i<<' ';}cout<<endl;*/for(int i=1;i<=q;i++){cin>>l>>r;int idx=upper_bound(res.begin(),res.end(),r)-res.begin();//找大于r的最小索引(顺序)if(idx-1>=0 && res[idx-1]>=l){cout<<res[idx-1]<<endl;}else{cout<<"-1"<<endl;}}return 0;
}

Python代码：

在python中：Python3二分查找库函数bisect(), bisect_left()和bisect_right()介绍-CSDN博客

import math
import bisect
res=[]
a,b=map(int,input().split())
q=int(input())
temp=math.gcd(a,b)
j=1
while j*j<=temp:if temp%j==0:res.append(j)res.append(int(temp/j))j+=1
res.sort()
for i in range(1,q+1):l,r=map(int,input().split())idx=bisect.bisect_right(res,r)if idx-1>=0 and res[idx-1]>=l:print(res[idx-1])else:print("-1")