2020暑期牛客多校训练营第七场（G）Topo Counting

Topo Counting

原题请看这里

题目描述：

给定一种有向无环图———的晒肉架图$(DRG)$，由唯一参数$N$控制。$DRG$包含$N$组顶点。第$i$组$V^i$包含$2N$个顶点：$V_1^i，V_2^i，\cdots ，V_{2N}^i$
在$DRG$中有两种类型的边：组内边（每组内的边）和内部-组边（组之间的边）。
组内边缘：对于第$i$组，存在以下组内边缘：
$(V_j^i，V_{j+N}^i)$，对于所有整数$j$使得$1\le j\le N$
$(V_j^i，V_{j+1}^i)$，对于所有整数$j$使得$1\le j\le N-1$或$N+1\le j\le 2N-1$
组间边缘：组边缘存在：
$(V_{i+N}^1，V_1^{i+1})$，对于所有整数$i$使得$1\le i\le N-1$；
$(V_i^1，V_{1+N}^i)$，对于所有整数$i$使得$2\le i\le N$
现在我们想知道以$N$为参数的$DRG$的拓扑序数。
有向图$G=(V，E)$的拓扑顺序是$V(G)$的所有顶点的排列$v_{p_1}，v_{p_2}，\cdots ，v_{p_{|V(G)|}})$使得对于所有$i<j$，$(v_{p_j}，v_{p_i})\notin E(G)$
为了避免计算巨大的整数，请对答案取模$M$。

输入描述：

输入仅包含两个整数$N，M(1\le N\le 3000，N\ast N\ast 2<M\le 2^{30})$，并且$M$是质数。

输出描述：

输出一个整数，表示答案。

样例：

样例输入1：

2 1073741789

样例输出1：

31

样例输入2：

3 1073741789

样例输出2：

7954100

思路：

首先我们画一张图$(n=4)$：（爬ppt（官方）的图）

从中我们可以发现：所有的子图都连向了第一号子图，仔细观察就会发现这个图就像一个晒肉架子… $n$唯一确定了这张图。
看到这张图，我们忽略图与图之间的连边就可以分开考虑，然后$dp$求出所有子图$($肉$)$的可能性，最后汇总即可，如果是分开的几个子图，设两个子图上的节点数为$d_1$,$d_2$，拓扑序的可能性为$a$，$b$,那么两个子图的可能性就是$C_{d1+d2}^{d1}ab$，这是很好解决的，但是图中是连着的，所以我们要分几种情况考虑。
对于一个子图，进行拓扑时删掉的点的数量上面的肯定比下面的多，而拓扑限制只有与一号子图的第一条连边，不删掉这条边是无法得到下面这块“肉”的，所以我们的$dp$需要维护的东西比较多：当前第二排架子删到了第几个；第一排架子删到了第几个；维护的“肉”还剩多少。这时你会发现：我们写出的$dp$是一个三维的$dp$，超时。
哇我辛辛苦苦求出来的居然超时呜呜呜呜
于是我们又想到如何压缩$dp$的维度：
如果在一个完整的图中只删去了一个点，我们只需要一维：

如果只删去第一行的，那么只需要两维：

如果把第二行也删了，也只需要两维：

这样我们就成功的把$dp$压成了二维。

$AC$ $Code$:

#include 
#define ll long long
using namespace std;
const int MAXN = 3005;
const int MAXM = MAXN * MAXN * 2;
int n , mod , dp[ MAXN ][ MAXN ] , mul[ MAXM ] , inv[ MAXM ];
int ksm (int a, int b) {
    int ret = 1;
    while (b) {
        if(b & 1) ret = 1ll * ret * a % mod;
        a = 1ll * a * a % mod;
        b >>= 1;
    }
    return ret;
}
int C (int n, int m) {
    if( n < m || m < 0 ) return 0;
    return 1ll * mul[ n ] % mod * inv[ m ] % mod * inv[ n - m ] % mod;
}
int Cat (int n, int m) { return ( C(n * 2 - m, n ) - C( n * 2 - m, n - m - 1 ) + mod ) % mod; }
int main () {
    scanf( "%d%d" , &n, &mod );
    mul[ 0 ] = dp[ 0 ][ 0 ] = 1;
    for ( int i = 1 ; i < MAXM ; ++i ) mul[ i ] = 1ll * mul[ i - 1 ] * i % mod;
    inv[ MAXM - 1 ] = ksm( mul[ MAXM - 1 ] , mod - 2 );
    for ( int i = MAXM - 2 ; ~ i ; --i ) inv[ i ] = 1ll * inv[ i + 1 ] * ( i + 1 ) % mod;
    for ( int i = 0 ; i < n ; ++i )
        for ( int j = 0 ; j < n ; ++j ) {
            if ( i == n - 1 && j == n - 1 ) continue;
            int left = 2 * n * n - min( i , j ) * n * 2 - i - j - 2;
            if ( j == i + 1)
                for ( int k = 0 ; k <= n; ++k )
                    dp[ i + 1 ][ j ] = ( dp[ i + 1 ][ j ] + 1ll * dp[ i ][ j ] * C( left - k , n * 2 - k ) % mod * Cat( n , k ) % mod ) % mod;
            else {
                dp[ i + 1 ][ j ] = ( dp[ i + 1 ][ j ] + dp[ i ][ j ] ) % mod;
                if ( i == j ) dp[ i ][ j + 1 ] = ( dp[ i ][ j + 1 ] + dp[ i ][ j ] ) % mod;
                else dp[ i ][ j + 1 ] = ( dp[ i ][ j + 1 ] + 1ll * dp[ i ][ j ] * C( left , n * 2 ) % mod * Cat( n , 0) % mod ) % mod;
            }
        }
    printf( "%d\n" , dp[ n - 1 ][ n - 1 ] );
}