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\n", 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;

}