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

1
2
4
6

 

Sample Output

3
5
13
25

 

<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
);
return 0;
} </span>