DP matix

题面去内网找

f[i][j]表示到第i列，有j个1放到了右区间里。
这种神思路。。（我考试时想到是不是可以横向找。。但没想出转移方程）
那么如何转移呢？
设l[i]为到i列时结束了的左区间的个数，r[i]到i列开始的右区间个数
1.有j个1放在右区间，那么也就有i-j放在了左区间。又因为在这一列之前结束的左区间在之前已经放上，所以就剩了i-j-l[i-1]个1，可以放在在第i列结束的左区间里。（当然澄清一下，这些个1可以放在之前某一列里）那么对答案的贡献就要乘上A(i-j-l[i-1],l[i]-l[i-1])了。
2.到i列时不放新的1，f[i][j]=f[i-1][j];如果放了一个：f[i][j+1]+=f[i-1][j]*(r[i]-j)；解释一下，新放的那个1，能放在i列已开始的右区间里，但要去掉已经放了的j个区间。
如此转移即可，但注意，1转移一定要在2转移完成后再进行。

#pragma GCC optimize("O3")
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define N 3005
#define mod 998244353
#define ll long long
using namespace std;
int read()
{
    int sum=0,f=1;char x=getchar();
    while(x<'0'||x>'9'){if(x=='-')f=-1;x=getchar();}
    while(x>='0'&&x<='9'){sum=(sum<<1)+(sum<<3)+x-'0';x=getchar();}
    return sum*f;
}
int n,m;
ll ans,L[N],R[N],suml[N],sumr[N],a[N][N],f[N][N];
void init()
{
    for(ll i=1;i<=m;i++)
    {
        a[i][0]=1;
        for(ll j=1;j<=i;j++)
            a[i][j]=a[i][j-1]*(i-j+1)%mod;
    }
}
int main()
{
    n=read();m=read();
    for(int i=1;i<=n;i++)
    {
        L[i]=read(),R[i]=read();
        suml[L[i]]++;sumr[R[i]]++;
    }
    for(int i=1;i<=m;i++)suml[i]+=suml[i-1],sumr[i]+=sumr[i-1];
    init();f[1][0]=1;
    a[0][0]=1;
    for(int i=1;i<=m;i++)
        for(ll j=0;j<=n;j++)
        {
            int x1=i-j-suml[i-1],x2=suml[i]-suml[i-1];
            if(x1<x2)f[i][j]=0;
            else
            {
                f[i][j]=f[i][j]*a[i-j-suml[i-1]][suml[i]-suml[i-1]]%mod;
                f[i+1][j]+=f[i][j];f[i+1][j]%=mod;
                f[i+1][j+1]=(f[i+1][j+1]+f[i][j]*(sumr[i+1]-j)%mod)%mod;
            }
        }
    printf("%lld",f[m][n]);
}