【BZOJ】【1001】 【BJOI2006】狼抓兔子

平面图最小割->对偶图最短路


  平面图最小割转对偶图最短路= =

  想到了就比较好写了……

  可能是我对区域的标号方式比较奇特?反正我没有特判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 }
View Code

1001: [BeiJing2006]狼抓兔子

Time Limit: 15 Sec  Memory Limit: 162 MB
Submit: 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

3 4
5 6 4
4 3 1
7 5 3
5 6 7 8
8 7 6 5
5 5 5
6 6 6

Sample Output

14

HINT

Source

[Submit][Status][Discuss]
posted @ 2015-04-08 17:33  Tunix  阅读(289)  评论(0编辑  收藏  举报