bzoj1066: [SCOI2007]蜥蜴（最大流）

1066: [SCOI2007]蜥蜴

题目：传送门

---

题解：

　　 哇QTT大佬一眼秒算法...ORT

　　 其实很容易就可以看出来是一道最大流

　　 因为有边的使用限制，那么就可以直接当成是流量来处理嘛

　　 因为是对点进行操作，那么就可以拆点啊

　　 一开始现将有柱子的点自己把限制条件连上，就是对点x拆成x1和x2那么就x1-->x2流量为限制

　　 然后就是无脑乱连，对于判断是否出了边界一开始还傻逼逼的没想出来。。。

　　 枚举多一圈界外的点，那么如果算曼哈顿距离（据说欧几里得也OK）符合条件，就判断是不是在界内：

　　 在界内就直接连嘛，对于点x和点y，x2-->y1流量无限

　　 在界外那就很nice，直接向ed连无限流量。

　　 最后就是对于每个蜥蜴的操作了，直接st连蜥蜴流量为1（因为柱子上只可能有一只）

　　

---

代码：

 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<cmath>
 #include<algorithm>
 #define inf 999999999
 using namespace std;
 struct node
 {
     int x,y,c,next,other;
 }a[];int len,last[];
 void ins(int x,int y,int c)
 {
     int k1,k2;
     k1=++len;
     a[len].x=x;a[len].y=y;a[len].c=c;
     a[len].next=last[x];last[x]=len;

     k2=++len;
     a[len].x=y;a[len].y=x;a[len].c=;
     a[len].next=last[y];last[y]=len;

     a[k1].other=k2;
     a[k2].other=k1;
 }
 int n,m,D,head,tail,st,ed;
 int list[],h[];
 bool bt_h()
 {
     memset(h,,sizeof(h));h[st]=;
     list[]=st;head=;tail=;
     while(head!=tail)
     {
         int x=list[head];
         for(int k=last[x];k;k=a[k].next)
         {
             int y=a[k].y;
             if(h[y]== && a[k].c>)
             {
                 h[y]=h[x]+;
                 list[tail++]=y;
             }
         }
         head++;
     }
     if(h[ed]>)return true;
     return false;
 }
 int find_flow(int x,int flow)
 {
     if(x==ed)return flow;
     int s=,t;
     for(int k=last[x];k;k=a[k].next)
     {
         int y=a[k].y;
         if(h[y]==h[x]+ && a[k].c> && s<flow)
         {
             s+=t=find_flow(y,min(a[k].c,flow-s));
             a[k].c-=t;a[a[k].other].c+=t;
         }
     }
     if(s==)h[x]=;
     return s;
 }
 char A[][],B[][];
 int d[][];
 int main()
 {
     scanf("%d%d%d",&n,&m,&D);int ss=;
     for(int i=;i<=n;i++)for(int j=;j<=m;j++)d[i][j]=++ss;
     for(int i=;i<=n;i++)scanf("%s",A[i]+);
     for(int i=;i<=n;i++)scanf("%s",B[i]+);
     len=;memset(last,,sizeof(last));
     st=*n*m+,ed=st+;
     for(int i=;i<=n;i++)for(int j=;j<=m;j++)if(A[i][j]!='')ins(d[i][j],d[i][j]+n*m,A[i][j]-'');
     for(int i=;i<=n;i++)
         for(int j=;j<=m;j++)
             if(A[i][j]!='')
                 for(int tx=;tx<=n+;tx++)
                     for(int ty=;ty<=m+;ty++)
                         if(abs(i-tx)+abs(j-ty)<=D)
                             if(tx>= && tx<=n && ty>= && ty<=m)ins(d[i][j]+n*m,d[tx][ty],inf);
                             else ins(d[i][j]+n*m,ed,inf);
     int sum=;
     for(int i=;i<=n;i++)
         for(int j=;j<=m;j++)
             if(B[i][j]=='L')
                 ins(st,d[i][j],),sum++;
     int ans=;
     while(bt_h())ans+=find_flow(st,inf);
     printf("%d\n",sum-ans);
     return ;
 }