bzoj1205: [HNOI2005]星际贸易

题目链接

bzoj1205: [HNOI2005]星际贸易

题解

辣鸡题面，毁我青春
辣鸡题面，毁我青
辣鸡题面，毁我
辣鸡题面，毁
第一问，背包dp
第二问
问题转化为在一个序列上经过好多点走到终点，
符合各项条件的情况下费用最小，其中有某些点必须经过。
g[u][fu]表示将要离开点u而还没有离开的时候有fu个燃料。
得到g[u][fu]=min(g[u][fu-1]+p[u],g[v][fu+2]) 、
复杂度n^3
单调队列维护fu相同时g[fu][v]的最小值
n^2

代码

#include<bits/stdc++.h> 
using namespace std;
inline int read() { 
	int x = 0,f = 1; 
	char c = getchar(); 
	while(c < '0' || c > '9'){if(c == '-')f = - 1; c = getchar();} 
 	while(c <= '9' && c >= '0')x = x * 10 + c - '0',c=  getchar();  
	return x * f; 
} 
const int maxn =  2007; 
int n,m,r,l0; //m max afford // r : 最大燃料量，//l0 不做维护飞的最大距离 
int a[maxn],l[maxn],p[maxn],c[maxn],f[maxn][maxn * 2]; 
int g[maxn][maxn * 2];  
bool ch[maxn][maxn]; 
int q[maxn ][maxn] ,h[maxn],t[maxn]; 
const int INF = 1e9; 
bool must[maxn]; 

int solve() { 
	memset(g,0x3f,sizeof g); 
	memset(t,-1,sizeof t); 
	q[r][t[r] = 0] = g[0][r] = 0 ; 
	for(int u = 1;u <= n;++ u) 
		for(int fu = 0;fu <= r;++ fu) { 
			if(fu > 0 && p[u]) g[u][fu] = min(g[u][fu],g[u][fu- 1] + p[u]); 
			int k = fu + 2; 
			if(k <= r) { 
				
				while(h[k] <= t[k] && l[u] - l[q[k][h[k]]] > l0) h[k] ++; 
				if(h[k] <= t[k]) g[u][fu] = std::min(g[u][fu],g[q[k][h[k]]][k] + c[u]); 
			} 
			while(h[fu] <= t[fu] && g[u][fu] < g[q[fu][t[fu]]][fu]) t[fu] --; 
			if(must[u]) t[fu] = h[fu] - 1; 
			q[fu][++ t[fu]] = u; 
		} 
	int ret = INF; 
	for(int i = 0;i <= r;++ i) ret = std::min(ret,g[n][i]); 
	return ret; 
} 
void get(int i,int j) { 
	if(! i ) return ; 
	if(ch[i][j]) must[i] = true , get(i - 1,j - a[i]); 
	else get(i - 1,j) ;  
} 
int main() { 
	n = read(),m = read(),r = read(),l0 = read(); 
	r = std::min(r,2 * n); 
	for(int b,i = 1;i <= n;++ i) { 
		a[i] = read(),b = read(),l[i] = read(),p[i] = read(),c[i] = read(); 
		for(int j = m;j >= 0;-- j) {  
			f[i][j] = f[i - 1][j]; 
			if(j >= a[i] && f[i - 1][j - a[i]] + b > f[i][j]) { 
				f[i][j] = f[i - 1][j - a[i]] + b; 
				ch[i][j] = 1; 
			} 
		} 
	}	
	get(n,m); 
	int ned = solve(); 
	if(ned >= INF) puts("Poor Coke!"); 
	else printf("%d %d
",f[n][m],f[n][m] - ned); 
	return 0; 
}