脑子不清醒的时候千万别写题。写题写不下去了千万别死扛,重构才是你唯一的出路QAQ

  昨天很想快点写道题,思路没有很清晰的时候就写了,结果……今天一怒之下决定重整思路重构代码,其实不过是半个小时的事情……

  提示很明显,总点数限制了范围。建立出虚树,在虚树上面 dp。 在虚树上面我们dp两遍,两遍一起处理出每一个(虚树上的)节点被谁管辖,在连接两者的路径上找到管理的分界点即可。处理两点之间的距离一定要上 ST 表啊……

#include <bits/stdc++.h>
using namespace std;
#define maxn 350000
#define CNST 20
int n, dep[maxn], gra[maxn][CNST];
int size[maxn], hson[maxn], ans[maxn];
int top, S[maxn], a[maxn], b[maxn];
int timer, dfn[maxn], sum[maxn];
bool mark[maxn], vis[maxn];
int T, Log[maxn * 2], bit[21], ST[maxn * 2][CNST], pos[maxn]; int read()
{
    int x = 0, k = 1;
    char c;
    c = getchar();
    while(c < '0' || c > '9') { if(c == '-') k = -1; c = getchar(); }
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * k;
} struct edge
{
    int cnp, to[maxn * 2], last[maxn * 2], head[maxn];
    edge() { cnp = 1; }
    void add(int u, int v)
    {
        if(u == v) return;
        to[cnp] = v, last[cnp] = head[u], head[u] = cnp ++;
        to[cnp] = u, last[cnp] = head[v], head[v] = cnp ++;
    }
}E1, E2; bool cmp(int x, int y) { return dfn[x] < dfn[y]; } void dfs(int u, int fa)
{
    size[u] = 1; dfn[u] = ++ timer;
    gra[u][0] = fa; dep[u] = dep[fa] + 1; ST[++ T][0] = dep[u]; pos[u] = T;
    for(int i = 1; i < CNST; i ++) gra[u][i] = gra[gra[u][i - 1]][i - 1];
    for(int i = E1.head[u]; i; i = E1.last[i])
    {
        int v = E1.to[i];
        if(v == fa) return;
        dfs(v, u); size[u] += size[v];
        ST[++ T][0] = dep[u];
    }
} int LCA(int x, int y)
{
    if(dep[x] < dep[y]) swap(x, y);
    for(int i = CNST - 1; ~i; i --)
        if(dep[gra[x][i]] >= dep[y]) x = gra[x][i];
    for(int i = CNST - 1; ~i; i --)
        if(gra[x][i] != gra[y][i]) x = gra[x][i], y = gra[y][i];
    return (x == y) ? x : gra[x][0];
} int Get_Dis(int x, int y)
{
    int ret = dep[x] + dep[y];
    if(!x) return dep[y]; else if(!y) return dep[x];
    x = pos[x], y = pos[y];
    if(x > y) swap(x, y);
    int k = Log[y - x + 1];
    ret = ret - 2 * min(ST[x][k], ST[y - bit[k] + 1][k]);
    return ret;
} int Get_Gra(int x, int y)
{
    for(int i = CNST - 1; ~i; i --)
        if(dep[gra[x][i]] > dep[y]) x = gra[x][i];
    return x;
} int Work(int u, int v, int k)
{
    bool flag = 0;
    int x = u;
    for(int i = CNST - 1; i >= 0; i --)
    {
        int tem = gra[x][i]; if(dep[tem] < dep[k]) continue;
        int d1 = Get_Dis(u, tem), d2 = Get_Dis(v, tem);
        if(d1 < d2) x = gra[x][i];
        else if(d1 == d2 && u < v) x = gra[x][i];
    }
     return x;
} void DP(int u, int fa)
{
    hson[u] = 0; sum[u] = 0;
    for(int i = E2.head[u]; i; i = E2.last[i])
    {
        int v = E2.to[i]; if(!v) continue;
        if(v == fa) continue;
        DP(v, u);
        int t = hson[v];
        int dis = Get_Dis(u, t), dis2 = Get_Dis(u, hson[u]);
        if(dis < dis2) hson[u] = t;
        else if(dis == dis2 && t < hson[u]) hson[u] = t;
        if(!hson[u]) hson[u] = t;
    }
    if(mark[u]) hson[u] = u;
} void DP2(int u, int fa)
{
    int tot = 0;
    if(hson[fa])
    {
        int t = hson[fa];
        int dis = Get_Dis(u, t), dis2 = Get_Dis(u, hson[u]);
        if(dis < dis2) hson[u] = t;
        else if(dis == dis2 && t < hson[u]) hson[u] = t;
        if(!hson[u]) hson[u] = t;
    }
    for(int i = E2.head[u]; i; i = E2.last[i])
    {
        int v = E2.to[i];
        if(v == fa) continue; DP2(v, u);
        int y = Get_Gra(v, u);
        int x = Work(hson[v], hson[u], u);
        int tem = size[x] - size[v];
        ans[hson[v]] += tem, ans[hson[u]] += size[y] - size[x];
        tot += size[y];
    }
    ans[hson[u]] += size[u] - tot;
    E2.head[u] = 0;
} void Solve()
{
    int K = read();
    for(int i = 1; i <= K; i ++) b[i] = a[i] = read(), mark[a[i]] = 1;
    sort(a + 1, a + 1 + K, cmp); S[top = 1] = 1; E2.cnp = 1; ans[0] = 0;
    for(int i = 1; i <= K; i ++)
    {
        int lca = LCA(S[top], a[i]);
        while(2333)
        {
            if(dep[S[top - 1]] <= dep[lca])
            {
                E2.add(lca, S[top --]);
                if(lca != S[top]) S[++ top] = lca;
                break;
            }
            else E2.add(S[top], S[top - 1]), top --;
        }
        S[++ top] = a[i];
    }
    while(top > 1) E2.add(S[top], S[top - 1]), top --;
    sum[1] = 0; DP(1, 0); DP2(1, 0);
    for(int i = 1; i <= K; i ++) printf("%d ", ans[b[i]]), mark[b[i]] = 0, ans[b[i]] = 0;
    puts("");
} void init()
{
    bit[0] = 1; for(int i = 1; i <= 20; i ++) bit[i] = bit[i - 1] << 1;
    Log[0] = -1; for(int i = 1; i < maxn * 2; i ++) Log[i] = Log[i >> 1] + 1;
} int main()
{
    init(); n = read();
    for(int i = 1; i < n; i ++)
    {
        int x = read(), y = read();
        E1.add(x, y);
    }
    dfs(1, 0);
    for(int i = 1; i <= Log[T]; i ++)
        for(int j = 1; j + bit[i - 1] <= T; j ++)
            ST[j][i] = min(ST[j][i - 1], ST[j + bit[i - 1]][i - 1]);
    int q = read();
    for(int i = 1; i <= q; i ++) Solve();
    return 0;
}

【题解】HNOI2014世界树的更多相关文章

  1. [题解] [HNOI2014] 世界树

    题面 [HNOI2014]世界树 题解 从数据范围很容易看出是个虚树DP(可惜看出来了也还是不会做) 虚树大家应该都会, 不会的话自己去搜吧, 我懒得讲了, 我们在这里只需要考虑如何DP即可 首先我们 ...

  2. [BZOJ3572][Hnoi2014]世界树

    [BZOJ3572][Hnoi2014]世界树 试题描述 世界树是一棵无比巨大的树,它伸出的枝干构成了整个世界.在这里,生存着各种各样的种族和生灵,他们共同信奉着绝对公正公平的女神艾莉森,在他们的信条 ...

  3. 【BZOJ3572】[Hnoi2014]世界树 虚树

    [BZOJ3572][Hnoi2014]世界树 Description 世界树是一棵无比巨大的树,它伸出的枝干构成了整个世界.在这里,生存着各种各样的种族和生灵,他们共同信奉着绝对公正公平的女神艾莉森 ...

  4. bzoj3572[Hnoi2014] 世界树 虚树+dp+倍增

    [Hnoi2014]世界树 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 1921  Solved: 1019[Submit][Status][Dis ...

  5. BZOJ 3572: [Hnoi2014]世界树

    BZOJ 3572: [Hnoi2014]世界树 标签(空格分隔): OI-BZOJ OI-虚数 OI-树形dp OI-倍增 Time Limit: 20 Sec Memory Limit: 512 ...

  6. bzoj 3572: [Hnoi2014]世界树 虚树 && AC500

    3572: [Hnoi2014]世界树 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 520  Solved: 300[Submit][Status] ...

  7. bzoj 3572 [Hnoi2014]世界树(虚树+DP)

    3572: [Hnoi2014]世界树 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 645  Solved: 362[Submit][Status] ...

  8. <虚树+树型DP> HNOI2014世界树

    <虚树+树型DP> HNOI2014世界树 #include <iostream> #include <cstdio> #include <cstring&g ...

  9. BZOJ3572:[HNOI2014]世界树——题解

    +++++++++++++++++++++++++++++++++++++++++++ +本文作者:luyouqi233. + +欢迎访问我的博客:http://www.cnblogs.com/luy ...

随机推荐

  1. netty之粘包分包的处理

    1.netty在进行字节数组传输的时候,会出现粘包和分包的情况.当个数据还好,如果数据量很大.并且不间断的发送给服务器,这个时候就会出现粘包和分包的情况. 2.简单来说:channelBuffer在接 ...

  2. Unity Container中的几种注册方式与示例

    1.实例注册 最简单的注册方式就是实例注册,Unity 容器负责维护对一个类型的单例引用,比如: 有如下的实际类型: namespace ConsoleSample { public class Sa ...

  3. 05-JVM对象探秘

    一.对象的内存布局         以Hotspot虚拟机为例,对象在内存中的结构可以分为三部分:对象头(header).实例数据(instance data).对齐填充(padding). 1.1. ...

  4. MySQL☞数值处理函数

    1.round():四舍五入函数 round(数值,参数):如果参数的值为正数,表示保留几位小数,如果参数的值为0,则只保留正数部分们如果参数的值为负数,表示对小数点前第几位进行四舍五入. Eg:(1 ...

  5. 通过 Python_Faker 生成测试数据

    通过 Python_Faker 生成测试数据 一.介绍 在软件需求.开发.测试过程中,有时候需要使用一些测试数据,针对这种情况,我们一般要么使用已有的系统数据,你不可能通过手工来生成(最傻的方法)可能 ...

  6. Anyproxy抓包工具

    1.安装Nodejs和AnyProxy以及安装模拟器(移动端抓包工具) 第一步:nodejs下载地址:http://nodejs.cn/download/ 下载Windows版本,直接运行安装即可, ...

  7. lintcode735. Replace With Greatest From Right

    Given an array of integers, replace every element with the next greatest element (greatest element o ...

  8. Spring Boot - Filter实现简单的Http Basic认证

    Copy自http://blog.csdn.net/sun_t89/article/details/51916834 @SpringBootApplicationpublic class Spring ...

  9. 如何让thinkpad X1C 用U盘 安装上专业版win10

    1 BIOS内置了文件 会导致win10 iso默认装家庭版 2 给iso 的resouse 目录中增加文件ei.cfg 3 内容如下 [EditionID]Professional[Channel] ...

  10. POJ 3525/UVA 1396 Most Distant Point from the Sea(二分+半平面交)

    Description The main land of Japan called Honshu is an island surrounded by the sea. In such an isla ...