1565: [NOI2009]植物大战僵尸 - BZOJ

Description

Input

Output
仅包含一个整数，表示可以获得的最大能源收入。注意，你也可以选择不进行任何攻击，这样能源收入为0。
Sample Input
3 2
10 0
20 0
-10 0
-5 1 0 0
100 1 2 1
100 0
Sample Output
25
HINT

在样例中, 植物P1,1可以攻击位置(0,0), P2, 0可以攻击位置(2,1)。
一个方案为，首先进攻P1,1, P0,1，此时可以攻击P0,0 。共得到能源收益为(-5)+20+10 = 25。注意, 位置(2,1)被植物P2,0保护，所以无法攻击第2行中的任何植物。
【大致数据规模】
约20%的数据满足1 ≤ N, M ≤ 5；
约40%的数据满足1 ≤ N, M ≤ 10；
约100%的数据满足1 ≤ N ≤ 20，1 ≤ M ≤ 30，-10000 ≤ Score ≤ 10000 。

 

 

详情请见胡伯涛《最小割模型在信息学竞赛中的应用》中最大权闭合图

但是有一点不同，因为有些点不能取，所以我们首先拓扑一下，把有用的点选出来，然后再跑最小割

const
    maxn=22;
    maxm=33;
    inf=10000000;
var
    d,a,first,q:array[0..maxn*maxm]of longint;
    next,last:array[0..maxn*maxm*maxn*maxm]of longint;
    map:array[0..maxn*maxm,0..maxn*maxm]of longint;
    n,m,tot,sum,cnt:longint;
 
function calc(i,j:longint):longint;
begin
    exit(i*m+j+1);
end;
 
procedure insert(x,y:longint);
begin
    inc(tot);
    last[tot]:=y;
    next[tot]:=first[x];
    first[x]:=tot;
    inc(d[y]);
end;
 
procedure init;
var
    i,j,k,l,r,x,y:longint;
begin
    read(n,m);
    for i:=0 to n-1 do
        for j:=0 to m-1 do
            begin
                read(a[calc(i,j)]);
                read(k);
                if j>0 then insert(calc(i,j),calc(i,j-1));
                for l:=1 to k do
                    begin
                        read(x,y);
                        insert(calc(i,j),calc(x,y));
                    end;
            end;
    l:=1;r:=0;
    for i:=0 to n -1do
        if d[calc(i,m-1)]=0 then
        begin
            inc(r);
            q[r]:=calc(i,m-1);
        end;
    while l<=r do
        begin
            if a[q[l]]>0 then inc(sum,a[q[l]]);
            if a[q[l]]>0 then inc(map[0,q[l]],a[q[l]]);
            if a[q[l]]<0 then inc(map[q[l],n*m+1],-a[q[l]]);
            i:=first[q[l]];
            while i<>0 do
                begin
                    dec(d[last[i]]);inc(map[last[i],q[l]],inf);
                    if d[last[i]]=0 then
                    begin
                        inc(r);
                        q[r]:=last[i];
                    end;
                    i:=next[i];
                end;
            inc(l);
        end;
    cnt:=r;
end;
 
var
    dis,vh,his,pre:array[0..maxn*maxm]of longint;
    flow:longint;
 
procedure sap;
var
    i,j,aug,min:longint;
    flag:boolean;
begin
    vh[0]:=cnt+2;
    i:=0;aug:=inf;
    while dis[i]<n*m+2 do
        begin
            his[i]:=aug;
            flag:=false;
            for j:=0 to n*m+1 do
                if (map[i,j]>0) and (dis[i]=dis[j]+1) then
                begin
                    flag:=true;
                    if aug>map[i,j] then aug:=map[i,j];
                    pre[j]:=i;
                    i:=j;
                    if i=n*m+1 then
                    begin
                        inc(flow,aug);
                        while i<>0 do
                            begin
                                inc(map[i,pre[i]],aug);
                                dec(map[pre[i],i],aug);
                                i:=pre[i];
                            end;
                        aug:=inf;
                    end;
                    break;
                end;
            if flag then continue;
            min:=n*m+1;
            for j:=0 to n*m+1 do
                if (map[i,j]>0) and (dis[j]<min) then min:=dis[j];
            dec(vh[dis[i]]);
            if vh[dis[i]]=0 then break;
            dis[i]:=min+1;
            inc(vh[min+1]);
            if i<>0 then
            begin
                i:=pre[i];
                aug:=his[i];
            end;
        end;
    writeln(sum-flow);
end;
 
begin
    init;
    sap;
end.