結果

問題 No.2513 Power Eraser
ユーザー 👑 amentorimaruamentorimaru
提出日時 2023-06-16 12:19:46
言語 C++17(clang)
(17.0.6 + boost 1.83.0)
結果
AC  
実行時間 1,700 ms / 6,000 ms
コード長 1,867 bytes
コンパイル時間 3,424 ms
コンパイル使用メモリ 144,772 KB
実行使用メモリ 54,524 KB
最終ジャッジ日時 2024-09-20 17:00:37
合計ジャッジ時間 54,264 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 5 ms
5,376 KB
testcase_04 AC 26 ms
5,376 KB
testcase_05 AC 6 ms
5,376 KB
testcase_06 AC 39 ms
5,632 KB
testcase_07 AC 18 ms
5,376 KB
testcase_08 AC 31 ms
5,376 KB
testcase_09 AC 44 ms
5,888 KB
testcase_10 AC 29 ms
5,376 KB
testcase_11 AC 6 ms
5,376 KB
testcase_12 AC 24 ms
5,376 KB
testcase_13 AC 544 ms
25,880 KB
testcase_14 AC 1,178 ms
48,496 KB
testcase_15 AC 1,259 ms
49,720 KB
testcase_16 AC 804 ms
28,100 KB
testcase_17 AC 938 ms
30,652 KB
testcase_18 AC 1,156 ms
48,576 KB
testcase_19 AC 1,640 ms
54,364 KB
testcase_20 AC 607 ms
26,680 KB
testcase_21 AC 1,217 ms
49,620 KB
testcase_22 AC 1,217 ms
49,480 KB
testcase_23 AC 1,666 ms
54,396 KB
testcase_24 AC 1,634 ms
54,396 KB
testcase_25 AC 1,658 ms
54,396 KB
testcase_26 AC 1,618 ms
54,240 KB
testcase_27 AC 1,700 ms
54,524 KB
testcase_28 AC 1,682 ms
54,396 KB
testcase_29 AC 1,646 ms
54,268 KB
testcase_30 AC 1,686 ms
54,392 KB
testcase_31 AC 1,681 ms
54,396 KB
testcase_32 AC 1,676 ms
54,248 KB
testcase_33 AC 2 ms
5,376 KB
testcase_34 AC 1 ms
5,376 KB
testcase_35 AC 1,680 ms
54,396 KB
testcase_36 AC 1,639 ms
54,400 KB
7_evil_case_1.txt AC 3,839 ms
109,088 KB
7_evil_case_2.txt AC 3,748 ms
108,960 KB
7_evil_case_3.txt AC 3,753 ms
109,084 KB
7_evil_case_4.txt AC 3,737 ms
108,956 KB
7_evil_case_5.txt AC 3,762 ms
109,080 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<vector>
#include<algorithm>
#include <atcoder/all>
using namespace std;
using ll = long long;
using namespace atcoder;
using mint2 = modint998244353;

vector<mint2> rev(vector<mint2>& f, ll mod) {
  vector<mint2> res = { mint2(1) / f[0] };
  ll t = 1;
  while (t < mod) {
    t *= 2;
    auto m = convolution(convolution(res, res), f);
    vector<mint2> nres(t);
    for (ll i = 0; i < t; i++) {      
      if (i < res.size())
        nres[i] += 2 * res[i];
      if (i < m.size())
        nres[i] += -m[i];
    }
    res = nres;
  }
  return res;
}

vector<mint2> div(vector<mint2>& f, vector<mint2>& g) {
  auto rf = f;
  auto rg = g;
  reverse(rf.begin(), rf.end());
  reverse(rg.begin(), rg.end());
  ll siz = f.size() - g.size() + 1;
  if (siz <= 0)
    return f;
  auto rv = rev(rg, siz);
  auto q = convolution(rf, rv);
  if (q.size() > siz)
    q.resize(siz);
  reverse(q.begin(), q.end());
  auto gq = convolution(g, q);
  auto res = f;
  for (ll i = 0; i < res.size(); i++) {    
    if (i < gq.size()) 
      res[i] -= gq[i];
  }
  while (res.back() == 0 && res.size() > 1)
    res.pop_back();
  return res;
}

int main() {
  ll n;
  cin >> n;
  vector<ll> a(n);
  for (ll i = 0; i < n; i++)
    cin >> a[i];
  ll m = 1;
  while (m < n)
    m *= 2;
  vector<mint2> def = { 1 };
  vector<vector<mint2>> muls(m * 2, def);
  for (ll i = 0; i < n; i++)
    muls[i + m] = { a[i],-1 };
  for (ll i = m - 1; i > 0; i--)
    muls[i] = convolution(muls[i * 2], muls[i * 2 + 1]);
  
  mint2 ans = 1;
  vector<vector<mint2>> divs(m * 2, def);
  for (ll i = 2; i < m + n; i++) {
    queue<ll> que;
    divs[i] = convolution(divs[i], divs[i / 2]);
    if (i % 2)
      divs[i] = convolution(muls[i - 1], divs[i]);
    divs[i] = div(divs[i], muls[i]);
    if (m <= i)
      ans *= divs[i][0];
  }
  cout << ans.val() << endl;
  return 0;
}
0