这题好久之前就被学长安利了...一直没写珍藏在收藏夹一个不为人知的角落233

  这题怎么做...我们来数形结合,横坐标为$t_i$被加的次数(可看作时间$t$),纵坐标为$v_i$,那么$t_i$实际上就是阶梯图形的面积。

  上图是父亲节点和子节点合并后的样子,阴影部分为答案即$t_i$,显然可以通过$deltat*deltav-deltatv$来得到。

  信息下传的时候儿子的$deltatv$加上父亲的$deltatv$后还要加上一个长方体的面积即$deltat_{son}*deltav_{fa}$,然后儿子的$deltat$和$deltav$都加上父亲的就行了。

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn=, inf=1e9;
struct poi{int too, pre;}e[maxn<<];
struct tjm{ll deltat, deltav, deltatv;}tree[maxn<<];
int n, m, tott, tot, ty, x, y, z;
int fa[maxn], d[maxn], size[maxn], top[maxn], w[maxn], last[maxn], son[maxn];
char buf[], *ptr=buf-;
inline void read(int &k)
{
int f=; k=; char c=*++ptr;
while(c<'' || c>'') c=='-'&&(f=-), c=*++ptr;
while(c<='' && c>='') k=k*+c-'', c=*++ptr;
k*=f;
}
inline void add(int x, int y){e[++tot]=(poi){y, last[x]}; last[x]=tot;}
void dfs1(int x)
{
size[x]=;
for(int i=last[x], too;i;i=e[i].pre)
{
d[too=e[i].too]=d[x]+;
dfs1(too);
size[x]+=size[too];
if(size[too]>size[son[x]]) son[x]=too;
}
}
void dfs2(int x, int tp)
{
top[x]=tp; w[x]=++tott;
if(son[x]) dfs2(son[x], tp);
for(int i=last[x], too;i;i=e[i].pre)
if((too=e[i].too)!=son[x]) dfs2(too, too);
}
inline void addone(int x, ll deltat, ll deltav, ll deltatv)
{
tree[x].deltatv+=deltatv+deltav*tree[x].deltat;
tree[x].deltat+=deltat;
tree[x].deltav+=deltav;
}
inline void down(int x)
{
addone(x<<, tree[x].deltat, tree[x].deltav, tree[x].deltatv);
addone(x<<|, tree[x].deltat, tree[x].deltav, tree[x].deltatv);
tree[x].deltat=tree[x].deltav=tree[x].deltatv=;
}
void update(int x, int l, int r, int cl, int cr, int delta, int ty)
{
if(cl<=l && r<=cr)
{
if(ty) tree[x].deltat+=delta;
else tree[x].deltatv+=tree[x].deltat*delta, tree[x].deltav+=delta;
return;
}
down(x);
int mid=(l+r)>>;
if(cl<=mid) update(x<<, l, mid, cl, cr, delta, ty);
if(cr>mid) update(x<<|, mid+, r, cl, cr, delta, ty);
}
ll query(int x, int l, int r, int cx)
{
if(l==r) return tree[x].deltat*tree[x].deltav-tree[x].deltatv;
down(x);
int mid=(l+r)>>;
if(cx<=mid) return query(x<<, l, mid, cx);
else return query(x<<|, mid+, r, cx);
}
inline void solve(int x, int y, int delta, int ty)
{
int f1=top[x], f2=top[y];
while(f1!=f2)
{
if(d[f1]<d[f2]) swap(x, y), swap(f1, f2);
update(, , n, w[f1], w[x], delta, ty);
x=fa[f1]; f1=top[x];
}
if(d[x]<d[y]) swap(x, y);
update(, , n, w[y], w[x], delta, ty);
}
int main()
{
fread(buf, , sizeof(buf), stdin); read(n);
for(int i=;i<=n;i++) read(fa[i]), add(fa[i], i);
dfs1(); dfs2(, ); read(m);
for(int i=;i<=m;i++) read(ty), read(x), read(y), solve(, x, y, ty-);
for(int i=;i<=n;i++) printf("%lld\n", query(, , n, w[i]));
}

  如果这题要求区间和呢?因为一个区间内可能有很多个阶梯形,但是上传的时候都是接上同样的一段父亲的阶梯,所以这个贡献还是可以计算的,只需要多记录一下区间内的$sumv,sumt$就很好求了。

最新文章

  1. IIS配置域用户自动登录
  2. Ubuntu14.04安装搜狗拼音输入法
  3. 区分LocalStorage和偏好数据
  4. 【随记】Hello World小记
  5. 实用技巧:使用 jQuery 异步加载 JavaScript 脚本
  6. 【转载】Java的四种引用
  7. abstract class与interface的区别与联系
  8. 微信小程序跳转页面
  9. placeholder改变颜色
  10. Blend在WPF开发过程中的作用
  11. mac下安装windows系统
  12. Java开发笔记(八)五种算术运算符
  13. C#.NET 大型通用信息化系统集成快速开发平台 4.1 版本 - 树形选择项目的标准例子
  14. Java二叉树的实现与特点
  15. 13. The Impact of New Technology on Business 新科技对商务的影响
  16. CSS外观属性
  17. WDA基础十三:常用模板管理
  18. 使用Visual Studio2013打开Visual Studio2015项目
  19. day40 数据结构-算法(二)
  20. 命名空间p方式的属性注入

热门文章

  1. vlc播放rtsp之怪事
  2. 大数据入门第十四天——Hbase详解(一)入门与安装配置
  3. Android病毒家族及行为(一)
  4. linux &amp; windows下重启oracle
  5. Python基础(条件判断和循环) if elif else for while break continue;
  6. Linux 设置core dump
  7. Selenium-ActionChainsApi接口详解
  8. 实验吧CTF天网管理系统
  9. 互评Final版本——二次元梦之队——“I Do”
  10. Leetcode题库——39.组合总和