仪仗队(欧拉函数模板)

 1 long long euler_phi(int n)
 2 {
 3     int m=(int)sqrt(n+0.5);
 4     int ans=n;
 5     for(int i=2;i<=m;i++)
 6 
 7         if(n%i==0)
 8         {
 9             ans=ans/i*(i-1);
10             while(n%i==0)
11                 n/=i;
12         }
13     if(n>1)ans=ans/n*(n-1);
14     return ans;
15 }

Description

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

Input

  共一个数N。

Output

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

Sample Input

  4

Sample Output

  9


Hint

【数据规模和约定】   对于 100% 的数据,1 ≤ N ≤ 40000

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <cstring>
 4 #include <algorithm>
 5 #include<cmath>
 6 #include<sstream>
 7 #include<string>
 8 using namespace std;
 9 long long euler_phi(int n)
10 {
11     int m=(int)sqrt(n+0.5);
12     int ans=n;
13     for(int i=2;i<=m;i++)
14 
15         if(n%i==0)
16         {
17             ans=ans/i*(i-1);
18             while(n%i==0)
19                 n/=i;
20         }
21     if(n>1)ans=ans/n*(n-1);
22     return ans;
23 }
24 
25 int main()
26 {
27     int n;
28     long long sum=0;
29     scanf("%d",&n);
30     for(int i=1;i<n;i++)
31        sum+=euler_phi(i);
32        sum=sum*2+1;
33        printf("%lld\n",sum);
34     return 0;
35 }

 

posted @ 2015-07-19 18:25  JoneZP  阅读(187)  评论(0编辑  收藏  举报