bzoj [POI2015]Myjnie

[POI2015]Myjnie

Time Limit: 40 Sec Memory Limit: 256 MBSec Special Judge

Description

有n家洗车店从左往右排成一排，每家店都有一个正整数价格p[i]。
有m个人要来消费，第i个人会驶过第a[i]个开始一直到第b[i]个洗车店，且会选择这些店中最便宜的一个进行一次消费。但是如果这个最便宜的价格大于c[i]，那么这个人就不洗车了。
请给每家店指定一个价格，使得所有人花的钱的总和最大。

Input

第一行包含两个正整数n,m(1<=n<=50，1<=m<=4000)。
接下来m行，每行包含三个正整数a[i],b[i],ci

Output

第一行输出一个正整数，即消费总额的最大值。
第二行输出n个正整数，依次表示每家洗车店的价格p[i]，要求1<=p[i]<=500000。
若有多组最优解，输出任意一组。

Sample Input

7 5

1 4 7

3 7 13

5 6 20

6 7 1

1 2 5

Sample Output

43

5 5 13 13 20 20 13

这个区间dp有点东西。。。 离散化一下 f[l][r][t] 表示完全属于区间的人中最小值刚好是t的最优贡献。 g[l][r][t] 表示完全属于区间的人中最小值大于等于t的最优贡献。 然后你的转移就成了枚举最小值是哪个店。 f[l][r][t]=max(g[l][i-1][t]+g[i+1][r][t]+cost(i)) 然而这个cost(i)有点尴尬。。。你可能还有有个预处理h[i][j]表示在处理当前这个dp区间的情况下，包含i这个位置的人且在最小值为j的情况下能给钱的人有多少个。。。 写吧~
1 mol 乱七八糟的错误。。。。。调一万年啊。。。。
```c++

include<bits/stdc++.h>

using namespace std;
const int maxm = 4004;
struct lpl{
int a, b, c;
}ini[maxm];
struct lpd{
int g, f, pos;
}dp[55][55][maxm];
int n, m, k, ans[55], pl[maxm], nam[maxm], mark[500005], h[55][maxm];

inline bool cmp(lpl A, lpl B) {return A.b == B.b ? A.a < B.a : A.b < B.b;}

inline void putit()
{
scanf("%d%d", &n, &m);
for(int i = 1; i <= m; ++i) scanf("%d%d%d", &ini[i].a, &ini[i].b, &ini[i].c), pl[i] = ini[i].c;
sort(ini + 1, ini + m + 1, cmp);
sort(pl + 1, pl + m + 1); k = 0;
for(int i = 1; i <= m; ++i)
if(pl[i] != pl[i - 1]){
k++; nam[k] = pl[i]; mark[pl[i]] = k;
}
for(int i = 1; i <= m; ++i) ini[i].c = mark[ini[i].c];
}

inline void workk()
{
for(int l = n; l >= 1; --l)
for(int r = l; r <= n; ++r){
for(int i = l; i <= r; ++i)
for(int j = 1; j <= k; ++j)
h[i][j] = 0;
for(int i = 1; i <= m && ini[i].b <= r; ++i){
if(ini[i].a < l) continue;
for(int j = ini[i].a; j <= ini[i].b; ++j) h[j][ini[i].c]++;
}
for(int i = l; i <= r; ++i)
for(int j = k - 1; j >= 1; --j)
h[i][j] += h[i][j + 1];
for(int i = l; i <= r; ++i)
for(int j = k; j >= 1; --j){
int tmp = dp[l][i - 1][j].g + dp[i + 1][r][j].g + h[i][j] * nam[j];
if(dp[l][r][j].f < tmp){
dp[l][r][j].f = tmp; dp[l][r][j].pos = i;
}
dp[l][r][j].g = max(dp[l][r][j].f, dp[l][r][j + 1].g);
}
}
}

void dfs(int l, int r, int t)
{
if(l > r) return;
if(dp[l][r][t].g == 0){
for(int i = l; i <= r; ++i) ans[i] = nam[t];
return;
}
for(int i = t; i <= k; ++i){
if(dp[l][r][t].g == dp[l][r][i].f){
int now = dp[l][r][i].pos;
ans[now] = nam[i];
dfs(l, now - 1, i); dfs(now + 1, r, i);
break;
}
}
}

inline void print()
{
printf("%d ", dp[1][n][1].g);
dfs(1, n, 1);
for(int i = 1; i <= n; ++i) printf("%d ", ans[i]);
}

int main()
{
putit();
workk();
print();
return 0;
}