可以看出来是最小割,然后你就去求最大流了
这么大的范围就是让你用网络流卡的?咋想的啊=。=???
建议还是老老实实用 平面图最小割等于其对偶图最短路 这个东西来做吧,虽然这个东西跑的也挺慢的,最后一个点跑了$2s$
对偶图就是被边分割出来的每个区域当成一个点,然后两个区域有公共边就连边,起点和终点的问题就在源汇点中间连一条边然后就能分出来了
#include<queue>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=,M=;
struct a{int node,dist;};
bool operator < (a x,a y)
{
return x.dist>y.dist;
}
priority_queue<a> hp;
int p[N],noww[*M],goal[*M],val[*M],dis[N],vis[N];
int n,m,rd,st,ed,bs,t1,t2,cnt;
void Read(int &x)
{
x=; char ch=getchar();
while(!isdigit(ch))
ch=getchar();
while(isdigit(ch))
x=(x<<)+(x<<)+(ch^),ch=getchar();
}
void Link(int f,int t,int v)
{
noww[++cnt]=p[f],p[f]=cnt;
goal[cnt]=t,val[cnt]=v;
}
int ID(int a,int b,int c)
{
if(a>n||!b) return st;
if(!a||b>m) return ed;
return (a-)*m+b+c*n*m;
}
void Dijkstra()
{
memset(dis,0x3f,sizeof dis);
dis[st]=,hp.push((a){st,});
while(!hp.empty())
{
a tt=hp.top(); hp.pop(); int tn=tt.node;
if(vis[tn]) continue; vis[tn]=true;
for(int i=p[tn];i;i=noww[i])
if(dis[goal[i]]>dis[tn]+val[i])
{
dis[goal[i]]=dis[tn]+val[i];
hp.push((a){goal[i],dis[goal[i]]});
}
}
}
int main()
{
register int i,j;
Read(n),Read(m);
n--,m--,st=*n*m+,ed=st+;
for(i=;i<=n+;i++)
for(j=;j<=m;j++)
{
Read(rd),t1=ID(i,j,),t2=ID(i-,j,);
Link(t1,t2,rd),Link(t2,t1,rd);
}
for(i=;i<=n;i++)
for(j=;j<=m+;j++)
{
Read(rd),t1=ID(i,j,),t2=ID(i,j-,);
Link(t1,t2,rd),Link(t2,t1,rd);
}
for(i=;i<=n;i++)
for(j=;j<=m;j++)
{
Read(rd),t1=ID(i,j,),t2=ID(i,j,);
Link(t1,t2,rd),Link(t2,t1,rd);
}
Dijkstra(); printf("%d",dis[ed]);
return ;
}