#define _CRT_SECURE_NO_WARNINGS
#include <bits/stdc++.h>
//#include "temp.cpp"
#include <cstdio>
using namespace std;
#ifndef ONLINE_JUDGE
#define dbg(x) cerr << #x <<" "; print(x); cerr << endl;
#else
#define dbg(x)
#endif
#define sz(x) (int((x).size()))
#define len(x) (int)x.length()
#define all(x) (x).begin(), (x).end()
#define rall(x) (x).rbegin(), (x).rend()
#define clr(x) (x).clear()
#define uniq(x) x.resize(unique(all(x)) - x.begin());
#define blt(x) __builtin_popcount(x)
#define pb push_back
#define popf pop_front
#define popb pop_back
void print(long long t) {cerr << t;}
void print(int t) {cerr << t;}
void print(string t) {cerr << t;}
void print(char t) {cerr << t;}
void print(double t) {cerr << t;}
void print(long double t) {cerr << t;}
void print(unsigned long long t) {cerr << t;}
template <class T, class V> void print(pair <T, V> p);
template <class T> void print(vector <T> v);
template <class T> void print(set <T> v);
template <class T, class V> void print(map <T, V> v);
template <class T> void print(multiset <T> v);
template <class T, class V> void print(T v[],V n) {cerr << "["; for(int i = 0; i < n; i++) {print(v[i]); cerr << " "; } cerr << "]";}
template <class T, class V> void print(pair <T, V> p) {cerr << "{"; print(p.first); cerr << ","; print(p.second); cerr << "}";}
template <class T> void print(vector <T> v) {cerr << "[ "; for (T i : v) {print(i); cerr << " ";} cerr << "]";}
template <class T> void print(deque <T> v) {cerr << "[ "; for (T i : v) {print(i); cerr << " ";} cerr << "]";}
template <class T> void print(set <T> v) {cerr << "[ "; for (T i : v) {print(i); cerr << " ";} cerr << "]";}
template <class T> void print(multiset <T> v) {cerr << "[ "; for (T i : v) {print(i); cerr << " ";} cerr << "]";}
template <class T, class V> void print(map <T, V> v) {cerr << "[ "; for (auto i : v) {print(i); cerr << " ";} cerr << "]";}
#include <ext/pb_ds/assoc_container.hpp>
using namespace __gnu_pbds;
#define nl '\n'
// for grid problems
int dx[8] = {-1,0,1,0,1,-1,1,-1};
int dy[8] = {0,1,0,-1,1,1,-1,-1};
// lowest / (1 << 17) >= 1e5 / (1 << 18) >= 2e5 / (1 << 21) >= 1e6
void fastIO() {
ios_base::sync_with_stdio(false);
cin.tie(nullptr); cout.tie(nullptr);
}
// file in/out
void setIO(string str = "") {
fastIO();
if (str != "") {
freopen((str + ".in").c_str(), "r", stdin);
freopen((str + ".out").c_str(), "w", stdout);
}
}
// Indexed Set
template <class T> using Tree = tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>;
const long double inf = 1e15;
const int N = 1e2 + 10;
long double dp[5 * N][5 * N][5 * N];
int n, k;
void solve_() {
cin >> n >> k;
vector<pair<long double, long double>> vec;
for(int i = 1; i <= n; i++) {
long double a, b; cin >> a >> b;
vec.emplace_back(b, a);
}
sort(all(vec));
if(n <= 100 || k == n) {
long double pi = inf;
for(int qq = 0; qq <= n; qq++) {
long double ci = qq + 1;
for(int i = 0; i <= n; i++) {
for(int j = 0; j <= n; j++) {
if(k == n) {
dp[i][j][i - j] = inf;
} else {
for(int e = 0; e <= n; e++) {
dp[i][j][e] = inf;
}
}
}
}
dp[0][0][0] = 0;
for(int i = 0; i < n; i++) {
for(int j = 0; j <= qq; j++) {
if(k == n) {
int e = i - j;
if(j + 1 <= qq && vec[i].first != -1) {
long double qn = j + 1;
dp[i + 1][j + 1][e] = min(dp[i + 1][j + 1][e], dp[i][j][e] + vec[i].first / qn);
}
dp[i + 1][j][e + 1] = min(dp[i + 1][j][e + 1], dp[i][j][e] + vec[i].second / ci);
dp[i + 1][j][e] = min(dp[i + 1][j][e], dp[i][j][e]);
} else {
for(int e = 0; e <= n && e + j <= i; e++) {
if(j + 1 <= qq && vec[i].first != -1) {
long double qn = j + 1;
dp[i + 1][j + 1][e] = min(dp[i + 1][j + 1][e], dp[i][j][e] + vec[i].first / qn);
}
dp[i + 1][j][e + 1] = min(dp[i + 1][j][e + 1], dp[i][j][e] + vec[i].second / ci);
dp[i + 1][j][e] = min(dp[i + 1][j][e], dp[i][j][e]);
}
}
}
}
pi = min(pi, dp[n][qq][max(0, k - qq)]);
}
cout << setprecision(15) << pi << '\n';
return;
}
long double p = inf;
for(int i = 0; i <= n; i++) {
long double cnt = 1, ans = 0;
vector<long double> vi;
int cc = 0;
for(int j = 0; j < i; j++) {
if(vec[j].first == -1) {
vi.push_back(vec[j].second);
continue;
}
cc++;
ans += vec[j].first / cnt;
cnt += 1;
}
for(int j = i; j < n; j++) {
vi.push_back(vec[j].second);
}
sort(all(vi));
int it = 0, mm = k - cc;
while(mm > 0) {
ans += vi[it] / cnt;
it++, mm--;
}
p = min(p, ans);
}
cout << setprecision(15) << p << '\n';
}
int main() {
setIO();
auto solve = [&](int test_case)-> void {
while(test_case--) {
solve_();
}
};
int test_cases = 1;
// cin >> test_cases;
solve(test_cases);
return 0;
}
Compilation message
Main.cpp: In function 'void setIO(std::string)':
Main.cpp:64:12: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
64 | freopen((str + ".in").c_str(), "r", stdin);
| ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Main.cpp:65:12: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
65 | freopen((str + ".out").c_str(), "w", stdout);
| ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
/tmp/ccrwesdu.o: in function `print(long long)':
Main.cpp:(.text+0x7da): relocation truncated to fit: R_X86_64_PC32 against symbol `std::cerr' defined in .bss._ZSt4cerr section in /usr/lib/gcc/x86_64-linux-gnu/10/libstdc++.a(globals_io.o)
/tmp/ccrwesdu.o: in function `print(int)':
Main.cpp:(.text+0x7f9): relocation truncated to fit: R_X86_64_PC32 against symbol `std::cerr' defined in .bss._ZSt4cerr section in /usr/lib/gcc/x86_64-linux-gnu/10/libstdc++.a(globals_io.o)
/tmp/ccrwesdu.o: in function `print(std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> >)':
Main.cpp:(.text+0x81e): relocation truncated to fit: R_X86_64_PC32 against symbol `std::cerr' defined in .bss._ZSt4cerr section in /usr/lib/gcc/x86_64-linux-gnu/10/libstdc++.a(globals_io.o)
/tmp/ccrwesdu.o: in function `print(char)':
Main.cpp:(.text+0x85a): relocation truncated to fit: R_X86_64_PC32 against symbol `std::cerr' defined in .bss._ZSt4cerr section in /usr/lib/gcc/x86_64-linux-gnu/10/libstdc++.a(globals_io.o)
/tmp/ccrwesdu.o: in function `print(double)':
Main.cpp:(.text+0x887): relocation truncated to fit: R_X86_64_PC32 against symbol `std::cerr' defined in .bss._ZSt4cerr section in /usr/lib/gcc/x86_64-linux-gnu/10/libstdc++.a(globals_io.o)
/tmp/ccrwesdu.o: in function `print(long double)':
Main.cpp:(.text+0x897): relocation truncated to fit: R_X86_64_PC32 against symbol `std::cerr' defined in .bss._ZSt4cerr section in /usr/lib/gcc/x86_64-linux-gnu/10/libstdc++.a(globals_io.o)
/tmp/ccrwesdu.o: in function `print(unsigned long long)':
Main.cpp:(.text+0x8aa): relocation truncated to fit: R_X86_64_PC32 against symbol `std::cerr' defined in .bss._ZSt4cerr section in /usr/lib/gcc/x86_64-linux-gnu/10/libstdc++.a(globals_io.o)
/tmp/ccrwesdu.o: in function `fastIO()':
Main.cpp:(.text+0x8d2): relocation truncated to fit: R_X86_64_PC32 against symbol `std::cin' defined in .bss._ZSt3cin section in /usr/lib/gcc/x86_64-linux-gnu/10/libstdc++.a(globals_io.o)
Main.cpp:(.text+0x8dd): relocation truncated to fit: R_X86_64_PC32 against symbol `std::cout' defined in .bss._ZSt4cout section in /usr/lib/gcc/x86_64-linux-gnu/10/libstdc++.a(globals_io.o)
/tmp/ccrwesdu.o: in function `setIO(std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> >)':
Main.cpp:(.text+0x925): relocation truncated to fit: R_X86_64_PC32 against symbol `std::cin' defined in .bss._ZSt3cin section in /usr/lib/gcc/x86_64-linux-gnu/10/libstdc++.a(globals_io.o)
Main.cpp:(.text+0x930): additional relocation overflows omitted from the output
/usr/bin/ld: failed to convert GOTPCREL relocation; relink with --no-relax
collect2: error: ld returned 1 exit status
