結果
問題 | No.2457 Stampaholic (Easy) |
ユーザー |
![]() |
提出日時 | 2023-09-01 21:53:52 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 144 ms / 4,000 ms |
コード長 | 983 bytes |
コンパイル時間 | 2,046 ms |
コンパイル使用メモリ | 196,664 KB |
最終ジャッジ日時 | 2025-02-16 16:33:09 |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 19 |
ソースコード
#include <bits/stdc++.h> using namespace std; const long long MOD = 998244353; long long modpow(long long a, long long b){ long long ans = 1; while (b > 0){ if (b % 2 == 1){ ans *= a; ans %= MOD; } a *= a; a %= MOD; b /= 2; } return ans; } long long modinv(long long a){ return modpow(a, MOD - 2); } vector<long long> get(int H, int K){ vector<long long> ans(K + 1, 0); int p = min(H - K + 1, K); for (int i = 1; i < p; i++){ ans[i] = 2; } ans[p] = H - 2 * (p - 1); return ans; } int main(){ int H, W, N, K; cin >> H >> W >> N >> K; vector<long long> A = get(H, K); vector<long long> B = get(W, K); long long ans = 0; long long C = modinv(H - K + 1) * modinv(W - K + 1) % MOD; for (int i = 1; i <= K; i++){ for (int j = 1; j <= K; j++){ long long P = i * j * C % MOD; ans += (MOD + 1 - modpow(MOD + 1 - P, N)) % MOD * A[i] % MOD * B[j] % MOD; ans %= MOD; } } cout << ans << endl; }