首页 > 其他分享 >UVA1674 闪电的能量 Lightning Energy Report 题解

UVA1674 闪电的能量 Lightning Energy Report 题解

时间:2024-05-26 11:34:50浏览次数:24  
标签:rt lazy int 题解 Energy tree Lightning 100010 top

题目传送门

前置知识

树链剖分 | 线段树

解法

树链剖分后维护一个支持区间修改,单点查询的线段树即可。

也可以树上差分,同 146. DFS 序 3,树上差分 1 的 \(1,2\) 操作,时间复杂度比树链剖分更优。

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long 
#define ull unsigned long long
#define sort stable_sort 
#define endl '\n'
struct node
{
    int nxt,to;
}e[100010];
int head[100010],c[100010],cc[100010],siz[100010],fa[100010],dep[100010],son[100010],top[100010],dfn[100010],cnt=0,tot=0;
struct SMT
{
    struct SegmentTree
    {
        int l,r,sum,lazy;
    }tree[200010];
    int lson(int x)
    {
        return x*2;
    }
    int rson(int x)
    {
        return x*2+1;
    }
    void pushup(int rt)
    {
        tree[rt].sum=tree[lson(rt)].sum+tree[rson(rt)].sum;
    }
    void build(int rt,int l,int r)
    {
        tree[rt].l=l;
        tree[rt].r=r;
        tree[rt].lazy=0;
        if(l==r)
        {
            tree[rt].sum=cc[l];
            return;
        }
        int mid=(l+r)/2;
        build(lson(rt),l,mid);
        build(rson(rt),mid+1,r);
        pushup(rt);
    }
    void pushdown(int rt)
    {
        if(tree[rt].lazy!=0)
        {
            tree[lson(rt)].sum+=tree[rt].lazy*(tree[lson(rt)].r-tree[lson(rt)].l+1);
            tree[rson(rt)].sum+=tree[rt].lazy*(tree[rson(rt)].r-tree[rson(rt)].l+1);
            tree[lson(rt)].lazy+=tree[rt].lazy;
            tree[rson(rt)].lazy+=tree[rt].lazy;
            tree[rt].lazy=0;
        }
    }
    void update(int rt,int x,int y,int val)
    {
        if(x<=tree[rt].l&&tree[rt].r<=y)
        {
            tree[rt].sum+=val*(tree[rt].r-tree[rt].l+1);
            tree[rt].lazy+=val;
            return;
        }
        pushdown(rt);
        int mid=(tree[rt].l+tree[rt].r)/2;
        if(x<=mid)
        {
            update(lson(rt),x,y,val);
        }
        if(y>mid)
        {
            update(rson(rt),x,y,val);
        }
    }
    int query(int rt,int x,int y)
    {
        if(x<=tree[rt].l&&tree[rt].r<=y)
        {
            return tree[rt].sum;
        }
        pushdown(rt);
        int mid=(tree[rt].l+tree[rt].r)/2,ans=0;
        if(x<=mid)
        {
            ans+=query(lson(rt),x,y);
        }
        if(y>mid)
        {
            ans+=query(rson(rt),x,y);
        }
        return ans;
    }
}T;
void add(int u,int v)
{
    cnt++;
    e[cnt].nxt=head[u];
    e[cnt].to=v;
    head[u]=cnt;
}
void dfs1(int x,int father)
{
    siz[x]=1;
    fa[x]=father;
    dep[x]=dep[father]+1;
    for(int i=head[x];i!=0;i=e[i].nxt)
    {
        if(e[i].to!=father)
        {
            dfs1(e[i].to,x);
            siz[x]+=siz[e[i].to];
            son[x]=(siz[e[i].to]>siz[son[x]])?e[i].to:son[x];
        }
    }
}
void dfs2(int x,int father,int id)
{
    top[x]=id;
    tot++;
    dfn[x]=tot;
    cc[tot]=c[x];
    if(son[x]!=0)
    {
        dfs2(son[x],x,id);
        for(int i=head[x];i!=0;i=e[i].nxt)
        {
            if(e[i].to!=father&&e[i].to!=son[x])
            {
                dfs2(e[i].to,x,e[i].to);
            }
        }
    }
}
void update1(int u,int v,int val)
{
    while(top[u]!=top[v])
    {
        if(dep[top[u]]>dep[top[v]])
        {
            T.update(1,dfn[top[u]],dfn[u],val);
            u=fa[top[u]];
        }
        else
        {
            T.update(1,dfn[top[v]],dfn[v],val);
            v=fa[top[v]];
        }
    }
    if(dep[u]<dep[v])
    {
        T.update(1,dfn[u],dfn[v],val);
    }
    else
    {
        T.update(1,dfn[v],dfn[u],val);
    }
}
int query1(int pos)
{
    return T.query(1,dfn[pos],dfn[pos]);
}
int main()
{
    int t,n,m,u,v,w,i,j;
    cin>>t;
    for(j=1;j<=t;j++)
    {
        cnt=tot=0;
        memset(e,0,sizeof(e));
        memset(head,0,sizeof(head));
        memset(son,0,sizeof(son));
        cin>>n;
        for(i=1;i<=n-1;i++)
        {
            cin>>u>>v;
            u++;
            v++;
            add(u,v);
            add(v,u);
        }
        dfs1(1,0);
        dfs2(1,0,1);
        T.build(1,1,n);
        cin>>m;
        for(i=1;i<=m;i++)
        {
            cin>>u>>v>>w;
            u++;
            v++;
            update1(u,v,w);
        }
        cout<<"Case #"<<j<<":"<<endl;
        for(i=1;i<=n;i++)
        {
            cout<<query1(i)<<endl;
        }
    }
    return 0;
}

标签:rt,lazy,int,题解,Energy,tree,Lightning,100010,top
From: https://www.cnblogs.com/The-Shadow-Dragon/p/18213457

相关文章

  • CF1141E Superhero Battle 题解
    题目传送门简化题意给定\(H,n\)和一个长度为\(n\)的序列\(d\),求一个最小的\(m\)使得\(H+\sum\limits_{i=1}^{m}d_{(i-1)\bmodn+1}\le0\)。解法将式子移项后得到\(\sum\limits_{i=1}^{m}-d_{(i-1)\bmodn+1}\geH\)。将\(\sum\limits_{i=1}^{m}-d_{(i-1)\bmo......
  • SP14943 DRTREE - Dynamically-Rooted Tree 题解
    题目传送门前置知识树链剖分|线段树解法树剖换根,子树查询板子。类似换根DP的思路,我们发现换根后仅有祖先、子树、深度等会随祖先的变化而变化。设\(rt_{i}\)表示第\(i\)次操作的树根,\(x_{i}\)表示第\(i\)次操作的节点。接着进行大力分讨。当\(rt_{i}=x_{i}\)......
  • 图的dfs遍历题解
    本蒟蒻又来做题啦!看看今天做啥题。。。。题目描述时间:1s 空间:256M题目描述:给出一个无向图和一个起点s,输出这个图从s结点开始的DFS遍历序列。规定:节点邻居按照输入的顺序遍历输入格式:共M+1行。第11行包含33个正整数N,M,s,表示有N个点,M条边,起点为s。第2~M+1行包含2个用......
  • 【题解】 [USACO 2009 Mar] Cow Frisbee Team S
    题目描述题意分析从\(N\)个整数中取若干个(不能不选),且取的数之和为\(F\)的倍数的总方案数对\(10^8\)取余的值。思路这道题是一道二维线性DP。那么根据线性DP的解题方法。首先,找出题目的阶段性。这道题显而易见方案数随着选择的个数而变化,那么这道题就以选择的个数作为阶段......
  • 【赛题解析】【网络建设与运维】2023年全国职业院校技能大赛中职组“网络建设与运维”
    在此之前,欢迎关注波比网络波比网络官方公众号:blbinet波比网络工作室官方公众号:blbistudio技能大赛各赛项交流群:https://www.blbi.cn/threads/40/更多正式赛题源文件访问:https://www.blbi.cn获取技术支持访问:https://www.blbi.cn/form/1/selectNISP、CIPS、PTE证书可......
  • 【赛题解析】【网络建设与运维】2023年全国职业院校技能大赛中职组“网络建设与运维”
    在此之前,欢迎关注波比网络波比网络官方公众号:blbinet波比网络工作室官方公众号:blbistudio技能大赛各赛项交流群:https://www.blbi.cn/threads/40/更多正式赛题源文件访问:https://www.blbi.cn获取技术支持访问:https://www.blbi.cn/form/1/selectNISP、CIPS、PTE证书可......
  • Tokio Marine & Nichido Fire Insurance Programming Contest 2024(AtCoder Beginner C
    A-WhoAtetheCake?题意:有三个嫌疑犯(1,2,3(号码))现在有两个证人他们指出谁不是嫌疑犯,你可以找到确定的那个罪人吗?找到输出这个人的号码没找到输出-1思路:如果两人指出的人是一个人则输出-1不是则输出6-a-b,因为1+2+3=6(sum)减去a,b肯定可以到达......
  • 题解:P8267 [USACO22OPEN] Counting Liars B & U208878 晴天
    其实,这个题,只需要最简单的枚举,加上最简单的二分查找即可~\(1\leN\le1000\)?枚举吧~咋枚举?显然,最好状态下Bessie的位置一定是某个\(p_i\),否则差一个就会导致有个奶牛要说谎。所以我们枚举(理论来讲要先去个重,这样快一点,不过貌似数据没有重的~)\(p_i\),每次遍历这帮奶牛看看有......
  • 2024XJTUPC西安交通大学校赛VP题解
    每次vp都自闭,已成习惯,时间还是分配的不合理,debug时做太多无用功。一键做题A.交小西的礼物输出a+b+2c+3d即可#pragmaGCCoptimize(3,"Ofast","inline")#include<bits/stdc++.h>#defineinf0x3f3f3f3f3fusingnamespacestd;usingll=longlong;usingpii=......
  • 题解【[SCOI2008] 配对】
    题目链接如果没有“配对数字不相同”的限制,将\(a,b\)数组排序后一一配对就能得到最小值。回到原题,考虑一种极端情况,\(\foralli\in[1,n],a_i=b_i\)即排序后全等。若\(n\)为偶数,一种显然的构造方法是:123456214365即分成两个两个一组,然后组内交换,这样跨越幅......