BZOJ4584: [Apio2016]赛艇

原题可以转化为将n个人按顺序放在若干个区间中，可以不放

将区间离散化到4n个区间（一直不懂泥萌是怎么离散化到2n个区间的我可能只能弄到4n个qwq
f[i][j]表示前i个区间，分布j个人的方案数
为了去重，转移时枚举第i个区间有k个人，强制第k个人一定在这个区间，其他人可以不选，那么方案是$C_{L+k-1}^k$，可重组合，选重复的代表不选
dp时可以原地倒着dp省掉一维也不需要滚动

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;

const int maxn = 2100;
const int mod = 1e9+7;
int add(int &a,const int &b){a+=b;if(a>=mod) a-=mod;}

int s[maxn],inv[maxn];
int pw(int x,int k)
{
    int re=1;
    for(;k;k>>=1,x=(ll)x*x%mod) if(k&1)
        re=(ll)re*x%mod;
    return re;
}
void pre()
{
    s[0]=1; for(int i=1;i<maxn;i++) s[i]=(ll)s[i-1]*i%mod;
    inv[maxn-1]=pw(s[maxn-1],mod-2);
    for(int i=maxn-2;i>=0;i--) inv[i]=(ll)inv[i+1]*(i+1)%mod;
    for(int i=1;i<maxn;i++) inv[i]=(ll)inv[i]*s[i-1]%mod;
}

int n,m;
struct inl{int p,i;}t[maxn]; int tp;
inline bool cmp(const inl x,const inl y){return x.p==y.p?x.i<y.i:x.p<y.p;}
struct node{int l,r;}a[maxn],c[maxn]; int cn;
int ci[maxn][maxn];
int f[maxn];

int main()
{
    pre();

    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d%d",&a[i].l,&a[i].r);
        t[++tp]=(inl){a[i].l,0},t[++tp]=(inl){a[i].r,1};
    }
    sort(t+1,t+tp+1,cmp); t[0].p=t[1].p-1;
    for(int i=1;i<=tp;i++) if(t[i].p!=t[i-1].p)
    {
        bool flag=false;
        if(t[i-1].p+1<t[i].p) 
        {
            if(t[i-1].i==0) c[cn].r=t[i].p-1;
            else if(t[i].i==1) flag=true,c[++cn]=(node){t[i-1].p+1,t[i].p};
            else c[++cn]=(node){t[i-1].p+1,t[i].p-1};
        }
        if(!flag) c[++cn]=(node){t[i].p,t[i].p};
    }
    for(int i=1;i<=cn;i++)
    {
        int l=c[i].r-c[i].l+1,u=2*n;
        ci[i][0]=1;
        for(int j=1;j<=u;j++)
            ci[i][j]=(ll)ci[i][j-1]*inv[j]%mod*(l+j-1)%mod;
    }
    for(int i=1;i<=n;i++)
    {
        int j=1;for(;c[j].r<a[i].l;j++); a[i].l=j;
        j=cn;for(;a[i].r<c[j].l;j--); a[i].r=j;
    }

    f[0]=1;
    for(int i=1;i<=cn;i++)
    {
        for(int j=n;j>=1;j--) if(a[j].l<=i&&i<=a[j].r)
        {
            int nowi=0;
            for(int k=1;k<=j;k++)
            {
                if(a[j-k+1].l<=i&&i<=a[j-k+1].r) nowi++;
                add(f[j],(ll)f[j-k]*ci[i][nowi]%mod);
            }
        }
    }
    int re=0;
    for(int i=1;i<=n;i++) add(re,f[i]);
    printf("%d\n",re);

    return 0;
}