結果

問題 No.2994 べき内積
ユーザー Kude
提出日時 2024-12-19 05:06:41
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 471 ms / 2,000 ms
コード長 2,066 bytes
コンパイル時間 4,092 ms
コンパイル使用メモリ 287,456 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-12-19 05:06:52
合計ジャッジ時間 9,904 ms
ジャッジサーバーID
(参考情報)
judge5 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<bits/stdc++.h>
namespace {
#pragma GCC diagnostic ignored "-Wunused-function"
#include<atcoder/all>
#pragma GCC diagnostic warning "-Wunused-function"
using namespace std;
using namespace atcoder;
#define rep(i,n) for(int i = 0; i < (int)(n); i++)
#define rrep(i,n) for(int i = (int)(n) - 1; i >= 0; i--)
#define all(x) begin(x), end(x)
#define rall(x) rbegin(x), rend(x)
template<class T> bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } else return false; }
template<class T> bool chmin(T& a, const T& b) { if (b < a) { a = b; return true; } else return false; }
using ll = long long;
using P = pair<int,int>;
using VI = vector<int>;
using VVI = vector<VI>;
using VL = vector<ll>;
using VVL = vector<VL>;
constexpr int p = 1009;
using mint = static_modint<p>;

} int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  int m, n;
  cin >> m >> n;
  m++, n++;
  VI k(m);
  vector<mint> a(n);
  rep(i, m) cin >> k[i];
  rep(i, n) { int x; cin >> x; a[i] = x; }
  int mx = 0;
  rrep(i, m) mx = min(mx * p + k[i], n - 1);
  VI knmx = k;
  knmx[0] -= mx;
  rep(i, m) {
    if (knmx[i] >= 0) break;
    int c = (-knmx[i] + p - 1) / p;
    knmx[i] += c * p;
    knmx[i+1] -= c;
  }
  vector<mint> coeff(mx + 1);
  {
    mint x = 1;
    rrep(i, m) x = x.pow(p) * a[0].pow(knmx[i]);
    coeff[mx] = x;
    rrep(i, mx) coeff[i] = coeff[i+1] * a[0];
  }
  {
    mint c = 1;
    int e = 0;
    rep(i, mx + 1) {
      coeff[i] *= e == 0 ? c : mint();
      if (i == mx) break;
      if (i + 1 == p) e--;
      else c /= i + 1;
      int idx = 0;
      while (k[idx] == 0) idx++;
      e += idx;
      c *= k[idx]--;
      rep(i, idx) k[i] = p - 1;
    }
  }
  using m9 = static_modint<2032156673>;
  vector<m9> f(n), fk(n);
  rep(i, n) if (i) f[i] = a[i].val();
  fk[0] = 1;
  vector<mint> ans(n);
  rep(k, mx + 1) {
    rep(i, n) ans[i] += 1ULL * fk[i].val() * coeff[k].val();
    if (k == mx) break;
    fk = convolution(f, fk);
    fk.resize(n);
    rep(i, n) fk[i] = fk[i].val() % p;
  }
  rep(i, n) cout << ans[i].val() << " \n"[i + 1 == n];
}
0