【PKUWC2018】随机算法【状压dp】【组合计数】

题意:一张nnn个点mmm条边的无向无权图,求以下算法计算最大独立集的正确率:随机一个排列,依次考虑排列中每一个点,如果不与任何一个独立集中的点相邻则将其加入独立集。模998244353998244353998244353

n≤20n\leq20n20

显然这是个计数问题。

显然是个状压dp。

f(S,i)f(S,i)f(S,i)表示当前已经考虑完了SSS,算出最大独立集大小为iii的方案数。

枚举接下来考虑的点jjj,设jjj及与其相邻的点中未被考虑的点的集合为TTT

因为相邻的关系是相互的,所以jjj一定是TTT中第一个被考虑的(废话)

在考虑完jjj后,剩下的∣T∣−1|T|-1T1个点可以在后面随便放,并不影响后面的计算,所以认为TTT中的点都考虑了。

f(S∣T,i+1)=∑f(S,i)An−∣S∣−1∣T∣−1f(S| T,i+1)=\sum f(S,i)A_{n-|S|-1}^{|T|-1}f(ST,i+1)=f(S,i)AnS1T1

最终答案为f(2n−1,mx)/n!f(2^n-1,mx)/n!f(2n1,mx)/n!,其中mxmxmx为真正的最大独立集大小

复杂度O(2nn2)O(2^nn^2)O(2nn2)

如果被卡常了可以在dp的时候顺便把mxmxmx算出来,不用单独再算

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
using namespace std;
const int MOD=998244353;
typedef long long ll;
inline int qpow(int a,int p)
{int ans=1;while (p){if (p&1) ans=(ll)ans*a%MOD;a=(ll)a*a%MOD;p>>=1;}return ans;
}
int fac[25],finv[25];
int g[25][25];
int T[25],cnt[1<<20],dp[25][1<<20];
int main()
{for (int i=1;i<(1<<20);i++) cnt[i]=cnt[i^(i&-i)]+1;int n,m;scanf("%d%d",&n,&m);fac[0]=1;for (int i=1;i<=n;i++) fac[i]=(ll)fac[i-1]*i%MOD;finv[n]=qpow(fac[n],MOD-2);for (int i=n-1;i>=0;i--) finv[i]=(ll)finv[i+1]*(i+1)%MOD;for (int i=1;i<=m;i++){int u,v;scanf("%d%d",&u,&v);--u,--v;g[u][v]=g[v][u]=1;}for (int i=0;i<n;i++){T[i]=1<<i;for (int j=0;j<n;j++)T[i]|=g[i][j]<<j;}dp[0][0]=1;int ans=0;for (int S=0;S<(1<<n);S++)for (int i=0;i<=cnt[S];i++){if (!dp[i][S]) continue;ans=max(ans,i);for (int j=0;j<n;j++)if (!(S&(1<<j)))dp[i+1][S|T[j]]=(dp[i+1][S|T[j]]+(ll)dp[i][S]*fac[n-cnt[S]-1]%MOD*finv[n-cnt[S|T[j]]])%MOD;}printf("%lld\n",(ll)dp[ans][(1<<n)-1]*finv[n]%MOD);return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/315010.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

架构杂谈《一》

架构杂谈《一》从传统单体架构到服务化架构的发展历程典型的单体架构分为三个层级&#xff0c;Web层、业务逻辑层和数据存储层&#xff0c;每个层的指责分别如下&#xff1a;Web 层&#xff1a;负责与用户交互或者对外提供接口业务逻辑层&#xff1a;为了实现业务逻辑而设计的流…

hdu 3308 LCIS 线段树 + 区间合并

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 日常水一篇题解。 带修改的求区间连续的递增序列&#xff0c;我们考虑用线段树维护。 直接维护mlenmlenmlen是区间最长的递增序列&#xff0c;lslsls是从左端点开始的最长递增序列&#xf…

【ARC103D】Distance Sums【树论】【构造】

题意&#xff1a;一棵nnn个点的无权树&#xff0c;给定每个点iii到其他所有点的距离之和did_idi​&#xff0c;保证did_idi​两两不同。构造或判断无法构造一棵满足条件的树。 n≤105n\leq 10^5n≤105 首先对于非根结点uuu&#xff0c;有 dudfau−sizu(n−sizu)d_ud_{fa_u}-si…

从零开始实现ASP.NET Core MVC的插件式开发(二) - 如何创建项目模板

标题&#xff1a;从零开始实现ASP.NET Core MVC的插件式开发(二) - 如何创建项目模板作者&#xff1a;Lamond Lu地址&#xff1a;https://www.cnblogs.com/lwqlun/p/11155666.html源代码&#xff1a;https://github.com/lamondlu/DynamicPlugins在上一篇中&#xff0c;我们介绍…

Codeforces Round #112 (Div. 2) E. Compatible Numbers sosdp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 以下假设all1<<22all1<<22all1<<22。 转化问题&#xff0c;对于每个aia_iai​我们都计算xaixorallxa_i\ \ xor \ \ allxai​ xor all&#xff0c;对于所有xxx的子集都…

【CF1349D】Slime and Biscuits【概率期望】【解方程】

题意&#xff1a;nnn个人共有mmm个饼干&#xff0c;每轮随机选一个饼干随机给一个另外的人&#xff0c;所有饼干都在一个人手里时游戏结束&#xff0c;求期望进行次数。模998244353998244353998244353。 n≤105,m≤3105n\leq10^5,m\leq3\times10^5n≤105,m≤3105 首先肯定是每…

多租户通用权限设计(基于 casbin)

所谓权限控制, 概念并不复杂, 就是确认某个操作是否能做, 本质上仅仅就是个bool判断.权限几乎是每个系统必不可少的功能, 和具体业务结合之后, 在系统中往往表现的非常复杂和难于控制, 很大部分原因是把权限和具体业务结合的太过紧密, 把业务的复杂度也加入到权限控制中来了.一…

Codeforces Round #225 (Div. 1) E. Vowels 容斥 + sosdp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个长度为333的串&#xff0c;串的每个字母都在a−za-za−z范围内&#xff0c;定义一个串合法当且仅当这个串中含有至少一个元音字母。现在他忘记了元音字母都有那几个&#xff0c;显然有2242^{24}224…

【CF1344D】Résumé Review【数学】【二分】

题意&#xff1a;给定长度为nnn的自然数序列aaa和整数kkk&#xff0c;要求构造长度为nnn的自然数序列bbb&#xff0c;使得0≤bi≤ai,∑bik0\leq b_i\leq a_i,\sum b_ik0≤bi​≤ai​,∑bi​k&#xff0c;且∑bi(ai−bi2)\sum b_i(a_i-b_i^2)∑bi​(ai​−bi2​)最大化。 n≤105…

容器化单页面应用中RESTful API的访问

最近在工作中&#xff0c;需要让运行在容器中的单页面应用程序能够访问外部的RESTful API。这个需求看起来并不困难&#xff0c;不过实现起来还是有些曲折的。在此&#xff0c;我就将这部分内容总结一下。在入正题之前&#xff0c;有个一问题&#xff0c;就是为什么要将单页面应…

2020牛客暑期多校训练营(第六场)H.Harmony Pairs 数位dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个nnn&#xff0c;问你有多少对0≤a≤b≤n0\le a\le b\le n0≤a≤b≤n满足S(a)>S(b)S(a)>S(b)S(a)>S(b)。其中S(a)S(a)S(a)表示aaa的十进制各位数之和。 思路&#xff1a; 显然的数位dpdpd…

【PKUWC2018】随机游走【Min-Max容斥】【树形dp】【FWT】

题意&#xff1a;给一棵nnn个点的无权树和xxx&#xff0c;qqq次询问&#xff0c;每次给定一个点集SSS,询问从xxx开始每次随机走一步&#xff0c;SSS中的每个点至少被经过一次的期望步数。 n≤18,q≤5000n\leq 18,q\leq 5000n≤18,q≤5000 题目求的相当于是SSS中的所有点 第一次…

LogoSharp:Logo语言的C#实现

上周在朋友中问了一圈&#xff0c;发现没有几个人知道Logo语言的&#xff0c;或许&#xff0c;这门古老的语言现在用的人已经非常少了&#xff0c;除了国外有少量的用户将Logo语言用于教育外&#xff0c;估计国内不会有人使用这门语言。其实&#xff0c;Logo语言本来也就是以编…

hdu 6899 Xor 数位dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 问你有多少对x,yx,yx,y满足以下条件&#xff1a; (1)x∈[0,A],y∈[0,B](1)x \in [0,A],y\in [0,B](1)x∈[0,A],y∈[0,B] (2)∣x−y∣≤K(2)|x-y|\le K(2)∣x−y∣≤K (3)xxory≤W(3)x \ \ xor \ \ y\le W(3)…

【PKUWC2018】Minimax【线段树合并】

题意&#xff1a;给定一棵nnn个点的二叉树&#xff0c;叶子的权值输入给定且互不相同&#xff0c;非叶子结点iii的权值有pip_ipi​的概率为儿子结点权值最大值&#xff0c;1−pi1-p_i1−pi​的概率为最小值。求根结点取每种值的概率。模998244353998244353998244353。 n≤3105n…

对数位dp的一些拙见

这不是一篇介绍数位dpdpdp的文章&#xff0c;只是我思考后的一些记录&#xff0c;怕以后就忘记了。 由于博主太菜不会组合数学&#xff0c;以下数位dpdpdp均采用记忆化搜索的方式。 首先最重要的就是状态设计了&#xff0c;正常来说数位dpdpdp的状态设计需要包含数的结构和状态…

Hello Kubernetes快速交互实验手册

“ K8S在线实验室提供了一个交互实验环境&#xff0c;现将这部分文章精简并翻译为中文&#xff0c;希望对那些不了解K8S的童鞋有帮助&#xff0c;能够快速玩起来有个感性认识...”原文&#xff1a;https://kubernetes.io/docs/tutorials翻译&#xff1a;Edison ZhouHello Kuber…

【ZJOI2018】历史【结论】【LCT思想】

题意&#xff1a;一棵nnn个点的树&#xff0c;每个点有权值aia_iai​&#xff0c;初始时给定。维护&#xff1a; 单点权值加上一个正数。 询问每个点恰好执行aia_iai​次 access 操作&#xff0c;任意安排顺序的条件下切换轻重链总次数的最大值。 n≤4105n\leq 4\times 10^5n…

2020 ICPC 上海 Sum of Log 数位dp + 状态剪枝

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 观察可以发现&#xff0c;由于i&j0i \And j0i&j0&#xff0c;所以log2(ij)log_2(ij)log2​(ij)表示的应该是ijijij的111的最高位。 一个显然的dpdpdp状态就出来了f[pos][flag1][fl…

从MVP到微软产品经理的几点心得

昨天&#xff0c;我应邀在微软2019年暑期实习生的入职培训活动上做了一个简短的分享 —— 以我自身的职业发展过程为例&#xff0c;给大家讲了一讲我是怎样从一个技术爱好者&#xff0c;到获得第一个MVP&#xff08;微软最有价值专家&#xff09;的荣誉奖项&#xff0c;以及后来…