平面图最小割->对偶图最短路
比较经典的一道平面图最小割转对偶图最短路的题。
推荐这个
我的对偶图标号方法:
注意特判n==1||m==1的情况
代码:
#include<queue>
#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1005
#define K N*N<<1
#define M N*N*3<<1
using namespace std;
struct edge{ int next,to,dis; }ed[M];
int n,m,k,h[K],s,t,dis[K];
bool f[K];
queue <int> que;//卡了十几发空间都没过,就直接搞STL了
//读优不知道哪里出了问题,就没打了
inline void addedge(int x,int y,int z){
ed[++k].next=h[x],ed[k].to=y,ed[k].dis=z,h[x]=k;
}
inline int spf(){//spfa
memset(dis,0x3f,sizeof(dis));
memset(f,false,sizeof(f));
while (!que.empty()) que.pop();
que.push(s),dis[s]=0;
while (!que.empty()){
int x=que.front(); que.pop(); f[x]=false;
for (int i=h[x],v;i;i=ed[i].next)
if (dis[x]+ed[i].dis<dis[v=ed[i].to]){
dis[v]=dis[x]+ed[i].dis;
if (!f[v]) que.push(v),f[v]=true;
}
}
return dis[t];
}
inline int id(int x,int y){//对偶图的标号
return (x-1)*(m-1)+y<<1;
}
inline void calc(){
int x,ans=0x7fffffff;
while (~scanf("%d",&x)) ans=min(ans,x);
printf("%d\n",ans);
}
int main(){
scanf("%d%d",&n,&m),s=0,t=((n-1)*(m-1)<<1)+1;
if (n==1||m==1) return calc(),0;//一定要特判
for (int i=1;i<=n;i++)
for (int j=1;j<m;j++){
int x,u=id(i,j),v=id(i-1,j)-1; scanf("%d",&x);
u=min(u,t),v=max(v,s),addedge(u,v,x),addedge(v,u,x);
}
for (int i=1;i<n;i++)
for (int j=1;j<=m;j++){
int x,u=id(i,j)-1,v=u-1; scanf("%d",&x);
if (j==1) v=t; if (j==m) u=s;
addedge(u,v,x),addedge(v,u,x);
}
for (int i=1;i<n;i++)
for (int j=1;j<m;j++){
int x,u=id(i,j),v=u-1; scanf("%d",&x);
addedge(u,v,x),addedge(v,u,x);
}
return printf("%d\n",spf()),0;
}