bzoj 2406: 矩阵 （有上下界的网络流）

2406: 矩阵

Time Limit: 10 Sec   Memory Limit: 128 MB
Submit: 229   Solved: 90
[ Submit][ Status][ Discuss]

Description

Input

第一行两个数n、m，表示矩阵的大小。

接下来n行，每行m列，描述矩阵A。

最后一行两个数L，R。

Output

第一行，输出最小的答案；

Sample Input

2 2

0 1

2 1

0 1

Sample Output

1

HINT

对于100%的数据满足N,M<=200,0<=L<=R<=1000,0<=Aij<=1000

Source

[ Submit][ Status][ Discuss] 
HOME   Back

题解：有上下界的网络流

直接二分要求的最小答案x，然后用有源汇有上下界的可行流判断。

s->i 容量[sumi[i]-x,sumi[i]+x] sumi[i]表示第i行的总和

j->t 容量[sumj[j]-x,sumj[j]+x] sumj[j]表示第j行的总和

i->j 容量[l,r]

然后按照有源汇有上下界的可行流建模即可。

#include
#include
#include
#include
#include
#include
#define N 200000
using namespace std;
int in[N],out[N],point[N],v[N],remain[N],nxt[N];
int deep[N],num[N],cur[N],ls,rs,n,m,tot;
int a[203][203],sumi[N],sumj[N],last[N],ck[N];
void add(int x,int y,int z)
{
	if (!z) return;
	tot++; nxt[tot]=point[x]; point[x]=tot; v[tot]=y; remain[tot]=z;
	tot++; nxt[tot]=point[y]; point[y]=tot; v[tot]=x; remain[tot]=0;
	//cout< p; p.push(t);
	while (!p.empty()) {
		int now=p.front(); p.pop();
		for (int i=point[now];i!=-1;i=nxt[i])
		 if (deep[v[i]]==t&&remain[i^1])
		  deep[v[i]]=deep[now]+1,p.push(v[i]);
	}
}
void isap(int s,int t)
{
	bfs(s,t);
	int ans=0,now=s;
	for (int i=1;i<=t;i++) cur[i]=point[i];
	for (int i=1;i<=t;i++) num[deep[i]]++;
	while (deep[s]=0) add(i,tt,k);
		else add(ss,i,-k);
		ck[i]=tot;
	}
	add(t,s,1e9);
	isap(ss,tt);
	//for (int i=1;i<=n+m+1;i++)
	 //cout<