全网唯一一篇容斥题解

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的挑战——全网唯一 一篇容斥题解相关推荐

  1. [BZOJ3523][Poi2014]KLO-Bricks——全网唯一 一篇O(n)题解+bzoj最优解

    Description 有n种颜色的砖块,第i种颜色的砖块有a[i]个,你需要把他们放成一排,使得相邻两个砖块的颜色不相同,限定第一个砖块的颜色是start,最后一个砖块的颜色是end,请构造出一种合 ...

  2. bzoj千题计划207:bzoj1879: [Sdoi2009]Bill的挑战

    http://www.lydsy.com/JudgeOnline/problem.php?id=1879 f[i][j] 表示匹配了i个字符,匹配字符串的状态为j的方案数 枚举下一个字符是什么 计算加 ...

  3. P2167 [SDOI2009]Bill的挑战

    P2167 [SDOI2009]Bill的挑战 题意: 有n个长度一样的字符串,字符串的每一位是?或者确定的字母,,求与这 N 个串中的刚好 K 个串匹配的字符串 T 的个数 1<=N<= ...

  4. 洛谷 P2167 [SDOI2009]Bill的挑战

    题目描述 输入输出格式 输入格式: 本题包含多组数据. 第一行:一个整数T,表示数据的个数. 对于每组数据: 第一行:两个整数,N和K(含义如题目表述). 接下来N行:每行一个字符串. 输出格式: 1 ...

  5. 【SDOI2009】Bill的挑战

    Description Sheng bill不仅有惊人的心算能力,还可以轻松地完成各种统计.在昨天的比赛中,你凭借优秀的程序与他打成了平局,这导致Sheng bill极度的不满.于是他再次挑战你.这次 ...

  6. FPGA----ZCU106基于axi-hp通道的pl与ps数据交互(全网唯一最详)

    1.大家好,今天给大家带来的内容是,基于AXI4协议的采用AXI-HP通道完成PL侧数据发送至PS侧(PS侧数据发送至PL侧并没有实现,但是保留了PL读取PS测数据的接口) 2.如果大家用到SoC这种 ...

  7. FPGA-ZCU106-PL侧读写ddr4(全网唯一)

    1.由于一直在PL侧做算法,外设接口接触的比较少,目前只做了sfp的UDP传输,但是由于课题的原因需要将一部分PL计算数据存储,而RAM存储空间比较小,因此本次给大带来了ZCU106的PL侧读写ddr ...

  8. 挑战自我的1000+篇文章总结

    挑战自我的1000+篇文章总结 本文收集和总结了有关挑战自我的1000+篇文章,由于篇幅有限只能总结近期的内容,想了解更多内容可以访问:http://www.ai2news.com/, 其分享了有关A ...

  9. 全网唯一的Microsoft?NETFramework3?5SP1简体中文完整版离线安装包安装netframewo

    简体中文完整版~~[很多软件都需要它支持]~<<全网络唯一可离线安装的>>!~我自己组合成功的!大家需要吗!加精! 安装netframework4.0装这个版本,目前.net程 ...

最新文章

  1. 2022-2028年中国集装箱涂料行业市场研究及前瞻分析报告
  2. 转载 http://blog.csdn.net/dengta_snowwhite/article/details/6418384
  3. java模拟数据库压测_java应用的优化【转】
  4. error: mod_deflate has been requested but can not be built due to prerequisite failures
  5. putty network error: connection refused
  6. python if语句和循环语句
  7. springboot jar中没有主清单属性_Spring Boot 常见错误及解决方法
  8. 有经验的面试官都是如何快速判断程序员能力的?
  9. mysql中使用case when语句
  10. Android——最全的系统对话框(AlertDialog)详解
  11. 【日期工具类】DateUtils
  12. win11壁纸|windows11桌面壁纸
  13. 超像素池化弱监督语义分割
  14. Win 10 版NVIDIA GeForce GTX 1060显卡驱动的下载及飞桨(Paddle)的安装
  15. FSR402电阻式薄膜压力传感器
  16. hpux 内存 够不够_hpux查看内存大小
  17. 关于word2016中mathtype无法使用以及“由于宏安全设置,无法找到宏或宏已被禁用”的解决方案
  18. 用STM32F407ZET6的HAL库写一个串口接收,发送代码,支持ringbuff
  19. JS键盘对应Code
  20. VUE经典面试题2022年,面试常用题

热门文章

  1. v-charts加载动画_加载动画-用户体验写作练习
  2. swiftui_SwiftUI的混合包
  3. 2021 年最值得了解的 Node.js 工具
  4. 知识付费不热了,得到们接下来故事怎么讲?
  5. lenos快速开发脚手架
  6. MQ消息队列之MSMQ
  7. RDS for MySQL Mysqldump常见问题及处理
  8. java中try 与catch的使用
  9. mysql内连接的自连接_mysql 内连接、外连接、自连接
  10. 猴子排圈求最后编号问题