結果
| 問題 |
No.1426 Got a Covered OR
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2021-01-08 11:00:47 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 45 ms / 2,000 ms |
| コード長 | 2,061 bytes |
| コンパイル時間 | 787 ms |
| コンパイル使用メモリ | 72,940 KB |
| 最終ジャッジ日時 | 2025-01-17 10:35:45 |
|
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 24 |
ソースコード
#include <iostream>
#include <vector>
using namespace std;
using ll = long long;
const ll MOD = 1000000007;
ll modpow(ll a, ll n) {
ll res = 1;
while (n > 0) {
if (n & 1) res = res * a % MOD;
a = a * a % MOD;
n >>= 1;
}
return res;
}
ll modinv(ll i) {
return modpow(i, MOD - 2);
}
ll fact(ll n) {
ll ans = 1;
for (int i = 0; i < n; i++) {
ans = ans * (i + 1) % MOD;
}
return ans;
}
ll comb(ll n, ll a) {
ll ans = 1, denom = 1;
for (int i = 0; i < a; i++) {
ans = ans * (n - i) % MOD;
denom = denom * (i + 1) % MOD;
}
ans = ans * modinv(denom) % MOD;
return ans;
}
ll calc(ll a, ll b, ll c) {
ll ret = 0;
for (int i = 0; i <= b; i++) {
ll v = comb(b, i) * modpow((1 << (b + c - i)) - 1, a) % MOD;
if (i % 2 == 0) {
ret += v;
ret %= MOD;
}
else {
ret += MOD - v;
ret %= MOD;
}
}
return ret;
}
int main()
{
ios::sync_with_stdio(false);
cin.tie(nullptr);
int i, n;
cin >> n;
vector<ll> b(n, 0);
for (i = 0; i < n; i++) {
cin >> b[i];
}
ll prev = 0;
ll streak = 0;
ll ans = 1;
for (i = 0; i < n; i++) {
if (b[i] == -1)
{
streak++;
}
else {
ll c01 = 0;
ll c11 = 0;
for (int j = 0; j < 32; j++) {
ll x = (prev >> j) % 2;
ll y = (b[i] >> j) % 2;
if (x == 0 && y == 1)
{
c01++;
}
else if (x == 1 && y == 1)
{
c11++;
}
else if (x == 1 && y == 0)
{
cout << 0 << "\n";
return 0;
}
}
ans *= calc(streak + 1, c01, c11);
ans %= MOD;
prev = b[i];
streak = 0;
}
}
cout << ans << "\n";
return 0;
}