BZOJ3232: 圈地游戏

比较经典的最小割
先二分V/C的最大值，$V/C=k -> V-kC ≥ 0$，将边权×k做一遍最小割，如果减完大于0就可行，因为要保留三位小数，所以可以将点权和边权都×1w后二分整数

---

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int maxr = 60;
const int maxn = 5100;
const int maxm = 110000;
struct edge
{
    int y,nex;
    ll c;
    edge(){}
    edge(int _y,ll _c,int _nex){y=_y;c=_c;nex=_nex;}
}a[maxm<<1]; int len,fir[maxn],fi[maxn];
int n,m,st,ed; ll sum;
ll s[maxr][maxr],e1[maxr][maxr],e2[maxr][maxr];
void ins(int x,int y,ll c)
{
    a[++len]=edge(y,c,fir[x]);fir[x]=len;
    a[++len]=edge(x,0,fir[y]);fir[y]=len;
}
int h[maxn];
queue<int>q;
bool bfs()
{
    memset(h,0,sizeof h); h[st]=1;
    q.push(st);
    while(!q.empty())
    {
        int x=q.front(); q.pop();
        for(int k=fir[x];k;k=a[k].nex)
        {
            int y=a[k].y;
            if(a[k].c&&!h[y]) h[y]=h[x]+1,q.push(y);
        }
    }
    return h[ed]>0;
}
ll dfs(int x,ll flow)
{
    if(x==ed) return flow;
    ll delta=0;
    for(int &k=fi[x];k;k=a[k].nex)
    {
        int y=a[k].y;
        if(h[y]==h[x]+1&&a[k].c)
        {
            ll minx=dfs(y,min(flow-delta,a[k].c));
            a[k].c-=minx; a[k^1].c+=minx;
            delta+=minx;
        }
        if(delta==flow) return delta;
    }
    if(!delta) h[x]=0;
    return delta;
}
ll Flow() 
{ 
    ll ret=0; 
    while(bfs()) 
    {
        for(int i=1;i<=ed;i++) fi[i]=fir[i];
        ret+=dfs(st,LLONG_MAX); 
    }
    return ret; 
}
bool judge(ll k)
{
    memset(fir,0,sizeof fir); len=1;
    ll r=sum*10000; 
    st=n*m+1,ed=st+1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++) 
            ins((i-1)*m+j,ed,s[i][j]*10000);
    for(int i=1;i<=m;i++)
    {
        ins(st,i,e1[0][i]*k);
        ins(st,(n-1)*m+i,e1[n][i]*k);
    }
    for(int i=1;i<n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            int tx=(i-1)*m+j,ty=tx+m;
            ins(tx,ty,e1[i][j]*k);
            ins(ty,tx,e1[i][j]*k);
        }
    }
    for(int i=1;i<=n;i++)
    {
        ins(st,(i-1)*m+1,e2[i][0]*k);
        ins(st,i*m,e2[i][m]*k);
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<m;j++)
        {
            int tx=(i-1)*m+j,ty=tx+1;
            ins(tx,ty,e2[i][j]*k);
            ins(ty,tx,e2[i][j]*k);
        }
    }
    r-=Flow();
    return r>0;
}
ll solve(ll l,ll r)
{
    while(l<=r)
    {
        ll mid=(l+r)>>1;
        if(judge(mid)) l=mid+1;
        else r=mid-1;
    }
    return l-1;
}

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++) scanf("%lld",&s[i][j]),sum+=s[i][j];
    for(int i=0;i<=n;i++)
        for(int j=1;j<=m;j++) scanf("%lld",&e1[i][j]);
    for(int i=1;i<=n;i++)
        for(int j=0;j<=m;j++) scanf("%lld",&e2[i][j]);
    printf("%.3lf\n",(double)(solve(1,1000000))/10000);

    return 0;
}