hdu 4983 Goffi and GCD(欧拉函数）

最新推荐文章于 2021-02-04 00:16:09 发布

huanzhizun

最新推荐文章于 2021-02-04 00:16:09 发布

阅读量523

点赞数

分类专栏：数学文章标签： ACM 数学欧拉函数

本文链接：https://blog.csdn.net/u013665921/article/details/38816527

版权

数学专栏收录该内容

25 篇文章 0 订阅

订阅专栏

k!=1时情况很简单，下面仅讨论k=1时情况。x=gcd(n-a,n),则n/x=gcd(n-b,n),因为n-a可以取到0...n-1也就是1....n，所以完全可以去掉n-这个限制条件，即gcd(a,n)=x、gcd(b,n)=n/x时个数，因为a<=n，所以gcd(a,n)的个数=u[n/x],u是欧拉函数。所以原式等于sigma(u[n/x]*u[x])其中x是n的约数。

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
#define mod 1000000007
#define LL __int64
using namespace std;
int prim[81000],tot;
int u[1000005];
bool vis[1000005];
void init(void){
    int i;
    for(i=1;i<=1000000;i++) u[i]=i;
    for(i=2;i<=1000000;i++){
        if(!vis[i]){
            for(int j=i;j<=1000000;j+=i){
                u[j]=u[j]/i*(i-1);
                vis[j]=1;
            }
            prim[tot++]=i;
        }
    }
}
int get(int n){
    int i,s,q;
    q=n;
    s=n;
    for(i=0;i<tot&&(LL)prim[i]*prim[i]<=q;i++){
        if(n%prim[i]==0){
            s=s/prim[i]*(prim[i]-1);
            while(n%prim[i]==0) n=n/prim[i];
        }
        if(n==1) break;
    }
    if(n!=1){
        s=s/n*(n-1);
    }
    return s;
}
int main()
{
    int i,k;
    int n;
    LL s=0,q;
    init();
    while(scanf("%d%d",&n,&k)!=EOF){
        if(n==1){
            printf("1\n");
            continue;
        }
        s=0;
        if(k>2){ printf("0\n");
            continue;
        }
        if(k==2){
            printf("1\n");
            continue;
        }
        for(i=1;(LL)i*i<=n;i++){
            if(n%i==0){
                if(i<=1000000)
                q=u[i];
                else q=get(i);
                if(n/i<=1000000) q=q*u[n/i];
                else q=q*get(n/i);
                q=q%mod;
                if((LL)i*i!=n)
                q=(q*2)%mod;
                s+=q;
                s%=mod;
            }
        }
        printf("%d\n",s);
    }
}