hdu 4983 Goffi and GCD

参考博客：http://blog.csdn.net/yanghuaqings/article/details/47167875

首先，gcd(n-a,n)=gcd(a,n);
    假设  gcd(n-a,n)=x    即存在k1,k2,k1!=k2&&gcd(k1,k2)=1使得
    n-a=k1*x,n=k2*x;
    a=n-k1*x=k2*x-k1*x=(k2-k1)*x;
    gcd(n-a,n)=gcd(k1*x,k2*x)=gcd((k2-k1)*x,k2*x)=x;(因为k1!=k2&&gcd(k1,k2)=1);


    gcd(n-a,n)*gcd(n-b,n)=n^k;
    如果k>2,不存在
    如果k==2,只有一种情况a=n,b=n;

  如果n==1,只有一种情况
    其余情况为k==1的时候


    假设 gcd(n-a,n)==x ,那么(n-1)/x 与 n/x 互质,其中 x 与 n 是已知的,即 n/x 是已知的,转化为求小于求与n/x互质的数的个数,
    即是欧拉函数 ma=ouler(n/x)
    同理可得对 gcd(n-b,n) == n/x; mb=ouler(x);
    因为a与b可以转换 2*ma*mb
    但是当 x*x==n, ma*mb

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;
const long long mod=1e9+7;
long long ouler(long long p)
{
    long long ans=p;
    long long x=sqrt(p);
    for(int i=2;i<=x;i++)
    {
        if(p%i==0)
        {
            while(p%i==0)
                p=p/i;
            ans=ans/i*(i-1);
        }
    }
    if(p>1)
        ans=ans*(p-1)/p;
    return ans;
}
int main()
{
    long long n,k;
    while(~scanf("%I64d%I64d",&n,&k))
    {
        if(n==1)
        {
            printf("1\n");
            continue;
        }
        if(k==2)
        {
            printf("1\n");
            continue;
        }
        if(k>2)
        {
            printf("0\n");
            continue;
        }
        long long ans=0;
        for(int i=1;i*i<=n;i++)
        {
            if(n%i==0)
            {
                if(i*i==n)
                    ans+=ouler(i)*ouler(n/i);
                else ans+=2*ouler(i)*ouler(n/i);
                ans%=mod;
            }
        }
        printf("%I64d\n",ans);
    }
    return 0;
}