結果

問題 No.2688 Cell Proliferation (Hard)
ユーザー てんぷらてんぷら
提出日時 2024-03-20 23:19:19
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 734 ms / 4,000 ms
コード長 1,852 bytes
コンパイル時間 6,033 ms
コンパイル使用メモリ 317,024 KB
実行使用メモリ 37,064 KB
最終ジャッジ日時 2024-03-20 23:19:39
合計ジャッジ時間 19,072 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,548 KB
testcase_01 AC 2 ms
6,548 KB
testcase_02 AC 2 ms
6,548 KB
testcase_03 AC 70 ms
7,004 KB
testcase_04 AC 687 ms
34,704 KB
testcase_05 AC 367 ms
20,128 KB
testcase_06 AC 162 ms
11,564 KB
testcase_07 AC 163 ms
11,364 KB
testcase_08 AC 639 ms
33,608 KB
testcase_09 AC 641 ms
32,448 KB
testcase_10 AC 683 ms
35,668 KB
testcase_11 AC 659 ms
34,108 KB
testcase_12 AC 176 ms
12,312 KB
testcase_13 AC 706 ms
37,048 KB
testcase_14 AC 598 ms
30,644 KB
testcase_15 AC 734 ms
37,064 KB
testcase_16 AC 572 ms
30,288 KB
testcase_17 AC 362 ms
21,272 KB
testcase_18 AC 696 ms
35,980 KB
testcase_19 AC 364 ms
21,376 KB
testcase_20 AC 284 ms
17,272 KB
testcase_21 AC 289 ms
17,252 KB
testcase_22 AC 355 ms
20,932 KB
testcase_23 AC 299 ms
18,056 KB
testcase_24 AC 698 ms
35,424 KB
testcase_25 AC 179 ms
12,512 KB
testcase_26 AC 713 ms
36,188 KB
testcase_27 AC 367 ms
21,512 KB
testcase_28 AC 578 ms
30,404 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <atcoder/all>
#include <bits/stdc++.h>
using ll = long long;
using ull = unsigned long long;
#define rep(i, n) for(int i = 0; i < (int)(n); i++)
#define REP(i, m, n) for(int i = (int)(m); i < (int)(n); i++)
using namespace std;
using namespace atcoder;
using mint = modint998244353;
const int inf = 1000000007;
const ll longinf = 1ll << 60;

namespace fps {
using P = vector<mint>;

void dft(P &f) {
    atcoder::internal::butterfly(f);
}

void idft(P &f) {
    atcoder::internal::butterfly_inv(f);
    mint c = mint(f.size()).inv();
    for(auto &e : f)
        e *= c;
}
P inv(P &f, int size = -1) {
    assert(f[0] != 0);
    if(size == -1) {
        size = f.size();
    }
    P g(size);
    g[0] = f[0].inv();
    for(int d = 1; d < size; d <<= 1) {
        P gg(2 * d), ff(2 * d);
        for(int j = 0; j < min((int)f.size(), 2 * d); j++)
            ff[j] = f[j];
        for(int j = 0; j < d; j++) {
            gg[j] = g[j];
        }

        dft(ff);
        dft(gg);
        for(int j = 0; j < 2 * d; j++)
            ff[j] *= gg[j];
        idft(ff);

        for(int j = 0; j < d; j++) {
            ff[j] = 0;
        }
        dft(ff);
        for(int j = 0; j < 2 * d; j++)
            ff[j] *= gg[j];
        idft(ff);

        for(int j = d; j < min(2 * d, size); j++)
            g[j] = -ff[j];
    }
    return g;
}
} // namespace fps
int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int p1, p2, q1, q2, t;
    cin >> p1 >> p2 >> q1 >> q2 >> t;
    mint p = mint(p1) / p2;
    mint q = mint(q1) / q2;
    vector<mint> bumbo(t + 1), bunshi(t + 1);
    ll c = 0;
    rep(i, t + 1) {
        bumbo[i] = -p * q.pow(c);
        c += i;
        bunshi[i] = q.pow(c);
    }
    bumbo[0] = 1;
    auto ret = convolution(bunshi, fps::inv(bumbo));
    cout << ret[t].val() << endl;
    return 0;
}
0