[SDOI2008]仪仗队

最新推荐文章于 2022-05-17 11:22:49 发布

原创最新推荐文章于 2022-05-17 11:22:49 发布 · 364 阅读

0 ·

CC 4.0 BY-SA版权

数学专栏收录该内容

97 篇文章

订阅专栏

本文详细解析了SDOI2008仪仗队问题，通过应用欧拉函数解决队伍整齐时可见学生人数的问题。文章介绍了如何筛选前n-1个欧拉函数的值并进行累加，最终得出正确答案。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

[SDOI2008]仪仗队

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

Input
　　共一个数N。

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

Sample Input
　　4

题解：欧拉函数的应用

我们知道除了 $(1, 0) (1, 1) (0, 1)$ 这三个特殊点，我们其他的点要满足条件必须要满足x,y坐标互质，所以我们这里就寻找前n-1个欧拉函数的值然后累加起来，因为我们这个答案只是y=x为对称轴的一半，所以我们还需要乘以2并且加上这个1。至于为什么筛1-(n-1)而且只加上1看这个图就明白了。
在这里插入图片描述

#include<bits/stdc++.h>
#define int long long
using namespace std;
typedef long long ll;
const int maxn=1e6+10;
bool ok[maxn];
int prime[maxn],phi[maxn],mu[maxn],cnt;
void sieve()
{
	phi[1]=mu[1]=1;
	for(ll i=2;i<maxn;++i)
	{
		if(!ok[i])
		{
			prime[cnt++]=i;
			phi[i]=i-1;
			mu[i]=-1;
		}
		for(int j=0;j<cnt;++j)
		{
			ll k=i*prime[j];
			if(k>=maxn)break;
			ok[k]=1;
			if(i%prime[j]==0)
			{
				phi[k]=phi[i]*prime[j];
				mu[k]=0;
				break; 
			}
			else
			{
			 phi[k]=phi[i]*(prime[j]-1);
			 mu[k]=-mu[i];
		    }
		}
	}
}
signed main()
{
	sieve();
    	int n; cin>>n;
	int res=0;
	for(int i=1;i<n;i++) res+=phi[i];
	cout<<(res<<1)+1<<endl;
}