仪仗队_SDOI2008_洛谷2158_数论

题目描述

---

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

输入格式：

---

共一个数N

输出格式：

---

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

数据规模和约定

---

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

Analysis

---

终于搞到数论专题lor 我才没跳过dp呢
首先说一个定理，只有xy坐标互质的点才能在原点被看到，这个可以用相似或者什么鬼的证明，不说了
在这题里面我们可以发现可以看到的人是沿矩阵对角线对称的，那么我们可以只考虑一半的规模。接着就要求$\sum_{i=0}^{i < n}\phi(i)$
至于为什么要求欧拉函数，是因为它的定义即是比它小且与它互质的数的数量，累计求和就行了。筛一个欧拉函数的表是线性的

果然我还是不适合vim吗

Code

---

#include <stdio.h>
#include <string.h>
#define rep(i, a, b) for (int i = a; i <= b; i ++)
#define fill(x, t) memset(x, t , sizeof(x))
#define N 80001
using namespace std;
int prime[N], phi[N], s[N];
inline int read(){
    int x = 0, v = 1;
    char ch = getchar();
    while (ch < '0' || ch > '9'){
        if (ch == '-'){
            v = -1;
        }
        ch = getchar();
    }
    while (ch <= '9' && ch >= '0'){
        x = (x << 1) + (x << 3) + ch - '0';
        ch = getchar();
    }
    return x * v;
}
inline int getPhi(int n){
    fill(prime, 0);
    fill(phi, 0);
    fill(s, 0);
    int cnt = 0;
    phi[1] = 1;
    s[1] = 1;
    rep(i, 2, n){
        if (!s[i]){
            prime[++ cnt] = i;
            phi[i] = i - 1;
        }
        rep(j, 1, cnt){
            if (i * prime[j] > n){
                break;
            }
            s[i * prime[j]] = 1;
            if (i % prime[j] == 0){
                phi[i * prime[j]] = phi[i] * prime[j];
                break;
            }else{
                phi[i * prime[j]] = phi[i] * (prime[j] - 1);
            }
        }
        // printf("%d\n", phi[i]);
    }
}
int main(void){
    int n = read();
    getPhi(n);
    int ans = 0;
    rep(i, 1, n - 1){
        ans += phi[i];
    }
    printf("%d\n", ans * 2 + 1);
    return 0;
}