[BZOJ1040] [ZJOI2008]骑士 解题报告

Description

　　Z国的骑士团是一个很有势力的组织，帮会中汇聚了来自各地的精英。他们劫富济贫，惩恶扬善，受到社会各界的赞扬。最近发生了一件可怕的事情，邪恶的Y国发动了一场针对Z国的侵略战争。战火绵延五百里，在和平环境中安逸了数百年的Z国又怎能抵挡的住Y国的军队。于是人们把所有的希望都寄托在了骑士团的身上，就像期待有一个真龙天子的降生，带领正义打败邪恶。骑士团是肯定具有打败邪恶势力的能力的，但是骑士们互相之间往往有一些矛盾。每个骑士都有且仅有一个自己最厌恶的骑士（当然不是他自己），他是绝对不会与自己最厌恶的人一同出征的。战火绵延，人民生灵涂炭，组织起一个骑士军团加入战斗刻不容缓！国王交给了你一个艰巨的任务，从所有的骑士中选出一个骑士军团，使得军团内没有矛盾的两人（不存在一个骑士与他最痛恨的人一同被选入骑士军团的情况），并且，使得这支骑士军团最具有战斗力。为了描述战斗力，我们将骑士按照1至N编号，给每名骑士一个战斗力的估计，一个军团的战斗力为所有骑士的战斗力总和。

 

　　调了很久的样子..首先分析这道题的隐藏条件

　　虽然看起来这张图杂乱无章，但想对于一个有n个骑士的联通块，必定有n条边，也就是说一定是一个存在且仅存在一个环的图

　　据说叫做内向树

　　由于只存在一个环，做法很简单，只需要将环上任意两个点之间的边断开，然后做两次树形DP

　　当控制v不取的时候，u可取可不取。当控制v取的时候，u一定不能取。

　　这样大体的思路就构建起来了。具体实现的时候要注意的地方很多。

　　首先我是用tarjan找的环，与传统的tarjan不同，这里我们要忽略单点，而且要同时保存两个环上的节点

　　找到了之后把边断开的操作可以在树形DP的时候特判一下

　　但是发现WA了之后手造了一组数据就发现tarjan出了问题

　　不管怎样的图好像都成了一个环

　　突然发现tarjan是应用于有向图的...而这道题的边恰好是无向边

　　是不是tarjan就不可以用了呢？实际上不会，因为每个节点最多只出去一条边

　　所以对于环上的任意一个节点只需要保存它到它讨厌的那个骑士的边就可以了

　　读入数据的时候我就既造了有向边又造了无向边

　　还是WA，后来手造一组数据又发现了问题...

　　有一个特殊情况就是可能环的size=2，这个时候实际上和树没什么区别

　　如果把边断开很可能造成的后果就是其中一个点连出去的点都无法到达了

　　而当size>2的时候每个点都会与两个点相连，断了左边的路右边还是可以到的

　　所以当环大小为2的时候需要特判，我直接将v赋为0就解决了这个问题~

　　还是一道非常锻炼思维和代码能力的题啊。

　　最近的状态比较好，有耐心一点一点敲代码，一点一点改

　　果然最后还是AC了呢ww

　　（好像这道题的题解代码都不长...

program bzoj1040;
const maxn=1000010;INF=1000000000001;
var n,time,top,flag,u,v,ans,head,tail,size:int64;
    i,j,j1:longint;
    a,w,fx:array[-1..maxn]of int64;
    ins,vis,used:array[-1..maxn]of boolean;
    low,dfn,link,link1,stack,opt:array[-1..maxn]of int64;
    fa,next,fa1,next1:array[-1..2*maxn]of int64;
    f:array[-1..maxn,0..1]of int64;

function min(a,b:int64):int64;
begin
    if a<b then exit(a) else exit(b);
end;

function max(a,b:int64):int64;
begin
    if a>b then exit(a) else exit(b);
end;

procedure tarjan(p:int64);
var j:longint;
begin
    used[p]:=false;
    inc(top);stack[top]:=p;ins[p]:=true;
    inc(time);dfn[p]:=time;low[p]:=time;
    j:=link1[p];
    while j<>0 do
    begin
        if used[fa1[j]] then
        begin
            tarjan(fa1[j]);
            low[p]:=min(low[p],low[fa1[j]]);
        end else if ins[fa1[j]] then low[p]:=min(low[p],dfn[fa1[j]]);
        j:=next1[j];
    end;
    if low[p]=dfn[p] then
    begin
        if stack[top]=p then
        begin
            ins[p]:=false;dec(top);
        end else
        begin
            v:=stack[top];size:=1;
            ins[v]:=false;
            repeat
                dec(top);
                u:=v;v:=stack[top];
                ins[v]:=false;
                                inc(size);
            until stack[top] = p;
            dec(top);
        end;
    end;
end;

procedure dfs0(p:int64);
var j:longint;
begin
        vis[p]:=false;
    j:=link[p];f[p,0]:=0;
    while j<>0 do
    begin
        if ((p=u)and(fa[j]=v))or((p=v)and(fa[j]=u)) then
        begin
            j:=next[j];continue;
        end;
        if vis[fa[j]] then
        begin
            dfs0(fa[j]);
            inc(f[p,0],max(f[fa[j],0],f[fa[j],1]));
            inc(f[p,1],f[fa[j],0]);
        end;
        j:=next[j];
    end;
    if p=v then f[p,1]:=-INF else inc(f[p,1],w[p]);
end;

procedure dfs1(p:int64);
var j:longint;
begin
        vis[p]:=false;
    j:=link[p];f[p,0]:=0;f[p,1]:=0;
    while j<>0 do
    begin
        if ((p=u)and(fa[j]=v))or((p=v)and(fa[j]=u)) then
        begin
            j:=next[j];continue;
        end;
        if vis[fa[j]] then
        begin
            dfs1(fa[j]);
            inc(f[p,0],max(f[fa[j],0],f[fa[j],1]));
            inc(f[p,1],f[fa[j],0]);
        end;
        j:=next[j];
    end;
    if p=v then f[p,0]:=-INF;
    inc(f[p,1],w[p]);
end;

procedure solve;
var tem:int64;
begin
    fillchar(vis,sizeof(vis),true);
    dfs0(u);
    tem:=max(f[u,0],f[u,1]);
        fillchar(vis,sizeof(vis),true);
    dfs1(u);
    tem:=max(f[u,0],tem);
        inc(ans,tem);
end;

procedure add(x,y:int64);
begin
    inc(j);fa[j]:=y;next[j]:=link[x];link[x]:=j;
    inc(j);fa[j]:=x;next[j]:=link[y];link[y]:=j;
        inc(j1);fa1[j1]:=y;next1[j1]:=link1[x];link1[x]:=j1;
end;

procedure bfs(p:longint);
var j:longint;
begin
    head:=0;tail:=1;opt[1]:=p;fx[p]:=flag;
    while head<>tail do
    begin
                inc(head);
        j:=link[opt[head]];
        while j<>0 do
        begin
            if fx[fa[j]]<>flag then
            begin
                fx[fa[j]]:=flag;
                inc(tail);opt[tail]:=fa[j];
            end;
            j:=next[j];
        end;
    end;
end;

begin
    readln(n);
    j:=0;j1:=0;
    for i:=1 to n do
    begin
        readln(w[i],a[i]);
        add(i,a[i]);
    end;
        ans:=0;
    fillchar(fx,sizeof(fx),0);
    fillchar(used,sizeof(used),true);
    fillchar(ins,sizeof(ins),false);
    flag:=0;time:=0;top:=0;
        for i:=1 to n do if fx[i]=0 then
    begin
                inc(flag);
                bfs(i);
                u:=0;v:=0;
            for j:=1 to tail do if used[opt[j]] then tarjan(opt[j]);
                if size=2 then v:=0;
            solve;
    end;
        writeln(ans);
end.