P3324 [SDOI2015]星际战争二分答案+网络流

直接二分答案用网络流判断就行了

#include<cstdio>
#define ll long long
#include<cstring> 
#include<queue>
using namespace std;
const int M=10000,inf=0x7fffffff;
const ll INF=1ll<<60;
int S,T,n,m;
int nex[M*2],to[M*2],head[M],tot,a[M],b[M],f[60][60];
int dep[M];ll maxn,cap[M*2];
void add(int x,int y,ll z){
    nex[++tot]=head[x];cap[tot]=z;
    to[tot]=y;head[x]=tot;
}
void init(){
    memset(head,0,sizeof head);tot=1;
}
struct dinic{
    queue<int>q;
    int bfs(){memset(dep,0,sizeof dep);
        while(!q.empty())q.pop();
        q.push(S);dep[S]=1;
        while(!q.empty()){
            int x=q.front();q.pop();
            for(int i=head[x],tmp;i<=tot&&i>=2;i=nex[i]){

                if(!dep[tmp=to[i]]&&cap[i])
                dep[tmp]=dep[x]+1,q.push(tmp);

            }
        } 
        return dep[T];
    }
    ll dfs(int x,ll minf){
        if(x==T) return minf;ll f;
        for(int i=head[x],tmp;i<=tot&&i>=2;i=nex[i]){
            if(dep[tmp=to[i]]==dep[x]+1&&cap[i]&&(f=dfs(tmp,min(cap[i],minf)))){
                cap[i]-=f;
                cap[i^1]+=f;
                return f;
            }
        }
        return 0;
    }
    ll danic(){ll ans=0,f;
        while(bfs()) {f=1;
            while(f)ans+=(f=dfs(S,INF));
        }

        return ans;
    }
}A;
bool check(ll x){
    init();
    for(int i=1;i<=n;i++) add(m+i,T,1ll*a[i]*10000),add(T,m+i,0);
    for(int i=1;i<=m;i++) add(S,i,1ll*b[i]*x),add(i,S,0);
     for(int i=1;i<=m;i++)for(int j=1;j<=n;j++)
     if(f[i][j]) add(i,j+m,INF),add(j+m,i,0); 

    return A.danic()==maxn;

}
int main(){
    scanf("%d%d",&n,&m);S=0,T=n+m+1;init();
    for(int i=1;i<=n;i++) scanf("%d",&a[i]),add(m+i,T,a[i]*10000),add(T,m+i,0);
    for(int i=1;i<=m;i++) scanf("%d",&b[i]),add(S,i,INF),add(i,S,0);
    for(int i=1;i<=m;i++)for(int j=1;j<=n;j++) {scanf("%d",&f[i][j]);
        if(f[i][j]==1) add(i,j+m,INF),add(j+m,i,0);
    }
    maxn=A.danic();

    int l=1,r=100000000;
    while(l<=r){int mid=(l+r)>>1;
        if(check(mid)) r=mid-1;else l=mid+1;
    }
    printf("%.4lf",(double)l/10000.0);
}