bzoj1001[狼爪兔子]最小割转最短路

我们都知道这个定理，最小割=最大流。

所以，一看到最小割的问题果断转化成最大流写，然后这道题就T了- -

然后就不知所措了，看了网上的题解才知道最小割可以转成最短路- -

 

bzoj1001:http://www.lydsy.com/JudgeOnline/problem.php?id=1001

 

题意就是用最少的狼封闭所有道路(即切断起点和终点)，典型的最小割模型。

推荐参考周东的<<两级相通--->>的论文

http://wenku.baidu.com/link?url=87F10nBWauMdSF-PaKHoG-3fZj0jFE63P6pHSeX6ZiguQqXOQxm41iLWW5IdZCp2MWFQ8JghamfeI68PtLqEv_JSWapGp5z415gNoYb031u

 

简单点讲，每个平面图都对应一个对偶图，对偶图中的最小环就是原图的最小割，如果删去对偶图中s-t这条边，就是相当于求最短路了！

如下图，蓝色的是原图，红色的是对偶图，绿色的是对偶图最短路（最小环删去了s-t的边）（也是原图最小割）。

 

 

把原图中每个点在对偶图中标号，重新建图，在新图中跑最短路就行了，效率比网络流快。

 

实在不想写dijistra,写spfa写惯了，就附上spfa的代码：

 

(PS:我的建图可能有点奇葩，参照了如右图片)：

 

至于前几次交上去WA了，是因为这种建图有一个Bug，就是n=1的时候，横向建了2次，纵向差不多

所以后来加上了特判！

 

 

#include<cstdio>
#include<iostream>
#include<queue>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdlib>
using namespace std;
const int inf=0x3f3f3f3f;
int n,m;
int head[2000000+20];
int dis[2000000+20];
int inq[2000000+20];
struct node
{
	int v,w,next;
}e[6000000+20];
int k;
int S,T;
void add(int u,int v,int w)
{
	e[k].v=v;
	e[k].w=w;
	e[k].next=head[u];
	head[u]=k++;
	
	e[k].v=u;
	e[k].w=w;
	e[k].next=head[v];
	head[v]=k++;
}
void spfa()
{
	memset(dis,inf,sizeof(dis));
	memset(inq,0,sizeof(inq));
	queue<int>q;
	dis[S]=0;
	inq[S]=1;
	q.push(S);
	while(!q.empty())
	{
		int u=q.front();
		q.pop();
		inq[u]=0;
		for(int i=head[u];~i;i=e[i].next)
		{
			int v=e[i].v;
			if(dis[u]+e[i].w<dis[v])
			{
				dis[v]=dis[u]+e[i].w;
				if(!inq[v])
				{
					inq[v]=1;
					q.push(v);
				}
			}
		}
	}
	printf("%d\n",dis[T]);
}
int main()
{
	memset(head,-1,sizeof(head));
	k=0;
	scanf("%d%d",&n,&m);
	if (n == 1 || m == 1)
    {
        if (n > m) swap(n, m);
        int ans = inf;
        for (int i = 1; i < m; ++i)
        {
        	int x;
            scanf("%d",&x);
            if (x < ans) ans = x;
        }
        printf("%d\n", ans);
        exit(0);
    }
	S=0,T=2*(m-1)*(n-1)+1;
	int w;
	for(int j=1;j<m;j++)
	{
		int id=j*2-1;
		scanf("%d",&w);
		add(S,id,w);
	}
	for(int i=1;i<n-1;i++)
	{
		for(int j=1;j<m;j++)
		{
			int t1=2*((i-1)*(m-1)+j);
			int t2=2*(i*(m-1)+j)-1;
			scanf("%d",&w);
			add(t1,t2,w);
		}
	}
	for(int j=1;j<m;j++)
	{
		int id=2*((n-2)*(m-1)+j);
		scanf("%d",&w);
		add(id,T,w);
	}
	//横向边
	
	for(int i=0;i<n-1;i++)
	{
		for(int j=1;j<=m;j++)
		{
			scanf("%d",&w);
			if(j==1)//与终点相连 
			{
				int id=2*(i*(m-1)+j);
				add(id,T,w);
			}
			else if(j==m)
			{
				int id=2*((i+1)*(m-1))-1;
				add(S,id,w); 
			}
			else
			{
				int t1=2*(i*(m-1)+j-1)-1;
				int t2=2*(i*(m-1)+j);
				add(t1,t2,w);
			}
		}
	} 
	//纵向
	
	for(int i=0;i<n-1;i++)
	{
		for(int j=1;j<m;j++)
		{
			scanf("%d",&w);
			int t1=2*(i*(m-1)+j)-1;
			int t2=2*(i*(m-1)+j);
			add(t1,t2,w);
		}
	} 
	//斜向
	
	spfa(); 
	return 0;
}