社区讨论
求助!WA30pts
P3175[HAOI2015] 按位或参与者 1已保存回复 1
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 1 条
- 当前快照
- 1 份
- 快照标识符
- @mhjt3kk7
- 此快照首次捕获于
- 2025/11/04 08:04 4 个月前
- 此快照最后确认于
- 2025/11/04 08:04 4 个月前
借鉴的是这篇题解https://www.luogu.com.cn/article/kdxxouxc
CPP#include <bits/stdc++.h>
using namespace std;
const int N = (1 << 20) + 100;
int n;
double res, f[2][N], p[N];
void sos(double *g) {
for(int s = 0; s < n; ++s)
for(int i = 0; i < (1 << n); ++i)
if(i >> s & 1)
g[i] += g[i ^ (1 << s)];
return;
}
void isos(double *g) {
for(int s = 0; s < n; ++s)
for(int i = 0; i < (1 << n); ++i)
if(i >> s & 1)
g[i] -= g[i ^ (1 << s)];
return;
}
signed main() {
ios :: sync_with_stdio(false);
cin.tie(0), cout.tie(0);
f[0][0] = 1;
cin >> n;
for(int i = 0; i < (1 << n); ++i) cin >> p[i];
sos(p);
for(int i = 1; ; ++i) {
if(1. * clock() / CLOCKS_PER_SEC > 0.95) break;
sos(f[i - 1 & 1]);
for(int j = 0; j < (1 << n); ++j) f[i & 1][j] = f[i - 1 & 1][j] * p[j];
isos(f[i & 1]);
res += i * f[i & 1][(1 << n) - 1];
f[i & 1][(1 << n) - 1] = 0;
}
if(res == 0) cout << "INF\n";
else cout << fixed << setprecision(10) << res << '\n';
return 0;
}
回复
共 1 条回复,欢迎继续交流。
正在加载回复...