欧拉函数应用

例题一：

hdu 2588

题意：给出N和M，在 [ 1, N ] 中找出X，满足 gcd (X,N)>=M，问：有多少的X ？

题解：欧拉函数的简单应用，依题意可知： gcd (X,N)>=M ，设X=s*a ，N=s*b （a与b互质），易推出s为X和N的最大公约数，          又因为 X<=N，可得 a<=b（只有M为1时，a与b才可能相等），显然题目就变成的求与b互素的个数，直接枚举s就行了。

 

///hdu 2588

#include<cstdio>
#include<algorithm>
#include<cstring>

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()
{
//
//    for(LL t=1;t<=12;t++)
//        printf("t=%lld,%lld\n",t,Euler(t));
    LL n,m,ans,T;
    scanf("%d",&T);
    while(T--)
    {
        ans=0;
        scanf("%lld%lld",&n,&m);

        for(int i=1;i*i<=n;i++)///因为n太大了，故不能直接从1开始枚举到n，故要优化
        {   ///例如n为12，i只需枚举到3就行了，不用枚举4之后的了，因为在枚举3时，
            ///不但处理了s为3时的情况，也处理了s为4的情况，故只需枚举前sqrt（n）就行了

            if(n%i==0){ ///模为0，表示存在一个数b与i相乘为n
                if(i>=m) ans+=Euler(n/i);///此时i就为s，n/i就为b
                if(i*i!=n&&n/i>=m) ans+=Euler(i);///此时你n/i就为s，i就为b
                ///为什么要判断i*i！=n，因为i*i==n时，那么即是s和b相同，只需算一次就行了
            }
        }

        printf("%lld\n",ans);
    }
    return 0;
}

 

例题二：

bzoj 2818

这题跟上一题有点类似，Gcd(x,y)=S（S为素数），依题意，，设x<y，x=S*a，y=S*b，即得 a<b（a与b互质），即是求与b互素的个数，这里因为y不同上一题（要y=n），所以这里b<=y/s，故结果是前b位数的欧拉之和，处理下前缀和就行了。

 

///bzoj 2818

#include<cstdio>
#include<algorithm>
#include<cstring>

using namespace std;
typedef long long LL;

const int MAXN=1e7+10;
const int MAX=1e6+10;

int tot,n,pri[MAXN],phi[MAXN];
LL sum[MAXN],ans;
bool mark[MAXN];

void getphi() ///线性筛欧拉函数模板
{
    phi[1]=1;

    for(int i=2;i<=n;i++)
    {
        if(!mark[i]){
            phi[i]=i-1;
            pri[++tot]=i;
        }

        for(int j=1;j<=tot;j++)
        {
            int x=pri[j];

            if(i*x>n) break;
            mark[i*x]=1;

            if(i%x==0) {
                phi[i*x]=phi[i]*x;break;
            }
            else phi[i*x]=phi[i]*phi[x];
        }
    }
}
int main()
{
    while(~scanf("%d",&n))
    {
        tot=0;
        getphi();

        for(int i=1;i<=n;i++) sum[i]=sum[i-1]+phi[i]; ///处理前缀和


        for(int i=1;i<=tot;i++) ///每组小于等于n的素数
            ans+=sum[n/pri[i]]*2-1;///依题意（2,4）和（4,2）为不同结果，故结果为两倍减1，
            ///这里每次减的是（2,2），（3,3）多出来的一组
            ///例如n=4，当枚举到第一个素数2时，b=n/pri[i](4/2)=2,sum[2]=2
            ///说明有（1,2），（1,1）互质的对，即是（x，y）->（1*b,2*b）=(2,4),(2,2),故乘于2会多出一对（2,2）
        printf("%d\n",ans);
    }
    return 0;
}