hdu 2588 GCD

对于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;
}