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

原创 2015年11月18日 20:41:34

题目

最近房地产商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.
版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

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

二元关系 这是两道题,大体相同,不过有一点点不同 我在题目描述中称为题1和题2Description最近房地产商GDOI(Group of Dumbbells Or Idiots)从NOI(Nut...

【2011集训队出题】Crash的数字表格

Description求∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j) n,m

【2011集训队出题】【BZOJ 2154】【JZOJ 1938】Crash的数字表格

Description 莫比乌斯反演 数论 求∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^n\sum_{j=1}^mlcm(i,j) n,m≤107n,m\leq 10^7Analys...

莫队算法,JZOJ1902,【2010集训队出题】小Z的袜子

题目很久之前就经常听到这个算法,但是一直都没有时间打一道莫队的题目,今天终于有时间了QAQ Description  作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿...

【2010集训队出题】小Z的袜子

Description 作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命……   具体来说,小Z把...

BZOJ2132 圈地计划-最小割

题目链接:右转进入题目 题目大意:自行参考原题 解:非常简单,限制有:不能同时有A和B,相邻的点同时A(或B)不能相连 所以显然要黑白染色(就是把相邻的点染成不同的两种颜色,这样检出两派点来...
  • Mys_C_K
  • Mys_C_K
  • 2017年05月11日 15:19
  • 130

【bzoj 2132】圈地计划(最小割)

长夜孤冷,魂梦无归,叹落花流水亦无踪

BZOJ P2132 圈地计划

同happiness dfs写错调了半天的很弱的我 #include #include #include #include #include using namespace std; const ...

【BZOJ 2132】 圈地计划

黑白染色+最小割~

BZOJ2132: 圈地计划

又是不得了的建图姿势。。。。 #include #include #include #define T 10001 const int INF=1...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:【2011集训队出题】圈地计划
举报原因:
原因补充:

(最多只允许输入30个字)