結果

問題 No.2475 Distance Permutation
ユーザー akakimidoriakakimidori
提出日時 2023-09-24 14:51:08
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,712 ms / 5,000 ms
コード長 1,885 bytes
コンパイル時間 4,941 ms
コンパイル使用メモリ 193,848 KB
実行使用メモリ 45,532 KB
最終ジャッジ日時 2023-09-24 14:52:03
合計ジャッジ時間 52,227 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 721 ms
45,324 KB
testcase_01 AC 723 ms
45,292 KB
testcase_02 AC 1,675 ms
45,244 KB
testcase_03 AC 1,676 ms
45,240 KB
testcase_04 AC 1,674 ms
45,332 KB
testcase_05 AC 1,697 ms
45,352 KB
testcase_06 AC 1,682 ms
45,280 KB
testcase_07 AC 1,712 ms
45,244 KB
testcase_08 AC 1,676 ms
45,276 KB
testcase_09 AC 1,642 ms
45,300 KB
testcase_10 AC 1,677 ms
45,256 KB
testcase_11 AC 1,672 ms
45,300 KB
testcase_12 AC 1,700 ms
45,364 KB
testcase_13 AC 1,246 ms
45,276 KB
testcase_14 AC 1,225 ms
45,312 KB
testcase_15 AC 1,237 ms
45,532 KB
testcase_16 AC 1,234 ms
45,296 KB
testcase_17 AC 1,663 ms
45,300 KB
testcase_18 AC 1,639 ms
45,332 KB
testcase_19 AC 1,632 ms
45,340 KB
testcase_20 AC 1,643 ms
45,324 KB
testcase_21 AC 1,633 ms
45,328 KB
testcase_22 AC 1,656 ms
45,284 KB
testcase_23 AC 1,633 ms
45,324 KB
testcase_24 AC 1,630 ms
45,332 KB
testcase_25 AC 1,630 ms
45,236 KB
testcase_26 AC 1,631 ms
45,368 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <atcoder/all>
#include <cmath>
#include <cstdint>
#include <iostream>
#include <map>

using namespace std;
using namespace atcoder;

using mint = modint998244353;
using i32 = int32_t;

#define REP(i, s, t) for (i32 i = (s); i < (t); ++i)
#define PER(i, s, t) for (i32 i = (t)-1; i >= (s); --i)

vector<mint> fact, ifact, inv;

const i32 m = 100000;

void init(void) {
    fact.resize(m + 1, 1);
    ifact.resize(m + 1, 1);
    inv.resize(m + 1, 1);
    REP(i, 1, m + 1) fact.at(i) = i * fact.at(i - 1);
    ifact.at(m) = fact.at(m).inv();
    PER(i, 1, m) {
        inv.at(i + 1) = ifact.at(i + 1) * fact.at(i);
        ifact.at(i) = ifact.at(i + 1) * (i + 1);
    }
}

i32 read_int(void) {
    i32 n;
    cin >> n;
    return n;
}

int main(void) {
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    init();
    const i32 k = read_int();
    const i32 q = read_int();
    vector<mint> dp(m + 1);
    dp[0] = 1;
    mint sum = 1;
    REP(i, 1, m + 1) {
        if (i > k) sum -= dp[i - k - 1];
        dp[i] = sum * inv[i];
        sum += dp[i];
    }
    const i32 batch = 2000;
    vector<vector<mint>> memo;
    REP(i, 0, m / batch) {
        vector<mint> a(dp.begin(), dp.begin() + i * batch);
        vector<mint> b = convolution(a, dp);
        if (b.empty()) {
            b.resize(m + 1);
        }
        memo.emplace_back(b);
    }
    auto calc = [&](const i32 n, const i32 k) -> mint {
        assert(k <= n);
        if (k < 0) return 0;
        const i32 p = k / batch;
        mint res = memo[p][n];
        REP(i, p * batch, k + 1) res += dp[i] * dp[n - i];
        return res;
    };
    REP(i, 0, q) {
        const i32 n = read_int();
        const i32 l = read_int();
        const i32 r = read_int();
        mint ans = calc(n - 1, r - 1) - calc(n - 1, l - 2);
        ans *= fact[n - 1];
        cout << ans.val() << "\n";
    }
    return 0;
}
0