bzoj 1408: [Noi2002]Robot （DP+欧拉函数）

最新推荐文章于 2018-09-20 11:04:53 发布

原创最新推荐文章于 2018-09-20 11:04:53 发布 · 416 阅读

0 ·

CC 4.0 BY-SA版权

动态规划同时被 2 个专栏收录

250 篇文章

订阅专栏

数论

158 篇文章

订阅专栏

本文针对NOI2002 Robot题目提供了一种解决方案，采用DP加欧拉函数的方法，通过计算phi值得到政客、军人和学者的数量。文章详细介绍了算法思路并提供了完整的代码实现。

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

1408: [Noi2002]Robot

Time Limit: 5 Sec Memory Limit: 64 MB
Submit: 559 Solved: 375
[ Submit][ Status][ Discuss]

Description

Input

Output

Sample Input

3
2 1
3 2
5 1

Sample Output

8
6
75

HINT

90号机器人有10个老师，加上它自己共11个。其中政客只有15号；军人有3号和5号；学者有8个，它们的编号分别是：2,6,9,10,18,30,45,90。

Source

[ Submit][ Status][ Discuss]

题解：DP+欧拉函数

这里的独立数其实就和phi。

欧拉函数phi是积性函数，那么如果i,j互质，那么phi(i*j)=phi(i)*phi(j)

那么我们用f[i][j]表示选到第i个数一共选了j个数的phi的和。

f[i][j]=f[i-1][j]+f[i-1][j-1]*(p[i]-1)

这样子可以分开计算奇偶，得到政客和军人，然后用总数-政客-军人得到学者。

那么总数就是sigma(d|n) phi(d)=n

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define N 1003
#define mod 10000
using namespace std;
int p[N],n,m,f[N][N];
int quickpow(int num,int x)
{
	int ans=1; int base=num%mod;
	while (x) {
		if (x&1) ans=base*ans%mod;
		x>>=1;
		base=base*base%mod;
	}
	return ans;
}
int main()
{
	freopen("a.in","r",stdin);
	scanf("%d",&n); m=1;
	int cnt=0; bool pd=true;
	for (int i=1;i<=n;i++) {
		cnt++;
		scanf("%d",&p[cnt]);
		int x; scanf("%d",&x);
		m=m*quickpow(p[cnt],x)%mod;
		if (x!=1) pd=false;
		if (p[cnt]==2) cnt--; 
	}
	m--;
	//cout<<m<<endl;
	//cout<<cnt<<endl;
	for (int i=1;i<=cnt;i++) 
	 f[i][1]=(f[i-1][1]+p[i]-1)%mod;
	for (int i=2;i<=cnt;i++)
	 for (int j=2;j<=i;j++) {
	 	f[i][j]=(f[i-1][j]+f[i-1][j-1]*(p[i]-1)%mod)%mod;
	 }
	int ans1=0,ans2=0;
	for (int i=2;i<=cnt;i+=2) ans1+=f[cnt][i],ans1%=mod;
	for (int i=1;i<=cnt;i+=2) ans2+=f[cnt][i],ans2%=mod;
	printf("%d\n%d\n",ans1,ans2);
	printf("%d\n",(m-ans1-ans2+mod+mod)%mod);
}