[BZOJ 2190][SDOI2008]仪仗队：欧拉函数

点击这里查看原题

显然，对于每个（x,y）（x,y均不为1），（x,y）能被看到的条件是gcd(x-1,y-1)=1，因此我们需要求gcd(x-1,y-1)=1的对数，由于第一排、第一列分别有一个可以看到的，因此答案还需要加2。

对每一排，x<=y中，满足条件的x的个数等于phi(y)；对每一列同理（但是（2,2）被算了两次），因此答案=2*sigma(phi(i))-1+2（1<=i<=n）

/*
User:Small
Language:C++
Problem No.:2705
*/
#include<bits/stdc++.h>
#define ll long long
#define inf 999999999
using namespace std;
const int M=4e5+5;
int phi[M],prime[M],cnt,n,ans;
bool not_prime[M];
int main(){
    freopen("data.in","r",stdin);//
    scanf("%d",&n);
    phi[1]=1;
    for(int i=2;i<=n;i++){
        if(!not_prime[i]){
            prime[++cnt]=i;
            phi[i]=i-1;
        }
        for(int j=1;j<=cnt&&prime[j]*i<=n;j++){
            not_prime[i*prime[j]]=1;
            if(i%prime[j]==0){
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            phi[i*prime[j]]=phi[i]*phi[prime[j]];
        }
    }
    for(int i=1;i<n;i++) ans+=phi[i];
    printf("%d\n",2*ans+1);
    return 0;
}