結果

問題 No.2513 Power Eraser
ユーザー 👑 amentorimaruamentorimaru
提出日時 2023-06-16 12:21:26
言語 C++17(clang)
(17.0.6 + boost 1.83.0)
結果
AC  
実行時間 1,704 ms / 6,000 ms
コード長 2,089 bytes
コンパイル時間 3,564 ms
コンパイル使用メモリ 144,340 KB
実行使用メモリ 54,512 KB
最終ジャッジ日時 2024-09-20 17:01:33
合計ジャッジ時間 55,608 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,944 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 5 ms
6,940 KB
testcase_04 AC 28 ms
6,944 KB
testcase_05 AC 7 ms
6,944 KB
testcase_06 AC 39 ms
6,940 KB
testcase_07 AC 18 ms
6,940 KB
testcase_08 AC 32 ms
6,944 KB
testcase_09 AC 45 ms
6,940 KB
testcase_10 AC 30 ms
6,940 KB
testcase_11 AC 6 ms
6,944 KB
testcase_12 AC 24 ms
6,944 KB
testcase_13 AC 559 ms
25,760 KB
testcase_14 AC 1,201 ms
48,724 KB
testcase_15 AC 1,281 ms
49,772 KB
testcase_16 AC 730 ms
27,912 KB
testcase_17 AC 917 ms
30,528 KB
testcase_18 AC 1,198 ms
48,420 KB
testcase_19 AC 1,678 ms
54,360 KB
testcase_20 AC 623 ms
26,672 KB
testcase_21 AC 1,261 ms
49,564 KB
testcase_22 AC 1,260 ms
49,348 KB
testcase_23 AC 1,684 ms
54,512 KB
testcase_24 AC 1,681 ms
54,388 KB
testcase_25 AC 1,681 ms
54,392 KB
testcase_26 AC 1,679 ms
54,388 KB
testcase_27 AC 1,685 ms
54,384 KB
testcase_28 AC 1,704 ms
54,264 KB
testcase_29 AC 1,682 ms
54,384 KB
testcase_30 AC 1,682 ms
54,260 KB
testcase_31 AC 1,696 ms
54,388 KB
testcase_32 AC 1,681 ms
54,388 KB
testcase_33 AC 2 ms
6,944 KB
testcase_34 AC 1 ms
6,940 KB
testcase_35 AC 1,674 ms
54,260 KB
testcase_36 AC 1,676 ms
54,388 KB
7_evil_case_1.txt AC 3,855 ms
108,956 KB
7_evil_case_2.txt AC 3,853 ms
108,960 KB
7_evil_case_3.txt AC 3,856 ms
109,084 KB
7_evil_case_4.txt AC 3,881 ms
108,956 KB
7_evil_case_5.txt AC 3,859 ms
108,952 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;
  vector<mint2> ff;
  ff.reserve(f.size());
  while (t < mod) {
    t *= 2;
    while (ff.size() < t && ff.size() < f.size()) {
      ff.push_back(f[ff.size()]);
    }
    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);
  if (rf.size() > siz)
    rf.resize(siz);
  if (rv.size() > siz)
    rv.resize(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