Codeforces 639E

显然当$c>0$时，分数最大化的方案是将题目按照$\frac{p_i}{t_i}$从大到小排序，$\frac{p_i}{t_i}$相同的可以任意交换。那么我们容易算出第$i$道题最早和最晚可能完成的时间$L_i$和$R_i$。
考虑对于两道题目$i$，$j$（$p_i<p_j$），合法的条件为$p_i(1-\frac{cLi}{T})\le p_j(1-\frac{c{R}_i}{T})$，也即$\frac{p_j{R}_j-p_i{L}_i}{p_j-p_i}\le \frac{T}{c}$，于是我们需要算出所有限制中最大的一个。我们将所有题目按$p_i$排序后，限制是一个斜率的式子，对于题目$j$，我们需要求出所有满足$p_i<p_j$的$i$中，点$(p_i,p_i{L}_i)$到$(p_j,p_j{R}_j)$的最大斜率，这个在下凸壳上三分即可。因为我们已经按$p_i$排了序，横坐标单调，直接维护凸壳就行了。
时间复杂度$\mathcal{O}(n\log n)$。

#include <bits/stdc++.h>
#define eps 1e-10
 
using namespace std;
 
typedef long double ldb;
typedef long long ll;
 
struct Point {
  ll x,y;
  Point() {}
  Point(ll a,ll b):x(a),y(b) {}
  Point operator - (Point b) {return Point(x-b.x,y-b.y);}
  bool operator < (const Point & b) const {return (x*b.y!=y*b.x)?x*b.y<y*b.x:y<b.y;}
};
 
inline ldb cross(Point x,Point y) {
  return (ldb)x.x*y.y-(ldb)x.y*y.x;
}
 
inline ldb rat(Point x,Point y) {
  return (ldb)(x.y-y.y)/(x.x-y.x);
}
 
Point st[150005];
int top;
 
void insert(Point x) {
  while (top>1&&cross(x-st[top],st[top]-st[top-1])>=0) top--;
  st[++top]=x;
}
 
ldb query(Point x) {
  int l=1,r=top;
  while (r-l>2) {
  	int m1=l+(r-l)/3,m2=r-(r-l)/3;
  	if (rat(x,st[m1])>rat(x,st[m2])) r=m2; else l=m1;
  }
  ldb ans=0;
  for(int i=l;i<=r;i++) ans=max(ans,rat(x,st[i]));
  return ans;
} 
 
Point p[150005];
ll lpos[150005],rpos[150005];
 
int id[150005];
 
bool cmp(int x,int y) {
  return p[x].y<p[y].y;
}
 
int main() {
  int n;
  scanf("%d",&n);
  for(int i=1;i<=n;i++) scanf("%lld",&p[i].y);
  for(int i=1;i<=n;i++) scanf("%lld",&p[i].x);
  sort(p+1,p+n+1);
  ll s=0;
  for(int i=1;i<=n;i++) {
  	if (i>1&&p[i].x*p[i-1].y==p[i].y*p[i-1].x) lpos[i]=lpos[i-1];
  	else lpos[i]=s;
  	s+=p[i].x;
  }
  for(int i=1;i<=n;i++) lpos[i]+=p[i].x;
  for(int i=n;i>0;i--) {
  	if (i<n&&p[i].x*p[i+1].y==p[i].y*p[i+1].x) rpos[i]=rpos[i+1];
  	else rpos[i]=s;
  	s-=p[i].x;
  }
  for(int i=1;i<=n;i++) id[i]=i;
  sort(id+1,id+n+1,cmp);
  ldb ans=eps;
  for(int i=1,j=0;i<=n;i=j+1) {
    while (j<n&&p[id[j+1]].y==p[id[i]].y) j++;
    for(int k=i;k<=j;k++) 
	  ans=max(ans,query(Point(p[id[k]].y,p[id[k]].y*rpos[id[k]])));
	for(int k=i;k<=j;k++)
	  insert(Point(p[id[k]].y,p[id[k]].y*lpos[id[k]]));
  }
  for(int i=1;i<=n;i++) s+=p[i].x;
  printf("%.10f\n",min(1.0,(double)(s/ans)));
  return 0;
}