Description
作为体育委员,C君负责这次运动会仪仗队的训练。仪仗队是由学生组成的N * N的方阵,为了保证队伍在行进中整齐划一,C君会跟在仪仗队的左后方,根据其视线所及的学生人数来判断队伍是否整齐(如下图)。 现在,C君希望你告诉他队伍整齐时能看到的学生人数。
Input
共一个数N。
Output
共一个数,即C君应看到的学生人数。
Sample Input
4
Sample Output
9
好多大佬都是用欧拉函数做的,吊的一批
然而我只会反演
把去掉一行和一列后,可以看出gcd(行数,列数)=1的点才可被看见
然后还要加上观测点上边和右边的点
代码:
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #define M 50010 5 using namespace std; 6 7 int n,tot; 8 int P[M],mul[M],sum[M]; 9 bool prime[M]; 10 11 void get_mul() 12 { 13 mul[1]=1; 14 for(int i=2;i<=50000;i++) 15 { 16 if(!prime[i]) 17 { 18 P[++tot]=i; 19 mul[i]=-1; 20 } 21 for(int j=1;j<=tot&&P[j]*i<=50000;j++) 22 { 23 prime[P[j]*i]=true; 24 if(i%P[j]==0) 25 { 26 mul[P[j]*i]=0; 27 break; 28 } 29 else mul[P[j]*i]=-mul[i]; 30 } 31 } 32 for(int i=1;i<=50000;i++) sum[i]=sum[i-1]+mul[i]; 33 } 34 35 int cal(int n,int m,int k) 36 { 37 int ans=0; 38 n/=k; m/=k; 39 for(int l=1,r;l<=n&&l<=m;l=r+1) 40 { 41 r=min(n/(n/l),m/(m/l)); 42 ans+=(n/l)*(m/l)*(sum[r]-sum[l-1]); 43 } 44 return ans; 45 } 46 47 int main() 48 { 49 get_mul(); 50 scanf("%d",&n); n--; 51 if(!n) cout<<0; 52 else cout<<cal(n,n,1)+2; 53 return 0; 54 }