Codeforces Global Round 14 E. Phoenix and Computers

刚开始自己想的是区间dp，写完之后发现一点问题调不出来，然后就搜了搜博客。
大佬使用线性dp写的，很容易理解，下面贴了链接我就不多写了。

大佬博客

#include <bits/stdc++.h>
using namespace std;
typedef  long long  ll;
#define int long long
const int N = 110000;
int dp[600][600];
int mod;
ll F[N], invn[N], invF[N], bit[N];

void init() {
    bit[1] = F[0] = F[1] = invn[0] = invn[1] = invF[0] = invF[1] = 1;
    for(int i = 2;i < N; i++){
        F[i] = F[i - 1] * i % mod;
        invn[i] = (mod - mod / i) * invn[mod % i] % mod;
        invF[i] = invF[i - 1] * invn[i] % mod;
        bit[i] = bit[i - 1] * 2 % mod;
    }
}

ll C(ll m, ll n) {
    if(m < 0 || n < 0 || n > m)
        return 0;
    ll ans = F[m];
    ans = ans * invF[n] % mod;
    ans = ans * invF[m - n] % mod;
    return ans;
}


signed main()
{

    int n;
    cin >> n >> mod;
    init();
    for (int i = 1; i <= n; ++i) {
        dp[i][i] = bit[i];
        for (int j = 0; j <= i; ++j) {
            for (int k = 1; k + 1 + i <= n; ++k) {
                dp[i+1+k][j+k] = (dp[i+1+k][j+k]+ dp[i][j]%mod * bit[k]%mod * C(j+k,k))%mod;
            }
        }
    }
    int ans =0 ;
    for (int l = 0; l <= n; ++l) {
        ans += dp[n][l];
    }
    cout << ans%mod  << endl;
}