QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#227359 | #6801. Blackjack | jzh | WA | 0ms | 3788kb | C++14 | 1.9kb | 2023-10-27 13:07:15 | 2023-10-27 13:07:15 |
Judging History
answer
#include <algorithm>
#include<bits/stdc++.h>
using namespace std;
using db = long double;
const int maxn = 510;
void solve() {
int n, a, b;
cin >> n >> a >> b;
vector<int> vec(n + 10);
for (int i = 0; i < n; i++) cin >> vec[i];
vector<vector<db>> dp(n + 10, vector<db>(b + 10, 0));
auto ndp = dp;
dp[0][0] = 1;
for (int i = 0; i < n; i++) {
for (auto &v: ndp) fill(begin(v), end(v), 0);
for (int x = 0; x <= i; x++) {
for (int y = 0; y <= b; y++) {
ndp[x][y] += dp[x][y];
}
for (int y = 0; y <= b - vec[i]; y++) {
ndp[x + 1][y + vec[i]] += dp[x][y] * (x + 1) / (n - x);
}
}
dp = ndp;
}
for (auto &v: ndp) fill(begin(v), end(v), 0);
db ans = 0;
for (int i = 0; i < n; i++) {
ndp = dp;
for (int x = 0; x < n; x++) {
for (int y = 0; y <= b - vec[i]; y++) {
assert(y + vec[i] <= b);
ndp[x + 1][y + vec[i]] -= dp[x][y] * (x + 1) / (n - x);
}
}
db temp = 0;
for (int x = 0; x < n; x++) {
db t1 = 0;
for (int y = max(0, a - vec[i] + 1); y <= min(a, b - vec[i]); y++) {
t1 += ndp[x][y];
if (ndp[x][y] / (n - x) > 1) {
cout << " error " << ndp[x][y] << " in " << x << ',' << y << endl;
}
}
temp += t1 / (n - x + 1e-9);
if (temp > 1) {
cout << " cur temp " << temp << " in " << i << "," << x << ',' << " t1 " << t1 << endl;
return;
}
}
ans += temp;
}
cout << fixed << setprecision(30) << (long double) ans << endl;
}
int main() {
ios::sync_with_stdio(false);
cin.tie(nullptr);
solve();
return 0;
}
详细
Test #1:
score: 0
Wrong Answer
time: 0ms
memory: 3788kb
input:
5 2 4 1 1 1 5 5
output:
0.000000000000000000000000000000
result:
wrong answer 1st numbers differ - expected: '0.1000000', found: '0.0000000', error = '0.1000000'