题目链接:GCD
题目大意:给你一个N,要求你在这个N里面找出有多少个x使得gcd(x,N) >= m
题目思路:朴素的暴力肯定T,所以我们可以想到进行优化,gcd(x,N) >= m,则题目一定有大于等于m的最大公约数s,则s*a = x,s*b = N,而且a,b必须得互质,否则,最大公约数就不止s了,因为N > x,则b > a且a与b互质,我们枚举了符合条件的b之后要去找比b小且互质的数,也就是找b的欧拉值,那么这道题就可以做出来了
#include <map>
#include <set>
#include <cmath>
#include <vector>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
ll Euler(ll n){
ll res = n;
for(ll i = 2;i*i <= n;i++){
if(n%i == 0) res = res/i*(i-1);
while(n%i == 0) n /= i;
}
if(n > 1) res = res/n*(n-1);
return res;
}
int main(){
ll t,n,m,ans;
scanf("%lld",&t);
while(t--){
ans = 0;
scanf("%lld%lld",&n,&m);
for(int i = 1;i*i <= n;i++){//i为枚举的s
if(n%i == 0){
if(i >= m)ans += Euler(n/i);
if(i*i != n&&n/i >= m) ans += Euler(i);
}
}
printf("%lld\n",ans);
}
return 0;
}