hdoj 4556 Stern-Brocot Tree(欧拉函数)

Stern-Brocot Tree

Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)
Total Submission(s): 518    Accepted Submission(s): 277

Problem Description

　　
　　
　　上图是一棵Stern-Brocot树，其生成规则如下：
　　从第1行到第n行，每行相邻两数a/b和c/d，产生中间数(a+c)/(b+d)，置于下一行中。将一行的分数（包括0/1,1/0），进行约分简化，则每一行（包括0/1,1/0,1/1），不会出现两个相同的分数。若分子或者分母大于n，则去掉该分数，将剩下的分数，从小到大排序，得到数列F。
　　现在请您编程计算第n行的数列F的个数。

 

Input

　　输入包含多组测试用例，每组输入数据是一个正整数n（n<=1000000）。

 

Output

　　对于每组的测试数据n，请输出第n行的数列F的个数。

 

Sample Input

1246

 

Sample Output

351325

 

<span style="font-size:18px;">#include<cstdio>#include<cstring>#include<cstdlib>#include<cmath>#include<algorithm>using namespace std;#define size 1000010__int64 euler[size];__int64 ans[size];void init(){  memset(euler,0,sizeof(euler));  euler[1]=1;  for(int i=2;i<size;i++)   if(!euler[i])    for(int j=i;j<size;j+=i)     {     if(!euler[j])      euler[j]=j;      euler[j]=euler[j]/i*(i-1); }}void fun(){  __int64 sum=0;  for(int i=1;i<size;i++)  {    sum+=euler[i];ans[i]=sum*2+1;  }}int main(){   init();   fun();   int n;   while(~scanf("%d",&n))   printf("%I64d\n",ans[n]);  return 0;} </span>