CodeVS 1022-覆盖

原题

题目描述 Description

有一个N×M的单位方格中,其中有些方格是水塘,其他方格是陆地.如果要用1×2的矩阵区覆盖(覆盖过程不容许有任何部分重叠)这个陆地,那么最多可以覆盖多少陆地面积.

输入描述 Input Description

输入文件的第一行是两个整数N,M  (1<=N，M<=100)，第二行为一个整数K( K<=50)，接下来的K行，每行两个整数X，Y表示K个水塘的行列位置。(1<=X<=N，1<=Y<=M)。

 

输出描述 Output Description

输出所覆盖的最大面积块（1×2面积算一块）。

样例输入 Sample Input

4 4

6

1 1

1 4

2 2

4 1

4 2

4 4

样例输出 Sample Output

4

题解

这道题一开始看以为是什么状压DP就一直没敢做,但最近周围各位神犇都在写这道题,就认真地再看了一次.

由于是求1*2或者2*1的块的个数,那就可以用二分图来做.

建图的时候就把相邻两个点连接起来.

每个点最多只能在二分图的左边出现一次,如果某个点在二分图的右边出现过一次,那么它同样不能加入二分图的左边.

由于我的算法是Dinic,因此要建一个原点与二分图左边的每个点一一相连.

同样的,二分图右边的每个点也要与汇点建一条边,但值得注意的是,二分图右边的点最多只能和汇点连一条边.

因此,我们就要建立两个数组usex/usey分别表示某一个点是否在二分图的左边或是右边出现,是为True,否为False.

建图之后直接套Dinic就OK了.

Code:

const fx:array[1..4] of longint=(1,-1,0,0);
const fy:array[1..4] of longint=(0,0,1,-1);
var ne,b,c,q:array[0..100000] of longint;
var fir,dis:array[0..1001] of longint;
var usex,usey:array[0..1001] of boolean;
var a:array[0..100,0..100] of boolean;
var i,j,n,m,tot,x,y,k,ans:longint;
function v(x,y:longint):longint; begin exit((x-1)*m+y); end;
procedure add(x,y:longint);
var ty,tx,i,xx,yy:longint;
begin
  tx:=v(x,y);
  usex[tx]:=true;
  for i:=1 to 4 do
  begin
    xx:=x+fx[i];yy:=y+fy[i];
    ty:=v(xx,yy);
    if (not a[xx,yy])and(not usex[ty])and(xx<=n)and(xx>0)and(yy<=m)and(yy>0) then
    begin
      inc(tot);b[tot]:=tx;c[tot]:=1;ne[tot]:=fir[0];fir[0]:=tot;//左边的点与原点建边
      inc(tot);b[tot]:=0;c[tot]:=0;ne[tot]:=fir[tx];fir[tx]:=tot;
      inc(tot);b[tot]:=ty;c[tot]:=1;ne[tot]:=fir[tx];fir[tx]:=tot;//两点之间建边
      inc(tot);b[tot]:=tx;c[tot]:=0;ne[tot]:=fir[ty];fir[ty]:=tot;
      if not usey[ty] then//只能和汇点建一条边!!!只能和汇点建一条边!!!只能和汇点建一条边!!!
      begin
        inc(tot);b[tot]:=m*n+1;c[tot]:=1;ne[tot]:=fir[ty];fir[ty]:=tot;
        inc(tot);b[tot]:=ty;c[tot]:=0;ne[tot]:=fir[n*m+1];fir[n*m+1]:=tot;
      end;
      usey[ty]:=true;
    end;
  end;
end;
procedure bfs;
var i,h,t,now:longint;
begin
  for i:=0 to n*m+1 do dis[i]:=-1;
  h:=1;t:=1;dis[0]:=0;q[1]:=0;
  while h<=t do
  begin
    now:=q[h];j:=fir[now];
    while j>0 do
    begin
      if (c[j]>0)and(dis[b[j]]=-1) then
      begin
        dis[b[j]]:=dis[now]+1;
        inc(t);q[t]:=b[j];
      end;j:=ne[j];
    end;inc(h);
  end;
end;
function dfs(x,mx:longint):longint;
var k,j:longint;
begin
  if x=n*m+1 then exit(mx);j:=fir[x];
  while j>0 do
  begin
    if (c[j]>0)and(dis[b[j]]=dis[x]+1) then
    begin
      if mx<c[j] then k:=dfs(b[j],mx) else k:=dfs(b[j],c[j]);
      dec(c[j],k);
      if odd(j) then inc(c[j+1],k) else inc(c[j-1],k);
      if k>0 then exit(k);
    end;j:=ne[j];
  end;
  exit(0);
end;
begin
  readln(n,m);
  readln(k);
  for i:=1 to k do begin readln(x,y);a[x,y]:=true; end;
  fillchar(fir,sizeof(fir),255);
  for i:=1 to n do for j:=1 to m do if (not a[i,j])and(not usey[v(i,j)]) then add(i,j);
  while true do
  begin
    bfs;//建分层图
    if dis[n*m+1]=-1 then break;//如果到不了汇点则跳出循环
    ans:=ans+dfs(0,maxlongint);
  end;
  writeln(ans);
end.

最后我发现我似乎是这道题唯一0ms过的人233333......
测试点#cover0.in 结果:AC 内存使用量: 256kB 时间使用量: 0ms 
测试点#cover1.in 结果:AC 内存使用量: 256kB 时间使用量: 0ms 
测试点#cover2.in 结果:AC 内存使用量: 256kB 时间使用量: 0ms 
测试点#cover3.in 结果:AC 内存使用量: 256kB 时间使用量: 0ms 
测试点#cover4.in 结果:AC 内存使用量: 256kB 时间使用量: 0ms 
测试点#cover5.in 结果:AC 内存使用量: 256kB 时间使用量: 0ms 
测试点#cover6.in 结果:AC 内存使用量: 256kB 时间使用量: 0ms 
测试点#cover7.in 结果:AC 内存使用量: 256kB 时间使用量: 0ms 
测试点#cover8.in 结果:AC 内存使用量: 256kB 时间使用量: 0ms 
测试点#cover9.in 结果:AC 内存使用量: 256kB 时间使用量: 0ms