HDU 4374 - One hundred layer

 

题目地址： http://acm.hdu.edu.cn/showproblem.php?pid=4374

 

2012年 第8场多校 ， 1005 。

 

方法： 单调队列。

 

先处理第一行，所有起点能到达的地方。

 

接下去每一行，先用优先队列往左扫一遍，再用优先队列往右扫一遍。取最佳值。

 

Trick ：

     1 .  处理时，每走一层，需要往左和往右延伸T的长度，直到往左到1、往右到m。 因为每一层最多朝某一方向走T单位。

     2 .  因为数据有负数，所以比较时设定的最小值需要足够小。

 

 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<map>
#include<vector>

#define ll (v<<1)
#define rr (v<<1|1)
#define tmid (l+r>>1)

using namespace std;

const int inf = 0x7fffffff;
const __int64 inf64 = (__int64)inf*inf;

__int64 mp[102][10010];

__int64 mx[10010],res[10010];
int rec[10010];
int ss,st;

int main(){
	int n,m,x,t,i,j,k,e,d,ex,dx;
	__int64 tmp,sav;
	while(~scanf("%d%d%d%d",&n,&m,&x,&t)){
		memset(res,0,sizeof(res));
		for(i=1;i<=n;i++)
			for(j=1;j<=m;j++)
				scanf("%I64d",&mp[i][j]);
		for(i=x-1,k=0;i>0 && k<t;i--,k++)
			mp[1][i]+=mp[1][i+1];
		for(i=x-t-1;i>0;i--)
			mp[1][i]=0;
		for(i=x+1,k=0;i<=m && k<t;i++,k++)
			mp[1][i]+=mp[1][i-1];
		for(i=x+t+1;i<=m;i++)
			mp[1][i]=0;
		e=d=x;
		for(i=2;i<=n;i++){
			e=e-t>0?e-t:1;
			d=d+t<m?d+t:m;
			sav=0;
			mx[ss=st=0]=mp[i-1][e];
			rec[0]=e;
			for(j=e;j<=d+t && j<=m;j++){ // 单调队列从左往右扫
				while(rec[ss]+t<j) ss++;
				if(j<=d){
					for(k=ss;k<=st;k++)
						if(mp[i-1][j]>=mx[k]+sav) break;
						else mx[k]+=sav;
					mx[st=k]=mp[i-1][j];
					rec[k]=j;
					sav=mp[i][j];
				}
				else sav+=mp[i][j];
				res[j]=mx[ss]+sav;
			}
			sav=0;
			mx[ss=st=0]=mp[i-1][d];
			rec[0]=d;
			for(j=d;j>=e-t && j>0;j--){  // 单调队列从右往左扫
				while(rec[ss]-t>j) ss++;
				if(j>=e){
					for(k=ss;k<=st;k++)
						if(mp[i-1][j]>=mx[k]+sav) break;
						else mx[k]+=sav;
					mx[st=k]=mp[i-1][j];
					rec[k]=j;
					sav=mp[i][j];
				}
				else sav+=mp[i][j];
				res[j]=max(mx[ss]+sav,res[j]);
			}
			for(j=1;j<=m;j++)
				mp[i][j]=res[j];
		}
		__int64 s=-inf64;
		for(int i=1;i<=m;i++)
			if(s<mp[n][i]) s=mp[n][i];
		printf("%I64d\n",s);
	}
	return 0;
}