[SDOI2009]Bill的挑战——全网唯一 一篇容斥题解
全网唯一一篇容斥题解
Description
Solution
看到这个题,大部分人想的是状压dp
但是我是个蒟蒻没想到,就用容斥切掉了。
并且复杂度比一般状压低,
(其实这个容斥的算法,提出来源于ywy_c_asm)
(然而我知道了这个算法,竟然和他写的不一样,而且比他跑的快)
进入正题:
我们需要统计恰好满足匹配k个的情况。
那么,我们可以先找出来,恰好满足n个,n-1,n-2。。。k个的情况。
分别记为ans[i]
ans[i]怎么算呢?
先给出公式:
ans[i]=cal(i)-∑C(j,i)×ans[j] 其中,i+1<=j<=n
cal(i)表示,从n个中任意选择i个,对于所有选择的情况,的方案数的和。
cal(i)可以dfs暴力C(n,i)枚举,每次统计答案。计入tot
void dfs(int x,int has){if(x==n+1){if(has!=up) return;ll lp=1;for(int j=1;j<=len;j++){las=-1;for(int i=1;i<=up;i++){if(a[mem[i]][j]!='?'){if(las==-1){las=a[mem[i]][j]-'a';}else if(las!=a[mem[i]][j]-'a') return;}}if(las==-1)lp=(lp*26)%mod;}(tot+=lp)%=mod;return;}if(has<up) {mem[++cnt]=x;dfs(x+1,has+1);mem[cnt--]=0;}if(n-x>=up-has) dfs(x+1,has); }
至于后面减去的部分。就是容斥的内容了。
大家可以自己画一个韦恩图理解一下。
这里有一个例子:n=4
现在我们要算ans[2],也就是恰好匹配2个的T的方案数
就是黄色的部分。
红色的数字是这个区域被算cal(i)的次数。
可见,三个点的重复区域,由于有C(3,2)种方法选到,所以会被算C(3,2)次。
所以减去所有的ans[3]即可。
其他情况同理。
最后输出ans[1]
组合数打表。
理论复杂度:
O(n×len×2^15)
Code
#include<bits/stdc++.h> using namespace std; typedef long long ll; const int N=20; const int M=52; const int mod=1000003; char a[N][M]; int len; int n,t,k; int mem[N],cnt; ll ans[N]; ll c[N][N]; ll sum; ll tot;//tot measures int up;//choose int las; void dfs(int x,int has){//dfs计算tot if(x==n+1){if(has!=up) return;ll lp=1;for(int j=1;j<=len;j++){las=-1;for(int i=1;i<=up;i++){if(a[mem[i]][j]!='?'){if(las==-1){las=a[mem[i]][j]-'a';}else if(las!=a[mem[i]][j]-'a') return;//两个字符不一样,无合法方案 }}if(las==-1)lp=(lp*26)%mod;//如果都是‘?’可以随便填,否则只有一种 }(tot+=lp)%=mod;return;}if(has<up) {mem[++cnt]=x;dfs(x+1,has+1);mem[cnt--]=0;}if(n-x>=up-has) dfs(x+1,has); }void clear(){memset(ans,0,sizeof ans);sum=0;len=0; } int main() {for(int i=0;i<=N-1;i++){c[i][0]=1;for(int j=1;j<=i;j++){c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;}}scanf("%d",&t);while(t--){clear();//清空数组,其实没有必要 scanf("%d%d",&n,&k);for(int i=1;i<=n;i++){scanf("%s",a[i]+1);}len=strlen(a[1]+1);//长度 for(int i=n;i>=k;i--){//ans[i]计算 tot=0;up=i;dfs(1,0);sum=0;for(int j=i+1;j<=n;j++){//容斥的处理 (sum+=c[j][i]*ans[j])%=mod;}ans[i]=(tot-sum+mod)%mod;}printf("%lld\n",ans[k]);}return 0; }
转载于:https://www.cnblogs.com/Miracevin/p/9585609.html
[SDOI2009]Bill的挑战——全网唯一 一篇容斥题解相关推荐
- [BZOJ3523][Poi2014]KLO-Bricks——全网唯一 一篇O(n)题解+bzoj最优解
Description 有n种颜色的砖块,第i种颜色的砖块有a[i]个,你需要把他们放成一排,使得相邻两个砖块的颜色不相同,限定第一个砖块的颜色是start,最后一个砖块的颜色是end,请构造出一种合 ...
- bzoj千题计划207:bzoj1879: [Sdoi2009]Bill的挑战
http://www.lydsy.com/JudgeOnline/problem.php?id=1879 f[i][j] 表示匹配了i个字符,匹配字符串的状态为j的方案数 枚举下一个字符是什么 计算加 ...
- P2167 [SDOI2009]Bill的挑战
P2167 [SDOI2009]Bill的挑战 题意: 有n个长度一样的字符串,字符串的每一位是?或者确定的字母,,求与这 N 个串中的刚好 K 个串匹配的字符串 T 的个数 1<=N<= ...
- 洛谷 P2167 [SDOI2009]Bill的挑战
题目描述 输入输出格式 输入格式: 本题包含多组数据. 第一行:一个整数T,表示数据的个数. 对于每组数据: 第一行:两个整数,N和K(含义如题目表述). 接下来N行:每行一个字符串. 输出格式: 1 ...
- 【SDOI2009】Bill的挑战
Description Sheng bill不仅有惊人的心算能力,还可以轻松地完成各种统计.在昨天的比赛中,你凭借优秀的程序与他打成了平局,这导致Sheng bill极度的不满.于是他再次挑战你.这次 ...
- FPGA----ZCU106基于axi-hp通道的pl与ps数据交互(全网唯一最详)
1.大家好,今天给大家带来的内容是,基于AXI4协议的采用AXI-HP通道完成PL侧数据发送至PS侧(PS侧数据发送至PL侧并没有实现,但是保留了PL读取PS测数据的接口) 2.如果大家用到SoC这种 ...
- FPGA-ZCU106-PL侧读写ddr4(全网唯一)
1.由于一直在PL侧做算法,外设接口接触的比较少,目前只做了sfp的UDP传输,但是由于课题的原因需要将一部分PL计算数据存储,而RAM存储空间比较小,因此本次给大带来了ZCU106的PL侧读写ddr ...
- 挑战自我的1000+篇文章总结
挑战自我的1000+篇文章总结 本文收集和总结了有关挑战自我的1000+篇文章,由于篇幅有限只能总结近期的内容,想了解更多内容可以访问:http://www.ai2news.com/, 其分享了有关A ...
- 全网唯一的Microsoft?NETFramework3?5SP1简体中文完整版离线安装包安装netframewo
简体中文完整版~~[很多软件都需要它支持]~<<全网络唯一可离线安装的>>!~我自己组合成功的!大家需要吗!加精! 安装netframework4.0装这个版本,目前.net程 ...
最新文章
- 2022-2028年中国集装箱涂料行业市场研究及前瞻分析报告
- 转载 http://blog.csdn.net/dengta_snowwhite/article/details/6418384
- java模拟数据库压测_java应用的优化【转】
- error: mod_deflate has been requested but can not be built due to prerequisite failures
- putty network error: connection refused
- python if语句和循环语句
- springboot jar中没有主清单属性_Spring Boot 常见错误及解决方法
- 有经验的面试官都是如何快速判断程序员能力的?
- mysql中使用case when语句
- Android——最全的系统对话框(AlertDialog)详解
- 【日期工具类】DateUtils
- win11壁纸|windows11桌面壁纸
- 超像素池化弱监督语义分割
- Win 10 版NVIDIA GeForce GTX 1060显卡驱动的下载及飞桨(Paddle)的安装
- FSR402电阻式薄膜压力传感器
- hpux 内存 够不够_hpux查看内存大小
- 关于word2016中mathtype无法使用以及“由于宏安全设置,无法找到宏或宏已被禁用”的解决方案
- 用STM32F407ZET6的HAL库写一个串口接收,发送代码,支持ringbuff
- JS键盘对应Code
- VUE经典面试题2022年,面试常用题