仪仗队 [Bzoj 2190，SDOI2008]

题目地址——

---

仪仗队

---

【题目描述】

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

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

---

【输入描述】

共一个数 $n$。

---

【输出描述】

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

---

【输入样例】

4

---

【输出样例】

9

---

【Solution】

题目转化：本题等价于求 $(x,y)=1$ 且 $x<=n$ 且 $y<=n$ 的数对个数。
就把 $\varphi(i)$ 加起来就好了。

---

【Code】

#include <iostream>
#include <cstdio>
#include <cmath>

#define LL long long

using namespace std;

LL n,ans;
bool no_prime[10000010];
int prime[10000010];
int phi[10000010];
int first[10000010];

int main(){

    scanf("%lld",&n);n--;

    phi[1]=1;
    no_prime[1]=true;
    phi[0]=0;

    for(LL i=2;i<=n;i++){
        if(!no_prime[i]){
            prime[++prime[0]]=i;
            phi[i]=i-1; 
        } 
        for(LL j=1;prime[j]*i<=n;j++){
            no_prime[i*prime[j]]=true;
            first[i*prime[j]]=i;
            if(i%prime[j]==0){
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            else phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
    }

    for(LL i=1;i<=n;i++)ans+=phi[i]*2;

    printf("%lld\n",ans+1);

    return 0;
}