1. 程式人生 > >Fence[單調佇列+DP]

Fence[單調佇列+DP]

傳送門

f[i][j] 表示第i個工匠做到j的最大價值

f[i][j]=max(f[i-1][j],f[i][j-1]) 

f[i][j]=max(f[i-1][x]+(j-x)*p[i])(j-l[i]<=x<=s[i]-1)

顯然j應滿足j>=s[i]

我們將f[i][j]變形

f[i][j]=max(f[i-1][x]-x*p[i])+j*p[i] 

單調佇列維護f[i-1][x]-x*p[i]

每次先將f[i-1][...]合法的(s[i]-l[i]<=x<=s[i]-1)中的放入單調佇列 , 每次更新前除去不合法狀態(x<=j-l[i])


#include<cstdio>
#include<algorithm>
#define K 150
#define N 16050
#define val(i,j) (f[i-1][j] - j * x[i].p)
using namespace std;
int f[K][N],n,k,q[N];
struct Node{int l,p,s;}x[K];
bool cmp(Node a,Node b){return a.s<b.s;}
int main(){
	scanf("%d%d",&n,&k);
	for(int i=1;i<=k;i++){
		scanf("%d%d%d",&x[i].l,&x[i].p,&x[i].s);
	}
	sort(x+1,x+k+1,cmp);
	for(int i=1;i<=k;i++){
		int l=1,r=0;
		for(int j=max(0,x[i].s-x[i].l);j<=x[i].s-1;j++){
			while(l<=r && val(i,q[r]) <= val(i,j)) r--;
			q[++r] = j;
		}
		for(int j=1;j<=n;j++){
			f[i][j]=max(f[i-1][j] , f[i][j-1]);
			if(j>=x[i].s){
				while(l<=r && q[l] < j-x[i].l) l++;
				if(l<=r) f[i][j] = max(f[i][j] , val(i,q[l]) + j*x[i].p);
			}
		}
	}printf("%d",f[k][n]); return 0;
}