歐拉函數~~

第一次寫歐拉函數，以前只知道其定義，以為編碼按其過程一步一步的實現應該就可以，結果發現不太一樣。

超時的想法：首先篩選素數，剛開始是篩選了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）的素數。

但是這樣之后還是超時了。。

之后看了百科上面的算法，他是對每一個素數，枚舉他的倍數，這樣每一步都是有效的，效率就提高了。

View Code 1 # include<stdio.h> 2 # include<string.h> 3 # include<stdlib.h> 4 # define N 3000005 5 bool vis[N]; 6 int phj[N]; 7 //__int64 sum[N]; 8 int main() 9 { 10 int i,j,a,b; 11 __int64 ans; 12 memset(vis,0,sizeof(vis)); 13 for(i=2;i*i<3000000;i++) 14 { 15 if(vis[i]) continue; 16 for(j=i+i;j<3000000;j+=i) 17 vis[j]=1; 18 } 19 for(i=3;i<3000000;i++) 20 phj[i]=i; 21 for(i=2;i<3000000;i++) 22 { 23 if(vis[i]==0) 24 { 25 for(j=i;j<3000000;j+=i) 26 phj[j]=phj[j]-phj[j]/i; 27 } 28 } 29 /*sum[2]=0; 30 for(i=3;i<=3000000;i++) 31 sum[i]=sum[i-1]+phj[i];*/ 32 while(scanf("%d%d",&a,&b)!=EOF) 33 { 34 ans=0; 35 for(i=a;i<=b;i++) 36 ans+=phj[i]; 37 printf("%I64d\n",ans); 38 //printf("%I64d\n",sum[b]-sum[a-1]); 39 } 40 return 0; 41 }

?

?

?

轉載于:https://www.cnblogs.com/183zyz/archive/2012/09/20/2687357.html

總結

以上是生活随笔為你收集整理的hdu 2824 The Euler function的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。