首页 > 代码库 > bzoj 1066 蜥蜴

bzoj 1066 蜥蜴

网络流。

建图:首先将每根柱子拆成两个点。

每根柱子的入点向出点连一条容量为柱子高度的边。

每根柱子的出点向可以到达的柱子的入点连一条容量为正无穷的边。

源点向每根初始有蜥蜴的柱子的入点连一条容量为一的边。

每根可以跳出地图的柱子的出点向汇点连一条容量为正无穷的边。

跑一遍最大流就是最多能逃出的蜥蜴数。

 1 #include<cstring>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<iostream>
 5 #include<cmath>
 6 #include<queue>
 7 using namespace std;
 8 const int dian=8005;
 9 const int bian=160005;
10 const int INF=0x3f3f3f3f;
11 int h[dian],nxt[bian],ver[bian],val[bian],ch[dian];
12 char mapp[25][25],map[25][25];
13 char lala[25];
14 int n,m,d,tot,ans;
15 int S,T;
16 void add(int aa,int bb,int cc){
17     tot++;ver[tot]=bb;nxt[tot]=h[aa];val[tot]=cc;h[aa]=tot;
18     tot++;ver[tot]=aa;nxt[tot]=h[bb];val[tot]=0;h[bb]=tot;
19 }
20 int bh(int aa,int bb,int cc){
21     return (aa-1)*m+bb+cc*n*m;
22 }
23 bool tell(){
24     memset(ch,-1,sizeof(ch));
25     queue<int>q;
26     q.push(S);
27     ch[S]=0;
28     while(!q.empty()){
29         int t=q.front();
30         q.pop();
31         for(int i=h[t];i;i=nxt[i])
32             if(ch[ver[i]]==-1&&val[i]){
33                 q.push(ver[i]);
34                 ch[ver[i]]=ch[t]+1;
35             }
36     }
37     return ch[T]!=-1;
38 }
39 int zeng(int a,int b){
40     if(a==T)
41         return b;
42     int r=0;
43     for(int i=h[a];i&&b>r;i=nxt[i])
44         if(ch[a]+1==ch[ver[i]]&&val[i]){
45             int t=zeng(ver[i],min(b-r,val[i]));
46             val[i]-=t,r+=t,val[i^1]+=t;
47         }
48     if(!r)
49         ch[a]=-1;
50     return r;
51 }
52 int dinic(){
53     int r=0,t;
54     while(tell())
55         while(t=zeng(S,INF))
56             r+=t;
57     return r;
58 }
59 int main(){
60     tot=1;
61     scanf("%d%d%d",&n,&m,&d);
62     S=2*n*m+1,T=2*n*m+2;
63     for(int i=1;i<=n;i++)
64         scanf("%s",mapp[i]+1);
65     for(int i=1;i<=n;i++)
66         for(int j=1;j<=m;j++)
67             map[i][j]=mapp[i][j]-0;
68     for(int i=1;i<=n;i++)
69         for(int j=1;j<=m;j++)
70             if(map[i][j]){
71                 add(bh(i,j,0),bh(i,j,1),map[i][j]);
72                 for(int k=i-d;k<=i+d;k++)
73                     for(int l=j-d;l<=j+d;l++)
74                         if(k>=1&&k<=n&&l>=1&&l<=m&&(k!=i||l!=j)&&abs(k-i)*abs(k-i)+abs(l-j)*abs(l-j)<=d*d&&map[k][l])
75                             add(bh(i,j,1),bh(k,l,0),INF);
76             }
77     for(int i=1;i<=n;i++)
78         for(int j=1;j<=m;j++)
79             if(min(min(i,n+1-i),min(j,m+1-j))<=d&&map[i][j])
80                 add(bh(i,j,1),T,INF);
81     for(int i=1;i<=n;i++){
82         scanf("%s",lala+1);
83         for(int j=1;j<=m;j++)
84             if(lala[j]==L){
85                 ans++;
86                 add(S,bh(i,j,0),1);
87             }
88     }
89     printf("%d",ans-dinic());
90     return 0;
91 }

注意:

数组大小要好好算。

读入均为字符串形式。

距离为欧几里得距离。

 

bzoj 1066 蜥蜴