社区讨论

关于dsu on tree过本题

P4556【模板】线段树合并 / [Vani 有约会] 雨天的尾巴参与者 2已保存回复 2

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
2 条
当前快照
1 份
快照标识符
@lo8dbkvx
此快照首次捕获于
2023/10/27 16:45
2 年前
此快照最后确认于
2023/10/27 16:45
2 年前
查看原帖
虽然能过,还很快,但是应该是能hack的
因为做完树上差分后,每个点带的数字或多或少,万一轻儿子子树中带的数字比较多,修改操作会进行多次,可能会卡掉
如果把子树数字个数作为size判断轻重儿子,复杂度能正确吗?
CPP
//dsu on tree 做法
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#include<cmath>
#include<map>
#include<bitset>
using namespace std;
typedef long long ll;
int n,m;
int dep[100010],top[100010],fa[100010],son[100010],sz[100010];
int tot=0,head[100010];
vector<int> vec[100010];
int ans[100010];
struct EDGE
{
    int nxt,to;
}edge[200010];
void add(int u,int v)
{
    edge[++tot].nxt=head[u];
    edge[tot].to=v;
    head[u]=tot;
}
void dfs1(int x,int f)
{
    sz[x]=1;
    fa[x]=f;
    dep[x]=dep[f]+1;
    for(int i=head[x];i;i=edge[i].nxt)
    {
        int y=edge[i].to;
        if(y==f) continue;
        dfs1(y,x);
        sz[x]+=sz[y];
        if(sz[y]>sz[son[x]]) son[x]=y;
    }
}
void dfs2(int x,int t)
{
    top[x]=t;
    if(son[x]) dfs2(son[x],t);
    for(int i=head[x];i;i=edge[i].nxt)
    {
        int y=edge[i].to;
        if(y==fa[x]||y==son[x]) continue;
        dfs2(y,y);
    }
}
int lca(int x,int y)
{
    int fx=top[x],fy=top[y];
    while(fx!=fy)
    {
        if(dep[fx]>dep[fy]) swap(x,y),swap(fx,fy);
        y=fa[fy],fy=top[y];
    }
    if(dep[x]>dep[y]) swap(x,y);
    return x;
}
int l,r,maxn[400010],pnt[400010];
void pushup(int k)
{
    if(maxn[k<<1]>=maxn[k<<1|1])
        maxn[k]=maxn[k<<1],pnt[k]=pnt[k<<1];
    else
        maxn[k]=maxn[k<<1|1],pnt[k]=pnt[k<<1|1];
}
void build(int k,int l,int r)
{
    if(l==r)
    {
        maxn[k]=0,pnt[k]=l;
        return;
    }
    int mid=(l+r)>>1;
    build(k<<1,l,mid);
    build(k<<1|1,mid+1,r);
    pushup(k);
}
void ask_add(int k,int l,int r,int x,int d)
{
    if(l==r)
    {
        maxn[k]+=d;
        return;
    }
    int mid=(l+r)>>1;
    if(x<=mid) ask_add(k<<1,l,mid,x,d);
    else ask_add(k<<1|1,mid+1,r,x,d);
    pushup(k);
}
void seg_clear(int k,int l,int r)
{
    maxn[k]=0,pnt[k]=l;
    if(l==r) return;
    int mid=(l+r)>>1;
    if(maxn[k<<1]) seg_clear(k<<1,l,mid);
    if(maxn[k<<1|1]) seg_clear(k<<1|1,mid+1,r);
}
void get_ans(int x,int s)
{
    for(int i=0,mn=vec[x].size();i<mn;++i)
    {
        int z=vec[x][i];
        if(z>0)
            ask_add(1,l,r,z,1);
        else
            ask_add(1,l,r,-z,-1);
    }
    for(int i=head[x];i;i=edge[i].nxt)
    {
        int y=edge[i].to;
        if(y==fa[x]||y==s) continue;
        get_ans(y,s);
    }
}
void dfs(int x)
{
    for(int i=head[x];i;i=edge[i].nxt)
    {
        int y=edge[i].to;
        if(y==fa[x]||y==son[x]) continue;
        dfs(y);
    }
    if(son[x]) dfs(son[x]);
    get_ans(x,son[x]);
    if(maxn[1]) ans[x]=pnt[1];
    if(x==son[fa[x]]) return;
    seg_clear(1,l,r);
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<n;++i)
    {
        int u,v;scanf("%d%d",&u,&v);
        add(u,v),add(v,u);
    }
    dfs1(1,0),dfs2(1,1);
    l=100000,r=1;
    for(int i=1;i<=m;++i)
    {
        int u,v,z;scanf("%d%d%d",&u,&v,&z);
        vec[u].push_back(z),vec[v].push_back(z);
        int anc=lca(u,v);
        vec[anc].push_back(-z),vec[fa[anc]].push_back(-z);
        r=max(r,z),l=min(l,z);
    }
    build(1,l,r);
    dfs(1);
    for(int i=1;i<=n;++i)
    {
        printf("%d\n",ans[i]);
    }
    return 0;
}

回复

2 条回复,欢迎继续交流。

正在加载回复...