【poj3718】 Facer's Chocolate Dream
http://poj.org/problem?id=3718 (题目链接)
题意
给出${2}$个长度为${n}$的${01}$串,问是否存在${m}$个长度为${n}$的有三个位置为${1}$的$01$串它们的异或和正好等于那两个$01$串的异或和。求方案数。
Solution
我们把那两个$01$串异或起来得到一个$S$,那么问题就变成了是否存在$m$个长度为$n$的有且仅有三个位置为$1$的$01$串的异或和等于$S$。
我们这样进行考虑,假设$S$有$v$个位置上为$1$,先求出异或出来的串有$v$个位置为$1$的方案数,因为$v$个$1$可以任意组合,我们再除以一个组合数$C_{n}^{v}$就是答案。
再添加一个条件,我们假设选出的$m$个数不同顺序算是不同的方案,求出来以后再除以一个$m!$就可以了。
这样我们就可以dp啦。$f[i][j]$表示已经选了$i$个数,异或起来有$j$个$1$的方案数,那么最后答案就是$$ans=\frac{f[m][v]}{m!C_n^v}$$
考虑转移,有这几种情况:
\begin{aligned} f[i][j]&=f[i-1][j-3]\binom{n-j+3}{3} \\ &+f[i-1][j-1]\binom{n-j+1}{2}(j-1) \\ &+f[i-1][j+1](n-j-1)\binom{j+1}{2} \\ &+f[i-1][j+3]\binom{j+3}{3} \end{aligned}
这应该很好理解吧,就是新加入的第$i$个数,它会对之前$i-1$个数已经异或出来的数造成的影响有$4$种:将三个$0$变成$1$;将两个$0$变成$1$,一个$1$变成$0$;将一个$0$变成$1$,两个$1$变成$0$;将三个$1$变成$0$。
考虑到这样转移,并没有考虑第$i$个数与之前的$i-1$个数重复的情况。如果第$i$个数与之前某个数重复了,那么就意味着剩下的$i-2$个数异或起来还是会有$j$个$1$。所以我们还要容斥一下:
\begin{aligned} f[i][j]-=f[i-2][j]\binom{n}{3}(i-1) \end{aligned}
表示第$i$个数的取值有$\binom{n}{3}$种,它可以与之前$i-1$个数中的任意一个相等。
那么就做完了,最后转移的时候注意判一下边界。
细节
LL
代码
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#define LL long long
#define inf (1ll<<30)
#define MOD 10007
#define Pi acos(-1.0)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std; const int maxn=1010;
int C[maxn][maxn],f[maxn][maxn],fac[maxn];
int n,m,v;
char s1[maxn],s2[maxn]; int power(int a,int b) {
int res=1;
while (b) {
if (b&1) res=(LL)res*a%MOD;
b>>=1;a=(LL)a*a%MOD;
}
return res;
}
int main() {
for (int i=0;i<=1000;i++) C[i][0]=1;
for (int i=1;i<=1000;i++)
for (int j=1;j<=i;j++) C[i][j]=(C[i-1][j]+C[i-1][j-1])%MOD;
fac[0]=1;for (int i=1;i<=1000;i++) fac[i]=(LL)fac[i-1]*i%MOD;
for (int i=1;i<=1000;i++) fac[i]=power(fac[i],MOD-2);
while (scanf("%d%d",&n,&m)!=EOF && n && m) {
scanf("%s%s",s1+1,s2+1);v=0;
for (int i=1;i<=n;i++) v+=(s1[i]!=s2[i]);
memset(f,0,sizeof(f));f[0][0]=1;
for (int i=1;i<=m;i++)
for (int j=0;j<=n;j++) {
if (j>=3) (f[i][j]+=(LL)C[n-j+3][3]*f[i-1][j-3]%MOD)%=MOD;
if (j>=1) (f[i][j]+=(LL)C[n-j+1][2]*(j-1)%MOD*f[i-1][j-1]%MOD)%=MOD;
if (j+1<=n) (f[i][j]+=(LL)C[j+1][2]*(n-j-1)%MOD*f[i-1][j+1]%MOD)%=MOD;
if (j+3<=n) (f[i][j]+=(LL)C[j+3][3]*f[i-1][j+3]%MOD)%=MOD;
if (i>=2) (f[i][j]+=MOD-(LL)f[i-2][j]*(C[n][3]-i+2)%MOD*(i-1)%MOD)%=MOD;
}
printf("%lld\n",(LL)f[m][v]*fac[m]%MOD*power(C[n][v],MOD-2)%MOD);
}
return 0;
}
【poj3718】 Facer's Chocolate Dream的更多相关文章
- poj3718 Facer's Chocolate Dream
题目链接 正解:组合数+$dp$. 今天考试的题,考试的时候感觉自己有点脑残过头了.. 似乎发现了所有$1$其实都是一样的,然后不知道怎么强制每种物品只选一个.. 然后就写了一个所有物品可以选任意个的 ...
- 【BZOJ】2016: [Usaco2010]Chocolate Eating(二分)
http://www.lydsy.com/JudgeOnline/problem.php?id=2016 这些最大最小显然是二分. 但是二分细节挺多的...这里注意二分的区间,可以累计所有的可能,然后 ...
- poj 2411 Mondriaan's Dream 【dp】
题目:id=2411" target="_blank">poj 2411 Mondriaan's Dream 题意:给出一个n*m的矩阵,让你用1*2的矩阵铺满,然 ...
- 【POJ2411】Mondriaan's Dream(轮廓线DP)
[POJ2411]Mondriaan's Dream(轮廓线DP) 题面 Vjudge 题解 这题我会大力状压!!! 时间复杂度大概是\(O(2^{2n}n^2)\),设\(f[i][S]\)表示当前 ...
- 【poj2411】Mondriaan's Dream 状态压缩dp
AC传送门:http://vjudge.net/problem/POJ-2411 [题目大意] 有一个W行H列的广场,需要用1*2小砖铺盖,小砖之间互相不能重叠,问有多少种不同的铺法? [题解] 对于 ...
- 【Codeforces 449A】Jzzhu and Chocolate
[链接] 我是链接,点我呀:) [题意] 题意 [题解] 设最后行分成了x行,列分成了y列. 那么答案就是floor(n/x)floor(n/y) 然后x+y-2=k //即平均分配x行.y列 我们可 ...
- 【Demo】QQ,github,微博第三方社交登录
本文主要讲解 集成 第三方社交账号登录 为什么会有这个需求? 主要是因为目前互联网的网站数量太多,如果在各个站点都注册一个账号 用户非常不容易记住每个账号的用户名和密码,并且非常难保证每个账号的密码足 ...
- 分享50款 Android 移动应用程序图标【下篇】
在这个移动程序流行的时代,持续增长的应用程序经济充满了商业机遇.任何对应用程序设计感兴趣的人,将会喜欢上这里的50个独特的 Android 应用程序图标.这些例子中的图标能够让应用程序的设计更具吸引力 ...
- 【转】KM匹配题集
转自:http://blog.csdn.net/shahdza/article/details/7779324 [HDU]2255 奔小康赚大钱 模板题★1533 Going Home 模板题★242 ...
随机推荐
- Oracle-归档日志详解(运行模式、分类)
一.Oracle日志分类 分三大类: Alert log files--警报日志,Trace files--跟踪日志(用户和进程)和 redo log 重做日志(记录数据库的更改 ...
- test zhenai
web.Document.InvokeScript("eval",new string[]{"document.getElementById('passwordbt'). ...
- Python 语言简介
Python是一种计算机程序设计语言.你可能已经听说过很多种流行的编程语言,比如非常难学的C语言,非常流行的Java语言,适合初学者的Basic语言,适合网页编程的JavaScript语言等等. 那P ...
- anaconda安装opencv3
opencv是C和C++语言编写的,很多教程都是基于C++语言进行学习的,可是机器学习最多的库是python写的,所以还是学学python怎么安装opencv3, 面向学习的大都是使用了anacond ...
- Js_增删改Cookie的值
//获得cookie 的值function cookie(name) { var cookieArray = document.cookie.split("; "); //得到分割 ...
- Mysql基础命令(二)select查询操作
条件查询 使用Where进行数据筛选结果为True的会出现在结果集里面 select 字段 from 表名 where 条件; # 例: select * from test_table where ...
- Unity XLua 官方教程学习
一.Lua 文件加载 1. 执行字符串 using UnityEngine; using XLua; public class ByString : MonoBehaviour { LuaEnv lu ...
- 微软职位内部推荐-Senior NLP Scientist
微软近期Open的职位: Job Title: Senior NLP Scientist Location: Suzhou, China Suzhou, one of the most vibrant ...
- Scrum立会报告+燃尽图(十二月十一日总第四十二次):贡献分配和收集用户报告
此作业要求参见:https://edu.cnblogs.com/campus/nenu/2018fall/homework/2484 项目地址:https://git.coding.net/zhang ...
- About The Algorithm Simplification
For mode 1, you have to ergod all the data in the files. So the key point to solve this problem is t ...