題意：給n個任務 每個任務有兩個值$a,b$ 現有許多機器 每臺最多可以執行兩次任務 若存在第二次任務則滿足$a_{second}<a_{first}$ 定義代價$val = \frac { \sum_{i \in S } a[i]} { \sum_{i \in S} b[i] }$ 其中$S$為當做第一次來執行的任務的集合 求$val$的最小值
$n \leq 60,a_i \leq 10^8,b_i \leq 100$

很容易想到二分最小值並且聯想到經典的$01$分數規劃 即
$\frac {\sum_{i=1}^{n}a[i]}{\sum_{i=1}^{n}b[i]} \leq x$
${\sum_{i=1}^{n}a[i]} \leq x \times {\sum_{i=1}^{n}b[i]}$

順著這個思路 我們來設計$dp$方程
先將任務按$a$降序排序
$f[i][j][k]$表示做到第$i$個有$j$個大於$a_i$的任務是當做第一次來做且沒有被安排第二次有$k$個等於$a_i$的任務是當做第一次來做且沒有被安排第二次
轉移式為
$f[i][j][k]=min(f[i+1][j-1][k]),f[i+1][j+1][k]+a_i-x \times b_i (a_i==a_{i+1})$

有個坑點是題目要求輸出$ceil(ans \times 1000)$ 我用$google$翻譯出來卻是把$ans \times 1000$四舍五入 以後還是讀英文題面好了

#include<bits/stdc++.h>
using namespace std;
#define mod 1000000007
#define ll long long
#define cl(x) memset(x,0,sizeof x)
 

#ifdef Poi
#define bug(x) cout<<(#x)<<" "<<(x)<<endl
#define debug(...) fprintf(stderr, __VA_ARGS__)
#else
#define bug(x)
#define debug(...)
#endif
const int INF = 0x7fffffff;
const int N=55;
inline int read(){
int x=0,rev=0,ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')rev=1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return rev?-x:x;
}
struct data{
double a,b;
}t[N];
double X,f[N][N][N];
bool vis[N][N][N];
int n;
ll ans;

bool cmp(data i,data j){
return i.a>j.a;
}
double dfs(int pos,int d,int g){
if(pos==n) return 0;
if(vis[pos][d][g]) return f[pos][d][g];
vis[pos][d][g]=1;
double mn=INF;
if(pos<n-1&&t[pos].a==t[pos+1].a) { 
if(d) mn=min(mn,dfs(pos+1,d-1,g));
mn=min(mn,dfs(pos+1,d,g+1)+t[pos].a-X*t[pos].b);
}
else if(t[pos].a!=t[pos+1].a||pos==n-1){
if(d) mn=min(mn,dfs(pos+1,d+g-1,0));
mn=min(mn,dfs(pos+1,d+g+1,0)+t[pos].a-X*t[pos].b);  
}
return f[pos][d][g]=mn;
}
bool judge(double k){
X=k,cl(vis),cl(f);
return dfs(0,0,0)<=0;
}

int main(){
#ifdef Poi 
freopen("in.txt","r",stdin);
#endif
n=read();
for(int i=0;i<n;i++) t[i].a=read();
//  bug(t[n-1].a);
for(int i=0;i<n;i++) t[i].b=read();
if(t[0].a==99999991) {
//  for(int i=0;i<n;i++) cout<<t[i].b<<endl;
}
sort(t,t+n,cmp);
double l=0,r=1e8;
for(int T=1;T<=100;T++){
double mid=(l+r)/2.0;
//  bug(mid);
if(judge(mid)) r=mid;
else l=mid;
}
//  bug(l);
ans=(ceil)(l*1000);
//  printf("%.6lf\n",l);
cout<<ans<<endl;
}

Codeforces 994F Compute Power 二分+DP