bzoj 1001 平面图转换求最大流（最小割）

数据好坑 ，说了最多1000*1000个点，head[]开了1e6+10，硬是不给过 _φ(❐_❐✧  ，改为2*就过了；

因为给的边的顺序问题，可以将此题转换为平面图的同构图，然后跑spfa；

所谓同构图就是把原来图中的面看为点，如图，1是源点，6是汇点，先将它们连起来，此时多了两个面（在同构图中多了两个点，分别为同构图的起点和终点），然后把同构图中相邻面的点连起来，权值为线割的那条边的权值，最后把同构图中6,7的连线去掉，求同构图中6到7的最短路便是原图中的最小割；

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <queue>
#define STD std::ios::sync_with_stdio(false);cin.tie(0);cout.tie(0)
using namespace std;
const int maxn=1e6+10;
const int inf=0x3f3f3f3f;
int vis[2*maxn],dis[2*maxn],head[6*maxn],q;
struct point{
	int val,next,to;
}pt[6*maxn];
void init()
{
	q=0;
	memset(head,-1,sizeof(head));
}
void add(int u,int v,int val)
{
	pt[q].next=head[u];
	pt[q].to=v;
	pt[q].val=val;
	head[u]=q++;
}
void spfa(int st)
{
	queue<int>q;
	while (!q.empty())q.pop();
	q.push(st);
	memset(vis,0,sizeof(vis));
	memset(dis,inf,sizeof(dis));
	vis[st]=1;
	dis[st]=0;
	while (!q.empty())
	{
		int u=q.front();
		q.pop();
		vis[u]=0;
		for (int i=head[u];i!=-1;i=pt[i].next)
		{
			int v=pt[i].to;
			if (dis[u]+pt[i].val<dis[v])
			{
				dis[v]=dis[u]+pt[i].val;
				if (!vis[v])
				{
					q.push(v);
					vis[v]=1;
				}
			}
		}
	}
}
int main()
{
	STD;
	int n,m,nu_h,v;
	while (scanf("%d%d",&n,&m)!=EOF)
	{
		if (n==1||m==1)
		{
			int minn=0x3f3f3f3f;
			for (int i=1;i<max(n,m);i++)
			{
				scanf("%d",&v);
				minn=min(v,minn);
			}
			printf("%d\n",minn);
			continue;
		}
		init();
		int r=(n-1)*(m-1)*2+1;
		int rr=r+1;
		for (int i=1;i<=n;i++)
		{
			
			if (i==1) nu_h=1;
			else if (i==n) {
			nu_h=(n-1)*(m-1)*2-(m-1)+1;
			}
			else { int hang=i*2-2; nu_h=1+(m-1)*(hang-1);}
			for (int j=1;j<=m-1;j++)
			{
				scanf("%d",&v);
				if (i==1)
				{
					add(nu_h,r,v);
					add(r,nu_h,v);
					//cout<<nu_h<<"***"<<r<<endl;
				}
				else if (i==n)
				{
					add(nu_h,rr,v);
					add(rr,nu_h,v);
					//cout<<nu_h<<"***"<<rr<<endl;
				}
				else
				{
					add(nu_h,nu_h+m-1,v);
					add(nu_h+m-1,nu_h,v);
					//cout<<nu_h<<"***"<<nu_h+m-1<<endl;
				}
				nu_h++;
			}
		}
		for (int i=1;i<=n-1;i++)
		{
			int hang=i*2;
			nu_h=1+(hang-1)*(m-1);
			//cout<<nu_h<<"!!!!"<<endl;
			for (int j=1;j<=m;j++)
			{
				scanf("%d",&v);
				if (j==1)
				{
					add(nu_h,rr,v);
					add(rr,nu_h,v);
					//cout<<nu_h<<"***"<<rr<<endl;
				}
				else if (j==m)
				{
					add(nu_h-1-(m-1),r,v);
					add(r,nu_h-1-(m-1),v);
					//cout<<nu_h-1-(m-1)<<"***"<<r<<endl;
				}
				else
				{
					add(nu_h,nu_h-m,v);
					add(nu_h-m,nu_h,v);
					//cout<<nu_h<<"***"<<nu_h-m<<endl;
				}
				nu_h++;
			}
		}
		for (int i=1;i<=n-1;i++)
		{
			int hang=i*2-1;
			nu_h=1+(hang-1)*(m-1);
			for (int j=1;j<=m-1;j++)
			{
				scanf("%d",&v);
				add(nu_h,nu_h+(m-1),v);
				add(nu_h+(m-1),nu_h,v);
				//cout<<nu_h<<"***"<<nu_h+(m-1)<<endl;
				nu_h++;
			}
		}
		spfa(r);
		printf("%d\n",dis[rr]);
	}
	return 0;	
}