1475: 方格取数

1475: 方格取数

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 578  Solved: 309
[Submit][Status][Discuss]

Description

在一个n*n的方格里，每个格子里都有一个正整数。从中取出若干数，使得任意两个取出的数所在格子没有公共边，且取出的数的总和尽量大。

Input

第一行一个数n；（n<=30） 接下来n行每行n个数描述一个方阵

Output

仅一个数，即最大和

Sample Input

2
1 2
3 5

Sample Output

6

HINT

 

Source

题解：又是一个网络流= =，貌似学到了一个新的建图模式——将两个块的分居两端，分别与源和汇相连，然后中间直接相互冲突的（比如本题中相邻的）连上正无穷边，然后网络流后再用总和剪下= =，然后没了，感觉网络流这东西越学越神奇QAQ么么哒

/**************************************************************
    Problem: 1475
    User: HansBug
    Language: Pascal
    Result: Accepted
    Time:304 ms
    Memory:612 kb
****************************************************************/
 
type
    point=^node;
    node=record
               g,w:longint;
               anti,next:point;
    end;
var
   i,j,k,l,m,n,s,t,flow:longint;
   a,e:array[0..10000] of point;
   b,c,d:array[0..10000] of longint;
function min(x,y:longint):longint;
         begin
              if x<y then min:=x else min:=y;
         end;
procedure add(x,y,z:longint);
          var p:point;
          begin
               new(p);p^.g:=y;p^.w:=z;p^.next:=a[x];a[x]:=p;
               new(p);p^.g:=x;p^.w:=0;p^.next:=a[y];a[y]:=p;
               a[x]^.anti:=a[y];a[y]^.anti:=a[x];
          end;
function spfa:boolean;
         var f,r,i:longint;p:point;
          begin
               for i:=s to t do c[i]:=maxlongint;
               f:=1;r:=2;d[1]:=s;c[s]:=0;
               while f<r do
                     begin
                          p:=a[d[f]];
                          while p<>nil do
                                begin
                                     if (p^.w<>0) and (c[p^.g]=maxlongint) then
                                        begin
                                             e[p^.g]:=p;
                                             c[p^.g]:=c[d[f]]+1;
                                             d[r]:=p^.g;inc(r);
                                        end;
                                     p:=p^.next;
                                end;
                          inc(f);
                     end;
               exit(c[t]<>maxlongint);
          end;
procedure deal;
          begin
               i:=t;l:=maxlongint;
               while i<>s do
                     begin
                          l:=min(e[i]^.w,l);
                          i:=e[i]^.anti^.g;
                     end;
               dec(flow,l);i:=t;
               while i<>s do
                     begin
                          if e[i]^.w<>maxlongint then dec(e[i]^.w,l);
                          if e[i]^.anti^.w<>maxlongint then inc(e[i]^.anti^.w,l);
                          i:=e[i]^.anti^.g;
                     end;
          end;
begin
     readln(n);flow:=0;s:=1;t:=n*n+2;
     for i:=1 to n*n+2 do a[i]:=nil;
     for i:=1 to n do
         for j:=1 to n do
             begin
                  read(k);inc(flow,k);
                  if odd(i+j) then add(1,(i-1)*n+j+1,k) else add((i-1)*n+j+1,n*n+2,k);
                  if odd(i+j) then
                     begin
                          if i>1 then add((i-1)*n+j+1,(i-2)*n+j+1,maxlongint);
                          if i<n then add((i-1)*n+j+1,(i)*n+j+1,maxlongint);
                          if j>1 then add((i-1)*n+j+1,(i-1)*n+j,maxlongint);
                          if j<n then add((i-1)*n+j+1,(i-1)*n+j+2,maxlongint);
                     end;
                  if j=n then readln;
             end;
     while spfa do deal;
     writeln(flow);
     readln;
end.