HDU-4556-Stern-Brocot Tree（欧拉函数）

Stern-Brocot Tree

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

解题思路：

首先对于这张图，不难发现，左右是对称的。那么只要计算右边*2+1 就可以求得答案。
Stern-Brocot Tree是一种能将所有的最简分数都表示出来的结构。
重点就在于 [能将所有的最简分数都表示出来] 且 题目要求 [去掉分子或者分母 大于n的项]。
那么这就又变成了另一个东西。
法里级数(Farey serires)。法里级数是介于0和1之间分母不超过N的所有最简分数的集合，且按照递增的次序排序。
如图：

知道了这个 就很好求了。
很明显，最简分数就是 gcd(x,y) == 1
那么每一行增加的项数就是 φ(n)个 。 求一个前缀和，就解决了。

AC代码：

#include <cstdio>
#include <vector>
#include <queue>
#include <cstring>
#include <cmath>
#include <map>
#include <set>
#include <stack>
#include <string>
#include <iostream>
#include <algorithm>
#include <iomanip>
using namespace std;
#define sd(n) scanf("%d",&n)
#define sdd(n,m) scanf("%d%d",&n,&m)
#define sddd(n,m,k) scanf("%d%d%d",&n,&m,&k)
#define pd(n) printf("%d\n", n)
#define pc(n) printf("%c", n)
#define pdd(n,m) printf("%d %d", n, m)
#define pld(n) printf("%lld\n", n)
#define pldd(n,m) printf("%lld %lld\n", n, m)
#define sld(n) scanf("%lld",&n)
#define sldd(n,m) scanf("%lld%lld",&n,&m)
#define slddd(n,m,k) scanf("%lld%lld%lld",&n,&m,&k)
#define sf(n) scanf("%lf",&n)
#define sc(n) scanf("%c",&n)
#define sff(n,m) scanf("%lf%lf",&n,&m)
#define sfff(n,m,k) scanf("%lf%lf%lf",&n,&m,&k)
#define ss(str) scanf("%s",str)
#define rep(i,a,n) for(int i=a;i<=n;i++)
#define per(i,a,n) for(int i=n;i>=a;i--)
#define mem(a,n) memset(a, n, sizeof(a))
#define debug(x) cout << #x << ": " << x << endl
#define pb push_back
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define mod(x) ((x)%MOD)
#define gcd(a,b) __gcd(a,b)
#define lowbit(x) (x&-x)
#define pii map<int,int>
#define mk make_pair
#define rtl rt<<1
#define rtr rt<<1|1
#define Max(x,y) (x)>(y)?(x):(y)
#define int long long


typedef pair<int,int> PII;
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
const int MOD = 1e9 + 7;
const ll mod = 10007;
const double eps = 1e-9;
const ll INF = 0x3f3f3f3f3f3f3f3fll;
//const int inf = 0x3f3f3f3f;
inline int read(){int ret = 0, sgn = 1;char ch = getchar();
while(ch < '0' || ch > '9'){if(ch == '-')sgn = -1;ch = getchar();}
while (ch >= '0' && ch <= '9'){ret = ret*10 + ch - '0';ch = getchar();}
return ret*sgn;}
inline void Out(int a){if(a>9) Out(a/10);putchar(a%10+'0');}
int qpow(int m, int k, int mod){int res=1%mod,t=m%mod;while(k){if(k&1)res=res*t%mod;t=t*t%mod;k>>=1;}return res;}
ll gcd(ll a,ll b){if(b > a) swap(a,b); return b==0?a : gcd(b,a%b);}
ll lcm(ll a,ll b){return a/gcd(a,b)*b;}
ll inv(ll x,ll mod){return qpow(x,mod-2,mod)%mod;}

const int N = 1e6+15;
int primes[N], euler[N], pcnt;
bool st[N];
void getEulers(int n)
{ //欧拉筛的扩展
	euler[1] = 1;
    pcnt = 0;
	for (int i=2; i <= n; i++ ){
		if ( !st[i] ){
			primes[pcnt++] = i;
			euler[i] = i-1;
		}
		for ( int j=0; primes[j] <= n/i; j++ ){
			st[ primes[j]*i ] = true;
			if ( i % primes[j] == 0 ) {
				euler[ i*primes[j] ] = euler[i] * primes[j];
				break;
			}
			euler[ i*primes[j] ] = euler[i] * ( primes[j]-1 );
		}
	}
}


signed main()
{
    getEulers(N);
    int t = 1,cas = 1;
    for(int i = 1 ; i <= (int)1e6+10 ; i ++)
        euler[i] += euler[i-1];
    int n;
    //cin>>t;
    while(cin>>n)
    {
        int res = 0;
        cout<<euler[n]*2+1<<endl;

    }
}