BZOJ 1066 [SCOI2007]蜥蜴

PS：变量名引起的血案。。

代名词：dinic,建图，简单网络流，点的流量限制。

step1:当点有容量时把点一分为二，cap为点的容量

step2:加入源点 汇点

step3:每个点 家一个点 我们可以这样加，mp[i][j]=++cnt;cnt++;

     最后一个cnt++ 的cnt是加入的点

step4:源点到每一个石头add(s,y,1);

         每个石头到一个汇点add(y`,t,inf);

       然后我们还要判断该石头能否到外面。

PS:建图小心变量。

#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <iostream>
#include <string>
#include <queue>
#include <string.h>
#include <cmath>

using namespace std;

#define N 200240
#define inf 0x3f3f3f

struct edge
{
  int next,v,c;
}e[N<<1];

int tot,head[N],d[N];
int q[N<<2];

int min(int a,int b)
{
   return a<b?a:b;
}

void add(int u,int v,int c)
{
 e[tot].v=v;e[tot].c=c;
 e[tot].next=head[u];head[u]=tot++;

 e[tot].v=u;e[tot].c=0;
 e[tot].next=head[v];head[v]=tot++;
}

int  bfs(int s,int t)
{
   memset(d,-1,sizeof(d));
   int f=0,r=0;
   q[r++]=s;
   d[s]=0;

   while (f<r)
   {
         int u=q[f++];
         if (u==t) return 1;
         for (int i=head[u];i!=-1;i=e[i].next)
         {
           int v=e[i].v;
           if (d[v]==-1&&e[i].c>0)
           {
             d[v]=d[u]+1;
             q[r++]=v;
           }
         }
   }
   return 0;
}

int dfs(int s,int t,int b)
{
  int r=0;
  if (s==t) return b;
  for (int i=head[s];i!=-1&&r<b;i=e[i].next)
  {
     int v=e[i].v;
     if (e[i].c>0&&d[v]==d[s]+1)
     {
       int x=min(e[i].c,b-r);
       x=dfs(v,t,x);
       r+=x;
       e[i].c-=x;
       e[i^1].c+=x;
    }
  }
  if (!r) d[s]-=2;
  return r;
}

int dinic(int s,int t)
{
  int ans=0,tmp;
  while (bfs(s,t))
  {
     while (tmp=dfs(s,t,inf)) ans+=tmp;
  }
  return ans;
}
void init()
{
    tot=0;
    memset(head,-1,sizeof(head));
}

char s[123][123];
char ss[123][123];
int mp[123][123];

int dis(int x,int y,int xx,int yy)
{
    return (xx-x)*(xx-x)+(yy-y)*(yy-y);
}


int main()
{
   int r,c,d;
   int cnt=0;

   scanf("%d%d%d",&r,&c,&d);

   for (int i=1;i<=r;i++){
   scanf("%s",s[i]+1);
   for (int j=1;j<=c;j++)
   if (s[i][j]!='0'){
   mp[i][j]=++cnt;
   cnt++;
   }
   }

   init();

   for (int i=1;i<=r;i++)
   scanf("%s",ss[i]+1);

   int h=0,t=cnt+1;

   for (int i=1;i<=r;i++)
   for (int j=1;j<=c;j++)
   if (s[i][j]>'0')
   add(mp[i][j],mp[i][j]+1,s[i][j]-'0');

   int ans=0;
   for (int i=1;i<=r;i++)
   for (int j=1;j<=c;j++){
   if (ss[i][j]=='L')
   {
       ans++;
       add(0,mp[i][j],1);
   }
   if (s[i][j]>'0')  {
   if (i<=d||(r-i)<d||j<=d||(c-j)<d) add(mp[i][j]+1,t,inf);
   }
   }

   for (int i=1;i<=r;i++)
   for (int j=1;j<=c;j++)
   {
       for (int p=1;p<=r;p++)
       for (int k=1;k<=c;k++)
       if (!(i==p&&j==k))
       {
           if (s[i][j]>'0'&&s[p][k]>'0'&&dis(i,j,p,k)<=d*d)
           add(mp[i][j]+1,mp[p][k],inf);
       }
   }
   printf("%d
",ans-dinic(h,t));
   return 0;
}