結果

問題 No.2994 べき内積
ユーザー 2qbingxuan2qbingxuan
提出日時 2024-12-22 04:27:01
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,921 bytes
コンパイル時間 3,010 ms
コンパイル使用メモリ 247,324 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-12-22 04:27:06
合計ジャッジ時間 5,294 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 WA -
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

#define all(x) begin(x), end(x)
#ifdef CKISEKI
#define safe cerr << __PRETTY_FUNCTION__ << " line " << __LINE__ << " safe\n"
#define debug(a...) debug_(#a, a)
#define orange(a...) orange_(#a, a)
#include <experimental/iterator>
void debug_(auto s, auto ...a) {
  cerr << "\e[1;32m(" << s << ") = (";
  int f = 0;
  (..., (cerr << (f++ ? ", " : "") << a));
  cerr << ")\e[0m\n";
}
void orange_(auto s, auto L, auto R) {
  cerr << "\e[1;33m[ " << s << " ] = [ ";
  using namespace experimental;
  copy(L, R, make_ostream_joiner(cerr, ", "));
  cerr << " ]\e[0m\n";
}
#else
#define safe ((void)0)
#define debug(...) safe
#define orange(...) safe
#endif

using lld = int64_t;

int main() {
  constexpr int p = 1009;
  constexpr int p2 = p * p;
  constexpr lld mod = p * p * (p - 1);

  cin.tie(nullptr)->sync_with_stdio(false);
  int M, N;
  cin >> M >> N;

  const lld inf = 1e15;
  lld K_capped = 0, K_mod = 0;

  {
    vector<int> k(M + 1);
    for (int i = 0; i <= M; i++) {
      cin >> k[i];
    }
    for (int i = M; i >= 0; i--) {
      K_capped = min(inf, K_capped * p + k[i]);
      K_mod = (K_mod * p + k[i]) % mod;
    }
  }


  vector<int> e(N + 1), r(N + 1);
  r[0] = 1;

  for (int i = 0; i <= N; i++)
    e[i] = rand() % p;
    // cin >> e[i];

  auto mul = [&](const auto &a, const auto &b) {
    vector<int> c(N + 1);
    for (int i = 0; i <= N; i++) {
      for (int j = 0; i + j <= N; j++) {
        c[i + j] += a[i] * b[j];
        if (c[i + j] >= p2) c[i + j] -= p2;
      }
    }
    for (int i = 0; i <= N; i++)
      c[i] %= p;
    return c;
  };

  // lld K = K_mod + (K_capped - K_mod) / mod * mod;
  lld K = K_mod;
  while (K) {
    if (K & 1) r = mul(r, e);
    e = mul(e, e);
    K >>= 1;
  }
  /*
  for (int i = 0; i <= N; i++)
    if (r[i] != 0)
      debug(i);
      */

  for (int i = 0; i <= N; i++)
    cout << r[i] << (i==N ? '\n' : ' ');
}
0