bzoj1066

练了这么久，应该就成了裸的最大流了吧

不1Y对不起自己

const inf=1000007;
type node=record
       next,point,flow:longint;
     end;

var edge:array[0..200010] of node;
    a:array[0..30,0..30] of longint;
    p,x,y,numh,h,cur,pre,f:array[0..1010] of longint;
    len,ans,sum,t,tot,i,j,n,m,d:longint;
    c:char;

function min(a,b:longint):longint;
  begin
    if a>b then exit(b) else exit(a);
  end;

procedure add(x,y,f:longint);
  begin
    inc(len);
    edge[len].point:=y;
    edge[len].flow:=f;
    edge[len].next:=p[x];
    p[x]:=len;
  end;

procedure sap;
  var u,q,i,j,tmp,neck:longint;
  begin
    u:=0;
    numh[0]:=tot+1;
    while h[0]<tot+1 do
    begin
      if u=tot then
      begin
        i:=0;
        while i<>tot do
        begin
          j:=cur[i];
          dec(edge[j].flow);
          inc(edge[j xor 1].flow);
          i:=edge[j].point;
        end;
        u:=0;
        inc(ans);
      end;
      q:=-1;
      i:=p[u];
      while i<>-1 do
      begin
        j:=edge[i].point;
        if (edge[i].flow>0) and (h[u]=h[j]+1) then
        begin
          q:=i;
          break;
        end;
        i:=edge[i].next;
      end;
      if q<>-1 then
      begin
        cur[u]:=q;
        pre[j]:=u;
        u:=j;
      end
      else begin
        dec(numh[h[u]]);
        if numh[h[u]]=0 then exit;
        tmp:=tot+1;
        i:=p[u];
        while i<>-1 do
        begin
          j:=edge[i].point;
          if edge[i].flow>0 then tmp:=min(tmp,h[j]);
          i:=edge[i].next;
        end;
        h[u]:=tmp+1;
        inc(numh[h[u]]);
        if u<>0 then u:=pre[u];
      end;
    end;
  end;

begin
  readln(n,m,d);
  len:=-1;
  fillchar(p,sizeof(p),255);
  for i:=1 to n do
  begin
    for j:=1 to m do
    begin
      read(c);
      if c<>'0' then
      begin
        inc(t);
        f[t]:=ord(c)-48;
        x[t]:=i;
        y[t]:=j;
        a[i,j]:=t;
      end;
    end;
    readln;
  end;
  for i:=1 to n do
  begin
    for j:=1 to m do
    begin
      read(c);
      if c='L' then
      begin
        inc(sum);
        add(0,a[i,j],1);
        add(a[i,j],0,0);
      end;
    end;
    readln;
  end;
  tot:=t shl 1+1;
  for i:=1 to t do
  begin
    for j:=i+1 to t do
      if d*d>=sqr(x[i]-x[j])+sqr(y[i]-y[j]) then
      begin
        add(i+t,j,inf);
        add(j,i+t,0);
        add(j+t,i,inf);
        add(i,j+t,0);
      end;
    if (x[i]+d>n) or (x[i]-d<=0) or (y[i]+d>m) or (y[i]-d<=0) then
    begin
      add(i+t,tot,inf);
      add(tot,i+t,0);
    end;
    add(i,i+t,f[i]);
    add(i+t,i,0);
  end;
  sap;
  writeln(sum-ans);
end.