hdu4374One hundred layer （DP+单调队列）

http://acm.hdu.edu.cn/showproblem.php?pid=4374

去年多校的题 今年才做 不知道这一年都干嘛去了。。

DP的思路很好想 dp[i][j] = max(dp[i-1][g]+sum[i][j]-sum[i][g-1],dp[i][j]) abs(g-j)<=t  不过复杂度是相当高的 所以呢 就出现了个单调队列 把它优化下

所谓的单调队列其实也就是一队列 始终保持着队头是最大的 若满足不了距离的条件 队头+1 队尾始终保持更新 让满足的了距离而且比队里的更大的元素进来

。。因为一初始化错了 　WA了十多次啊 注意：有负数 

#include <iostream>
#include<cstdio>
#include<cstring>
#include<stdlib.h>
#include<algorithm>
using namespace std;
#define N 105
#define M 10005
#define LL __int64
int a[N][M],que[M];
int dp[N][M];
int s1[N][M],s2[N][M];
int main()
{
    int i,j,n,m,x,t,str,end;
    while(cin>>n>>m>>x>>t)
    {
        memset(dp,128,sizeof(dp));
        memset(s1,0,sizeof(s1));
        memset(s2,0,sizeof(s2));
        for(i = 1; i <= n; i++)
            for(j = 1; j <= m ; j++)
            {
                scanf("%d",&a[i][j]);
                s1[i][j] = s1[i][j-1]+a[i][j];
            }
        for(i = 1 ; i <= n ;i++)
            for(j = m ; j >= 1; j--)
            s2[i][j] = s2[i][j+1]+a[i][j];
        for(i = x ; i>=1&&i>=x-t ; i--)
        dp[1][i] = s1[1][x]-s1[1][i-1];
        for(i = x ; i <= m&&i<=x+t ; i++)
        dp[1][i] = s2[1][x]-s2[1][i+1];
        for(i = 2; i <= n ; i++)
        {
            str = 0;end=0;
            for(j = 1;j <= m ; j++)
            {
                while(str<end&&dp[i-1][j]-s1[i][j-1]>dp[i-1][que[end-1]]-s1[i][que[end-1]-1])
                end--;
                que[end++] = j;
                dp[i][j] = max(dp[i][j],dp[i-1][que[str]]-s1[i][que[str]-1]+s1[i][j]);
                if(j-que[str]>=t)
                str++;
            }
            str = 0;end=0;
            for(j = m ; j >= 1; j--)
            {
                while(str<end&&dp[i-1][j]-s2[i][j+1]>dp[i-1][que[end-1]]-s2[i][que[end-1]+1])
                end--;
                que[end++] = j;
                dp[i][j] = max(dp[i][j],dp[i-1][que[str]]-s2[i][que[str]+1]+s2[i][j]);
                if(que[str]-j>=t)
                str++;
            }
        }
        int maxz=dp[n][1];
        for(i = 1; i <= m ; i++)
        maxz = max(maxz,dp[n][i]);
        cout<<maxz<<endl;
    }
    return 0;
}