JZOJ5192. 【NOI2017模拟7.2】容器

题目大意

在一个长度为n的序列里面，一开始全部为0，
要选择k个非空的区间，对这个区间里面的每个位置加1，
最终序列上面的每一个位置都不能大于t，
求有多少种合法的方案。
注意：k个区间是有序的，即[1,1] , [2,2]与[2,2] , [1,1]算两种不同的方案。

题解

其中t=1的部分分给出了提示，
在这个部分分里面每个位置最多只能被一个区间所覆盖。
对于区间的顺序问题，
可以先按照左端点从小到大排序，然后再最后乘上k的阶乘就可以了。
那么设$f_{i,j}$表示当前在第i个位置，已经放了j个区间，而且强制i这个位置必须被最后一个区间覆盖。
根据这个思路，将它扩展到t的其他情况。
这样就会要多出一维，
设$f_{i,j,k}$表示当前在第i个位置，已经放了j个区间，其中有k个区间没有右端点，即可以延续。
转移的时候，就枚举一个inc和dec，表示在i这个位置，新增左端点的个数为inc，在这个位置的右端点有dec个，
那么就可以从$f_{i,j,k}$转移到$f_{i+1,j+inc,k+inc-dec}$。
现在就要考虑转移系数，首先i这个位置作为右端点，就可以从k个里面任意选择，$C ^{ dec}_k$，然后，新加入inc个区间可以插入到j个区间里面的任意位置，$C^{inc}_{inc+j}$
最后的答案就很简单了。

code

#include <queue>
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string.h>
#include <cmath>
#include <math.h>
#include <time.h>
#define ll long long
#define N 100003
#define M 43
#define db double
#define P putchar
#define G getchar
using namespace std;
char ch;
void read(int &n)
{
    n=0;
    ch=G();
    while((ch<'0' || ch>'9') && ch!='-')ch=G();
    ll w=1;
    if(ch=='-')w=-1,ch=G();
    while('0'<=ch && ch<='9')n=(n<<3)+(n<<1)+ch-'0',ch=G();
    n*=w;
}

int max(int a,int b){return a>b?a:b;}
int min(int a,int b){return a<b?a:b;}
ll abs(ll x){return x<0?-x:x;}
ll sqr(ll x){return x*x;}
void write(ll x){if(x>9) write(x/10);P(x%10+'0');}

const int mo=1011110011;
long long C[M][M],f[M][M][M];
int n,m,t,ans;

int main()
{
    for(register int i=1;i<M;i++)
    {
        C[i][1]=C[i][i]=1;
        for(register int j=2;j<i;j++)
            C[i][j]=(C[i-1][j]+C[i-1][j-1])%mo;
    }
    freopen("container.in","r",stdin);
    freopen("container.out","w",stdout);
    read(n);read(m);read(t);
    f[0][0][0]=1;
    for(register int i=0;i<n;i++)
        for(register int j=0;j<=m;j++)
            for(register int k=0;k<=min(t,j);k++)
                if(f[i][j][k])
                for(register int inc=0;inc<=min(t-k,m-j);inc++)
                    for(register int dec=max(0,k+inc-t);dec<=k+inc;dec++)
                        f[i+1][j+inc][k+inc-dec]=(f[i+1][j+inc][k+inc-dec]+
                            f[i][j][k]*C[j+inc+1][inc+1]%mo*C[k+inc+1][dec+1]%mo)%mo;
    printf("%lld",f[n][m][0]);
    return 0;
}