元のタイトルポータル
01スコア計画
∑ a ∑ b <= mid \ frac {\ sum a} {\ sum b} <= midΣbΣA<=m i d
∑ a <= mid ∗ ∑ b \ sum a <= mid * \ sum bΣa<=M I D∗Σb
wi = ai − mid ∗ bi w_i = a_i-mid * b_iw私=a私−M I D∗b私
∑ wi <= 0 \ sum w_i <= 0があるかどうか尋ねますΣw私<=0
aaによるとaは大から小にグループ化され、同じaaグループは、続くにBBグループB昇順
このような小さいこと満たす厳しい条件よりも
できDPI、j個のdp_ {I、Jd pI 、J前者を表すiiグループ jでjjが選択されていますj最初のタスク
、次にバックパックを作る
O(n 3 log)O(n ^ 3log) O (n3 log)
コード:
#include <bits/stdc++.h>
#define maxn 55
#define LL long long
using namespace std;
struct data{
LL a, b;
}a[maxn];
LL dp[maxn][maxn];
vector <int> object[maxn];
int n, m;
inline int read(){
int s = 0, w = 1;
char c = getchar();
for (; !isdigit(c); c = getchar()) if (c == '-') w = -1;
for (; isdigit(c); c= getchar()) s = (s << 1) + (s << 3) + (c ^ 48);
return s * w;
}
bool cmp(data x, data y){
return x.a == y.a ? x.b > y.b : x.a > y.a; }
bool check(LL mid){
memset(dp, 0x3f, sizeof(dp));
dp[0][0] = 0;
for (int i = 1, sum = 0; i <= m; ++i){
sum += object[i].size();
for (int j = 0; j <= sum; ++j){
LL x = 0;
for (int k = 0; k <= object[i].size() && k <= j && j - k >= sum - j; ++k){
if (k) x += 1LL * 1000 * a[object[i][k - 1]].a - mid * a[object[i][k - 1]].b;
dp[i][j] = min(dp[i][j], dp[i - 1][j - k] + x);
}
}
}
for (int i = 0; i <= n; ++i)
if (dp[m][i] <= 0) return 1;
return 0;
}
int main(){
n = read();
for (int i = 1; i <= n; ++i) a[i].a = read();
for (int i = 1; i <= n; ++i) a[i].b = read();
sort(a + 1, a + 1 + n, cmp);
a[0].a = a[1].a - 1;
for (int i = 1; i <= n; ++i){
if (a[i].a != a[i - 1].a) ++m;
object[m].push_back(i);
}
LL l = 0, r = 1e11, ans;
while (l <= r){
LL mid = (l + r) >> 1;
if (check(mid)) ans = mid, r = mid - 1; else l = mid + 1;
}
printf("%lld\n", ans);
return 0;
}