Problem
Mean
求平面图最小割。
Analysis
数据范围太大导致直接跑最大流会TLE(然而据说只要姿势不太丑都能水过去)。
正确做法是平面图转对偶图,就可以把最小割问题转化成求最短路(Orz,太神了……)。
注意要特判 n=1
与 m=1
的情况。
参考资料:周冬2008年集训队论文《两极相通——浅析最大—最小定理在信息学竞赛中的应用》
Code
#include<cstdio>
const int N=1005,INF=~0U>>2;
int n,m,e,ed,cur,x,h=1,t,i,j,ans=INF,r[N][N],c[N][N],d[N][N],v[N*N*6],f[N*N*6],nxt[N*N*6],g[N*N*2],dis[N*N*2],q[N*N*2];
bool in[N*N*2];
void addedge(int x,int y,int w){
v[++ed]=y,f[ed]=w;
nxt[ed]=g[x];
g[x]=ed;
}
void add(int x,int w){
if(w>dis[x]) return;
dis[x]=w;
if(!in[x]){
in[x]=1;
if(w<dis[q[h]]){
if(--h<1) h+=e;
q[h]=x;
}else{
if(++t>e) t-=e;
q[t]=x;
}
}
}
int main(){
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)
for(j=1;j<m;j++) scanf("%d",&r[i][j]);
for(i=1;i<n;i++)
for(j=1;j<=m;j++) scanf("%d",&c[i][j]);
for(i=1;i<n;i++)
for(j=1;j<m;j++) scanf("%d",&d[i][j]);
if(n==1 && m==1){printf("0");return 0;}
if(n==1){
for(int i=1;i<m;i++) if(r[1][i]<ans) ans=r[1][i];
printf("%d",ans);
return 0;
}
if(m==1){
for(int i=1;i<n;i++) if(c[i][1]<ans) ans=c[i][1];
printf("%d",ans);
return 0;
}
e=(n-1)*(m-1)*2+1;
for(i=1;i<n;i++)
for(j=1;j<m;j++){
cur=(i-1)*(m-1)*2+j;
if(i==1) addedge(cur,e,r[i][j]),addedge(e,cur,r[i][j]);
if(i==n-1) addedge(cur+m-1,0,r[i+1][j]),addedge(0,cur+m-1,r[i+1][j]);
if(j==1) addedge(cur+m-1,0,c[i][j]),addedge(0,cur+m-1,c[i][j]);
if(j==m-1) addedge(cur,e,c[i][j+1]),addedge(e,cur,c[i][j+1]);
if(i!=n-1) addedge(cur+(m-1)*2,cur+m-1,r[i+1][j]),addedge(cur+m-1,cur+(m-1)*2,r[i+1][j]);
if(j!=m-1) addedge(cur,cur+m,c[i][j+1]),addedge(cur+m,cur,c[i][j+1]);
addedge(cur,cur+m-1,d[i][j]),addedge(cur+m-1,cur,d[i][j]);
}
for(i=0;i<=e;i++) dis[i]=INF;
add(0,0);
while(h!=t+1){
int tmp=h;
if(++h>e) h-=e;
for(i=g[x=q[tmp]],in[x]=0;i;i=nxt[i]) add(v[i],dis[x]+f[i]);
}
printf("%d",dis[e]);
return 0;
}