题目描述
作为体育委员,C君负责这次运动会仪仗队的训练。仪仗队是由学生组成的N * N的方阵,为了保证队伍在行进中整齐划一,C君会跟在仪仗队的左后方,根据其视线所及的学生人数来判断队伍是否整齐(如下图)。 现在,C君希望你告诉他队伍整齐时能看到的学生人数。
输入格式:
共一个数N
输出格式:
共一个数,即C君应看到的学生人数。
数据规模和约定
对于 100% 的数据,1 ≤ N ≤ 40000
Analysis
终于搞到数论专题lor 我才没跳过dp呢
首先说一个定理,只有xy坐标互质的点才能在原点被看到,这个可以用相似或者什么鬼的证明,不说了
在这题里面我们可以发现可以看到的人是沿矩阵对角线对称的,那么我们可以只考虑一半的规模。接着就要求
∑i<ni=0ϕ(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;
}