【模拟8.01】matrix（DP杂题，思维题）

很神的题，感谢lnc大佬的指点。

先设1-LL[i]统称左区间，RR[i]-m为右区间

用L[i]统计从1-i列，出现的左区间端点的前缀和，R[i]是右区间....

f[i][j]中j表示当前在第i列，右区间的左端点（RR[i]）到i存在的1的个数，总体表示当前方案数。

所以，我们分几种情况

两种是直接转移的

     f[i+1][j]=f[i][j]表示i向右移动但j不变方案转移

     f[i+1][j+1]=f[i][j]*(r[i+1]-j)表明又选了一个，当前新的1可以从r[i+1]（即右区间的个数）-j转移

一种是用排列.......

     f[i][j]=f[i][j]*A(i-j-L[i-1],L[i]-L[i-1])

     表示当前与上步操作中多出的区间是L[i]-L[i-1],这是必须填上1的

     i-j-L[i-1]表示最多放的，又因为顺序不定，所以是排列......

#include<iostream>
#include<cstdio>
#include<string>
#include<algorithm>
#include<cmath>
#include<vector>
#include<map>
#include<cstring>
#define int long long
#define MAXN 3010
#define mod 998244353
using namespace std;
int jie[MAXN],ni[MAXN],ni_c[MAXN];
int f[MAXN][MAXN];
int A(int x,int y)
{
    if(y==0)return 1;
   // printf("x=%lld x-y=%lld %lld %lld
",x,x-y,jie[x],ni_c[x-y]);
    return jie[x]*ni_c[x-y]%mod;
}
int n,m;
int l[MAXN],r[MAXN];
signed main()
{
    scanf("%lld%lld",&n,&m);
    jie[0]=1; ni[0]=1; ni_c[0]=1;
    jie[1]=1; ni[1]=1; ni_c[1]=1;
    for(int i=2;i<=m;++i)
    {
        jie[i]=(jie[i-1]*i)%mod;
        //printf("jie[%lld]=%lld
",i,jie[i]);
        ni[i]=(mod-mod/i)*ni[mod%i]%mod;
        ni_c[i]=(ni_c[i-1]*ni[i])%mod;
    }
    for(int i=1;i<=n;++i)
    {
        int x,y;
        scanf("%lld%lld",&x,&y);
        l[x]++;r[y]++;
    }
    for(int i=1;i<=m;++i)
    {
        l[i]+=l[i-1];
        r[i]+=r[i-1];
    }
    f[1][0]=1;
    for(int i=1;i<=m;++i)
    {
        for(int j=0;j<=r[i];++j)
        {
             if(l[i]-l[i-1]>i-j-l[i-1])break;
             f[i][j]=(f[i][j]*A(i-j-l[i-1],l[i]-l[i-1]))%mod;
             f[i+1][j]+=f[i][j];
             f[i+1][j+1]+=f[i][j]*(r[i+1]-j)%mod;
             //printf("f[%lld][%lld]=%lld
",i,j,f[i][j]);
        }
    }
    printf("%lld
",f[m][n]);
}