设为首页 加入收藏

TOP

CodeForces #196(Div. 2) 337D Book of Evil (树形dp)
2014-11-23 20:10:29 来源: 作者: 【 】 浏览:4
Tags:CodeForces #196 Div. 337D Book Evil 树形

题目: 点击打开链接

题意
给一棵n个结点的树,任意两个节点的距离是指连接两点的最短的边数
在树上的某个结点有一个“恶魔之书”,这本书会让距离它d以内的节点都受到影响
已知有m个节点收到了影响,问最多有几个结点可能放着“恶魔之书”?

思路
要判断某个点是不是放着书,就要判断这个点的周围d距离以内是否包含所有受影响的m节点
而如果某个节点距离最远的那个受影响节点的距离是L,如果L <= d,那么说明所有受影响的m节点都在d以内,就可判断这个点可能放着书

那么,我们只要能够求出每个节点距离最远的影响节点是多少,就可以O(n)的时间求出答案了。

所以可以用树形dp求解:

f(u, 0): 表示u为顶点的子树中,距u最远的“受影响节点”的距离
f(u, 1): 表示整个树删去u为顶点的子树,但是依旧保留u点为顶点,这个树中距离u最远的“受影响节点”的距离

所有的f(u, 0)可以一次dfs搞定, O(n)
f(u, 1)可以由顶节点一直推下去

f(v, 1) = max{f[brother1][0], f[brother2][0]..., f[brother3][0], f[father][1] | brother是v的兄弟节点,fa是v的父节点} + 1


这一步可以再一次dfs解决,同样是O(n)


代码

**========================================== *
   This is a solution for ACM/ICPC problem * *  
 @source:CodeForces 337D Book of Evil *  
 @type: 树形dp *  

 @author: shuangde *   
@blog: blog.csdn.net/shuangde800 *  
 @email: zengshuangde@gmail.com *===========================================*
/#include #include 
#include
 #include 
#include
 #include 
#include 
#define MP make_pairusing namespace std;
typedef long long int64;t
ypedef pair PII;
const int INF = 0x3f3f3f3f;
const double PI  = acos(-1.0);const int MAXN = 1e5+10;
namespace Adj{    int head[MAXN], size;   
 struct Node{        int v, next;  
   }E[MAXN*2+100];  
  void initAdj() {        size = 0;         
memset(head, -1, sizeof(head));
    }    void addEdge(int u, int v) {        E[size].v = v;  
      E[size].next = head[u];     
   head[u] = size++;  
  }}using namespace Adj;int n, m, d;int f[MAXN][2];
bool vis[MAXN], p[MAXN];int ans;void dfs(int u) {    vis[u] = true;  
  int& ans = f[u][0] = (p[u]   0 : -1);  
  for (int e = head[u]; e != -1; 
e = E[e].next) {        int v = E[e].v;  
      if (vis[v]) continue;     
   dfs(v);   
     if (f[v][0] != -1)          
   ans = max(ans, f[v][0] + 1);   
 }}// 维护m1,m2保存第一大,第二大inline void update(int w, int v, 
PII& m1, PII& m2) {    if (w >= m1.first) {     
   m2 = m1; m1.first = w; m1.second = v;    
} else if (w >= m2.first) {        m2.first = w; 
m2.second = v;    }}void dp(int u) {    vis[u] = true;   
 PII m1 = MP(-1, 0), m2=MP(-1, 0);    
int tmp = max(f[u][0], f[u][1]);  
  if (tmp != -1 && tmp <= d) {        ++ans;    
 }    update(f[u][1], u, m1, m2);   
 for (int e = head[u]; e != -1; e = E[e].next) {        int v = E[e].v;     
   if (vis[v]) continue;    
     if (f[v][0] != -1) {            update(f[v][0]+1, v, m1, m2);  
      }    }    for (int e = head[u]; e != -1; e = E[e].next) {   
     int v = E[e].v;  
      if (vis[v]) continue;       
  f[v][1] = -1;     
   if (v!=m1.second && m1.first!=-1) {        
    f[v][1] = max(f[v][1], m1.first+1);    
    } else if (m2.first != -1) {            f[v][1] = max(f[v][1], m2.first+1);  
      }        dp(v); 
   }}int main(){    while (~scanf("%d%d%d"
, &n, &m, &d)) {        initAdj();    

    memset(p, 0, sizeof(p));       
 for (int i = 0; i < m; ++i) {            int x;   
         scanf("%d", &x);        
    p[x] = true;   
     }         for (int i = 0; i < n - 1; ++i) {            int u, v;     
       scanf("%d%d", &u, &v);        
    addEdge(u, v);       
     addEdge(v, u);   
     }        memset(vis, 0, sizeof(vis));     
   dfs(1);     
   f[1][1] = p[1] 0:-1;   
     ans = 0;   
     memset(vis, 0, sizeof(vis));   
     dp(1);        
printf("%d\n", ans);  
  }    return 0;} 
 

】【打印繁体】【投稿】【收藏】 【推荐】【举报】【评论】 【关闭】 【返回顶部
分享到: 
上一篇UVA10168 下一篇舒适的路线

评论

帐  号: 密码: (新用户注册)
验 证 码:
表  情:
内  容:

·微服务 Spring Boot (2025-12-26 18:20:10)
·如何调整 Redis 内存 (2025-12-26 18:20:07)
·MySQL 数据类型:从 (2025-12-26 18:20:03)
·Linux Shell脚本教程 (2025-12-26 17:51:10)
·Qt教程,Qt5编程入门 (2025-12-26 17:51:07)