对于x <= n n = p*d, x = q *d d = gcd(n,x) 则得出p q一定互质,否则p q中有公因子那么gcd(n,x) > d 。
所以只要对于每一个d >= m ,求出与p互质的数 的个数,也就是p的欧拉函数。
#include <iostream>
#include<stdio.h>
#include<cmath>
using namespace std;
int euler_phi(int n)
{
int m = (int)sqrt(n+0.5);
int ans = n;
for(int i = 2;i <= m;i++)
{
if(n % i == 0)
{
ans=ans/i*(i - 1);
while(n%i == 0) n/=i;
}
}
if(n > 1)ans = ans/n*(n - 1);
return ans;
}
int work(int n,int m)
{
int ans = 0;
for(int i = 1;i*i <= n;i++)
{
if(i*i == n && i >= m) ans+=euler_phi(i);
else if(n%i == 0 )
{
if(i >= m) {ans+=euler_phi(n/i);}
if(n/i >= m) {ans+=euler_phi(i);}
}
}
return ans;
}
int main()
{
int t,n,m;
scanf("%d",&t);
while(t--)
{
scanf("%d %d",&n,&m);
printf("%d\n",work(n,m));
}
return 0;
}