結果
問題 | No.1751 Fortune Nim |
ユーザー |
![]() |
提出日時 | 2021-11-19 22:20:58 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 86 ms / 2,000 ms |
コード長 | 1,421 bytes |
コンパイル時間 | 1,491 ms |
コンパイル使用メモリ | 167,392 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2025-01-01 18:52:49 |
合計ジャッジ時間 | 4,621 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 30 |
ソースコード
#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> mf = {1};vector<long long> mfi = {1};long long modfact(int n){if (mf.size() > n){return mf[n];} else {for (int i = mf.size(); i <= n; i++){long long next = mf.back() * i % MOD;mf.push_back(next);mfi.push_back(modinv(next));}return mf[n];}}long long modfactinv(int n){if (mfi.size() > n){return mfi[n];} else {return modinv(modfact(n));}}long long modbinom(int n, int k){if (n < 0 || k < 0 || k > n){return 0;} else {return modfact(n) * modfactinv(k) % MOD * modfactinv(n - k) % MOD;}}int main(){int N;cin >> N;vector<int> a(N);for (int i = 0; i < N; i++){cin >> a[i];}long long S = 0;int X = 0;for (int i = 0; i < N; i++){S += a[i];X ^= a[i];}long long cnt;if (X == 0){cnt = S;} else {cnt = S;for (int i = 0; i < N; i++){if ((a[i] ^ X) < a[i]){cnt = min(cnt, S - a[i] + (a[i] ^ X) + 1);}}}long long r = (MOD - 1) * modinv(3) % MOD;long long ans = (1 + MOD - modpow(r, cnt)) * modinv(2) % MOD;cout << ans << endl;}