QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#773096#8037. Gambler's RuinAndevikingWA 1125ms111824kbC++203.0kb2024-11-23 00:51:542024-11-23 00:51:55

Judging History

你现在查看的是最新测评结果

  • [2024-11-23 00:51:55]
  • 评测
  • 测评结果:WA
  • 用时:1125ms
  • 内存:111824kb
  • [2024-11-23 00:51:54]
  • 提交

answer

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using pii = pair<int, int>;
#define range(x) (x).begin(), (x).end()
vector<long double> p;
vector<int> c;
vector<long double> ck, rck;
vector<ll> sum, rsum;
int n;
const long double EPS = 1e-12L;
long double f(long double x, long double y)
{
    ll sx = 0, sy = 0;
    int pos = upper_bound(range(ck), x) - ck.begin();
    if (pos - 1 >= int(sum.size()))
        sx = sum.back();
    else if (pos)
        sx = sum[pos - 1];
    // cout << pos << '\n';
    // cout << rck.size() << '\n';
    pos = upper_bound(range(rck) - pos, y) - rck.begin();
    if (pos - 1 >= int(rsum.size()))
        sy = rsum.back();
    else if (pos)
        sy = rsum[pos - 1];
    // cout << pos << '\n';
    // cout << sx << ' ' << sy << ' ' << x << ' ' << y << ' ' << sx + sy - max(sx * x, sy * y) << '\n';
    // if (sx + sy - max(sx * x, sy * y) > 15)
    //     cout << x << ' ' << y << ' ' << sx << ' ' << sy << '\n';
    return sx + sy - max(sx * x, sy * y);
}
long double get(long double x)
{
    long double l = 0, r = 1000010;
    long double lans = f(x, l), rans = f(x, r);
    long double ans = max(lans, rans);
    for (int _ = 0; _ < 3000; ++_) {
        long double lmid = l + (r - l) / 3;
        long double rmid = r - (r - l) / 3;

        lans = f(x, lmid), rans = f(x, rmid);
        ans = max(lans, rans);
        if (lans - rans >= EPS)
            r = rmid;
        else
            l = lmid;
    }
    // cout << x << ' ' << l << ' ' << ans << '\n';
    return ans;
}
void solve()
{
    cin >> n;
    p.resize(n + 5);
    c.resize(n + 5);
    for (int i = 1; i <= n; ++i)
        cin >> p[i] >> c[i];
    vector<pair<long double, int>> ckk;
    for (int i = 1; i <= n; ++i)
        ckk.emplace_back(p[i], c[i]);
    sort(range(ckk));
    for (const auto &[a, b] : ckk) {
        rck.push_back(1.L / (1 - a));
        rsum.push_back(b);
        if (rsum.size() > 1)
            rsum.back() += rsum.end()[-2];
    }

    reverse(range(ckk));
    for (const auto &[a, b] : ckk) {
        ck.push_back(1.L / a);
        sum.push_back(b);
        if (sum.size() > 1)
            sum.back() += sum.end()[-2];
    }
    // for (const auto &c : ck)
    //     cout << c << ' ';
    // cout << '\n';
    long double l = 0, r = 1000010;
    long double lans = get(l), rans = get(r);
    long double ans = max(lans, rans);
    for (int _ = 0; _ < 1000; ++_) {
        long double lmid = l + (r - l) / 3;
        long double rmid = r - (r - l) / 3;

        lans = get(lmid), rans = get(rmid);
        ans = max(lans, rans);
        if (lans - rans >= EPS)
            r = rmid;
        else
            l = lmid;
    }
    // cout << l << '\n';
    // cout << f(1 / .6L, 1.9) << '\n';
    // cout << get(1 / .6L) << '\n';
    cout << fixed << setprecision(10) << ans << '\n';
    // cout << f(0, 0) << '\n';
}

int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    cout.tie(0);
    int T = 1;
    while (T--)
        solve();
    return 0;
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 207ms
memory: 3816kb

input:

2
1 15
0 10

output:

10.0000000000

result:

ok found '10.0000000', expected '10.0000000', error '0.0000000'

Test #2:

score: 0
Accepted
time: 221ms
memory: 3948kb

input:

3
0.4 100
0.5 100
0.6 100

output:

33.3333333333

result:

ok found '33.3333333', expected '33.3333333', error '0.0000000'

Test #3:

score: 0
Accepted
time: 190ms
memory: 4060kb

input:

1
0 1

output:

-0.0000000000

result:

ok found '-0.0000000', expected '0.0000000', error '0.0000000'

Test #4:

score: 0
Accepted
time: 176ms
memory: 3992kb

input:

2
1 0
1 100

output:

-0.0000000000

result:

ok found '-0.0000000', expected '0.0000000', error '0.0000000'

Test #5:

score: 0
Accepted
time: 182ms
memory: 4044kb

input:

1
0.5 100

output:

0.0000000000

result:

ok found '0.0000000', expected '0.0000000', error '-0.0000000'

Test #6:

score: 0
Accepted
time: 207ms
memory: 3888kb

input:

3
0.4 100
0.6 100
0.6 100

output:

0.0000000000

result:

ok found '0.0000000', expected '0.0000000', error '-0.0000000'

Test #7:

score: 0
Accepted
time: 205ms
memory: 3952kb

input:

3
0.2 100
0.2 100
0.8 100

output:

50.0000000000

result:

ok found '50.0000000', expected '50.0000000', error '0.0000000'

Test #8:

score: 0
Accepted
time: 201ms
memory: 3948kb

input:

2
0.999999 1000000
0.999998 2

output:

0.9999990000

result:

ok found '0.9999990', expected '0.9999990', error '0.0000000'

Test #9:

score: 0
Accepted
time: 205ms
memory: 3768kb

input:

2
0 100000
0.000001 1

output:

-0.0000000000

result:

ok found '-0.0000000', expected '0.0000000', error '0.0000000'

Test #10:

score: 0
Accepted
time: 205ms
memory: 3924kb

input:

2
0 1000000000
0.000001 1

output:

1.0000000000

result:

ok found '1.0000000', expected '1.0000000', error '0.0000000'

Test #11:

score: -100
Wrong Answer
time: 1125ms
memory: 111824kb

input:

1000000
0.375733 595197307
0.505261 377150668
0.517039 15795246
0.448099 228176467
0.529380 871983979
0.905546 876268308
0.095891 272104456
0.500302 916153337
0.128705 355768079
0.070600 78747362
0.444107 466118868
0.194987 298494965
0.462293 593292779
0.287909 838058266
0.237226 934603199
0.391909 ...

output:

85717936612193.0287170410

result:

wrong answer 1st numbers differ - expected: '85718080941203.0156250', found: '85717936612193.0312500', error = '0.0000017'