问题描述
ZPS经过长期的努力争取,终于成为了0901班的领操员,他要带领0901班参加广播操比赛。现在0901班的队伍可以看作是一个n*n的点阵,每个人都站在格点上。现在作为领操员的ZPS站(0,0)点,他想知道如果0901班的队伍站齐了,他能看到多少个人的脸(假设每个人的身高相同,体积相同)。
输入
一个正整数n。
输出
ZPS能看到多少个人的脸(当然他是看不到自己的脸的)。
样例输入
3
样例输出
5
算法讨论
显然当两个点的横纵坐标之差互质时才能看见,我们筛出phi后统计即可。因为是n*n的矩阵,所以只用统计对角线一侧的点数乘以2后再加上对角线上的一个就行了。
#include <cstdio>
#define MAX_N 100006
using namespace std;
int prime[MAX_N],v[MAX_N],phi[MAX_N],n,m;
long long s;
int main()
{
scanf("%d",&n);
phi[1]=1;
for (int i=2;i<=n;i++)
{
if (!v[i])
{
v[i]=i;
prime[++m]=i;
phi[i]=i-1;
}
for (int j=1;j<=m;j++)
if (prime[j]>v[i] || prime[j]>n / i)
break;
else
{
v[prime[j]*i]=prime[j];
if (i % prime[j]==0)
phi[i*prime[j]]=phi[i]*prime[j];
else
phi[i*prime[j]]=phi[i]*(prime[j]-1);
}
}
for (int i=1;i<n;i++)
s+=phi[i];
s=s*2+1;
if (n==1)
printf("0");
else
printf("%lld",s);
}