結果
問題 | No.2396 等差二項展開 |
ユーザー | SSRS |
提出日時 | 2023-07-28 21:44:10 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,254 ms / 6,000 ms |
コード長 | 872 bytes |
コンパイル時間 | 1,765 ms |
コンパイル使用メモリ | 171,056 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-10-06 18:12:01 |
合計ジャッジ時間 | 12,037 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 31 |
ソースコード
#include <bits/stdc++.h> using namespace std; int main(){ long long N, M; int L, K; long long B; cin >> N >> M >> L >> K >> B; M %= B; auto convolution = [&](vector<long long> a, vector<long long> b){ vector<long long> c(L, 0); for (int i = 0; i < L; i++){ for (int j = 0; j < L; j++){ if (i + j < L){ c[i + j] += a[i] * b[j]; c[i + j] %= B; } else { c[i + j - L] += a[i] * b[j] % B * M; c[i + j - L] %= B; } } } return c; }; vector<long long> f(L, 0); f[0] += 1; f[0] %= B; if (L == 1){ f[0] += M; f[0] %= B; } else { f[1] += 1; f[1] %= B; } vector<long long> g(L, 0); g[0] = 1; g[0] %= B; while (N > 0){ if (N % 2 == 1){ g = convolution(g, f); } f = convolution(f, f); N /= 2; } cout << g[K] << endl; }