当前位置:编程学习 > C/C++ >>

Zoj 3527 Shinryaku! Kero Musume (DP_章鱼图上的树形DP)

题目大意:一有向图图n个点,n条边,每个点有且只有一条出边。取某个点会有信仰值,同时某个点与它的后继结点同时取的话, 它的信仰值会改变一个值,问怎么取点,使得总信仰值最大 

解题思路:拓扑排序+树形DP.这一题的图很有特色,我第一次碰到,官方题解说着实章鱼图,因为每个点只有一条出边,那么要么肯定存在一个或多个环,其他的点依靠着这些环。这题攒了一天才AC,原因是一个int64写成了int,悲催啊...
    对于本题,我的解法是先处理章鱼的触须,也就是那些依靠环的边,然后再处理一个个的环。前者就是普通的树形dp,因为每个点只有建或不建两种状态,设dp[i][0]表示i点不建尼姑院获得的最大信仰值,dp[i][1]表示i点建尼姑院获得的最大信仰值,那么触须部分就可以用dp[next][0] = max(dp[i][0],dp[i][1]),dp[next][1] = max(dp[i][0],dp[i][1]-g[i]),或者要加个g[i]是因为同时建尼姑院要减去一定的信仰值。而这个转移的顺序是怎么实现的?我总不能每次都去找叶子然后再一次次往上更新吧?用拓扑排序一步一步往上更新,直到那些环为止,他们的入边永远不为0.
     处理完了触须,我们要开始解剖章鱼真身了。因为是个环,我们不能简单地像前面那样转移,那样的话从i点开始到i点结束的时候i点的状态会很混乱。做这题我想到一种解决环上问题的技巧--拆环。随便从某个点开始就将环拆成链。这题分两种情况进入环,当再碰到这个点时特判下,其他情况转移方程同上。

测试数据:
2
3 -1 2
2 -2 1

4
3 -2 2
4 -3 3
2 -1 1
5 -2 2

4
3 0 2
4 0 3
2 0 1
5 0 2

4
3 -2 2
4 -3 3
2 -2 1
5 -9 2

OutPut:
3
9
14
8

C艹代码:
[cpp]
#include <stdio.h> 
#include <string.h> 
#include <vector> 
using namespace std; 
#define MAX 110000 
#define int64 long long 
#define max(a,b) ((a)>(b)?(a):(b)) 
 
 
struct node { 
 
    int val,g,v; 
}arr[MAX]; 
int st[MAX],top,cnt[MAX],tot;            //拓朴排序用 
int mmap[MAX],n,flag[MAX];               //flag表示是否为环中的点,dp[i][0]表示不在i建,1表示在i建 
int64 dp[MAX][2],power[MAX][2],ans; 
 
 
void Initial() { 
 
    ans = top = 0; 
    memset(dp,0,sizeof(dp)); 
    memset(cnt,0,sizeof(cnt)); 
    memset(flag,0,sizeof(flag)); 
    memset(power,0,sizeof(power)); 

void Debug_Input(){ 
 
    //printf("top = %d\n",tot); 
    for (int i = 1; i <= n; ++i) 
        printf("power[%d][0]=%d power[%d][1]=%d\n",i,power[i][0],i,power[i][1]); 

void ToooPooo() { 
 
    int i,j,k,cur,next; 
 
 
    for (i = 1; i <= n; ++i) 
        if (!cnt[i]) st[++top] = i; 
 
 
    while (top) { 
 
        cur = st[top--]; 
        tot++; 
        flag[cur] = 1;                  //非环中的点 
        dp[cur][1] += arr[cur].val;     //此处建庙则不管后序的情况加上增加的信仰值 
 
 
        next = mmap[cur];           
        dp[next][0] += max(dp[cur][0],dp[cur][1]); 
        dp[next][1] += max(dp[cur][0],dp[cur][1]+arr[cur].g); 
 
 
        cnt[next]--;                 
        if (cnt[next] == 0) st[++top] = next; 
    } 

void Dfs(int s,int kind) { 
 
    flag[s] = 1; 
    int next = mmap[s]; 
    power[s][1] += arr[s].val;         //此处建庙则不管后序的情况加上增加的信仰值 
 
 
    if (!flag[next])  { 
 
        power[next][0] = dp[next][0]; 
        power[next][1] = dp[next][1]; 
        power[next][0] += max(power[s][0],power[s][1]); 
        power[next][1] += max(power[s][0],power[s][1]+arr[s].g); 
        Dfs(next,kind); 
    } 
    else { 
         
        if (kind == 0) power[next][0] = max(power[s][0],power[s][1]); 
        else  power[next][1] = max(power[s][0],power[s][1]+arr[s].g); 
    } 
    if (kind == 0) flag[s] = 0; 

 
 
int main() 

    int i,j,k; 
 
 
    while (scanf("%d",&n) != EOF) { 
 
        Initial(); 
        for (i = 1; i <= n; ++i) { 
 
            scanf("%d%d%d",&arr[i].val,&arr[i].g,&arr[i].v); 
            cnt[arr[i].v]++; 
            mmap[i] = arr[i].v; 
        } 
 
 
        ToooPooo(); 
        for (i = 1; i <= n; ++i) 
            if (flag[i] == 0) { 
 
                int next = mmap[i]; 
                dp[i][1] += arr[i].val; 
                int64 temp = max(dp[i][1],dp[i][0]); 
                power[i][0] = dp[i][0]; 
                power[i][1] = dp[i][1]; 
                power[next][0] = dp[next][0]; 
                power[next][1] = dp[next][1]; 
                power[next][0] += dp[i][0]; 
         &nb

补充:软件开发 , C++ ,
CopyRight © 2022 站长资源库 编程知识问答 zzzyk.com All Rights Reserved
部分文章来自网络,