洛谷P3676 小清新数据结构题（动态点分治+树链剖分）

传送门

　　感觉这题做下来心态有点崩……$RMQ$求$LCA$没有树剖快我可以理解为是常数太大……然而我明明用了自以为不会退化的点分然而为什么比会退化的点分跑得反而更慢啊啊啊啊~~~

　　先膜一波zsy大佬

　　讲讲做法。题目的要求是给定一个根$p$，求$\sum _{i=1}^ns_i^2$，其中$s_i$表示子树中的点权和

　　我们设$sum=\sum _{i=1}^n val_i$，即整棵树的点权和。先考虑一下$\sum _{i=1}^ns_i$怎么求。考虑一下每一个点的贡献，每一个点都会对被计算$dep_i+1$次（其中$dep_i$表示$dist(i,p)$），那么很显然$\sum _{i=1}^ns_i=\sum_{i=1}^nval_i*(dep_i+1)=\sum_{i=1}^nval_i*dep_i+sum$。然后考虑一下$val_i*dep_i$，如何动态维护？->幻想乡战略游戏……

　　简单来说，就是建好点分树，然后每一次及时修改和查询

　　然后我们令$calc(p)$以$p$为根时的$\sum _{i=1}^nval_i*dep_i$

　　然后考虑如下式子$$\sum_{i=1}^n\sum_{j=1}^nval_i*val_j*dist(i,j)$$

　　是不是可以理解为在所有的点对$(i,j)$之间的所有边上加上权值$val_i*val_j$（刚好有$dist(i,j)$条边），然后再求整棵树的权值？

　　然后我们考虑一下每条边的权值，肯定等于两侧的子树点权和的乘积。那么，不论是以哪一个点$p$为根，它的权值都等于$s_i*(sum-s_i)$，其中$s_i$表示这条边指向的儿子的子树的点权和

　　那么，上面的式子就可以变成这样$$\sum_{i=1}^n\sum_{j=1}^nval_i*val_j*dist(i,j)=\sum_{i=1}^ns_i*(sum-s_i)$$

　　又因为上式左边是不变的，所以不管选取哪一个$p$为根，右边都是不变的

　　令$W=\sum_{i=1}^ns_i*(sum-s_i)$，然后可以直接$O(n)dp$出$W$，然后考虑对点的修改对$W$造成的影响

　　$W=\sum_{i=1}^n\sum{j=1}^nval_i*val_j*dist(i,j)$，设点$u$的变化量为$Δv$，那么$ΔW=Δv*\sum_{j=1}^nval_j*dist(i,j)$，相当于$Δv*calc(i)$，然后可以考虑和一般的动态点分一样计算

　　然后最后询问的答案就是$$W=\sum_{i=1}^ns_i*(sum-s_i)$$

　　$$\sum_{i=1}^ns_i^2=\sum_{i=1}^ns_i*sum-W$$

　　$$\sum_{i=1}^ns_i^2=sum(calc(i)+sum)-W$$

 // luogu-judger-enable-o2
 //minamoto
 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #define ll long long
 using namespace std;
 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
 char buf[<<],*p1=buf,*p2=buf;
 template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,:;}
 inline int read(){
     #define num ch-'0'
     char ch;bool flag=;int res;
     while(!isdigit(ch=getc()))
     (ch=='-')&&(flag=true);
     for(res=num;isdigit(ch=getc());res=res*+num);
     (flag)&&(res=-res);
     #undef num
     return res;
 }
 char sr[<<],z[];int C=-,Z;
 inline void Ot(){fwrite(sr,,C+,stdout),C=-;}
 inline void print(ll x){
     if(C><<)Ot();if(x<)sr[++C]=,x=-x;
     while(z[++Z]=x%+,x/=);
     while(sr[++C]=z[Z],--Z);sr[++C]='\n';
 }
 const int N=;
 int ver[N<<],head[N],Next[N<<];
 int val[N],fa[N],pa[N],d[N],sz[N],son[N],top[N];
 int n,q,tot;
 void dfs1(int u,int fa){
     pa[u]=fa,d[u]=d[fa]+,sz[u]=;
     for(int i=head[u];i;i=Next[i]){
         int v=ver[i];if(v==fa) continue;
         dfs1(v,u);
         sz[u]+=sz[v];if(sz[v]>sz[son[u]]) son[u]=v;
     }
 }
 void dfs2(int u,int fa){
     top[u]=fa;
     if(son[u]) dfs2(son[u],fa);else return;
     for(int i=head[u];i;i=Next[i])
     if(ver[i]!=pa[u]&&ver[i]!=son[u])
     dfs2(ver[i],ver[i]);
 }
 int LCA(int u,int v){
     while(top[u]^top[v]){
         if(d[top[u]]<d[top[v]]) swap(u,v);
         u=pa[top[u]];
     }
     return d[u]<d[v]?u:v;
 }
 int dis(int u,int v){return d[u]+d[v]-(d[LCA(u,v)]<<);}
 int size,rt,vis[N];
 ll sum[N],sum1[N],sum2[N],sigma,omega,ans;
 void getrt(int u,int fa){
     sz[u]=,son[u]=;
     for(int i=head[u];i;i=Next[i]){
         int v=ver[i];if(v==fa||vis[v]) continue;
         getrt(v,u),sz[u]+=sz[v],cmax(son[u],sz[v]);
     }
     cmax(son[u],size-sz[u]);
     if(son[u]<son[rt]) rt=u;
 }
 void solve(int u,int f){
     fa[u]=f,vis[u]=;int totsz=size;
     for(int i=head[u];i;i=Next[i]){
         int v=ver[i];
         if(!vis[v]){
             size=sz[v]>sz[rt]?totsz-sz[rt]:sz[v];
             rt=;
             getrt(v,),solve(rt,u);
         }
     }
 }
 inline void modify(int u,int v){
     sum[u]+=v;
     for(int i=u;fa[i];i=fa[i]){
         int dist=dis(u,fa[i]);
         sum[fa[i]]+=v;
         sum1[fa[i]]+=dist*v;
         sum2[i]+=dist*v;
     }
 }
 inline ll calc(int u){
     ll res=sum1[u];
     for(int i=u;fa[i];i=fa[i]){
         int dist=dis(fa[i],u);
         res+=(ll)dist*(sum[fa[i]]-sum[i]);
         res+=sum1[fa[i]]-sum2[i];
     }
     return res;
 }
 void DP(int u,int fa){
     sz[u]=val[u];
     for(int i=head[u];i;i=Next[i]){
         int v=ver[i];
         if(v!=fa) DP(v,u),sz[u]+=sz[v];
     }
     omega+=1ll*sz[u]*(sigma-sz[u]);
 }
 int main(){
     n=read(),q=read();
     for(int i=;i<n;++i){
         int u=read(),v=read();
         ver[++tot]=v,Next[tot]=head[u],head[u]=tot;
         ver[++tot]=u,Next[tot]=head[v],head[v]=tot;
     }
     dfs1(,),dfs2(,);
     size=n,son[rt=]=n+;
     getrt(,),solve(rt,);
     for(int i=;i<=n;++i)
     val[i]=read(),modify(i,val[i]),sigma+=val[i];
     DP(,);
     while(q--){
         int opt=read(),x=read();
         if(opt&){
             int y=read();y-=val[x];
             modify(x,y),sigma+=y,omega+=y*calc(x);
             val[x]+=y;
         }
         else print((calc(x)+sigma)*sigma-omega);
     }
     Ot();
     return ;
 }