QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#354022#8137. 'Ello, and What Are You After, Then?SolitaryDreamWA 823ms4384kbC++172.4kb2024-03-14 20:39:272024-03-14 20:39:27

Judging History

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

  • [2024-03-14 20:39:27]
  • 评测
  • 测评结果:WA
  • 用时:823ms
  • 内存:4384kb
  • [2024-03-14 20:39:27]
  • 提交

answer

#include <bits/stdc++.h>
using namespace std;
const int N = 1e3 + 10;
int B, C, S;
int n, m[N];
struct Task {
    int f, t;
    double e;
};
vector<Task> npc[N];

inline double Calc(int o, double E, double k) {
    vector<double> A(m[o]), B(m[o]);
    vector<double> mx;
    for (int i = 0; i < m[o]; ++i) {
        auto [f, t, e] = npc[o][i];
        A[i] = f * (e - E) * t - C * k * f;
        B[i] = S * k * f;
        mx.push_back(max(A[i], B[i]));
    }
    int cnt = max(m[o] - ::B, 0);
    sort(mx.begin(), mx.end(), greater<double>());
    double ans = 0;
    for(int i = 0; i < cnt; i += 1)
        ans += mx[i];
    for(int i = cnt; i < m[o]; i += 1)
        if(mx[i] > 0)
            ans += mx[i];
    return ans;
}
inline bool Check(double E) {
    double maxL = -1e3, minR = 0;
    for (int i = 1; i <= n; ++i) {
        // npc i
        double L = -1e3, R = 0;
        while((R - L) / max(1., R) > 1e-6) {
            double m1 = L + (R - L) / 3;
            double m2 = R - (R - L) / 3;
            if (Calc(i, E, m1) > Calc(i, E, m2)) L = m1;
            else R = m2;
        }
        double p = L;
        // cout << L << ' ' << R << ' ' << Calc(i, E, p) << endl;
        // exit(0);
        if (Calc(i, E, p) > 0) return 1;
        L = maxL, R = p;
        while((R - L) / max(1., -L) > 1e-6) {
            double mid = (L + R) * 0.5;
            if (Calc(i, E, mid) >= 0) L = mid;
            else R = mid;
        }
        double pl = L;
        L = p, R = minR;
        while ((R - L) / max(1., -L) > 1e-6) {
            double mid = (L + R) * 0.5;
            if (Calc(i, E, mid) >= 0) R = mid;
            else L = mid;
        }
        double pr = L;
        if (pr < maxL || pl > minR) return 1;
        maxL = max(maxL, pl);
        minR = min(minR, pr);
    }
    return 0;
}
int main() {
    ios::sync_with_stdio(0);
    cin.tie(0);
    cin >> B >> C >> S;
    cin >> n;
    for (int i = 1; i <= n; ++i) {
        cin >> m[i];
        npc[i].resize(m[i]);
        for (auto &[f, t, e] : npc[i]) cin >> f >> t >> e;
        sort(npc[i].begin(), npc[i].end(), [](Task u, Task v) { return u.f < v.f; });
    }
    double L = 0, R = 1e4;
    while((R - L) / max(1., R) > 1e-6) {
        double mid = (L + R) * 0.5;
        if (Check(mid)) L = mid;
        else R = mid;
    }
    cout << fixed << setprecision(15);
    cout << L << endl;
    return 0;
}

详细

Test #1:

score: 100
Accepted
time: 1ms
memory: 3880kb

input:

0 1 6
2
1
1 1 1
2
1 10 1
1 10 10

output:

6.999997422099113

result:

ok found '6.9999974', expected '7.0000000', error '0.0000004'

Test #2:

score: 0
Accepted
time: 1ms
memory: 3832kb

input:

2 1 2
1
4
10 2 1
10 1 1
1 10 1
1 1 10

output:

5.909092724323273

result:

ok found '5.9090927', expected '5.9090909', error '0.0000003'

Test #3:

score: -100
Wrong Answer
time: 823ms
memory: 4384kb

input:

14 1000 1000
1000
30
113 80 1188
92 145 1074
130 56 1296
139 102 1142
60 76 1317
128 126 1208
73 120 1155
91 89 1197
115 64 979
80 118 592
110 97 556
83 105 578
94 51 848
98 134 757
107 138 1038
105 143 892
92 72 893
88 103 961
87 148 879
105 84 823
85 134 607
100 82 1084
199 58 801
138 85 743
214 1...

output:

1453.366279602050781

result:

wrong answer 1st numbers differ - expected: '1453.3645790', found: '1453.3662796', error = '0.0000012'