全网唯一一篇容斥题解

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. TNS-12560: TNS: 协议适配器错误 Oracle11g 创建数据库中问题处理(必须运行Netca以配置监听程序)
  2. hdu 2686(多线程dp)
  3. deferred对象和promise对象(二)---deferred对象
  4. Android 高级编程 RecyclerView 控件的使用
  5. mysql哪些字段适合索引_(转)Mysql哪些字段适合建立索引
  6. Git Flow分支管理
  7. 轨道坐标系_天文坐标系分类
  8. PREV-52 小数第n位 (大数)
  9. 计算机外接硬盘存储器,电脑存“片”究竟是用内置硬盘好,还是外置的好?
  10. POJ1067 取石子游戏 Wythoff's game C语言
  11. 杂项-数学软件:MATLAB
  12. EndNote X8参考文献样式下载与自定义
  13. 如何用Mockplus快速做一个手风琴菜单?
  14. Win7远程桌面连接不上问题解决方案
  15. 域名查询服务商的方法
  16. 复旦大学2017--2018学年第一学期(17级)高等代数I期末考试第八大题解答
  17. 一键智能抠图-原理与实现
  18. 正态分布以及推测统计简述
  19. 手把手教你使用NBS
  20. 电子招标采购平台源码php,电子招投标平台_企业招标管理系统_招标采购管理软件-高亚科技8manage SRM...

热门文章

  1. java uuid 效率_java uuid第一次性能
  2. python基础(文件、异常、模块、类、对象)
  3. enkey java_近期的Java项目(前端)
  4. 手把手教你写个小程序定时器管理库
  5. 面试体验:Facebook 篇(转)
  6. react 引入 mobx @babel/core: 7.2.2
  7. jmeter学习指南之聚合报告
  8. 近半年能力没进步原因分析与求助
  9. 原型模式——创建型模式
  10. [转]Cordova + Ionic in Visual Studio - 101 Tutorial [Part I]