BZOJ2436: [Noi2011]Noi嘉年华

因为时间的具体数值对安排活动没什么影响，可以先离散化时间
需要知道第i个活动$Li\backsim Ri$必须举办，活动较少的嘉年华活动数量的最大值，可以转化为$(i=1\backsim Li)\backsim(j=Ri\backsim T)$的活动必须在同一个嘉年华举办的答案
设$f[i][j]$表示$i\backsim j$的活动在同一个嘉年华举办的答案，需要几个辅助的状态
设$pre[i][j]$表示时间$1~i$，有$j$个活动在A举行，B活动数的最大值，
$suf[i][j]$表示后缀，状态定义类似$pre$
$num[i][j]$表示$i\backsim j$的活动数量

$pre[i][j]=max(pre[k][j]+num[k+1][i],pre[k][j-num[k+1][i]])$
$suf$ 类似

于是得

$$f[i][j]=max(min(k+l,pre[i-1][k]+num[i][j]+suf[j+1][l]))$$
直接暴力转移是 $O(n^4)$的
可以发现，当 $k$增加时，$l$单调不升
可以优化到 $O(n^3)$

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

inline void swap(int &x,int &y){x^=y;y^=x;x^=y;}
inline void up(int &x,const int &y){if(x<y)x=y;}
const int maxn = 410;

struct node
{
    int s,t;
}a[maxn];
int c[maxn];
int n,m,K;

map<int,int>lsh;

void pr()
{
    sort(c+1,c+m+1);
    c[0]=c[1]-1;
    K=0;
    for(int i=1;i<=m;i++) if(c[i]!=c[i-1])
        lsh[c[i]]=++K;

    for(int i=1;i<=n;i++) a[i].s=lsh[a[i].s],a[i].t=lsh[a[i].t];
}

int num[maxn][maxn],pre[maxn][maxn],suf[maxn][maxn];
int f[maxn][maxn];

inline int clac(int i,int j,int x,int y)
{
    return min(x+y+num[i][j],pre[i][x]+suf[j][y]);
}

int main()
{
    scanf("%d",&n); m=0;
    for(int i=1;i<=n;i++) 
    {
        scanf("%d%d",&a[i].s,&a[i].t);
        a[i].t+=a[i].s;
        c[++m]=a[i].s,c[++m]=a[i].t;
    }
    pr();
    for(int i=1;i<=n;i++)
    {
        for(int l=a[i].s;l>=1;l--) for(int r=a[i].t;r<=K;r++)
            num[l][r]++;
    }

    swap(n,K);
    memset(pre,-1,sizeof pre);
    pre[1][0]=0;
    for(int r=2;r<=n;r++) for(int k=num[1][r];k>=0;k--)
    {
        for(int l=1;l<=r;l++)
        {
            int x=num[l][r];
            if(k>=x) up(pre[r][k],pre[l][k-x]);
            if(pre[l][k]!=-1) up(pre[r][k],pre[l][k]+x);
        }
        up(pre[r][k],pre[r][k+1]);
    }
    memset(suf,-1,sizeof suf);
    suf[n][0]=0;
    for(int l=n-1;l>=1;l--) for(int k=num[l][n];k>=0;k--)
    {
        for(int r=l;r<=n;r++)
        {
            int x=num[l][r];
            if(k>=x) up(suf[l][k],suf[r][k-x]);
            if(suf[r][k]!=-1) up(suf[l][k],suf[r][k]+x);
        }
        up(suf[l][k],suf[l][k+1]);
    }

    int mx=0;
    for(int i=1;i<n;i++) for(int j=i+1;j<=n;j++)
    {
        int cc=num[i][j];
        int y=num[j][n];
        int t1,t2;
        for(int x=0;x<=num[1][i];x++)
        {
            t1=clac(i,j,x,y);
            while(y&&(t2=clac(i,j,x,y-1))>=t1) 
                y--,t1=t2;
            up(f[i][j],t1);
        }
        up(mx,f[i][j]);
    }

    printf("%d\n",mx);
    for(int i=1;i<=K;i++)
    {
        int re=0;
        for(int l=a[i].s;l>=1;l--) for(int r=a[i].t;r<=n;r++)
            up(re,f[l][r]);
        printf("%d\n",re);
    }

    return 0;
}