平面图最小割->对偶图最短路
平面图最小割转对偶图最短路= =
想到了就比较好写了……
可能是我对区域的标号方式比较奇特?反正我没有特判n==1||m==1也能过2333(机智吧~(滚开啦你个自恋狂!))
想当年我刚学网络流的时候第一次上B站就是写了这道题TAT顺利挂掉
时空复杂度比较丑,偷懒只写了spfa
![](https://i-blog.csdnimg.cn/blog_migrate/8f900a89c6347c561fdf2122f13be562.gif)
![](https://i-blog.csdnimg.cn/blog_migrate/961ddebeb323a10fe0623af514929fc1.gif)
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 }
1001: [BeiJing2006]狼抓兔子
Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 11872 Solved: 2773
[Submit][Status][Discuss]
Description
现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓兔子还是比较在行的,而且现在的兔子还比较笨,它们只有两个窝,现在你做为狼王,面对下面这样一个网格的地形:
左上角点为(1,1),右下角点为(N,M)(上图中N=4,M=5).有以下 三种类型的道路 1:(x,y)<==>(x+1,y) 2:(x,y)<==>(x,y+1) 3:(x,y)<==>(x+1,y+1) 道路上的权值表示这条路上最多能够通过的兔子数,道路是无向的. 左上角和右下角为兔子的两个窝,开始时所有的兔子都聚集在左上角(1,1)的窝里,现在它们要跑到右下解(N,M)的窝中去,狼王开始伏击这些兔子.当然 为了保险起见,如果一条道路上最多通过的兔子数为K,狼王需要安排同样数量的K只狼,才能完全封锁这条道路,你需要帮助狼王安排一个伏击方案,使得在将兔 子一网打尽的前提下,参与的狼的数量要最小。因为狼还要去找喜羊羊麻烦.
Input
第 一行为N,M.表示网格的大小,N,M均小于等于1000.接下来分三部分第一部分共N行,每行M-1个数,表示横向道路的权值. 第二部分共N-1行,每行M个数,表示纵向道路的权值. 第三部分共N-1行,每行M-1个数,表示斜向道路的权值. 输入文件保证不超过10M
Output
输出一个整数,表示参与伏击的狼的最小数量.
Sample Input
5 6 4
4 3 1
7 5 3
5 6 7 8
8 7 6 5
5 5 5
6 6 6