计数dp
计数dp
计数类的$dp$没做过几个,所以之前都放到"思维"标签下了,后来发现原来这属于一类问题啊...搬过来了.
管道取珠:https://www.lydsy.com/JudgeOnline/problem.php?id=1566
题意概述:
有这样的两根管道接在一起,上边有$n$个珠子,下边有$m$个,每次可以从上或下取出一个(向右滑出去),只会有两种颜色的珠子,求不同的取珠序列的总数的平方和.(只要两个取出序列的珠子颜色都相同就视为相同的序列,不用考虑上下).$n,m<=300$
$shallwe$:完全可以作为$NOIP$ $D2T1$.
既然要求平方和,是不是要先求每种序列的方案数再平方求和呢?然而这样做是做不出来的.
出题人怎么说就怎么做是没有前途的...这里有一个挺神奇的转化:视为两个人各拿一个这样的管道取珠,当两个人取到的方案完全相同时答案$+1$,考虑这样做为什么是对的.比如取出某个序列有$a$种方法,那么这第一个人从这$a$种中任取一种,第二个人就有$a$种方法取出与之相同的序列,也就是平方了.太奇妙了.
这个转移其实非常简单,考虑一个最普通的:$dp[i][j][k][z]$表示第一个人在上边选了$i$个,下面选了$j$个,第二个人在上边选了$k$个,下边选了$z$个的方案数.既然序列要相同,那两个人取的珠子总数应该是相同的,可以直接删掉任意一维.现在时间复杂度已经能过了,但是空间复杂度又过不了了...显然需要用一个滚动数组,但是这样的状态没有办法滚动,所以重新设计:$dp[s][i][k]$,为了方便理解直接沿用了上边的下标名称,表示共取了$s$个,第一个人在上边取了$i$个,第二个人在上边取了$j$个的方案数.此时可以发现$s$每次只会增加$1$,滚动掉.
# include <cstdio>
# include <iostream>
# include <queue>
# include <cstring>
# include <string>
# define R register int
# define ll long long
# define mod using namespace std; const int maxn=;
int n,m;
char s[maxn];
int a[maxn],b[maxn];
int dp[][maxn][maxn]; int ad (int a,int b)
{
a=a+b;
if(a>=mod) a-=mod;
return a;
} int main()
{
scanf("%d%d",&n,&m);
scanf("%s",s+);
for (R i=n;i>=;--i)
if(s[i]=='A') a[n-i+]=;
scanf("%s",s+);
for (R i=m;i>=;--i)
if(s[i]=='A') b[m-i+]=;
dp[][][]=;
for (R s=;s<=n+m;++s)
{
int las=s&;
int no=las^;
memset(dp[no],,sizeof(dp[no]));
for (R i=;i<=n;++i)
for (R k=;k<=n;++k)
{
int j=s-i,z=s-k,x=dp[las][i][k];
if(!x) continue;
if(a[i+]==a[k+]) dp[no][i+][k+]=ad(dp[no][i+][k+],x);
if(a[i+]==b[z+]) dp[no][i+][k]=ad(dp[no][i+][k],x);
if(b[j+]==a[k+]) dp[no][i][k+]=ad(dp[no][i][k+],x);
if(b[j+]==b[z+]) dp[no][i][k]=ad(dp[no][i][k],x);}
}
printf("%d",dp[(n+m)&][n][n]);
return ;
}
管道取珠
中国象棋:https://www.lydsy.com/JudgeOnline/problem.php?id=1801
题意概述:在一个$n \times m$的棋盘上放任意多个炮,使得它们不能相互攻击,求方案数.$n,m<=100$
前置知识:炮的攻击方式
一开始以为是炮可以攻击同一行和同一列,然而并不是..."马后炮"了解一下?
转化题目:每一行,每一列不能出现超过$2$个棋子的方案数.
一开始想到了状压,但是状压压不了$100$位,即使用$bitset$强行压起来了也没法转移.后来发现具体是哪一行有两个棋子或者一个棋子并不重要,我们只需要知道有多少列有两个、一个、没有棋子就可以转移了,因为每列的本质其实是一样的。$dp[i][j][k]$表示前$i$行中$j$列有一个棋子,$k$列有两个棋子的方案数,转移比较麻烦.
·如果这一行一个都不放,可以直接继承上面的状态;
·如果放一个,可以放到已有一个的地方,那么这一行数量加一,方案数为$(j+1) \times dp[i-1][j+1][k-1]$,也可以放到本来为空的地方$dp[i-1][j-1][k] \times (m-(j-1)-k)$;
·如果放两个,可以都放到已有一个的地方$dp[i-1][j+2][k-2] \times c[j+2][2]$,也可以都放到空位上$dp[i-1][j-2][k] \times c[m-(j-2)-k][2]$,也可以一个放到已有一个的地方,另一个放到空位上$dp[i-1][j][k-1] \times (m-j-(k-1)) \times j$。($c$数组中是组合数)
# include <cstdio>
# include <iostream>
# include <cstring>
# define ll long long
# define R register int
# define mod using namespace std; int n,m,ans;
long long dp[][][];
int c[][]; int ad (int i,int j,int k)
{
long long a=;
a=dp[i-][j][k]%mod;
if(j) a=(a+dp[i-][j-][k]*(m-(j-)-k))%mod;
if(k) a=(a+(dp[i-][j+][k-]*(j+)))%mod;
if(j>=) a=(a+(dp[i-][j-][k]*c[m-(j-)-k][]))%mod;
if(k) a=(a+(dp[i-][j][k-]*(m-j-(k-))%mod*j))%mod;
if(k>=) a=(a+(dp[i-][j+][k-]*c[j+][]))%mod;
return a;
} int main()
{
scanf("%d%d",&n,&m);
dp[][][]=;
c[][]=;
for (R i=;i<=m+;++i)
{
c[i][]=;
for (R j=;j<=;++j)
c[i][j]=(c[i-][j-]+c[i-][j])%mod;
}
for (R i=;i<=n;++i)
for (R j=;j<=m;++j)
for (R k=;k+j<=m;++k)
dp[i][j][k]=ad(i,j,k);
for (R i=;i<=m;++i)
for (R j=;j+i<=m;++j)
ans=(ans+dp[n][i][j])%mod;
printf("%d",ans);
return ;
}
中国象棋
计数dp的更多相关文章
- HDU5800 To My Girlfriend 背包计数dp
分析:首先定义状态dp[i][j][s1][s2]代表前i个物品中,选若干个物品,总价值为j 其中s1个物品时必选,s2物品必不选的方案数 那么转移的时候可以考虑,第i个物品是可选可可不选的 dp[i ...
- CodeForces 176B Word Cut (计数DP)
Word Cut Time Limit:2000MS Memory Limit:262144KB 64bit IO Format:%I64d & %I64u Submit St ...
- [DP之计数DP]
其实说实在 我在写这篇博客的时候 才刚刚草了一道这样类型的题 之前几乎没有接触过 接触过也是平时比赛的 没有系统的做过 可以说0基础 我所理解的计数dp就是想办法去达到它要的目的 而且一定要非常劲非常 ...
- HDU4815/计数DP
题目链接[http://acm.hdu.edu.cn/showproblem.php?pid=4815] 简单说一下题意: 有n道题,每到题答对得分为a[ i ],假如A不输给B的最小概率是P,那么A ...
- HDU 6377 度度熊看球赛 (计数DP)
度度熊看球赛 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Subm ...
- [SDOI2010]地精部落[计数dp]
题意 求有多少长度为 \(n\) 的排列满足 \(a_1< a_2> a_3 < a_4 \cdots\) 或者 $a_1> a_2 < a_3 > a_4\cdo ...
- 【BZOJ】2111: [ZJOI2010]Perm 排列计数 计数DP+排列组合+lucas
[题目]BZOJ 2111 [题意]求有多少1~n的排列,满足\(A_i>A_{\frac{i}{2}}\),输出对p取模的结果.\(n \leq 10^6,p \leq 10^9\),p是素数 ...
- 【AtCoder】AGC022 F - Leftmost Ball 计数DP
[题目]F - Leftmost Ball [题意]给定n种颜色的球各k个,每次以任意顺序排列所有球并将每种颜色最左端的球染成颜色0,求有多少种不同的颜色排列.n,k<=2000. [算法]计数 ...
- 【BZOJ】4559: [JLoi2016]成绩比较 计数DP+排列组合+拉格朗日插值
[题意]n位同学(其中一位是B神),m门必修课,每门必修课的分数是[1,Ui].B神碾压了k位同学(所有课分数<=B神),且第x门课有rx-1位同学的分数高于B神,求满足条件的分数情况数.当有一 ...
随机推荐
- Effective C++ 50条款
条款1:尽量用const和inline而不用#define 以const 行使常量折叠,用inline 代替常用操作的宏定义,而且库里面有很多常用函数可用.当然不能抛弃宏,宏还是很有用滴.偶最近才发现 ...
- Oracle XE快捷版(速成版)的限制
1.CPU上限:无论把数据库安装在多少核的服务器上,都只会提供一个CPU核心的运算能力 2.安装和执行限制:只能安装一个实例且只能运行一个实例 3.用户数据上限:最大11G的用户数据 4.内存使用上限 ...
- HDU3416(KB11-O spfa+最大流)
Marriage Match IV Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others ...
- Codeforces687C(SummerTrainingDay03-D DP)
C. The Values You Can Make time limit per test:2 seconds memory limit per test:256 megabytes input:s ...
- 超强、超详细Redis入门教程
(1)什么是redis? Redis 是一个基于内存的高性能key-value数据库. (有空再补充,有理解错误或不足欢迎指正) (2)Reids的特点 Redis本质上是一个Key-Value类型的 ...
- python-状态模式
源码地址:https://github.com/weilanhanf/PythonDesignPatterns 说明: 在软件开发过程中,各种应用程序可能会根据不同的情况做出不同的处理.最直接的方案就 ...
- Jquery中on绑定的一些小坑
---恢复内容开始--- 今天我们来说说关于JQuery中事件绑定中on绑定的一些小问题,直接上代码了,大家拷下去就可以用 <!DOCTYPE html> <html lang=&q ...
- html5+css3图片旋转特效
<!DOCTYPE html> <html lang="en"> <head><meta charset="UTF-8" ...
- @NotNull、@NotEmpty、@NotBlank的区别
Spring中@NotNull.@NotEmpty.@NotBlank的区别@NotNull:用于基本数据类型@NotEmpty:用于集合类@NotBlank:用于String上面
- Java基础笔记(1) 语言 JAVA的历史 Java的搭建环境
本文除了搭建是重点,其他的都当阅读小说一样去看就好了,不想看可以直接抓住重点,我会改变颜色勾出重点! 英语是人与人交流沟通的重要方式之一.JAVA:是人与计算机沟通交流重要方式之一.我们除了用java ...