BZOJ 1001 狼抓兔子 平面图最小割

题意:n*m地图 每个点可以走到其右边,下方,右下方.这3道路都有各自的流量限制.
源点s为左上角,终点t为右下角,
n,m<=1000,问最少需要多少个人堵住道路 使得s-t无路径存在.

堵住割上的边 则s-t无路径存在.花费最小也就是求最小割/最大流 n,m<=1e3 TLE...

本图是平面图,平面图的最小割等于其新图的最短路.

新图的点为原图的面.(建新图时对面合理的编号..)

新图边的权值等于:原图中分割这两个面的边的权值. 

#include <bits/stdc++.h>
using namespace std;
typedef pair<int,int> ii;
const int inf=0x3f3f3f3f;
const int N=4e6+20;
vector<ii> e[N];
int s,t;
int n,m;
void add(int u,int v,int w)
{
	e[u].push_back(ii(v,w));
	e[v].push_back(ii(u,w));
}
bool inq[N];
int dis[N];
queue<int> q;
void SPFA(int s,int t)
{
	memset(dis,inf,sizeof(dis));
	memset(inq,0,sizeof(inq));
	dis[s]=0,inq[s]=1,q.push(s);
	while(!q.empty())
	{
		int u=q.front();
		q.pop(),inq[u]=0;
		for(int i=0;i<e[u].size();i++)
		{
			int v=e[u][i].first,w=e[u][i].second;
			if(dis[v]>dis[u]+w)
			{
				dis[v]=dis[u]+w;
				if(!inq[v])
					inq[v]=1,q.push(v);
			}
		}
	}
	cout<<dis[t]<<endl;
}
int main()
{
	int k;
	scanf("%d%d",&n,&m);
	s=0,t=(n-1)*(m-1)*2+1;
	//ho 
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m-1;j++)
		{
			scanf("%d",&k);
			if(i==1)	add(j*2,t,k);
			else if(i==n)	add(s,(n-2)*(m-1)*2+j*2-1,k);
			else add((i-2)*(m-1)*2+j*2-1,(i-1)*(m-1)*2+j*2,k);
		//	int u=(i-2)*(m-1)*2+j*2-1,v=(i-1)*(m-1)*2+j*2;
		//	cout<<i<<' '<<j<<' '<<u<<' '<<v<<endl;
		}
	}
	//ver
	for(int i=1;i<=n-1;i++)
	{
		for(int j=1;j<=m;j++)
		{
			scanf("%d",&k);
            if (j==1) add((i-1)*(m-1)*2+1,0,k);  
            else if (j==m) add(t,i*(m-1)*2,k);  
            else add((i-1)*(m-1)*2+(j-1)*2,(i-1)*(m-1)*2+(j-1)*2+1,k);  
		}
	}
	       
    for(int i=1;i<=n-1;i++)
    {  
   		for(int j=1;j<=m-1;j++)
	 	{	          
    	    scanf("%d",&k);  
            add((i-1)*(m-1)*2+j*2-1,(i-1)*(m-1)*2+j*2,k);  
    	}
    }  
	SPFA(s,t);
	return 0;
}