数据好坑 ,说了最多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;
}