bzoj2406 矩阵 二分+有源汇上下界网络流

题意略。

%%ymw大神，考场上直接切。
先二分答案。
然后连边，先预处理一下横行纵行的前缀和，然后b的每行每列取值范围就是r[i]-mid,r[i]+mid，列同理。
那么我们把横行和纵行看做点，拆成两部分形成一个二分图，横行纵行之间连边，容量为(r[i]-mid,r[i]+mid)，列同理。
然后跑可行流来判断。
如果要求输出答案的话提前记录一下就好了。

#include<cstdio>#include<algorithm>#include<cstring>#include<map>#define fo(i,a,b) for(int i=a;i<=b;i++)#define fd(i,a,b) for(int i=a;i>=b;i--)using namespace std;const int N=4e5+5;int n,m;int head[N],go[N],next[N],cur[N],val[N];int tot=1;int num[300][300];const int inf=0x3f3f3f3f;int cnt,s,t,S,T,L,R;int dis[N],q[N],du[N];int row[N],col[N];inline void clear(){    tot=1;    memset(head,0,sizeof(head));    memset(du,0,sizeof(du));    s=0,t=n+m+1,S=n+m+2,T=n+m+3;}inline void ins(int x,int y,int z){    go[++tot]=y;    next[tot]=head[x];    val[tot]=z;    head[x]=tot;}inline void add(int x,int y,int z){    ins(x,y,z);    ins(y,x,0);}inline bool bfs(){    int t=0,w=1;    memset(dis,-1,sizeof(dis));    q[1]=S,dis[S]=0;    while (t<w)    {        int x=q[++t];        for(int i=head[x];i;i=next[i])        {            int v=go[i];            if (val[i]&&dis[v]==-1)            {                dis[v]=dis[x]+1;                q[++w]=v;                if (v==T)return 1;            }        }    }    return 0;}inline int dfs(int x,int f){    if (x==T||!f)return f;    int w,used=0;    for(int i=cur[x];i;i=next[i])    {        int v=go[i];        if (dis[v]==dis[x]+1)        {            w=f-used;            w=dfs(v,min(w,val[i]));            val[i]-=w;            val[i^1]+=w;            used+=w;            if (used==f)break;         }    }    return used;}inline int dinic(){    int ans=0;    while (bfs())    {        fo(i,0,T)cur[i]=head[i];        ans+=dfs(S,inf);    }    return ans;}bool pd(int mid){    clear();    add(t,s,inf);    fo(i,1,n)    fo(j,1,m)    add(i,j+n,R-L),du[i]-=L,du[j+n]+=L,num[i][j]=tot;    fo(i,1,n)    {        int p=max(0,row[i]-mid),q=mid+row[i];        add(s,i,q-p);        du[s]-=p,du[i]+=p;    }    fo(i,1,m)    {        int p=max(0,col[i]-mid),q=mid+col[i];        add(i+n,t,q-p);        du[i+n]-=p,du[t]+=p;    }    int cnt=0;    fo(i,s,t)    if (du[i]>0)add(S,i,du[i]),cnt+=du[i];    else if (du[i]<0)add(i,T,-du[i]);    return dinic()==cnt;}int main(){    scanf("%d%d",&n,&m);    fo(i,1,n)    {        fo(j,1,m)        {            int x;            scanf("%d",&x);            row[i]+=x,col[j]+=x;        }    }    scanf("%d%d",&L,&R);    int ans=-1;    int l=0,r=200000;    while (l<=r)    {        int mid=(l+r)>>1;        if (pd(mid))r=mid-1,ans=mid;        else l=mid+1;    }    printf("%d\n",ans);}