Time:2016.08.18
Author:xiaoyimi
转载注明出处谢谢
传送门
思路:
显然这个题是求这个图的最小割
在这个n*m个点,2*(n*(m-1)+(n-1)m+(n-1)(m-1))条边的图中直接跑最大流的结果是
那怎么办呢
关于平面图与对偶图,有一个比较牛逼的性质
如果让对偶图的每条边的权值对应平面图中的每条权值,那么对偶图中的最短路就是平面图中的最小割
啥,你不知道什么是平面图,什么是对偶图……
没事,之前的时候我也不知道……
总之这个性质非常牛逼就对了
你跑一遍SPFA最短路的时间总比你跑dinic/isap的时间要少的多了吧
注意:
建图略显蛋疼……
discuss里提示还要判断n=1/m=1的情况
代码:
#include<cstdio>
#include<iostream>
#include<cstring>
#include<queue>
#define M 1000002
using namespace std;
int n,m,s,t,cnt,tot;
int dis[M<<1],first[M<<1];
int a[1001][1001],b[1001][1001],c[1001][1001],id[1001][1001][2];
bool vis[M<<1];
queue<int> q;
struct edge{
int v,w,next;
}e[M*6];
int in()
{
int t=0;char ch=getchar();
while (ch<'0'||ch>'9') ch=getchar();
while (ch>='0'&&ch<='9') t=(t<<1)+(t<<3)+ch-48,ch=getchar();
return t;
}
void add(int x,int y,int z)
{
e[++tot]=(edge){y,z,first[x]};first[x]=tot;
e[++tot]=(edge){x,z,first[y]};first[y]=tot;
}
main()
{
n=in();m=in();
t=n*m+1;
for (int i=1;i<=n;++i)
for (int j=1;j<m;++j)
a[i][j]=in();
for (int i=1;i<n;++i)
for (int j=1;j<=m;++j)
b[i][j]=in();
if (n==1)
{
a[0][0]=1<<30;
for (int i=1;i<m;i++) a[0][0]=min(a[0][0],a[1][i]);
printf("%d\n",a[0][0]);
return 0;
}
if (m==1)
{
b[0][0]=1<<30;
for (int i=1;i<n;i++) b[0][0]=min(b[0][0],b[i][1]);
printf("%d\n",b[0][0]);
return 0;
}
for (int i=1;i<n;++i)
for (int j=1;j<m;++j)
c[i][j]=in(),
id[i][j][0]=++cnt,
id[i][j][1]=++cnt,
add(id[i][j][0],id[i][j][1],c[i][j]);
for (int i=1;i<=n-1;++i)
add(s,id[i][1][0],b[i][1]),
add(id[i][m-1][1],t,b[i][m]);
for (int i=1;i<n;++i)
for (int j=1;j<m-1;++j)
add(id[i][j][1],id[i][j+1][0],b[i][j+1]);
for (int i=1;i<=m-1;++i)
add(s,id[n-1][i][0],a[n][i]),
add(id[1][i][1],t,a[1][i]);
for (int i=1;i<n-1;++i)
for (int j=1;j<m;++j)
add(id[i][j][0],id[i+1][j][1],a[i+1][j]);
q.push(s);
memset(dis,63,sizeof(dis));
dis[s]=0;
for (;!q.empty();q.pop())
{
int x=q.front();
vis[x]=1;
for (int i=first[x];i;i=e[i].next)
if (e[i].w+dis[x]<dis[e[i].v])
{
dis[e[i].v]=dis[x]+e[i].w;
if (!vis[e[i].v])
vis[e[i].v]=1,
q.push(e[i].v);
}
vis[x]=0;
}
printf("%d\n",dis[t]);
}