洛谷 P2254 [NOI2005]瑰丽华尔兹(单调栈优化DP)

题目描述

不妨认为舞厅是一个N行M列的矩阵,矩阵中的某些方格上堆放了一些家具,其他的则是空地。钢琴可以在空地上滑动,但不能撞上家具或滑出舞厅,否则会损坏钢琴和家具,引来难缠的船长。每个时刻,钢琴都会随着船体倾斜的方向向相邻的方格滑动一格,相邻的方格可以是向东、向西、向南或向北的。而艾米丽可以选择施魔法或不施魔法:如果不施魔法,则钢琴会滑动;如果施魔法,则钢琴会原地不动。

艾米丽是个天使,她知道每段时间的船体的倾斜情况。她想使钢琴在舞厅里滑行的路程尽量长,这样1900 会非常高兴,同时也有利于治疗托尼的晕船。但艾米丽还太小,不会算,所以希望你能帮助她。

输入格式

输入文件的第一行包含5个数N, M, x, y和K。N和M描述舞厅的大小,x和y为钢琴的初始位置;我们对船体倾斜情况是按时间的区间来描述的,且从1开始计算时间,比如“在[1, 3]时间里向东倾斜,[4, 5]时间里向北倾斜”,因此这里的K表示区间的数目。

以下N行,每行M个字符,描述舞厅里的家具。第i 行第j 列的字符若为‘ . ’,则表示该位置是空地;若为‘ x ’,则表示有家具。

以下K行,顺序描述K个时间区间,格式为:si ti di(1 ≤ i ≤ K)。表示在时间区间[si, ti]内,船体都是向di方向倾斜的。di为1, 2, 3, 4中的一个,依次表示北、南、西、东(分别对应矩阵中的上、下、左、右)。输入保证区间是连续的,即

s1 = 1 ti = si-1 + 1 (1 < i ≤ K)

tK = T

输出格式

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

输入输出样例

输入 #1
4 5 4 1 3
..xx.
.....
...x.
.....
1 3 4
4 5 1
6 7 3
输出 #1
6

说明/提示

钢琴的滑行路线:

钢琴在“×”位置上时天使使用一次魔法,因此滑动总长度为6。

【数据范围】

50%的数据中,1≤N, M≤200,T≤200;

100%的数据中,1≤N, M≤200,K≤200,T≤40000。


题意:题目说的很清楚,就是给你一张图有空'.',有障碍'x',钢琴起始位置在sx,sy,给你K个输入,表示钢琴在时间s到t可以向d方向滑动,但是你可施加魔法,是的钢琴在某个时间不滑动,问你钢琴最多可以滑动多少格;

题解:首先考虑对于时间t来dp: f[t][i][j]表示在第t时刻在第i行第j列所能获得的最长距离。 转移方程:f[t][i][j]=max(f[t-1][i][j],f[t][i*][j*]+1)(i*,j*为上一个合理的位置) 但对于100%TLE且MLE。 所以必须优化,首先把时间t换成区间k, 令f[k][i][j]表示在第k段滑行区间中在位置i,j所能获得最长距离 注意到在第k段时间内只能向某个方向最多走x步(x为区间长度),得到转移方程 f[k][i][j]=max(f[k-1][i][j],f[k][i*][j*]+dis(i,j,i*,j*))(i*,j*为上一个合理的位置) 这个做法的时间复杂度是O(kn^3),会超时,需要进一步优化 用单调队列优化掉内层的一个n,就可以做到O(kn^2),可以AC,本代码中还使用了滚动数组优化 用单调递减队列求最大值时,遇到障碍清空整个队列即可,另外队列比较时需要加上偏移量dis。

参考代码:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define pii pair<int,int>
#define pil pair<int,ll>
#define mkp make_pair
#define fi first
#define se second
const int INF=0x3f3f3f3f;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}
const int maxn=210;
pii q[maxn];
int n,m,sx,tx,K,ans,dp[maxn][maxn];
char G[maxn][maxn];
int dx[5]={0,-1,1,0,0},dy[5]={0,0,0,-1,1};
inline void work(int x,int y,int len,int d)
{
    int l=0,r=0;
    for(int i=1;x>=1&&x<=n&&y>=1&&y<=m;++i,x+=dx[d],y+=dy[d])
    {
        if(G[x][y]=='x') l=r=0;
        else
        {
            while(l<r&&q[r-1].fi+i-q[r-1].se<dp[x][y]) --r;
            q[r++]=mkp(dp[x][y],i);
            if(q[r-1].se-q[l].se>len) ++l;
            dp[x][y]=q[l].fi+i-q[l].se;
            ans=max(ans,dp[x][y]);
        }
    }
}
int main()
{
    n=read();m=read();sx=read();tx=read();K=read();
    for(int i=1;i<=n;++i) scanf("%s",G[i]+1);
    memset(dp,-INF,sizeof dp);
    dp[sx][tx]=0;
    for(int j=1;j<=K;++j)
    {
        int s,t,d,len;
        s=read();t=read();d=read();
        len=t-s+1;
        if(d==1) for(int i=1;i<=m;++i) work(n,i,len,d);
        if(d==2) for(int i=1;i<=m;++i) work(1,i,len,d);
        if(d==3) for(int i=1;i<=n;++i) work(i,m,len,d);
        if(d==4) for(int i=1;i<=n;++i) work(i,1,len,d);
    }
    printf("%d
",ans);

    return 0;
}
View Code
原文地址:https://www.cnblogs.com/csushl/p/11343802.html