[题解]AT_abc281_g [ABC281G] Farthest City

思路

定义 \(dp_{i,j}\) 表示前若干层一共有 \(i\) 个点，且在这些层中最外层的点数为 \(j\) 的方案数。

那么，我们先求出最外层的方案数，在 \(n - (i - j) - 1\)（即减去前若干层中除最外层点的数量再减去 \(n\) 号点）中选出 \(j\) 个作为最外层的点。所以，贡献为 \(C_{n - (i - j) - 1}^j\)。

令 \(k\) 为前 \(i - j\) 个数的最后一层的点数。那么，最后一层的 \(j\) 个数可以与上一层的 \(k\) 个点两两连一条边（如果连接了上上层，则 \(d\) 会发生改变，所以不行），即对于每一个点都有 \(\sum_{p = 1}^{k}C_{k}^{p}\)，又根据二项式定理，化简为 \(2^k - 1\)。因此，对于所有的 \(j\) 个点，贡献为 \((2^k - 1)^j\)。

然后，对于那 \(j\) 个点，内部也可以两两连边，而最多建为一张完全图，边数为 \(\frac{j \times (j - 1)}{2}\)，记作 \(x\)。所以，方案数为 \(\sum_{p = 1}^xC_{x}^p\)，又由二项式定理得 \(2^x\)，也就是 \(2^{\frac{j \times (j - 1)}{2}}\)。

综上，可得出状态转移方程：

\[ dp_{i,j} = \sum_{k = 1}^{i - j}(dp_{i - j,k} \times C_{n - (i - j) - 1}^j \times (2^k - 1)^j \times 2^{\frac{j \times (j - 1)}{2}}) \]

答案为：\(dp_{n,1}\)。

注意要初始化 \(2^k,(2^k - 1)^j,C_n^m\)。

Code

#include <bits/stdc++.h>  
#define int long long  
#define re register  
  
using namespace std;  
  
const int N = 510;  
int n,mod;  
int pot[N * N];  
int C[N][N],q[N][N],dp[N][N];  
  
inline int read(){  
    int r = 0,w = 1;  
    char c = getchar();  
    while (c < '0' || c > '9'){  
        if (c == '-') w = -1;  
        c = getchar();  
    }  
    while (c >= '0' && c <= '9'){  
        r = (r << 3) + (r << 1) + (c ^ 48);  
        c = getchar();  
    }  
    return r * w;  
}  
  
inline int qmi(int a,int b){  
    int res = 1;  
    while (b){  
        if (b & 1) res = res * a % mod;  
        a = a * a % mod;  
        b >>= 1;  
    }  
    return res;  
}  
  
inline void init(){//初始化   
    C[0][0] = C[0][1] = pot[0] = 1;  
    for (re int i = 1;i <= n;i++){  
        C[i][0] = C[i][i] = 1;  
        for (re int j = 1;j < i;j++) C[i][j] = (C[i - 1][j] + C[i - 1][j - 1]) % mod;  
    }  
    for (re int i = 1;i <= n * (n - 1) / 2;i++) pot[i] = pot[i - 1] * 2 % mod;  
    for (re int i = 1;i <= n;i++){  
        for (re int j = 1;j <= n;j++) q[i][j] = qmi(pot[i] - 1,j);  
    }  
}  
  
signed main(){  
    n = read();  
    mod = read();  
    init();  
    dp[1][1] = 1;  
    for (re int i = 2;i <= n;i++){  
        for (re int j = 1;j < i;j++){  
            for (re int k = 1;k <= i - j;k++) dp[i][j] = (dp[i][j] + dp[i - j][k] * C[n - (i - j) - 1][j] % mod * q[k][j] % mod * pot[j * (j - 1) / 2] % mod) % mod;  
        }  
    }  
    printf("%lld",dp[n][1]);  
    return 0;  
}