BZOJ 2190([SDOI2008]仪仗队-O(n)线性筛欧拉函数)

2190: [SDOI2008]仪仗队

Time Limit: 10 Sec Memory Limit: 259 MB

Submit: 521 Solved: 331

[Submit][Status][Discuss]

Description

　　作为体育委员，C君负责这次运动会仪仗队的训练。仪仗队是由学生组成的N * N的方阵，为了保证队伍在行进中整齐划一，C君会跟在仪仗队的左后方，根据其视线所及的学生人数来判断队伍是否整齐(如下图)。 　　

　　现在，C君希望你告诉他队伍整齐时能看到的学生人数。

Input

　　共一个数N。

Output

　　共一个数，即C君应看到的学生人数。

Sample Input

　　4

Sample Output

　　9

HINT

【数据规模和约定】

　　对于 100% 的数据，1 ≤ N ≤ 40000

线性筛欧拉函数的方法相信大家都懂

这里不赘述

总而言之，这题就是求有多少(i,j)满足gcd(i,j)=1,i<n,j<n,正宗欧拉函数。。。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<functional>
#include<algorithm>
#include<functional>
using namespace std;
#define MAXN (40000+10)
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
int n,a[MAXN],size=0,phi[MAXN]={0};
bool b[MAXN]={0};
int get_phi()
{
phi[1]=1;
Fork(i,2,n)
{
if (!b[i])
{
a[++size]=i;phi[i]=i-1;
}
For(j,size)
{
if (i*a[j]>MAXN) break;
b[i*a[j]]=1;
if (i%a[j]==0)
{
phi[i*a[j]]=phi[i]*a[j];break;
}
else phi[i*a[j]]=phi[i]*(a[j]-1);
}
}
}
int main()
{
cin>>n;
get_phi();
long long sum=0;
For(i,n-1) sum+=phi[i];
sum=sum*2+(n>=2);
//	For(i,n) cout<<phi[i]<<' ';
cout<<max(sum,(long long)0)<<endl;

return 0;
}