链接
题面
有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。
若有多组最优解，输出任意一组。
题解
这道题……是一个正常的（？）区间DP。
先将所有c[i]离散化。
f[i][j][k]表示区间[i, j]内的店、最低价格为k，能取到的最大价值。这个最大价值包括所有“完全包含在[i, j]”内的顾客提供的价值。
g[i][j][k]表示区间[i, j]内的店、最低价格>=k，能取到的最大价值。
转移的时候，枚举最低价格k所在的位置p，采取类似分治的思想，设完全包含在[i, j]内、横跨位置p、且c[i] >= k的顾客有cnt[p][k]个，那么有：
\[f[i][j][k] = max_{p = i}^{j} g[i][p - 1][k] + g[p + 1][j][k] + k * cnt[p][k]\]
这个cnt数组在每次枚举[i, j]的时候都预处理一下。
因为要输出方案，所以额外记录一下最优解都是在哪里取到的，最后dfs求出方案。
总复杂度是\(O(n^3m)\)。
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
#define enter putchar('\n')
#define space putchar(' ')
template <class T>
void read(T &x){
char c;
bool op = 0;
while(c = getchar(), c < '0' || c > '9')
if(c == '-') op = 1;
x = c - '0';
while(c = getchar(), c >= '0' && c <= '9')
x = x * 10 + c - '0';
if(op) x = -x;
}
template <class T>
void write(T x){
if(x < 0) putchar('-'), x = -x;
if(x >= 10) write(x / 10);
putchar('0' + x % 10);
}
const int N = 52, M = 4002;
int n, m, a[M], b[M], c[M], cnt[N][M], ans[N], lst[M], idx, f[N][N][M], g[N][N][M], gk[N][N][M];
short fp[N][N][M];
void dfs(int l, int r, int k){
int val = gk[l][r][k], pos = fp[l][r][val];
ans[pos] = lst[val];
if(pos > l) dfs(l, pos - 1, val);
if(pos < r) dfs(pos + 1, r, val);
}
int main(){
read(n), read(m);
for(int i = 1; i <= m; i++)
read(a[i]), read(b[i]), read(c[i]), lst[i] = c[i];
sort(lst + 1, lst + m + 1);
idx = unique(lst + 1, lst + m + 1) - lst - 1;
for(int i = 1; i <= m; i++)
c[i] = lower_bound(lst + 1, lst + idx + 1, c[i]) - lst;
for(int i = n; i; i--)
for(int j = i; j <= n; j++){
for(int k = 1; k <= idx; k++)
for(int h = i; h <= j; h++)
cnt[h][k] = 0;
for(int k = 1; k <= m; k++)
for(int h = i; h <= j; h++)
if(a[k] >= i && a[k] <= h && b[k] >= h && b[k] <= j)
cnt[h][c[k]]++;
for(int h = i; h <= j; h++)
for(int k = idx; k; k--)
cnt[h][k] += cnt[h][k + 1];
for(int k = 1; k <= idx; k++)
for(int p = i; p <= j; p++)
if(g[i][p - 1][k] + g[p + 1][j][k] + lst[k] * cnt[p][k] >= f[i][j][k])
fp[i][j][k] = p, f[i][j][k] = g[i][p - 1][k] + g[p + 1][j][k] + lst[k] * cnt[p][k];
for(int k = idx; k; k--){
if(g[i][j][k + 1] > f[i][j][k])
gk[i][j][k] = gk[i][j][k + 1], g[i][j][k] = g[i][j][k + 1];
else
gk[i][j][k] = k, g[i][j][k] = f[i][j][k];
}
}
write(g[1][n][1]), enter;
dfs(1, n, 1);
for(int i = 1; i <= n; i++)
write(ans[i]), i == n ? enter: space;
return 0;
}