[题解]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;  
}