BZOJ1499: [NOI2005]瑰丽华尔兹
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
..xx.
.....
...x.
.....
1 3 4
4 5 1
6 7 3
Sample Output
本蒟蒻的第一道单调队列优化 DP 。
感觉啥都不会。。。
首先最显然的DP: dp[ i ][ x ][ y ] = max{ dp[ i-1 ][ x' ][ y' ] } + 1 。
(反正我是想不出来。。。)
ii 表示的是时间点, (x,y) 表示位置,而 (x',y') 表示上一个合法位置。
But !铁定 TLE。(NOI的题怎么可能这么水)
再看题目——时间段!
那我们换种表示试试:dp[ i ][ x ][ y ] = max{ dp[ i-1 ][ x' ][ y' ] } + dist( ( x , y ) ,( x' , y' ) )
i 表示的是第 i 个时间段结束后, (x,y) 这个位置最长的滑行距离。
注意 (x,y) 与 (x',y') 必定是在同一列或同一行上的,但不一定相邻。
但是还是超时。
那么回头再看状态转移方程,发现有要求 max ,而且求可能拓展的状态有线性关系(在同一列或同一行上)。
主角登场——
单调队列
首先按照读入的每个时间段的方向,穷举每个方向的初始点。
在单调队列里记录所在列或行某一个节点减去它到这一列或行初始位置的距离。
对于每个点,用步数(就是初始点到当前点的距离)加上队首。
煮个栗子:目前是第 22 短时间,向下,也就是向南。
dp[1][2][2] 已经求出来了,是 3 ,并且在队首。
那么队首记录的就是 3-2=1
那么 dp[2][3][2]=1+3=4
假设 dp[1][2][2] 还是队首;
那么 dp[2][4][2]=1+4=5
于是代码就很愉快的敲出来了。
附代码:
#include<iostream> #include<algorithm> #include<cstdio> #include<cstring> #define MAXN 210 #define MAXM 100010 #define MAX 999999999 using namespace std; const int fx[5]={0,-1,1,0,0},fy[5]={0,0,0,-1,1}; int n,m,x,y,k,ans=0,head,tail; int val[MAXN][MAXN],que[MAXM],pos[MAXM],dp[MAXN][MAXN][MAXN]; inline int read(){ int date=0,w=1;char c=0; while(c<'0'||c>'9'){if(c=='-')w=-1;c=getchar();} while(c>='0'&&c<='9'){date=date*10+c-'0';c=getchar();} return date*w; } inline void push(int now,int v,int x,int y){ if(v==-MAX)return; while(head<=tail&&v-now>=que[tail])tail--; que[++tail]=v-now; pos[tail]=now; } void solve(int p,int x,int y,int d,int t){ if(t<1)return; int now=1; head=1;tail=0; while(1<=x&&x<=n&&1<=y&&y<=m){ if(!val[x][y]){head=1;tail=0;} else push(now,dp[p-1][x][y],x,y); while(head<=tail&&now-pos[head]>t)head++; if(head<=tail)dp[p][x][y]=que[head]+now; else dp[p][x][y]=-MAX; ans=max(ans,dp[p][x][y]); x+=fx[d];y+=fy[d]; now++; } } void work(){ int s,t,d; for(int i=1;i<=k;i++){ s=read();t=read();d=read(); switch(d){ case 1:for(int j=1;j<=m;j++)solve(i,n,j,d,t-s+1);break; case 2:for(int j=1;j<=m;j++)solve(i,1,j,d,t-s+1);break; case 3:for(int j=1;j<=n;j++)solve(i,j,m,d,t-s+1);break; case 4:for(int j=1;j<=n;j++)solve(i,j,1,d,t-s+1);break; } } printf("%d\n",ans); } void init(){ char ch[MAXN]; n=read();m=read();x=read();y=read();k=read(); memset(dp,127,sizeof(dp)); for(int i=1;i<=n;i++){ scanf("%s",ch+1); for(int j=1;j<=m;j++){ val[i][j]=(ch[j]=='x'?0:1); dp[0][i][j]=-MAX; } } dp[0][x][y]=0; } int main(){ init(); work(); return 0; }