[bzoj 4380--POI2015]Myjnie

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

这道题考虑区间dp，f[l][r][k]表示从l到r的洗车店，最便宜的价格是k时，花的钱的max。价格可以离散化，dp方程很容易推吧。。。（那这道题就做完了）

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0' || ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9')x=x*10+ch-'0',ch=getchar();
    return x*f;
}
inline void write(int x)
{
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+'0');
}
struct node
{
    int l,r,d;
}a[4010];
bool cmp(node a,node b){return a.d<b.d;}
int s[4010],ex[4010],f[55][55][4010],ff[55][55][4010],g[55][55][4010],gg[55][55][4010],cnt[55][4010];
int ans[55];
inline void dfs(int l,int r,int i)
{
    if(l>r)return ;
    int p=ff[l][r][i];
    ans[p]=ex[i];
    dfs(l,p-1,gg[l][p-1][i]),dfs(p+1,r,gg[p+1][r][i]);
}
int main()
{
    //freopen("a.in","r",stdin);
    //freopen("a.out","w",stdout);
    int n=read(),m=read();
    for(int i=1;i<=m;i++)a[i].l=read(),a[i].r=read(),a[i].d=read();
    sort(a+1,a+m+1,cmp);
    int ul=0;
    for(int i=1;i<=m;i++)
    {
        if(a[i].d!=a[i-1].d)ul++,ex[ul]=a[i].d;
        s[i]=ul;
    }
    for(int d=1;d<=n;d++)
    {
        for(int l=1;l+d-1<=n;l++)
        {
            int r=l+d-1;
            for(int i=l;i<=r;i++)
            {
                for(int j=1;j<=ul;j++)cnt[i][j]=0;
            }
            for(int i=1;i<=m;i++)
            {
                if(a[i].l<l || a[i].r>r)continue;
                for(int j=a[i].l;j<=a[i].r;j++)cnt[j][s[i]]++;
            }
            for(int i=l;i<=r;i++)
            {
                for(int j=ul-1;j>=1;j--)cnt[i][j]+=cnt[i][j+1];
            }
            for(int i=ul;i>=1;i--)
            {
                for(int p=l;p<=r;p++)
                {
                    int wy=g[l][p-1][i]+g[p+1][r][i]+ex[i]*cnt[p][i];
                    if(f[l][r][i]<=wy)f[l][r][i]=wy,ff[l][r][i]=p;
                    g[l][r][i]=g[l][r][i+1];gg[l][r][i]=gg[l][r][i+1];
                    if(g[l][r][i]<=f[l][r][i])g[l][r][i]=f[l][r][i],gg[l][r][i]=i;
                }   
            }
        }
    }
    int p=1;
    for(int i=2;i<=ul;i++)if(f[1][n][p]<f[1][n][i])p=i;
    dfs(1,n,p);
    printf("%d\n",f[1][n][p]);
    for(int i=1;i<=n;i++)printf("%d ",ans[i]);
    return 0;
}