hdu 2824 The Euler function

欧拉函数~~

第一次写欧拉函数，以前只知道其定义，以为编码按其过程一步一步的实现应该就可以，结果发现不太一样。

超时的想法：首先筛选素数，刚开始是筛选了N以内的所有素数，然后对于每个数，从小到大枚举素数，

for(;;)//枚举素数

if(n%prime[i]==0)

{

   phj[n]=phj[n]-phj[n]/prime[i];

   while(n%prime[i]==0)

   {

        n/=prime[i];

    }

}

超时了之后，发现素数其实没必要枚举到N,只用枚举到sqrt(N)

因为对于每一个数n，for（）循环之后，如果n不为1，那n肯定是一个大于sqrt（N）的素数。

但是这样之后还是超时了。。

之后看了百科上面的算法，他是对每一个素数，枚举他的倍数，这样每一步都是有效的，效率就提高了。

# include<stdio.h>
# include<string.h>
# include<stdlib.h>
# define N 3000005
bool vis[N];
int phj[N];
//__int64 sum[N];
int main()
{
    int i,j,a,b;
    __int64 ans;
    memset(vis,0,sizeof(vis));
    for(i=2;i*i<3000000;i++)
    {
        if(vis[i]) continue;
        for(j=i+i;j<3000000;j+=i)
            vis[j]=1;
    }
    for(i=3;i<3000000;i++)
        phj[i]=i;
    for(i=2;i<3000000;i++)
    {
        if(vis[i]==0)
        {
            for(j=i;j<3000000;j+=i)
                phj[j]=phj[j]-phj[j]/i;
        }
    }
    /*sum[2]=0;
    for(i=3;i<=3000000;i++)
        sum[i]=sum[i-1]+phj[i];*/
    while(scanf("%d%d",&a,&b)!=EOF)
    {
        ans=0;
        for(i=a;i<=b;i++)
            ans+=phj[i];
        printf("%I64d\n",ans);
        //printf("%I64d\n",sum[b]-sum[a-1]);
    }
    return 0;
}