Stern-Brocot Tree HDU - 4556
上图是一棵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
题意:
法里数列长度f[n]=f[n-1]+φ(n)(这里的长度值0-1的长度)
因为这个题还有对称的另一半所以只需要欧拉函数乘2即可f[n]=f[n-1]+2*φ(n)
code:
#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
#define maxn 1000005
#define ll __int64
ll phi[maxn],ans[maxn];
void init()
{
phi[1]=1;
for(int i=2;i<maxn;i++)
{
if(!phi[i])
{
for(int j=i;j<maxn;j+=i)
{
if(!phi[j])
phi[j]=j;
phi[j]-=phi[j]/i;
}
}
}
ans[1]=3;
for(int i=2;i<maxn;i++)ans[i]=ans[i-1]+2*phi[i];
}
int main()
{
init();
int n;
while(~scanf("%d",&n))
{
printf("%I64d\n",ans[n]);
}
return 0;
}