欧拉函数，求素数 

  1. /*=======================================================*\ 
  2.  |               递推求欧拉函数phi(i)
  3.  欧拉函数\varphi(n)是小于或等于n的正整数中与n互质的数的数目
  4. \*=======================================================*/ 
  6. #define N 3000000
  7. __int64 phi[N + 100];
  9. void Euler()
  10. {
  11. 	int i, j;
  12. 	for(i = 1; i < N; i++)
  13. 	    phi[i] = i;
  14. 	
  15. 	for(i = 2; i < N; i++)
  16. 	    if(i == phi[i])								//若i为素数 第一次不是
  17. 	        for(j = i; j < N; j += i)				//让j累加i
  18. 	            phi[j] = (phi[j] / i) * (i - 1);	//求欧拉函数
  19. }
  23. /*=========================================================*\ 
  24. 	  筛选法求一个数的最大素因子在素数表中是第几个素数
  25. 	
  26.     利用筛选法得到素数，每次用一个素数更新它的所有倍数，
  27. 				 更新为这个素数的标号，
  28. 	由于大的素数会更新小的素数，所以最后得到的就是最大素因子
  29. \*=========================================================*/ 
  31. #include <string>
  32. #define maxn 1000000
  33. int prime[maxn];
  34. void f()
  35. {
  36. 	memset(prime, 0, sizeof(prime));
  37. 	for(int i=2, flag=0; i<maxn; i++)
  38. 	{
  39. 		if(!prime[i])
  40. 		{
  41. 			flag++;
  42. 			for(int j=i; j<maxn; j+=i)
  43. 				prime[j] = flag;
  44. 		}
  45. 	}
  46. }





附件列表

     

    posted @ 2015-01-29 11:32  sober_reflection  阅读(389)  评论(0)  编辑  收藏  举报