【2011集训队出题】圈地计划

题目

最近房地产商GDOI(Group of Dumbbells Or Idiots)从NOI(Nuts Old Idiots)手中得到了一块开发土地。据了解，这块土地是一块矩形的区域，可以纵横划分为N×M块小区域。GDOI要求将这些区域分为商业区和工业区来开发。根据不同的地形环境，每块小区域建造商业区和工业区能取得不同的经济价值。更具体点，对于第i行第j列的区域，建造商业区将得到Aij收益，建造工业区将得到Bij收益。另外不同的区域连在一起可以得到额外的收益，即如果区域(I,j)相邻（相邻是指两个格子有公共边）有K块（显然K不超过4）类型不同于(I,j)的区域，则这块区域能增加k×Cij收益。经过Tiger.S教授的勘察，收益矩阵A,B,C都已经知道了。你能帮GDOI求出一个收益最大的方案么？

分析

我们可以很明显的看到约束条件，二元关系，和它们选和不选的利润，所以自然而然想到了最小割。
那么怎么割呢，那么多条件。
我们再分析，我么可以把这些点分成两组：
商业区和工业区。
那么我们便可以根据每个点和它周围的关系连边，设两个点(x,y)(l,r)
若它们都是商业区，那么利润为Ax,y+Al,r
若它们都是工业区，那么利润为Bx,y+Bl,r
若它们不同，且x,y为商业区，那么利润为Ax,y+Bl,r+Cx,y+Cl,r
然后我们会想到这两个点有关系只能是相邻的点，故我们可以从这里下手：
黑白染色。
然后我们便这样染色：
源点向每个黑点连Ax,y，每个白点Bx,y
每个黑点连汇点Bx,y，每个白点连汇点Ax,y
于是每个黑点向每个白点连Cx,y+Cl,r双向边
至于为什么对，请读者自行思考。

var
    n,m,t,i,j,k,x,y,nu,ans,pp,an:longint;
    a1,b1,c,co:array[1..100,1..100] of longint;
    q:array[1..4,1..2] of longint=((0,1),(1,0),(-1,0),(0,-1));
    las,b,nex,f,a,dis:array[0..30000] of longint;
function min(l,r:longint):longint;
begin
   if l<r then exit(l);exit(r);
end;
procedure insert(x,y,z,k:longint);
begin
   inc(nu);b[nu]:=y;nex[nu]:=las[x];las[x]:=nu;f[nu]:=z;
   inc(nu);b[nu]:=x;nex[nu]:=las[y];las[y]:=nu;f[nu]:=k;
end;
function bfs:boolean;
var l,r,p:longint;
begin
    fillchar(dis,sizeof(dis),0);dis[0]:=1;l:=0;r:=1;a[1]:=0;
    while l<r do begin
       inc(l);p:=las[a[l]];
       while p<>0 do begin
          if (dis[b[p]]=0)and(f[p]>0) then begin
             dis[b[p]]:=dis[a[l]]+1;inc(r);a[r]:=b[p];
          end;p:=nex[p];
       end;
    end;
    exit(dis[t]>0);
end;
function dinic(x,y:longint):longint;
var p,o:longint;
begin
    if x=t then exit(y);
    dinic:=0;p:=las[x];
    while p<>0 do begin
       if (dis[b[p]]=dis[x]+1)and(f[p]>0) then begin
          o:=dinic(b[p],min(y,f[p]));
          if o>0 then begin
             dec(f[p],o);inc(f[p xor 1],o);inc(dinic,o);dec(y,o);
          end;
       end;p:=nex[p];
    end;
end;
begin
    readln(n,m);t:=n*m+1;nu:=1;
    for i:=1 to n do for j:=1 to m do begin read(a1[i,j]);ans:=ans+a1[i,j];end;
    for i:=1 to n do for j:=1 to m do begin read(b1[i,j]);ans:=ans+b1[i,j];end;
    for i:=1 to n do for j:=1 to m do read(c[i,j]);
    for i:=1 to n do for j:=1 to m do if (i+j)mod 2=1 then co[i,j]:=1;
    for i:=1 to n do
        for j:=1 to m do begin
        if co[i,j]=1 then begin insert(0,(i-1)*m+j,a1[i,j],0);insert((i-1)*m+j,t,b1[i,j],0);end
        else begin insert(0,(i-1)*m+j,b1[i,j],0);insert((i-1)*m+j,t,a1[i,j],0);end;
        for k:=1 to 4 do begin
           x:=i+q[k,1];y:=j+q[k,2];
           if (x<1)or(y<1)or(x>n)or(y>m) then continue;
           insert((i-1)*m+j,(x-1)*m+y,c[i,j]+c[x,y],0);pp:=pp+c[i,j]+c[x,y];
        end;
        end;
    while bfs do
    an:=an+dinic(0,maxlongint);
    writeln(ans+pp div 2-an);
end.