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")