[洛谷 P3787] 冰精冻西瓜

题目描述

琪露诺是拥有操纵冷气程度的能力的妖精，一天她发现了一片西瓜地。这里有n个西瓜，由n-1条西瓜蔓连接，形成一个有根树，琪露诺想要把它们冷冻起来慢慢吃。

这些西瓜蔓具有神奇的性质，可以将经过它的冷气的寒冷程度放大或缩小，每条西瓜蔓放大/缩小冷气寒冷程度的能力值为Wi，表示冷气经过它后，寒冷程度值x会变为x*wi。每个西瓜也有一个寒冷程度值，炎热的夏日，所有西瓜的寒冷程度值初始都为0。

琪露诺会做出两种动作:

①.对着西瓜i放出寒冷程度为x的冷气。这股冷气顺着西瓜蔓向“西瓜树”的叶子节点蔓延，冷气的寒冷程度会按照上面的规则变化。遇到一个西瓜连了多条西瓜蔓时，每条叶子节点方向的西瓜蔓均会获得与原先寒冷程度相等的冷气。途径的所有西瓜的寒冷程度值都会加上冷气的寒冷程度值。

⑨.向你询问西瓜i的寒冷程度值是多少。

等等，为什么会有⑨？因为笨蛋琪露诺自己也会忘记放了多少冰呢。

所以，帮她计算的任务就这么交给你啦。

输入输出格式

输入格式：

第一行一个整数n,表示西瓜的数量。

西瓜编号为1~n，1为这棵“西瓜树”的根。

接下来n-1行，每行有两个整数u,v和一个实数w，表示西瓜u和西瓜v之间连接有一条藤蔓，它放大/缩小冷气寒冷程度的能力值为w。

接下来一行一个整数m，表示操作的数量。

接下来m行，每行两个或三个整数。

第一个数只能是1或9。

如果为1，接下来一个整数i和一个实数x，表示对西瓜i放出寒冷程度为x的冷气。

如果为9，接下来一个整数i，表示询问编号为i的西瓜的寒冷程度值。

输出格式：

对于每个操作⑨，输出一行一个实数，表示对应西瓜的寒冷程度值。

输入输出样例

输入样例#1：

4
1 2 1.00000000
2 3 0.00000000
3 4 1.00000101
9
1 1 3.00000000
9 2
9 3
1 2 1.42856031
9 4
9 2
1 3 4.23333333
9 2
9 4

输出样例#1：

3.00000000
0.00000000
0.00000000
4.42856031
4.42856031
4.23333761

说明

子任务可能出现如下的特殊性质:

“西瓜树”退化为一条链

输入数据中的实数均保留8位小数，选手的答案被判作正确当且仅当输出与标准答案误差不超过10^-7。请特别注意浮点数精度问题。

实际数据中，冷气的寒冷程度x的范围为 [-0.1,0.1]

(样例中的冷气寒冷程度的范围为[1,5])

Hint:冷气只会向释放点的子树蔓延，不会向根方向蔓延

题目说的很清楚，这里就不解释了。（再说我也解释不来）

其实，这种题目第一眼看到就是这种类型的维护来维护去的。由于本题的特殊性，我们可以获得以下几点：

1.每当有1操作时，其实相当于从根节点出发，但是那冷气的初始寒冷值是要改变的。怎么改变？设从根节点到节点x路径上的所以w的积为dec，那么寒冷值就要变为原来的dec(倍)；

2.当出现w=0时，无法直接运用上述规律；

3.当w=0时，建一棵新树，使当前节点成为根节点，并与原来的树相对独立（可能有很多棵，请深入思考，怎样做到相对独立）；（Important！！！）

4.由于要运用到树状数组，所以我们要将有关联的节点放到一起来——怎么办？构造入栈时刻数组in和出栈时刻数组out。节点x影响到的将会是in[x]~out[x](in[x]<=out[x])，因此，我们在Update时，应该这么做：Update(in[x],dec),Update(out[x]+1,-dec)，相当于构造后缀和；同样，在Query时应该是Query(in[x])，但是，别忘了，这只是对于根节点来说的（尽管根节点不一定就得到那么多），我们还应该乘以w[x]（w[x]即为从根节点到节点x路径上的所以w的积）。

根据以上几点，不难构建出代码。

代码如下：

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<vector>
 #include<queue>
 #define M(Ar) memset(Ar,0,sizeof Ar);
 using namespace std;
 ;
 int n,Q,fa[maxn],cloc,incl[maxn],outcl[maxn];
 bool vis[maxn];
 double w[maxn],cold[maxn];
 vector <int> G[maxn];
 vector <double> A[maxn];
 queue <int> q;
 inline int read(){
     ; char ch=getchar();
     ') ch=getchar();
     +ch-',ch=getchar();
     return x;
 }
 inline void add(int x,int y,double z){
     G[x].push_back(y),A[x].push_back(z);
     G[y].push_back(x),A[y].push_back(z);
 }
 inline void DFS(int x){
     vis[x]=,incl[x]=++cloc; int sz=G[x].size();
     ; i<sz; i++)
     ) w[G[x][i]]=1.00000000*w[x]*A[x][i],DFS(G[x][i]);
     else if (!vis[G[x][i]]) q.push(G[x][i]);
     outcl[x]=cloc;
 }
 inline void Update(int x,double z){
     for (int i=x; i<=n; i+=i&(-i)) cold[i]+=z;
 }
 inline double Query(int x){
     ;
     for (int i=x; i; i-=i&(-i)) ret+=cold[i];
     return ret;
 }
 int main(){
     n=read();
     ; i<n; i++){
         int x=read(),y=read(); double z; scanf("%lf",&z);
         add(x,y,z);
     }
     q.push(),cloc=,M(vis);
     ,DFS(q.front()),q.pop();
     Q=read(),M(cold);
     for ( ; Q; Q--){
         int op=read();
         ){
             int x=read(); double z; scanf("%lf",&z);
             double dec=z/w[x];
             Update(incl[x],dec),Update(outcl[x]+,-dec);
         }else
         ){
             int x=read();
             printf("%.8lf\n",Query(incl[x])*w[x]);
         }
     }
     ;
 }