欧拉函数

模板

#include<stdio.h>
const int MAXN=3000010;
bool  vis[MAXN];
int prime[MAXN/10];
long long euler[MAXN];
int main()
{
    int i,j,temp,top;
    int  n;
    top=0;
    euler[1]=1;
    for (i=2; i<MAXN; i++)
    {
        if (!vis[i])
        {
            prime[top++]=i;
            euler[i]=i-1;///质数p的欧拉函数=p-1
        }
        for(j=0;prime[j]*i<MAXN;j++)
        {
            vis[prime[j]*i]=1;
            if(i%prime[j]==0)///遇到最小素因数退出循环
            {
                ///因为i去除过prime[j]的倍数，所以不需要再次去除
                euler[i*prime[j]]=euler[i]*prime[j];
                 break;
            }
            else
            {
                ///因为i不含有prime[j]，需要去除prime[j]的倍数
                ///因为i与prime[j]互质
                ///=> euler[i*prime[j]]=euler[i]*euler[prime[j]];
                euler[i*prime[j]]=euler[i]*(prime[j]-1);
            }
        }
    }
    for(i=2;i<MAXN;i++)
    {
        euler[i]+=euler[i-1];///对欧拉函数求前缀和
    }
    int a,b;
    while(scanf("%d%d",&a,&b)!=EOF)
    {
        printf("%I64d
", euler[b]- euler[a-1]);
    }
    return 0;
}