单调队列优化DP

瑰丽华尔兹

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

输入格式

输入文件的第一行包含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

输出格式

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

Sample Input

4 5 4 1 3
..xx.
.....

...x.

.....

1 3 4

4 5 1

6 7 3

Sample Output

6

Hint

【样例说明】
　　钢琴的滑行路线：

　　钢琴在“×”位置上时天使使用一次魔法，因此滑动总长度为6。
【数据范围】
　　50%的数据中，1≤N, M≤200，T≤200；
　　100%的数据中，1≤N, M≤200，K≤200，T≤40000。
　　

 题目出处：P2254 [NOI2005] 瑰丽华尔兹

最容易想出来的DP状态就是设f[t,i,j]表示第t时刻[i,j]位置的最长距离
现在的最优值要么是这一秒原地不动，要么是这一秒从上下左右某个地方滑过来的

so可以很轻松地写出DP方程f[t,i,j]=max(f[t−1,i,j],f[t−1,x,y]),[x,y]]为合法的位置
答案为max(f[T,i,j]),T=last t，时间复杂度O(TNM),明显爆炸

把DP状态改一改：设f[t,i,j]表示第t个时间段[i,j]位置的最长距离
那么DP方程:

f[t,i,j]=min(f[t−1,i,j],f[t−1,x±dis,y±dis])
dis通过枚举得到，分四种情况考虑转移即可

时间复杂度O(min(N,M)KNM)，按正常情况来说，爆炸
怎么加优化呢？四个优先队列分别记录东西南北的最优值，每次碰到障碍物就清空某个队列就可以了
时间复杂度O(KNM)
 

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=210;
int gra[maxn][maxn],f[maxn][maxn],head,rear,n,m,x,y,k,s,e,dir,ans,dx[]={0,-1,1,0,0},dy[]={0,0,0,-1,1};
//gra存图，f表示当在这个点走过的最大路径长度， dx,dy表示向上下左右的移动 
struct queuenode
{
	int fnub,cnt;//fnub记录队节点本身的路程值，cnt在该行中的位置 
}queue[maxn];//单进双出队列 
char c;
void update(int x,int y,int direction,int distance)//每次一个时间段的更新，这个函数只能更新图的某一行或某一列 
{
	head=rear=0;
	for(int i=1;((x>=1&&x<=n)&&(y>=1&&y<=m));i++)//保证不出界 
	{
		if(gra[x][y]==1)
		{
			head=rear=0;
			x=x+dx[dir];//移动 
			y=y+dy[dir];
			continue;
		}
		while(head<rear&&queue[rear-1].fnub+(i-queue[rear-1].cnt)<=f[x][y])//队尾行驶到该点的距离没有该点原本的大，那么队尾出队 
		rear--;
		queue[rear].cnt=i;//入队 
		queue[rear++].fnub=f[x][y];
		while(head<rear&&queue[head].cnt+distance<i)//队头到该点的距离超过了改时间段允许的最大距离时出队 
		head++;
		f[x][y]=queue[head].fnub+(i-queue[head].cnt);
		ans=max(ans,f[x][y]);//更新答案 
		x=x+dx[dir],y=y+dy[dir];//移动 
	}
	return;
}
int main()
{
	scanf("%d%d%d%d%d",&n,&m,&x,&y,&k);
	getchar();
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			scanf("%c",&c);
			if(c=='x')
			gra[i][j]=1;
		}
		getchar();
	}
	memset(f,-0x3f,sizeof(f));
	f[x][y]=0;
	while(k--)
	{
		scanf("%d%d%d",&s,&e,&dir);//在每次一个时间段将更新图的每个点O(m*n) 
		int dis=e-s+1;
		switch(dir)//每个方向的update是能更新一行或者一列，for拓展到二维遍历整个图 
		{
			case 1:  for(int i=1;i<=m;i++)
			         {
			         	update(n,i,dir,dis);
		             }
					 break;
			case 2:  for(int i=1;i<=m;i++)
			         {
			         	update(1,i,dir,dis);
		             }
			         break;
			case 3:  for(int i=1;i<=n;i++)
			         {
			         	update(i,m,dir,dis);
		             }
		           	 break;
			case 4:  for(int i=1;i<=n;i++)
			         {
			         	update(i,1,dir,dis);
		             }
			         break;
		}
	}
	printf("%d\n",ans);
    return 0;
}

2021.8.2