HDU4556(递推求欧拉函数)

 

HDU4556(递推求欧拉函数)

分类： 数论 2013-05-21 18:18  61人阅读  评论(0)  收藏  举报

题目：Stern-Brocot Tree

 

本题跟POJ3090和POJ2478基本一个意思，很简单。

[cpp]  view plain copy

1. #include <iostream>  
2. #define N 1000005  
3. typedef unsigned long long LL;  
4.   
5. LL p[N];  
6.   
7. void phi()  
8. {  
9.     LL i,j;  
10.     for(i=1;i<N;i++) p[i]=i;  
11.     for(i=2;i<N;i+=2) p[i]>>=1;  
12.     for(i=3;i<N;i+=2)  
13.     {  
14.         if(p[i]==i)  
15.         {  
16.             for(j=i;j<N;j+=i)  
17.             {  
18.                 p[j]=p[j]/i*(i-1);  
19.             }  
20.         }  
21.     }  
22. }  
23.   
24. int main()  
25. {  
26.     LL n;  
27.     phi();  
28.     for(LL i=2;i<N;i++)  
29.         p[i]+=p[i-1];  
30.     while(std::cin>>n)  
31.     {  
32.         std::cout<<2*p[n]+1<<std::endl;  
33.     }  
34.     return 0;  
35. }