分数规划。。
二分答案 建图 跑费用流看看最后答案是不是>=0就好了
为什么我的代码好像挺慢的QAQ
#include<bits/stdc++.h>
#define me(a,x) memset(a,x,sizeof a)
using namespace std;
const int N=210,M=21000;
const double inf=1e11;
inline int read(){
char ch=getchar(); int x=0,f=1;
while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0'; ch=getchar();}
return x*f;
}
struct E{int x,y,c;double d;int next;}a[M]; int len,fir[N];
void ins(int x,int y,int c,double d){
a[++len]=(E){x,y,c,d,fir[x]},fir[x]=len;
a[++len]=(E){y,x,0,-d,fir[y]},fir[y]=len;
}
double d[N]; bool v[N]; int l[N],t;
queue<int>q;
bool bfs(){
for(int i=0;i<=t;i++)l[i]=v[i]=0,d[i]=-inf;
q.push(0); d[0]=0,v[0]=1;
while(!q.empty()){
int x=q.front(); q.pop(),v[x]=0;
for(int k=fir[x];k;k=a[k].next){
int y=a[k].y;
if(d[y]<d[x]+a[k].d && a[k].c){
l[y]=k; d[y]=d[x]+a[k].d;
if(!v[y]){q.push(y); v[y]=1;}
}
}
}
return d[t]!=-inf;
}
int b[102][102],p[102][102],n;
double flow(double u){
len=1; me(fir,0);
for(int i=1;i<=n;i++)ins(0,i,1,0),ins(i+n,t,1,0);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)ins(i,j+n,1,b[i][j]-p[i][j]*u);
double s=0;
while(bfs()){
s+=d[t];
for(int k=l[t];k;k=l[a[k].x])
a[k].c--,a[k^1].c++;
}
return s;
}
int main()
{
n=read(),t=n+n+1;
for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)b[i][j]=read();
for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)p[i][j]=read();
double l=0,r=1e6,ans;
while(r-l>=1e-8){
double mid=(l+r)/2;
if(flow(mid)>=0)ans=mid,l=mid+1e-8;
else r=mid-1e-8;
}
printf("%.6lf\n",ans);
return 0;
}