• P4183 [USACO18JAN]Cow at Large P


    题目链接

    要考虑子树不要光顾着想对于一个叶子的情况啊!!

    暴力

    考虑贝西所在节点固定时怎么做。我们以贝西为根节点,如果某个农夫能从叶节点迅速控制某个子树的根节点,那么这棵子树就只用一个农夫,否则是其儿子的农夫数之和。用式子表示就是:设 (f_x) 表示 (x) 子树所需要农夫的最少数量,(g_x) 表示能到达 (x) 的最近叶节点的距离,那么有:

    [f_x = 1,dep_x ge g_x ]

    [f_x = sum_yf_y,dep_x < g_x ]

    树形DP即可。复杂度 (O(n^2))

    正解

    考虑如何优化暴力。感觉不是很好优化的样子。我们不去管 (f_x) 具体是多少了,我们只去管最终能对答案贡献 1 的子树。我们发现 (dep_x ge g_x) 这个条件有个特点,满足这个条件的一定是一棵棵子树。因此我们可以在子树的顶点上把这个 1 加上。

    发现还是不好做。考虑树形差分(好像叫这个),一棵子树(整棵树除外)的点的度数的和为 (2siz - 1),即:

    [sum_{i}d_i = 2siz - 1 ]

    移项:

    [sum_{i}2 - d_i = 1 ]

    这样,如果每个点的贡献是 (2 - d_i) 的话,整棵子树的贡献就是 1.(神奇)

    这样,问题就转化为了:

    [ans_p = sum_{x}(2 - d_x)[dis(x,p) ge g_x] ]

    发现这个可以看作点对 ((x, p)) 之间的关系,因此可以点分治做。加在 (g_x - dep_x),查 (le dep_p) 的贡献和。

    关键代码:

    int d[N];
    int g[N];
    void dfs_dp(int cur, int faa) {
      if (d[cur] == 1)  g[cur] = 0;
      else  g[cur] = n;
      for (register int i = head[cur]; i; i = e[i].nxt) {
        int to = e[i].to; if (to == faa)  continue;
        dfs_dp(to, cur); MIN(g[cur], g[to] + 1);
      }
    }
    void dp_dfs(int cur, int faa) {
      for (register int i = head[cur]; i; i = e[i].nxt) {
        int to = e[i].to; if (to == faa)  continue;
        MIN(g[to], g[cur] + 1); dp_dfs(to, cur);
      }
    }
    
    int siz[N], Siz, mxSiz, root;
    bool vis[N];
    void find_root(int cur, int faa)
    uint ans[N];
    void dfs_calc(int cur, int faa, int nwd) {
      ans[cur] += query(nwd);
      for (register int i = head[cur]; i; i = e[i].nxt) {
        int to = e[i].to; if (to == faa || vis[to]) continue;
        dfs_calc(to, cur, nwd + 1);
      }
    }
    void dfs_add(int cur, int faa, int nwd) {
      add(g[cur] - nwd, 2 - d[cur]);
      for (register int i = head[cur]; i; i = e[i].nxt) {
        int to = e[i].to; if (to == faa || vis[to]) continue;
        dfs_add(to, cur, nwd + 1);
      }
    }
    void dfs_del(int cur, int faa, int nwd) {
      del(g[cur] - nwd);
      for (register int i = head[cur]; i; i = e[i].nxt) {
        int to = e[i].to; if (to == faa || vis[to]) continue;
        dfs_del(to, cur, nwd + 1);
      }
    }
    
    int stk[N], stop;
    void dfs(int cur) {
      for (register int i = head[cur]; i; i = e[i].nxt) {
        int to = e[i].to; if (vis[to])  continue;
        dfs_calc(to, cur, 1);
        dfs_add(to, cur, 1); stk[++stop] = to;
      }
      ans[cur] += query(0);
      dfs_del(cur, 0, 0);
      add(g[cur], 2 - d[cur]);
      for (register int i = stop; i; --i) {
        int to = stk[i]; dfs_calc(to, cur, 1); dfs_add(to, cur, 1);
      }
      stop = 0;
      dfs_del(cur, 0, 0);
      vis[cur] = true;
      for (register int i= head[cur]; i; i = e[i].nxt) {
        int to = e[i].to; if (vis[to])  continue;
        Siz = siz[to], mxSiz = n;
        find_root(to, cur);
        dfs(root);
      }
    }
    
    int main() {
      read(n); nn = n + n;
      for (register int i = 1; i < n; ++i) {
        int u, v; read(u), read(v);
        addedge(u, v), addedge(v, u);
        ++d[u]; ++d[v];
      }
      dfs_dp(1, 0); dp_dfs(1, 0);
      Siz = mxSiz = n;
      find_root(1, 0);
      dfs(root);
      for (register int i = 1; i <= n; ++i) {
        if (d[i] == 1)  ans[i] = 1;
        Print(ans[i]);
      }
      return 0;
    }
    
  • 相关阅读:
    结对项目黄金点游戏(邓乐&曾亮)
    软件工程个人项目:一个能自动生成小学四则运算的程序
    IE7下面zindex失效的问题
    参加一战到底有感
    Javascript中的prototype
    浏览器的页面渲染
    使用document.domain实现ajax跨子域
    PHP模拟http请求
    六天带你了解活动营销之第一天营销活动形式多样化
    jQuery源码学习第一天jQuery框架学习
  • 原文地址:https://www.cnblogs.com/JiaZP/p/13586937.html
Copyright © 2020-2023  润新知