codeforces 1515 E - Phoenix and Computers

E. Phoenix and Computers

这场晚上脑子真的不行我处刑我自己以为能凭手速上分结果没注意细节被fst话说您的pretest也太弱了叭这都能让我过你过不过分呜呜

flag再补个F,和队友的差距越来越大了这可不太好呜呜呜

题意: 有一排n台电脑,当把第i-1台和第i+1台开起来的时候第i台会自动开,问把这n台电脑都开起来的方案数

dfs找出规律是开k台电脑完全手动开不让他自动开机的次数是 2 k − 1 2^{k - 1} 2k1, 刚开始的想法是如果能知道自动开机的电脑的位置,可以将n台电脑分成几段,每一段可以用完全手动开的序列独立考虑,合并的时候是组合数学的插空法(?(就是将一个长度为a的序列和一个长度为b的序列不改变各自的顺序合在一起大概是 C a + b a C_{a+ b}^{a} Ca+ba

但是我不会找自动开机的位置

这不是个dp就是个dfs,然后我一看数据量我觉得铁dfs,我是什么傻逼

关于那个 2 k − 1 2^{k-1} 2k1, 官方给了两种证明,一种是:先开1 ,然后我们就只能接着2345…这样开,是 C k − 1 0 C_{k - 1}^{0} Ck10

先开2,345… 还是得按这个顺序但是1可以插在任意位置,就是 C k − 1 1 C_{k - 1}^{1} Ck11

同理

总共就是 C k − 1 0 + C k − 1 1 + ⋯ + C k − 1 k − 1 = 2 k − 1 C_{k - 1}^{0} + C_{k - 1}^{1} + \dots + C_{k - 1}^{k - 1} = 2 ^{k - 1} Ck10+Ck11++Ck1k1=2k1

自动开机设为0,手动设为1,最后这排电脑的开机是1111101111110111111… 这种序列

d p [ l e n ] [ c n t ] dp[len] [cnt] dp[len][cnt] := 到第len个位置,前面len 个电脑全开,手动开了cnt次的方案数 (第len个是自动开)

d p [ l e n + k + 1 ] [ c n t + k ] + = d p [ l e n ] [ c n t ] × C c n t + k k × 2 k − 1 . dp[len+k+1] [cnt+k] += dp[len] [cnt] \times C_{cnt + k}^{k} \times 2^{k - 1}. dp[len+k+1][cnt+k]+=dp[len][cnt]×Ccnt+kk×2k1.

答案就是 ∑ i = 1 n d p [ n + 1 ] [ i ] \sum_{ i = 1}^{n} dp[n + 1][i] i=1ndp[n+1][i]

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

constexpr int N = 510;

ll mod;

ll dp[N][N], f[N][N], pow2[N];

inline void init() {
    for (int i = 0; i < N; ++i) {
        f[0][i] = 1;
    }
    for (int i = 1; i < N; ++i) {
        for (int j = 1; j <= i; ++j) {
            f[j][i] = (f[j][i - 1] + f[j - 1][i - 1]) % mod;
        }
    }

    pow2[0] = 1;
    for (int i = 1; i < N; ++i) {
        pow2[i] = pow2[i - 1] * 2 % mod;
    }
}

int main() {
    ll n;
    cin >> n >> mod;

    init();

    dp[0][0] = 1;
    for (int i = 0; i <= n; ++i) {
        for (int j = 0; j <= n; ++j) {
            for (int k = 1; i + k <= n; ++k) {
                dp[i + k + 1][j + k] = (dp[i + k + 1][j + k] +
                                        dp[i][j] * f[k][j + k] % mod * pow2[k - 1] % mod) % mod;
            }
        }
    }

    ll ans = 0;
    for(int i = 1; i <= n; ++ i) {
        ans =(ans + dp[n + 1][i]) % mod;
    }
    cout << ans << endl;

    return 0;
}
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值