[bzoj4003]城市攻占
倍增，对于每一个点计算他走到$2^i$次祖先所需要的攻击力以及最终会变成什么（一个一次函数），简单处理即可
（然而这样是错的，因为他只保证了骑士的攻击力可以存，并没有保证这个一次函数的系数可以存）
（其实还可以用科学记数法即pair<long double,int>来存即可，只要注意精度&常数）
正解是模拟，维护当前子树中骑士血量的左偏树（支持合并），然后考虑不断删除堆顶，修改可以用打标记来实现（因为乘的是正的，所以不改变顺序）
1 #include<bits/stdc++.h> 2 using namespace std; 3 #define N 300005 4 #define ll long long 5 struct ji{ 6 int nex,to,len; 7 }edge[N<<1]; 8 pair<ll,int>val[N]; 9 int E,n,m,x,r[N],dis[N],ls[N],rs[N],head[N],d[N],p[N],ans1[N],ans2[N]; 10 ll y,tc[N],tj[N],h[N],v[N]; 11 void upd(int k,ll x,ll y){ 12 tc[k]=tc[k]*x; 13 tj[k]=tj[k]*x+y; 14 val[k].first=val[k].first*x+y; 15 } 16 void down(int k){ 17 upd(ls[k],tc[k],tj[k]); 18 upd(rs[k],tc[k],tj[k]); 19 tc[k]=1; 20 tj[k]=0; 21 } 22 int merge(int x,int y){ 23 if ((!x)||(!y))return x+y; 24 down(x); 25 down(y); 26 if (val[x]>val[y])swap(x,y); 27 rs[x]=merge(rs[x],y); 28 if (dis[ls[x]]<dis[rs[x]])swap(ls[x],rs[x]); 29 dis[x]=dis[rs[x]]+1; 30 return x; 31 } 32 void add(int x,int y){ 33 edge[E].nex=head[x]; 34 edge[E].to=y; 35 head[x]=E++; 36 } 37 void dfs(int k,int sh){ 38 d[k]=sh; 39 for(int i=head[k];i!=-1;i=edge[i].nex){ 40 dfs(edge[i].to,sh+1); 41 r[k]=merge(r[k],r[edge[i].to]); 42 } 43 while ((r[k])&&(val[r[k]].first<h[k])){ 44 down(r[k]); 45 ans1[k]++; 46 ans2[r[k]]=d[val[r[k]].second]-d[k]; 47 r[k]=merge(ls[r[k]],rs[r[k]]); 48 } 49 if (!p[k])upd(r[k],1,v[k]); 50 else upd(r[k],v[k],0); 51 } 52 int main(){ 53 scanf("%d%d",&n,&m); 54 memset(head,-1,sizeof(head)); 55 for(int i=1;i<=n;i++)scanf("%lld",&h[i]); 56 for(int i=2;i<=n;i++){ 57 scanf("%d%d%lld",&x,&p[i],&v[i]); 58 add(x,i); 59 } 60 for(int i=1;i<=m;i++){ 61 scanf("%lld%d",&y,&x); 62 val[i]=make_pair(y,x); 63 tc[i]=1; 64 r[x]=merge(r[x],i); 65 } 66 dfs(1,0); 67 while (r[1]){ 68 ans2[r[1]]=d[val[r[1]].second]+1; 69 r[1]=merge(ls[r[1]],rs[r[1]]); 70 } 71 for(int i=1;i<=n;i++)printf("%d\n",ans1[i]); 72 for(int i=1;i<=m;i++)printf("%d\n",ans2[i]); 73 }
