matrix （dp）

10.4

思路：
神奇dp
f[i][j]表示处理到前i列（保证前i列合法：被完全包含的区间中已经放了1，没有完全包含的区间中有j个不放（有j列在i的右侧放1）（且没有被完全包含的左区间中一定不放1）
对于所有的左区间，我们只在扫到它的l时处理（加个1进去），而对于右区间我们可以有选择的处理

#include<cstdio>
#include<cstring>
#include<algorithm>
#define rep(i,a,b) for(i=a;i<=b;i++)
using namespace std;
const int N=3010;
const int maxn=3000;
const int mod=998244353;
int c[N][N],dp[N][N];
int mud[N];
int sl[N],sr[N],s[N];
int cl[N],cr[N];

int main(){
    freopen ("matrix.in", "r", stdin);
    freopen ("matrix.out", "w", stdout);
    int k,c1,c2;
    int l,r;
    int n,m;
    for(int i=0; i<=maxn; i++){
        c[i][0]=1;
        for(int j=1; j<=i; j++) 
            c[i][j] = (c[i-1][j] + c[i-1][j-1]) % mod;
    }//递推组合数 
    mud[0] = 1;
    for(int i=1; i<=maxn; i++)  mud[i] = 1ll * mud[i-1] * i % mod;//递推阶乘 
    scanf("%d%d", &n, &m);
    for(int i=1; i<=n; i++){
        scanf("%d%d", &l, &r);
        cl[l]++;
        cr[r]++;
    }
    for(int i=1; i<=m; i++) sl[i] = sl[i-1] + cl[i];
    for(int i=1; i<=m; i++) sr[i] = sr[i-1] + cr[i];
    for(int i=1; i<=m; i++) s[i] = s[i-1] + cl[i] + cr[i];
    dp[0][0] = 1;
    for(int i=1; i<=m; i++){
        for(int j=0; j<=sr[i-1]; j++){
            if( !dp[i-1][j] ) continue;
            //f[i][j]表示处理到前i列（保证前i列合法：
            //被完全包含的区间中已经放了1，没有完全包含的区间中有j个不放（有j列在i的右侧放1）（且没有被完全包含的左区间中一定不放1）
            //j等于没有被完全包含的左区间个数加上选择不放1的右区间个数 
            //并且保证每列最多一个1
            //对于所有的左区间，我们只在扫到它的l时处理（加个1进去），而对于右区间我们可以有选择的处理 
            k = i - 1 - s[i-1] + j;//i的左边还可以填多少个1 
            //i左边有i-1个位置，s[i-1]表示被包含的区间总数（其中有j个里面不放1，放在了i的后面）
            //所以左边已经放了s[i-1]-j个1 
            if(k < 0) continue;
            c1 = cl[i] <= k ? 1ll * c[k][cl[i]] * mud[cl[i]] % mod : 0;
            //不算此列，我们有k个空位，要放 cl[i]个1进去 （交换两个1的位置是不同的方案，所以是排列） 
            c2 = cl[i] <= k+1 ? 1ll * c[k+1][cl[i]] * mud[cl[i]] % mod : 0;//组合乘阶乘就是排列 
            //算上此列，我们有k+1个空位，要放 cl[i]个1进去 
            dp[i][j+cr[i]] = (dp[i][j+cr[i]] + 1ll * dp[i-1][j] * c2) % mod;
            //如果我们新增的右区间全部选择不放1（放在后面），那么我们留到之后处理的1就有j+cri个
            //方案数就等于上一个状态的方案数 * c2
            //当前这一列不处理新加入的区间，就是说i左边的k个空位加上i本身一共有k+1个空位，
            // 要处理掉cli个区间（跳出这个区间时一定要放一个1进去，只会跳出左区间），就是要加入cli个1
            if(j + cr[i]) dp[i][j+cr[i]-1] = (dp[i][j+cr[i]-1] + 1ll * dp[i-1][j] * c1 % mod * (j+cr[i])) % mod;
            //考虑在i这一位处理一个新加入的区间（在新加入的区间里选择一个放上1），然后就只剩下k个位置了
            //用这k个位置处理跳出的cli个区间，放cli个1到这k个空位中 
        }
    }
    printf("%d", dp[m][0]) % mod;//
    return 0;
}