結果

問題 No.2475 Distance Permutation
ユーザー nok0
提出日時 2023-09-09 00:51:59
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
CE  
(最新)
AC  
(最初)
実行時間 -
コード長 2,226 bytes
コンパイル時間 2,773 ms
コンパイル使用メモリ 202,436 KB
最終ジャッジ日時 2025-02-16 20:31:39
ジャッジサーバーID
(参考情報)
judge5 / judge4
このコードへのチャレンジ
(要ログイン)
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。(リジャッジ後のコンパイルエラーは公開されます)
ただし、clay言語の場合は開発者のデバッグのため、公開されます。

コンパイルメッセージ
In file included from /usr/include/c++/13/vector:63,
                 from /usr/include/atcoder/convolution.hpp:8,
                 from /usr/include/atcoder/convolution:1,
                 from /usr/include/atcoder/all:1,
                 from main.cpp:5:
/usr/include/c++/13/bits/allocator.h: In destructor ‘std::_Vector_base<char, std::allocator<char> >::_Vector_impl::~_Vector_impl()’:
/usr/include/c++/13/bits/allocator.h:184:7: error: inlining failed in call to ‘always_inline’ ‘std::allocator< <template-parameter-1-1> >::~allocator() noexcept [with _Tp = char]’: target specific option mismatch
  184 |       ~allocator() _GLIBCXX_NOTHROW { }
      |       ^
In file included from /usr/include/c++/13/vector:66:
/usr/include/c++/13/bits/stl_vector.h:133:14: note: called from here
  133 |       struct _Vector_impl
      |              ^~~~~~~~~~~~

ソースコード

diff #

#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")

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

#include "testlib.h"

using namespace std;
using namespace atcoder;

using mint = modint998244353;
using i32 = int32_t;
using u64 = uint64_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);
    }
}

int main(int argc, char **argv) {
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    registerValidation(argc, argv);
    init();
    const i32 k = inf.readInt(1, m);
    inf.readSpace();
    const i32 q = inf.readInt(1, m * 2);
    inf.readEoln();

    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 auto ep = convolution(dp, dp);
    REP(i, 0, q) {
        const i32 n = inf.readInt(1, m);
        inf.readSpace();
        const i32 l = inf.readInt(1, n);
        inf.readSpace();
        const i32 r = inf.readInt(l, n);
        inf.readEoln();

        __int128_t sum = 0;
        const i32 len = r - l + 1;
        if (len <= n - 1 - len) {
            REP(i, l, r + 1) {
                sum += (u64)dp[i - 1].val() * dp[n - i].val();
            }
        } else {
            sum = ep[n - 1].val();
            REP(i, 1, l) {
                sum += (u64)998244353 * 998244353 - (u64)dp[i - 1].val() * dp[n - i].val();
            }
            REP(i, r + 1, n + 1) {
                sum += (u64)998244353 * 998244353 - (u64)dp[i - 1].val() * dp[n - i].val();
            }
        }
        sum %= 998244353;
        mint ans = sum;
        ans *= fact[n - 1];
        cout << ans.val() << "\n";
    }
    inf.readEof();
    return 0;
}
0