【JZOJ 3927】【NOIP2014模拟11.6】可见点数 (欧拉函数)

问题描述
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);
}