BZOJ 2190 [SDOI2008]仪仗队 题解

2190: [SDOI2008]仪仗队

Time Limit: 10 Sec   Memory Limit: 259 MB
Submit: 735   Solved: 484
[ Submit][ Status]

Description

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

　　 　　

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

Input

　　共一个数N。

Output

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

Sample Input

　　4

Sample Output

　　9

HINT

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

Source

数论

 

【分析】：

      sigma{(a,b)=1} (0<整数a,b<n)转化为2 * sigma{φ(i)} - 1(0<整数i<n),用欧拉函数

【代码1】：<单求版>

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<iostream>
#include<vector>
#include<stack>
#include<queue>
using namespace std;
int N;
long long phi(int x)
{
      long long ans=x;
      for(int i=2;i*i<=x;i++)
      {
            if(!(x%i))
            {
                  ans=ans/i*(i-1);
                  while(!(x%i)) x/=i;
            }
      }
      if(x>1)   ans=ans/x*(x-1);
      return ans;
}
int main()
{
      //freopen("input.in","r",stdin);
	  //freopen("output.out","w",stdout); 
	  scanf("%d",&N);
	  long long ans=2;
	  for(int i=2;i<N;i++)
	        ans+=phi(i);
      printf("%lld\n",2*ans-1);
	  //system("pause");
      return 0;
}

【代码2】：<群求版>

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<iostream>
#include<vector>
#include<stack>
#include<queue>
using namespace std;
#define MAXN 40001 
int N,ans[MAXN],p[MAXN];
bool a[MAXN];
int main()
{
      //freopen("input.in","r",stdin);
	  //freopen("output.out","w",stdout); 
	  scanf("%d",&N);
	  if(N<=1)   return printf("0\n"),0;
	  int num=0;
	  for(int i=2;i<N;i++)
      {
            if(!a[i])
            {
                  p[num++]=i;
                  ans[i]=i-1;
            }
            for(int j=0;(j<num && i*p[j]<=N);j++)
            {
                  a[i*p[j]]=true;
                  if(!(i%p[j]))
                  {
                        ans[i*p[j]]=ans[i]*p[j];
                        break;
                  }
                  else
                        ans[i*p[j]]=ans[i]*(p[j]-1);
            }
      }
      long long ans1=2;
      for(int i=2;i<N;i++)
            ans1+=ans[i];
      cout<<2*ans1-1<<endl;
	  //system("pause");
      return 0;
}

转载注明出处： http://blog.csdn.net/u011400953