Bzoj - 1499 [Noi2005] 瑰丽华尔兹 [DP][单调队列优化]

1499: [NOI2005]瑰丽华尔兹

Time Limit: 3 Sec   Memory Limit: 64 MB
Submit: 1715   Solved: 1043
[ Submit][ Status][ Discuss]

Description

你跳过华尔兹吗？当音乐响起，当你随着旋律滑动舞步，是不是有一种漫步仙境的惬意？众所周知，跳华尔兹时，最重要的是有好的音乐。但是很少有几个人知道，世界上最伟大的钢琴家一生都漂泊在大海上，他的名字叫丹尼•布德曼•T.D.•柠檬•1900，朋友们都叫他1900。 1900在20世纪的第一年出生在往返于欧美的邮轮弗吉尼亚号上，很不幸他刚出生就被抛弃了，成了孤儿。1900孤独的成长在弗吉尼亚号上，从未离开过这个摇晃的世界。也许是对他命运的补偿，上帝派可爱的小天使艾米丽照顾他。可能是天使的点化，1900拥有不可思议的钢琴天赋：从未有人教，从没看过乐谱，但他却能凭着自己的感觉弹出最沁人心脾的旋律。当1900的音乐获得邮轮上所有人的欢迎时，他才8岁，而此时的他已经乘着海轮往返欧美大陆50余次了。虽说是钢琴奇才，但1900还是个孩子，他有着和一般男孩一样的好奇和调皮，只不过更多一层浪漫的色彩罢了：这是一个风雨交加的夜晚，海风卷起层层巨浪拍打着弗吉尼亚号，邮轮随着巨浪剧烈的摇摆。船上的新萨克斯手马克斯•托尼晕船了，1900招呼托尼和他一起坐上舞厅里的钢琴，然后松开了固定钢琴的闸，于是，钢琴随着海轮的倾斜滑动起来。准确的说，我们的主角1900、钢琴、邮轮随着1900的旋律一起跳起了华尔兹，随着“嘣嚓嚓”的节奏，托尼的晕船症也奇迹般的消失了。后来托尼在回忆录上写道：大海摇晃着我们使我们转来转去快速的掠过灯和家具我意识到我们正在和大海一起跳舞真是完美而疯狂的舞者晚上在金色的地板上快乐的跳着华尔兹是不是很惬意呢？也许，我们忘记了一个人，那就是艾米丽，她可没闲着：她必须在适当的时候施展魔法帮助1900，不让钢琴碰上舞厅里的家具。不妨认为舞厅是一个N行M列的矩阵，矩阵中的某些方格上堆放了一些家具，其他的则是空地。钢琴可以在空地上滑动，但不能撞上家具或滑出舞厅，否则会损坏钢琴和家具，引来难缠的船长。每个时刻，钢琴都会随着船体倾斜的方向向相邻的方格滑动一格，相邻的方格可以是向东、向西、向南或向北的。而艾米丽可以选择施魔法或不施魔法：如果不施魔法，则钢琴会滑动；如果施魔法，则钢琴会原地不动。艾米丽是个天使，她知道每段时间的船体的倾斜情况。她想使钢琴在舞厅里滑行路程尽量长，这样1900会非常高兴，同时也有利于治疗托尼的晕船。但艾米丽还太小，不会算，所以希望你能帮助她。

Input

输入文件的第一行包含5个数N, M, x, y和K。N和M描述舞厅的大小，x和y为钢琴的初始位置（x行y列）；我们对船体倾斜情况是按时间的区间来描述的，且从1开始计量时间，比如“在[1, 3]时间里向东倾斜，[4, 5]时间里向北倾斜”，因此这里的K表示区间的数目。以下N行，每行M个字符，描述舞厅里的家具。第i行第j列的字符若为‘ . ’，则表示该位置是空地；若为‘ x ’，则表示有家具。以下K行，顺序描述K个时间区间，格式为：si ti di。表示在时间区间[si, ti]内，船体都是向di方向倾斜的。di为1, 2, 3, 4中的一个，依次表示北、南、西、东（分别对应矩阵中的上、下、左、右）。输入保证区间是连续的，即 s1 = 1 si = ti-1 + 1 (1 < i ≤ K) tK = T

Output

输出文件仅有1行，包含一个整数，表示钢琴滑行的最长距离(即格子数)。

Sample Input

4 5 4 1 3
..xx.
.....
...x.
.....
1 3 4
4 5 1
6 7 3

Sample Output

6

这道题先从简单的方法考虑

暴力的方法是从时间角度DP，设f[t][x][y]为t时刻在x, y位置的最长距离。

转移方程显然

if(mp[x][y] == '.')

f[t][x][y] = max(f[t - 1][x][y], f[t - 1][x - mx[d]][y - my[d]] + 1)

mx[5] = {0, -1, 1, 0, 0} my[5] = {0, 0, 0, -1, 1} 表示移动的方向

即 从上一个合法位置加一 前提是该位置为'.'，否则f[t][x][y] = -INF

O(tnm)显然超时，而且数组开的大小不定，用滚动的话可以过7个点(Codevs)

现在想想优化

首先是状态定义，可以改为f[k][x][y]代表第k次移动后处于

x, y的最长距离（之所以那么定义是因为区间内的 转移是等价的，

同一时间段方向不变）

现在总的复杂度O(knm * 转移复杂度)

要解决转移的问题

以前的转移是相邻时刻间的转移，取决于区间长度即O(t - s + 1)

观察转移式(假设现在是x方向的转移)

f[k][x][y] = max(f[k - 1][x][y], f[k - 1][i][y] + x - i)(1 <= i <= n)

发现x可以提出来 所以就是我们要找所有前面的状态里最大的一个转移过来

所以我们将所有的f[k][x][y] - x(用x替换i 有点懵的话多看几遍这一段)

加入单调队列 维护一个单减队列 每次加入时去掉队尾的不优的答案，统计答案前

去掉 队首距离超过len的答案，取队首更新答案即可

写法有技巧 代码自己理解

#include <bits/stdc++.h>
#define INF (0x3f3f3f3f)
using namespace std;

int n, m, sx, sy, k, mx[5] = {0, -1, 1, 0, 0}, my[5] = {0, 0, 0, -1, 1};
int f[2][205][205], cur, que[205][2], h, t;
char mp[205][205];

void DP( int x, int y, int rg, int typ, int len )
{
	h = t = 0;
	for(int i = 1; i <= rg; ++i)
	{
		if(mp[x][y] == '.')
		{
			int tmp = f[cur ^ 1][x][y] - i;
			while(h < t && tmp > que[t][0]) t--;
			que[++t][0] = tmp;
			que[t][1] = i;
			while(h < t && i - que[h + 1][1] > len) h++;
			f[cur][x][y] = que[h + 1][0] + i;
		}
		else f[cur][x][y] = -INF, h = t = 0;
		x += mx[typ], y += my[typ];
	}
} 

int main() 
{
	while(~scanf( "%d%d%d%d%d", &n, &m, &sx, &sy, &k ))
	{
		for(int i = 1; i <= n; ++i)
			scanf( "%s", mp[i] + 1 );
		for(int i = 1; i <= n; ++i)
			fill(f[0][i] + 1, f[0][i] + m + 1, -INF);
		f[0][sx][sy] = 0;
		int s, t, d, ans = 0;
		for(int i = 1; i <= k; ++i)
		{
			cur ^= 1;
			scanf( "%d%d%d", &s, &t, &d );
			if(d == 1) for(int j = 1; j <= m; ++j) DP(n, j, n, 1, t - s + 1);
			if(d == 2) for(int j = 1; j <= m; ++j) DP(1, j, n, 2, t - s + 1);
			if(d == 3) for(int j = 1; j <= n; ++j) DP(j, m, m, 3, t - s + 1);
			if(d == 4) for(int j = 1; j <= n; ++j) DP(j, 1, m, 4, t - s + 1);
		}
		for(int i = 1; i <= n; ++i)
			for(int j = 1; j <= m; ++j)
				ans = max(ans, f[cur][i][j]);
		cout << ans << endl;
	}
	return 0;
}