学习到的内容:
1.gcd(a,b)=1,则gcd(a+b,b)=1 即gcd(a,b)=gcd(a+b,b)
2.a*b能整除c, gcd(a,c)=1 则问题转化为b能整除c;
3. 由1.2.得, a*b能整除a+c ac互质, 则问题转化为b能整除a+c
4.求两个数互质的方法1.打表,2.枚举然后gcd=1
5.有a,b,gcd(a,b)的式子,可设d=gcd(a,b),a=k1*d,b=k2*d,再观察式子的特征
题解:
设d=gcd(a,b) a=k1*d,b=k2*d题目要求满足b*d能整除a+b所有a,b的值
b*gcd(a,b) / a+b 为整数 转化为k2*d*d / k1*d+k2*d 即k2*d / k1+k2为整数
k1,k2必互质,由学习到的内容3进一步转化问题为 d/(k1+k2)为整数
列举出所有k1,k2互质的值, 由于d整除(k1+k2) 则d>k1 d>k2,
因为a=k1*d<=n b=k1*d<=m 并且d>k1 d>k2则k1<sqrt(n) k2<sqrt(m);
对于每组k1k2互质,里面满足题目的对应的a,b的对数为min(n/k1/(k1+k2),m/k2/(k1+k2))
解释: n/k1表示1~n范围内有多少个a的系数为k1,则满足k1*d<=n的d的个数就是a的个数即n/k1,满足的d为1~n/k1 1~n/k1中满足d整除(k1+k2)的d的个数为n/k1/(k1+k2) 对于k2同理,则满足题目的a,b的对数(即对于这组k1,k2,满足题目条件d的个数)
判断k1,k2是否互质1.枚举每个k1,k2对, 求gcd(k1,k2)是否为1 2.预处理
预处理的方法
for (ll i = 2; i <= min(n/i,m/i); i++) {
for (ll k1 = i; k1 <= n/k1; k1 += i) {
for (ll k2 = i; k2 <= sqm; k2 += i) {
bad[k1][k2] = true;
}
}
}
bad数组中被打为true的k1,k2值不互质,没被打的互质
ac代码
#include<iostream>
using namespace std;
const int N=2e6+10;
int gcd(int a,int b){
return b?gcd(b,a%b):a;
}
int main(){
int T;
cin>>T;
while(T--){
int n,m;
cin>>n>>m;
long long ans=0;
for(int k1=1;k1<=n/k1;k1++){
for(int k2=1;k2<=m/k2;k2++){
if(gcd(k1,k2)==1)ans+=min(n/k1/(k1+k2),m/k2/(k1+k2));
}
}
cout<<ans<<endl;
}
}