bzoj 1066 蜥蜴

最大流。

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

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

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

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

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

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

 #include<cstring>
 #include<cstdio>
 #include<algorithm>
 #include<iostream>
 #include<cmath>
 #include<queue>
 using namespace std;
 const int dian=;
 const int bian=;
 const int INF=0x3f3f3f3f;
 int h[dian],nxt[bian],ver[bian],val[bian],ch[dian];
 char mapp[][],map[][];
 char lala[];
 int n,m,d,tot,ans;
 int S,T;
 void add(int aa,int bb,int cc){
     tot++;ver[tot]=bb;nxt[tot]=h[aa];val[tot]=cc;h[aa]=tot;
     tot++;ver[tot]=aa;nxt[tot]=h[bb];val[tot]=;h[bb]=tot;
 }
 int bh(int aa,int bb,int cc){
     return (aa-)*m+bb+cc*n*m;
 }
 bool tell(){
     memset(ch,-,sizeof(ch));
     queue<int>q;
     q.push(S);
     ch[S]=;
     while(!q.empty()){
         int t=q.front();
         q.pop();
         for(int i=h[t];i;i=nxt[i])
             if(ch[ver[i]]==-&&val[i]){
                 q.push(ver[i]);
                 ch[ver[i]]=ch[t]+;
             }
     }
     return ch[T]!=-;
 }
 int zeng(int a,int b){
     if(a==T)
         return b;
     int r=;
     for(int i=h[a];i&&b>r;i=nxt[i])
         if(ch[a]+==ch[ver[i]]&&val[i]){
             int t=zeng(ver[i],min(b-r,val[i]));
             val[i]-=t,r+=t,val[i^]+=t;
         }
     if(!r)
         ch[a]=-;
     return r;
 }
 int dinic(){
     int r=,t;
     while(tell())
         while(t=zeng(S,INF))
             r+=t;
     return r;
 }
 int main(){
     tot=;
     scanf("%d%d%d",&n,&m,&d);
     S=*n*m+,T=*n*m+;
     for(int i=;i<=n;i++)
         scanf("%s",mapp[i]+);
     for(int i=;i<=n;i++)
         for(int j=;j<=m;j++)
             map[i][j]=mapp[i][j]-'';
     for(int i=;i<=n;i++)
         for(int j=;j<=m;j++)
             if(map[i][j]){
                 add(bh(i,j,),bh(i,j,),map[i][j]);
                 for(int k=i-d;k<=i+d;k++)
                     for(int l=j-d;l<=j+d;l++)
                         if(k>=&&k<=n&&l>=&&l<=m&&(k!=i||l!=j)&&abs(k-i)*abs(k-i)+abs(l-j)*abs(l-j)<=d*d&&map[k][l])
                             add(bh(i,j,),bh(k,l,),INF);
             }
     for(int i=;i<=n;i++)
         for(int j=;j<=m;j++)
             if(min(min(i,n+-i),min(j,m+-j))<=d&&map[i][j])
                 add(bh(i,j,),T,INF);
     for(int i=;i<=n;i++){
         scanf("%s",lala+);
         for(int j=;j<=m;j++)
             if(lala[j]=='L'){
                 ans++;
                 add(S,bh(i,j,),);
             }
     }
     printf("%d",ans-dinic());
     return ;
 }

注意：

数组大小要好好算。

读入均为字符串形式。

距离为欧几里得距离。