BZOJ 4819: [Sdoi2017]新生舞会

分数规划。。
二分答案 建图 跑费用流看看最后答案是不是>=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;
}