結果
問題 |
No.1189 Sum is XOR
|
ユーザー |
|
提出日時 | 2020-09-22 16:55:13 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 105 ms / 2,000 ms |
コード長 | 1,736 bytes |
コンパイル時間 | 1,639 ms |
コンパイル使用メモリ | 168,836 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-06-26 06:30:30 |
合計ジャッジ時間 | 4,062 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 21 |
ソースコード
#include <bits/stdc++.h> //#include "atcoder/all" typedef long long int ll; using namespace std; // using namespace atcoder; #define MAXN 1024 #define mod 998244353 ll Pow(ll a, ll n) { if (n == 0) return 1; if (n % 2) return (a * Pow(a, n - 1)) % mod; else { ll k = Pow(a, n / 2); return (k * k) % mod; } } ll gyakusuu[MAXN + 1]; ll cnt[MAXN]; ll dp[MAXN][11]; int main() { for (int i = 1; i <= MAXN; i++) { gyakusuu[i] = Pow(i, mod - 2); } ll n, K; cin >> n >> K; if (K >= 11) { cout << 0 << endl; return 0; } for (int i = 0; i < n; i++) { ll a; cin >> a; cnt[a]++; } memset(dp, 0, sizeof(dp)); dp[0][0] = 1; for (int i = 0; i < 1024; i++) { for (int j = 0; j < 1024; j++) { // p xor i = j とななるような p ll count = 1; ll p = 0; bool ok = true; ll i1 = i; ll j1 = j; for (int k = 1; k <= 10; k++) { if ((j1 % 2 == 0)&&(i1 % 2 == 1)) { ok = false; } else if ((j1 % 2 == 1)&&(i1 % 2 == 0)) { p += count; } count *= 2; i1 /= 2; j1 /= 2; } if (ok) { for (int k = 1; k <= 10; k++) { dp[j][k] = (dp[j][k] + dp[i][k - 1] * cnt[p]) % mod; } } } } ll ret = 0; for (int i = 1; i <= 1023; i++) { ret = (ret + dp[i][K]) % mod; } for (int i = 1; i <= K; i++) { ret = ret * gyakusuu[i] % mod; } cout << ret << endl; return 0; }