2016级算法第三次上机-G.Winter is coming
904 Winter is coming
思路
难题。首先简化问题, \(n\) 个0与 \(m\) 个1排成一列,连续的0不能超过x个,连续的1不能超过y个,求排列方法数。
显然会想到这是动态规划。最快想到的方法是 \(dp[i][j][x][y]\) 表示已经有i个北境兵j个野人参与排列,且末尾有x个连续北境士兵或y个连续野人士兵的方案数。这方法显然是正确的,但是光是 \(dp[200][200][10][10]\) 数组已经十分接近本题内存限制了,保证MLE。状态转移方法是大模拟,四层for循环,每次增加一个人放在最后,讨论各种情况。具体代码可见MLE参考代码,比较好理解。
不过这个方法已经和答案很接近了,只需要稍微优化一下。可以发现第三维和第四维很多空闲的空间被浪费了,我们没有必要用两个维度来分别记录有几个0或1,可以把第四维变成一个标志位,0代表北境军,1代表野人军,而第三维记录最后一段连续的人数,这样空间变成原来的1/6,算是简单的优化了一下,思想并没有变。具体可见优化代码,在此感谢孟尧提供。
本题还可以继续优化,换种方式,dp[i][j][k]:已经有i个北境兵j个野人参与排列,第三维k是标志位(0代表北境军,1代表野人军)的排列方法数。状态转移方程变为:: \(dp[i][j][0]=∑(dp[i-k][j][1])%MOD\) ;其中 \(k∈[1,min(i,x)]\) 。同理, \(dp[i][j][1]=∑(dp[i][j-k][0])%MOD\) ;其中 \(k∈[1,min(j,y)]\) 。相信你很快就能看懂,这里用 \(dp[i-k][j][1]\) 来代表最后有k个0,相当于同时把三四维合并了,巧妙至极。具体可见最优参考代码。
分析
本题不卡时间,卡的是内存。目的是让大家在解决问题的时候有优化的思想(实际的目的是把它从中等题变成难题)。
DP只能意会,不可言传。大家在做DP题的时候一定要理清思路,一般是先不管空间,毕竟以空间换时间,大多数题都是先卡时间再卡空间的。
以本题为例粗略讲解一下DP,以后不会再讲。记住DP具备的两个要素:最优子结构和子问题重叠,见《算法导论》225页。本问题明显备最优子结构,最少的排列数是由多个较短一点的最少排列数组成。DP的多层循环也是有规律的,因为子问题的重叠,你得先把子问题算出来,才能计算更深层的。这里i和j从小到大地计算,保证所加上的都是已经计算过的,才不会出现问题,如果这题加一个dp[i+1][j][1],那明显不对了,因为这个还没有计算过。状态转移方程有时候很微妙,需要一番数学推理。
最优参考代码
//
// Created by AlvinZH on 2017/10/24.
// Copyright (c) AlvinZH. All rights reserved.
//
#include <cstdio>
#include <cstring>
#include <iostream>
#define MOD 1000007
using namespace std;
int n, m, x, y;
int dp[205][205][2];
int main()
{
while(~scanf("%d %d %d %d", &n, &m, &x, &y))
{
memset(dp, 0, sizeof(dp));
for (int i = 0; i <= x; ++i)
dp[i][0][0] = 1;
for (int i = 0; i <= y; ++i)
dp[0][i][1] = 1;
for (int i = 1; i <= n; ++i) {
for (int j = 1; j <= m; ++j) {
for (int k = 1; k <= min(i,x); ++k)
dp[i][j][0] = (dp[i][j][0] + dp[i-k][j][1]) % MOD;
for (int k = 1; k <= min(j,y); ++k)
dp[i][j][1] = (dp[i][j][1] + dp[i][j-k][0]) % MOD;
}
}
printf("%d\n", (dp[n][m][0] + dp[n][m][1]) % MOD);
}
}
/* 把第三四维合并,因为我们只要在状态转移的时候保证最后连续一段不超过x或y就好了,第三维用来记录最后连续的是0还是1就好了。
* dp[i][j][k]:已经有i个北境兵j个野人参与排列,k为标志位(0代表北境军,1代表野人军)的排列方法数。
*/
优化参考代码
/*
Author: 孟爻(12593)
Result: AC Submission_id: 403884
Created at: Sun Nov 12 2017 23:05:10 GMT+0800 (CST)
Problem: 904 Time: 73 Memory: 11232
*/
#include <cstdio>
#include <cstring>
long f[205][205][15][2];
long M = 1000007;
long n,m,x,y;
int main() {
while(~scanf("%ld%ld%ld%ld",&n,&m,&x,&y))
{
memset(f,0,sizeof(f));
f[0][0][0][0]=1;
for(long i=0; i<=n; i++) {
for(long j=0; j<=m; j++) {
if(i) {
for(long k=1; k<=x; k++)
f[i][j][k][0]=(f[i][j][k][0]+f[i-1][j][k-1][0])%M;
for(long k=0; k<=y; k++)
f[i][j][1][0]=(f[i][j][1][0]+f[i-1][j][k][1])%M;
}
if(j) {
for(long k=1; k<=y; k++)
f[i][j][k][1]=(f[i][j][k][1]+f[i][j-1][k-1][1])%M;
for(long k=0; k<=x; k++)
f[i][j][1][1]=(f[i][j][1][1]+f[i][j-1][k][0])%M;
}
}
}
long ans=0;
for(long k=1; k<=x; k++)
ans=(ans+f[n][m][k][0])%M;
for(long k=1; k<=y; k++)
ans=(ans+f[n][m][k][1])%M;
printf("%ld\n",ans);
}
}
MLE代码
//
// Created by AlvinZH on 2017/10/24.
// Copyright (c) AlvinZH. All rights reserved.
//
#include <cstdio>
#include <cstring>
#define MOD 1000007
int n, m, x, y;
int dp[205][205][12][12];
int main()
{
while(~scanf("%d %d %d %d", &n, &m, &x, &y))
{
memset(dp, 0, sizeof(dp));
dp[0][0][0][0] = 1;
for (int i = 0; i <= n; ++i) {
for (int j = 0; j <= m; ++j) {
for (int k = 0; k <= x; ++k) {
for (int l = 0; l <= y; ++l) {
if(dp[i][j][k][l] == 0) continue;
if(i != n && k != x)//末尾是北境兵
{
dp[i+1][j][k+1][0] += dp[i][j][k][l];
dp[i+1][j][k+1][0] %= MOD;
}
if(j != m && l != y)//末尾是野人兵
{
dp[i][j+1][0][l+1] += dp[i][j][k][l];
dp[i][j+1][0][l+1] %= MOD;
}
}
}
}
}
int ans = 0;
for (int i = 1; i <= x; ++i)
ans = (ans + dp[n][m][i][0]) % MOD;
for (int i = 1; i <= y; ++i)
ans = (ans + dp[n][m][0][i]) % MOD;
printf("%d\n", ans);
}
}
/*
* dp[i][j][x][y]表示已经有i个北境兵j个野人参与排列,且末尾有x个连续北境士兵或y个连续野人士兵的方案数。
*/
2016级算法第三次上机-G.Winter is coming的更多相关文章
- 2016级算法第三次上机-C.AlvinZH的奇幻猜想——三次方
905 AlvinZH的奇幻猜想--三次方 思路 中等题.题意简单,题目说得简单,把一个数分成多个立方数的和,问最小立方数个数. 脑子转得快的马上想到贪心,从最近的三次方数往下减,反正有1^3在最后撑 ...
- 2016级算法第五次上机-G.ModricWang的撒币游戏
1062 ModricWang的撒币游戏 思路 此题为2017年ACM-ICPC亚洲区域赛乌鲁木齐赛区的A题,现场94个队中有38个队做出此题.在这里作为满分以外的题,是为了让大家看一下外面一些题的风 ...
- 2016级算法第三次上机-B.Bamboo和巧克力工厂
B Bamboo和巧克力工厂 分析 三条流水线的问题,依然是动态规划,但是涉及的切换种类比较多.比较易于拓展到n条流水线的方式是三层循环,外层是第k个机器手,里面两层代表可切换的流水线 核心dp语句: ...
- 2016级算法第六次上机-G.ModricWang likes geometry
1116 ModricWang likes geometry 思路 难题,非常考察几何知识,放在这里作为计算几何场次的最难的题. 原题地址 原版题解 代码
- 2016级算法第四次上机-G.ModricWang的序列问题 II
1021 ModricWang的序列问题II 思路 此题与上一题区别不是很大,只是增加了一个长度限制,当场通过的人数就少了很多. 大体解题过程与上一题相同.区别在于对\(f[]\) 的操作.没有长度限 ...
- 2016级算法第三次上机-F.ModricWang的导弹防御系统
936 ModricWang的导弹防御系统 思路 题意即为:给出一个长度为n的序列,求出其最长不降子序列. 考虑比较平凡的DP做法: 令\(nums[i]\) 表示这个序列,\(f[x]\) 表示以第 ...
- 2016级算法第三次上机-E.ModricWang's Polygons
930 ModricWang's Polygons 思路 首先要想明白,哪些多边形可能是格点正多边形? 分情况考虑: 三角形不可能,因为边长为有理数的正三角形的面积为无理数,而格点三角形的面积为有理数 ...
- 2016级算法第三次上机-D.双十一的抉择
915 双十一的抉择 思路 中等题.简化题目:一共n个数,分成两组,使得两组的差最接近0,就是说要使两组数都尽可能的接近sum/2. 思路还是很混乱的,不知道如何下手,暴力也挺难的,还不能保证对.想一 ...
- 2016级算法第三次上机-A.Bamboo的小吃街
A Bamboo的小吃街 分析 经典的两条流水线问题,题目描述基本类似于课件中的流水线调度,符合动态规划最优子结构性质 关键的动态规划式子为: dp[0][j] = min(dp[0][j - 1], ...
随机推荐
- HDFS Hbase
Hbase:Hbase是Hadoop的数据库,是bigtable的实现,基于HDFSHDFS:文件系统,是gfs的实现 Hive,Hbase,HDFS等之间的关系Hive:Hive不支持更改数据的操作 ...
- Linux下Maven的安装与使用
pache Maven,是一个软件(特别是Java软件)项目管理及自动构建工具,由Apache软件基金会所提供.基于项目对象模型(POM)概念,Maven利 用一个中央信息片断能管理一个项目的构建.报 ...
- Qt自定义插件编程小结
qt自定义组件开发步骤演示.以下所有步骤的前提是自己先编译Qtcreator源码,最好生成release版的QtCreator,否则自定义的插件嵌入QtCreator会失败!!!(这个网上教程很多) ...
- web.xml配置详解[转]
引文: 对于一个J2EE领域的程序员而言,基本上每天都会和web应用打交道. 什么是web应用?最简单的web应用什么样?给你一个web应用你该从何入手? 1.什么是web应用? web应用是一种可以 ...
- Excel 常用快捷键
Excel 常用快捷键 1. 移动整列 使用Shift快捷键可以快速移动整列:选中该列,当鼠标变成十字箭头时,按住Shift键,然后将该列移动到想要的位置. 2 绝对引用 使用F4快捷键可以快速设置绝 ...
- Sed命令n,N,d,D,p,P,h,H,g,G,x解析3
摘自:https://blog.csdn.net/WMSOK/article/details/78463199 Sed命令n,N,d,D,p,P,h,H,g,G,x解析 2017年11月06日 23: ...
- "软掩膜"和“硬掩膜”-智能IC卡
目录 一.“软掩膜”和“硬掩膜”... 2 二.EMV迁移进程... 3 三.PBOC规范和EMV规范对比... 3 四.总结... 5 五.关于SDA和DDA. 6 一.“软掩膜”和“硬掩膜” “软 ...
- 4.4.3 Java中的指针:Unsafe类
java多线程编程的无锁CAS底层都是通过 Unsafe进行操作的:源码如下 public final boolean compareAndSet(int expect, int update) { ...
- sql查询层级分类
先上个效果图吧 CTE递归查询里面用了一些小的技巧,查询出结果以后在前端用表格展示出来,层级视觉效果还是很明显的 with tree as(select [ID],[Name],[Address],[ ...
- PHP 5.2 5.3 5.4 5.5 memcache dll扩展
在windows下PHP5.2版本的memcache扩展dll文件好找,5.3的可能不是很好找,这里提供PHP5.2.5.3.5.4.5.5的php_memcache.dll扩展,需要的可以下载. 全 ...