我是靠谱客的博主 勤劳月饼,最近开发中收集的这篇文章主要介绍loj#2540. 「PKUWC2018」随机算法,觉得挺不错的,现在分享给大家,希望可以做个参考。

概述

题目
思路
f [ S ] f[S] f[S]表示已经覆盖了 S S S的状态的最大独立集大小
g [ S ] g[S] g[S]表示此时的算法正确的方案数(最后除以 n ! n! n!就是期望了)
设新加入的扩充最大独立集的点为 i i i,与i相连的点的状态为 t o i to_i toi,那么
g [ S ∣ t o i ] = ∑ g [ S ] ⋅ A n − ∣ S ∣ − 1 ∣ t o i − S ∣ − 1 g[S|to_i]=sum g[S]cdot A^{|to_i-S|-1}_{n-|S|-1} g[Stoi]=g[S]AnS1toiS1(集合的减这里表示差集,因为“” m a r k d o w n markdown markdown里输不进QAQ)
排列的意义:已经选了 i i i这个点,那么与 i i i相连的点也要覆盖掉,这些点按题目中的算法随意排列都是正确的

#include<cstdio>
const int N=1<<20,M=998244353;
int n,m,i,x,y,to[N],f[N],g[N],s,v,F[21],G[21],cnt[N];
#define A(n,m) (1ll*F[n]*G[(n)-(m)]%M)
int main(){
	scanf("%d%d",&n,&m);
	for (;m--;) scanf("%d%d",&x,&y),x--,y--,to[x]|=1<<y,to[y]|=1<<x;
	F[0]=G[0]=F[1]=G[1]=1;
	for (i=2;i<=n;i++) F[i]=1ll*F[i-1]*i%M,G[i]=1ll*(M-M/i)*G[M%i]%M;
	for (i=1;i<=n;i++) G[i]=1ll*G[i]*G[i-1]%M;
	for (i=0;i<n;i++) to[i]|=1<<i;
	for (s=0;s<1<<n;s++) cnt[s]=cnt[s>>1]+(s&1);
	g[0]=1;
	for (s=0;s<1<<n;s++) if (g[s]){
		for (i=0;i<n;i++) if (!(s&(1<<i))){
			v=s|to[i];
			if (f[v]<f[s]+1) f[v]=f[s]+1,g[v]=0;
			if (f[v]==f[s]+1) g[v]=(g[v]+1ll*g[s]*A(n-cnt[s]-1,cnt[v^s]-1))%M;
		}
	}
	printf("%lld",1ll*g[(1<<n)-1]*G[n]%M);
}

最后

以上就是勤劳月饼为你收集整理的loj#2540. 「PKUWC2018」随机算法的全部内容,希望文章能够帮你解决loj#2540. 「PKUWC2018」随机算法所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(52)

评论列表共有 0 条评论

立即
投稿
返回
顶部