#include <bits/stdc++.h>
#pragma GCC optimize ("Ofast", "unroll-loops")
using namespace std;
#define ll long long
#define int ll
#define FOR(i,a,b) for (int i = (a); i<(b); i++)
#define REP(i,n) FOR(i,0,n)
#define REP1(i,n) FOR(i,1,n+1)
#define RREP(i,n) for (int i=(n)-1; i>=0; i--)
#define f first
#define s second
#define pb push_back
#define ALL(x) x.begin(),x.end()
#define SZ(x) (int)(x.size())
#define SQ(x) (x)*(x)
#define pii pair<int, int>
#define pdd pair<double ,double>
#define pcc pair<char, char>
#define endl '\n'
//#define TOAD
#ifdef TOAD
#define bug(x) cerr<<__LINE__<<": "<<#x<<" is "<<x<<endl
#define IOS()
#else
#define bug(...)
#define IOS() ios::sync_with_stdio(0), cin.tie(0), cout.tie(0)
#endif
const ll inf = 1ll<<60;
const int iinf=2147483647;
const ll mod = 1e9+7;
const ll maxn=505;
const double PI=acos(-1);
ll pw(ll x, ll p, ll m=mod){
ll ret=1;
while (p>0){
if (p&1){
ret*=x;
ret%=m;
}
x*=x;
x%=m;
p>>=1;
}
return ret;
}
ll inv(ll a, ll m=mod){
return pw(a,m-2);
}
//=======================================================================================
int n, k;
vector<pii> vc;
bool cmp(pii a, pii b){
if (a.s==-1&&b.s==-1){
return a.f<b.f;
}
if (a.s==-1||b.s==-1){
return a.s>b.s;
}
return a.s<b.s;
}
long double dp[maxn][maxn], sf[maxn][maxn];
signed main (){
IOS();
cout<<fixed<<setprecision(5);
cin>>n>>k;
REP(i, n){
int x, y; cin>>x>>y;
vc.pb({x, y});
}
sort(ALL(vc), cmp);
RREP(i, n){
vector<int> cur;
FOR(j, i, n){
cur.pb(vc[j].f);
}
sort(ALL(cur));
REP(j, k+1){
if (j>SZ(cur)) break;
if (j==0) sf[i][j]=0.0;
else sf[i][j]=sf[i][j-1]+cur[j-1];
}
}
long double res=sf[0][k];
REP1(t, k){
REP(i, n){
REP(j, t+1){
dp[i][j]=inf;
}
}
REP(i, n){
if (i==0){
dp[i][0]=(0.0+vc[i].f)/(t+1);
if (vc[i].s!=-1) dp[i][1]=(0.0+vc[i].s);
continue;
}
REP(j, t+1){
if (j>i+1) break;
if (j==0){
dp[i][0]=dp[i-1][0]+(0.0+vc[i].f)/(t+1);
}
else dp[i][j]=min(dp[i-1][j]+(0.0+vc[i].f)/(t+1), ((vc[i].s!=-1)?(dp[i-1][j-1]+(0.0+vc[i].s)/(j)):inf));
}
}
REP(i, n){
if (i+1<t){
continue;
}
if (n-i-1<k-t){
continue;
}
res=min(res, dp[i][t]+sf[i+1][k-t]/(t+1));
}
//res=min(res, dp[t-1][t]+(0.0+sf[t][k-t])/(t+1)); // 0 base
}
cout<<res<<endl;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
0 ms |
340 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
0 ms |
340 KB |
Output is correct |
5 |
Correct |
22 ms |
5092 KB |
Output is correct |
6 |
Correct |
101 ms |
6220 KB |
Output is correct |
7 |
Correct |
342 ms |
7756 KB |
Output is correct |
8 |
Correct |
714 ms |
7760 KB |
Output is correct |
9 |
Correct |
1139 ms |
7804 KB |
Output is correct |
10 |
Correct |
628 ms |
7732 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
0 ms |
340 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
0 ms |
340 KB |
Output is correct |
5 |
Correct |
22 ms |
5092 KB |
Output is correct |
6 |
Correct |
101 ms |
6220 KB |
Output is correct |
7 |
Correct |
342 ms |
7756 KB |
Output is correct |
8 |
Correct |
714 ms |
7760 KB |
Output is correct |
9 |
Correct |
1139 ms |
7804 KB |
Output is correct |
10 |
Correct |
628 ms |
7732 KB |
Output is correct |
11 |
Correct |
0 ms |
212 KB |
Output is correct |
12 |
Correct |
155 ms |
6464 KB |
Output is correct |
13 |
Correct |
150 ms |
6428 KB |
Output is correct |
14 |
Correct |
146 ms |
6372 KB |
Output is correct |
15 |
Correct |
736 ms |
7724 KB |
Output is correct |
16 |
Correct |
665 ms |
7760 KB |
Output is correct |
17 |
Correct |
633 ms |
7768 KB |
Output is correct |
18 |
Correct |
1302 ms |
7808 KB |
Output is correct |
19 |
Correct |
1219 ms |
7772 KB |
Output is correct |
20 |
Correct |
1138 ms |
7764 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
0 ms |
340 KB |
Output is correct |
3 |
Correct |
0 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
0 ms |
340 KB |
Output is correct |
6 |
Correct |
0 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
340 KB |
Output is correct |
8 |
Correct |
0 ms |
340 KB |
Output is correct |
9 |
Incorrect |
0 ms |
340 KB |
Output isn't correct |
10 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
0 ms |
340 KB |
Output is correct |
3 |
Correct |
0 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
0 ms |
340 KB |
Output is correct |
6 |
Correct |
0 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
340 KB |
Output is correct |
8 |
Correct |
0 ms |
340 KB |
Output is correct |
9 |
Incorrect |
0 ms |
340 KB |
Output isn't correct |
10 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
0 ms |
340 KB |
Output is correct |
3 |
Correct |
0 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
0 ms |
340 KB |
Output is correct |
6 |
Correct |
0 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
340 KB |
Output is correct |
8 |
Correct |
0 ms |
340 KB |
Output is correct |
9 |
Incorrect |
0 ms |
340 KB |
Output isn't correct |
10 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
1353 ms |
7764 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
0 ms |
340 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
0 ms |
340 KB |
Output is correct |
5 |
Correct |
22 ms |
5092 KB |
Output is correct |
6 |
Correct |
101 ms |
6220 KB |
Output is correct |
7 |
Correct |
342 ms |
7756 KB |
Output is correct |
8 |
Correct |
714 ms |
7760 KB |
Output is correct |
9 |
Correct |
1139 ms |
7804 KB |
Output is correct |
10 |
Correct |
628 ms |
7732 KB |
Output is correct |
11 |
Correct |
0 ms |
212 KB |
Output is correct |
12 |
Correct |
155 ms |
6464 KB |
Output is correct |
13 |
Correct |
150 ms |
6428 KB |
Output is correct |
14 |
Correct |
146 ms |
6372 KB |
Output is correct |
15 |
Correct |
736 ms |
7724 KB |
Output is correct |
16 |
Correct |
665 ms |
7760 KB |
Output is correct |
17 |
Correct |
633 ms |
7768 KB |
Output is correct |
18 |
Correct |
1302 ms |
7808 KB |
Output is correct |
19 |
Correct |
1219 ms |
7772 KB |
Output is correct |
20 |
Correct |
1138 ms |
7764 KB |
Output is correct |
21 |
Correct |
1 ms |
340 KB |
Output is correct |
22 |
Correct |
0 ms |
340 KB |
Output is correct |
23 |
Correct |
0 ms |
340 KB |
Output is correct |
24 |
Correct |
1 ms |
340 KB |
Output is correct |
25 |
Correct |
0 ms |
340 KB |
Output is correct |
26 |
Correct |
0 ms |
340 KB |
Output is correct |
27 |
Correct |
1 ms |
340 KB |
Output is correct |
28 |
Correct |
0 ms |
340 KB |
Output is correct |
29 |
Incorrect |
0 ms |
340 KB |
Output isn't correct |
30 |
Halted |
0 ms |
0 KB |
- |