P5746 [NOI2002] 机器人M号

合集 - 题解(*•ω•)(20)

1.数塔2024-02-162.[洛谷P2014] 选课2024-02-173.玉蟾宫 （悬线DP）2024-03-014.MooFest POJ-1990 （树状数组）2024-02-195.情书密码 （树状数组）2024-02-186.[SCOI 2009] 迷路 (矩阵快速幂)2024-03-137.[BZOJ3306] 树2024-03-188.[NOI2002] 荒岛野人2024-03-289.SHUFFLE 洗牌 (扩展欧几里得+龟速乘)2024-03-2910.[SDOI2009] Bill的挑战 (状压DP)2024-04-0911.P4168 [Violet] 蒲公英 (莫队的强制在线)2024-04-2212.P3667 [USACO17OPEN] Bovine Genomics G （set容器+二分）2024-04-2413.P7903 兜心の顶（构造）2024-05-1214.[BZOJ2720 Violet 5]列队春游（概率期望+组合数学）2024-05-2915.划分大理石(多重背包)2024-06-14

16.P5746 [NOI2002] 机器人M号2024-06-24

17.苍与红的试炼(数位DP)2024-06-2518.【唐赛】高一小学期22024-07-0819.P10779 BZOJ4316 小 C 的独立集 (仙人掌DP)2024-07-2420.关于求合法括号子序列个数2024-09-09

收起

阅读目录

• P5746 [NOI2002] 机器人M号 题解

回到顶部

P5746 [NOI2002] 机器人M号 题解

头图

不会告诉你是因为洛谷提交题解被打回来好几次所以决定在博客上复制粘贴的

虽然...但是还是有一些思考量的...

分析

显然：

当 $i\ne 1$ 时

• $i$ 是 $j$ 的老师 $\iff i\mid j$

• 第 $i$ 号机器人的独立数为 $\phi (i)$

读题，设两个奇素数 $x,y$ ，分情况讨论：

• $x$ 为政客：则 $x\times y$ 为军人

• $x$ 为军人：则 $x\times y$ 为政客

所以政客和军人只能由不同的奇素数乘起来得到。

考虑DP。

设 $f(i,j)$ 为前 $i$ 个质数中，选奇数 $(j=1)$ 或偶数 $(j=0)$ 个质数的答案，则动态转移方程为

$$f(i,j)=f(i-1,j\mathrm{xor}1)\times \phi (p_i)+f(i-1,j)$$

学者的独立数只需快速幂求 $m$ 的值，再减去政客和军人，再减 $1$ (不明白减 $1$ 重新爬回去读题 怒)即可。

Elaina's Code

Elaina's Code

#include<bits/stdc++.h>
using namespace std;
const int N=1005;
#define rd read()
#define Elaina 0
inline int read(){
    int x=0,f=1;
    char ch=getchar();
    for(;!isdigit(ch);ch=getchar()) if(ch=='-') f=-1;
    for(;isdigit(ch);ch=getchar()) x=x*10+ch-'0';
    return x*f;
}

const int mod=10000;
int n,m,f[N][2];

int qpow(int a,int b){
	int ret=1;
	for(;b;b>>=1){
		if(b&1) ret=ret*a%mod;
		a=a*a%mod;
	}
	return ret;
}

signed main(){
	int k=rd;
	f[0][0]=1;
	int cnt=1;
	for(int i=1;i<=k;i++){
		int p=rd,e=rd;
		(cnt*=qpow(p,e))%=mod;
		f[i][1]=(f[i-1][0]*(p==2?0:(p-1))%mod+f[i-1][1])%mod;
		f[i][0]=(f[i-1][1]*(p==2?0:(p-1))%mod+f[i-1][0])%mod;
	}
	f[k][0]--;
	printf("%lld\n%lld\n%lld\n",(f[k][0]%mod+mod)%mod,f[k][1],((cnt-f[k][0]-f[k][1]-1)%mod+mod)%mod);
	
	return Elaina;
}

都看到这了，真的不点个赞吗(>ω<*)