[NOIP2016]天天爱跑步-题解

2019-10-08 08:46:47来源:博客园 阅读 ()

新老客户大回馈,云服务器低至5折

[NOIP2016]天天爱跑步-题解

题面传送门

解答

设第\(j\)号玩家在\(V_j\)时刻出发。

弱化问题:如果树退化成了一条链。则在\(j\)处的观察员能观察到的\(i\)号玩家当且仅当
\[ i玩家经过j,且 \begin{cases} dep_j - W_j = dep_{S_i} - V_j, &i向下跑步 \\ dep_j + W_j = dep_{S_i} + V_j, &i向上跑步 \end{cases} \]
一个点在树上的贡献是连续的,可以考虑把路径在LCA处(倍增找LCA)拆成两条链用离线+树上差分统计人经过每个点的时刻。

树上差分

如何只让一段树上的链被更新呢?(类比左闭右开区间,考虑下端为“闭”节点,上端为“开”节点的链)

考虑树上的一个点,只有链的一端在它的子树内部时才被链覆盖。
所以在dfs子树之前记下答案,在dfs之后更新链端点,真的答案就是答案的差。

#include <bits/stdc++.h>
using namespace std;

const int mx = 3e5+10, LOG = 20, inf = 1e9;

int n, m, ans[mx], s[mx], t[mx], w[mx];
int head[mx], nxt[mx<<1], to[mx<<1], tot;
int dep[mx] = {-1}, fa[mx][LOG];
int down[mx<<2], up[mx<<1];
vector<pair<int, int> > v[mx];

void add(int u, int v) {
    ++tot;
    nxt[tot] = head[u];
    to[tot] = v;
    head[u] = tot;
}

void dfs0(int u) {
    for (int i = 1; i < LOG; ++i) {
        fa[u][i] = fa[fa[u][i-1]][i-1];
    }
    for (int e = head[u], v; e; e = nxt[e]) {
        if ((v=to[e]) == fa[u][0]) continue;
        fa[v][0] = u;
        dep[v] = dep[u]+1;
        dfs0(v);
    }
}

int lca(int u, int v) {
    if (dep[u] < dep[v]) swap(u, v);
    for (int i = LOG-1; ~i; --i)
        if (dep[fa[u][i]] >= dep[v])
            u = fa[u][i];
    if (u == v) return u;
    for (int i = LOG-1; ~i; --i)
        if (fa[u][i] != fa[v][i])
            u = fa[u][i], v = fa[v][i];
    return fa[u][0];
}

#define pb push_back
#define mp make_pair
#define fi first
#define se second

void addr(int s, int t, int w=0) {
    assert(dep[s] != dep[t]);
    if (dep[s] < dep[t]) w += inf, swap(s, t);
    v[s].pb(mp(w, 1));
    v[t].pb(mp(w, -1));
}

void dfs(int u) {
    ans[u] -= up[dep[u]+w[u]] + down[(mx<<1) + dep[u]-w[u]];
    for (int e = head[u]; e; e = nxt[e]) {
        if (to[e] == fa[u][0]) continue;
        dfs(to[e]);
    }
    for (int i = 0; i < (int)v[u].size(); ++i) {
        if (v[u][i].fi > inf/2) down[(mx<<1) + v[u][i].fi-inf] += v[u][i].se;
        else up[v[u][i].fi] += v[u][i].se;
    }
    ans[u] += up[dep[u]+w[u]] + down[(mx<<1) + dep[u]-w[u]];
}

int main() {
    scanf("%d%d", &n, &m);
    int u, v;
    for (int i = 0; i < n-1; ++i) {
        scanf("%d%d", &u, &v);
        add(u, v), add(v, u);
    }
    dfs0(1);
    for (int i = 1; i <= n; ++i) scanf("%d", w+i);
    for (int i = 1; i <= m; ++i) {
        scanf("%d%d", s+i, t+i);
        u = lca(s[i], t[i]);
        if (u == s[i]) addr(fa[s[i]][0], t[i], dep[s[i]]);
        else if (u == t[i]) addr(s[i], fa[t[i]][0], dep[s[i]]);
        else {
            addr(s[i], fa[u][0], dep[s[i]]);
            addr(u, t[i], dep[u] - (dep[s[i]]-dep[u]));
        }
    }
    dfs(1);
    for (int i = 1; i <= n; ++i) printf("%d%c", ans[i], " \n"[i==n]);
}

原文链接:https://www.cnblogs.com/topsecret/p/noip2016-running.html
如有疑问请与原作者联系

标签:

版权申明:本站文章部分自网络,如有侵权,请联系:west999com@outlook.com
特别注意:本站所有转载文章言论不代表本站观点,本站所提供的摄影照片,插画,设计作品,如需使用,请与原作者联系,版权归原作者所有

上一篇:a++与++a,谈谈C++的参数传递

下一篇:fopen