平面图最小割转对偶图最短路= =
想到了就比较好写了……
可能是我对区域的标号方式比较奇特？反正我没有特判n==1||m==1也能过2333（机智吧～（滚开啦你个自恋狂！））
想当年我刚学网络流的时候第一次上B站就是写了这道题TAT顺利挂掉
时空复杂度比较丑，偷懒只写了spfa
1 /************************************************************** 2 Problem: 1001 3 User: Tunix 4 Language: C++ 5 Result: Accepted 6 Time:3012 ms 7 Memory:121388 kb 8 ****************************************************************/ 9 10 //BZOJ 1001 11 #include<vector> 12 #include<cstdio> 13 #include<cstring> 14 #include<cstdlib> 15 #include<iostream> 16 #include<algorithm> 17 #define rep(i,n) for(int i=0;i<n;++i) 18 #define F(i,j,n) for(int i=j;i<=n;++i) 19 #define D(i,j,n) for(int i=j;i>=n;--i) 20 #define pb push_back 21 using namespace std; 22 inline int getint(){ 23 int v=0,sign=1; char ch=getchar(); 24 while(ch<'0'||ch>'9'){ if (ch=='-') sign=-1; ch=getchar();} 25 while(ch>='0'&&ch<='9'){ v=v*10+ch-'0'; ch=getchar();} 26 return v*sign; 27 } 28 const int N=3e6+10,M=6e6+10,INF=~0u>>2; 29 typedef long long LL; 30 /******************tamplate*********************/ 31 int head[N],next[M],to[M],len[M],cnt; 32 void ins(int x,int y,int z){ 33 to[++cnt]=y; next[cnt]=head[x]; head[x]=cnt; len[cnt]=z; 34 } 35 void add(int x,int y,int z){ins(x,y,z);ins(y,x,z);} 36 /********************edge***********************/ 37 int n,m; 38 inline int pack(int i,int j){return i*(m<<1)+j;} 39 int Q[M],d[N]; 40 bool inq[N]; 41 int spfa(int S,int T){ 42 int l=0,r=0; 43 Q[++r]=S; 44 memset(d,0x7f,sizeof d); 45 memset(inq,0,sizeof inq); 46 d[S]=0; 47 while(l!=r){ 48 int x=Q[l++]; 49 if (l==M) l=0; 50 inq[x]=0; 51 for(int i=head[x];i;i=next[i]) 52 if (d[to[i]]>d[x]+len[i]){ 53 d[to[i]]=d[x]+len[i]; 54 if (!inq[to[i]]){ 55 r++; if(r==M) r=0; 56 Q[r]=to[i]; 57 inq[to[i]]=1; 58 } 59 } 60 } 61 return d[T]; 62 } 63 int main(){ 64 #ifndef ONLINE_JUDGE 65 freopen("1001.in","r",stdin); 66 freopen("1001.out","w",stdout); 67 #endif 68 n=getint(); m=getint(); 69 int x; 70 F(i,1,n) F(j,1,m-1){ 71 x=getint(); 72 add(pack(i-1,j<<1),pack(i,j<<1|1),x); 73 } 74 F(i,1,n-1) F(j,1,m){ 75 x=getint(); 76 add(pack(i,(j<<1)-1),pack(i,j<<1),x); 77 } 78 F(i,1,n-1) F(j,1,m-1){ 79 x=getint(); 80 add(pack(i,j<<1),pack(i,j<<1|1),x); 81 } 82 // if (n==1) {solve1(); return 0;} 83 int S=0,T=pack(n,m<<1); 84 F(j,1,m-1){ 85 add(S,pack(0,j<<1),0); 86 add(pack(n,j<<1|1),T,0); 87 } 88 F(i,1,n-1){ 89 add(pack(i,1),T,0); 90 add(S,pack(i,m<<1),0); 91 } 92 int ans=spfa(S,T); 93 // printf("%d\n",ans); 94 cout<<ans<<endl; 95 return 0; 96 }