2018-8-10 模拟赛T3(可持久化线段树)

出题人说：正解离线按DFS序排序线段维护区间和

但是对于树上每个点都有一个区间和一个值，两个点之间求1~m的区间和，这不就是用可持久化线段树吗。

只不过这个线段树需要区间修改，不过不需要标记下传，询问时加起来就好了。

对于每一个节点x，建一个1~m的线段树版本

询问时，先求出u和v的lca和lca的父亲flca

询问在{u+v}和{lca,flca}的差集中的区间和就好了

可持久化数据结构耗费的空间是巨大的^_^

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 typedef long long lnt;
 struct pnt{
     int hd;
     int a;
     int b;
     int oul;
     int fa;
     int dp;
     lnt w;
     int root;
 }p[];
 struct ent{
     int twd;
     int lst;
 }e[];
 struct tree{
     int l,r;
     lnt val;
     lnt laz;
 };
 struct Traos{
     tree tr[];
     int siz;
     int ks;
     lnt ansl;
     void build(int l,int r,int &spc)
     {
         if(!spc)
             spc=++siz;
         if(l==r)
             return ;
         int md=(l+r)/;
         build(l,md,tr[spc].l);
         build(md+,r,tr[spc].r);
     }
     void updte(int &spc,int last,int ll,int rr,int l,int r,lnt v)
     {
         if(l>rr||r<ll)
             return ;
         spc=++siz;
         tr[spc]=tr[last];
         if(ll<=l&&rr>=r)
         {
             tr[spc].val+=(lnt)(r-l+)*v;
             tr[spc].laz+=v;
             return ;
         }
         tr[spc].val+=(lnt)(min(rr,r)-max(l,ll)+)*v;
         int mid=(l+r)/;
         updte(tr[spc].l,tr[last].l,ll,rr,l,mid,v);
         updte(tr[spc].r,tr[last].r,ll,rr,mid+,r,v);
         return ;
     }
     lnt sumls(int spc1,int spc2,int ll,int rr,int l,int r)
     {
         if(l>rr||ll>r)
             return 0ll;
         if(ll<=l&&rr>=r)
             return (lnt)(tr[spc2].val-tr[spc1].val);
         int mid=(l+r)/;
         return (lnt)((lnt)(min(rr,r)-max(ll,l)+)*(tr[spc2].laz-tr[spc1].laz))+sumls(tr[spc1].l,tr[spc2].l,ll,rr,l,mid)+sumls(tr[spc1].r,tr[spc2].r,ll,rr,mid+,r);

     }
 }T;
 int cnt;
 int ont;
 int n,m,q;
 int ola[][];
 int rt[];
 int lg[];
 void ade(int f,int t)
 {
     cnt++;
     e[cnt].lst=p[f].hd;
     e[cnt].twd=t;
     p[f].hd=cnt;
 }
 void dfs(int x,int f)
 {
     p[x].fa=f;
     p[x].dp=p[f].dp+;
     p[x].oul=++ont;
     ola[][ont]=x;
     for(int i=p[x].hd;i;i=e[i].lst)
     {
         int to=e[i].twd;
         if(to!=f)
         {
             dfs(to,x);
             ola[][++ont]=x;
         }
     }
 }
 int mxs(int a,int b)
 {
     return p[a].dp<p[b].dp?a:b;
 }
 void kls()
 {
     for(int i=;i<=;i++)
     {
         for(int j=;j+(<<i)-<=ont;j++)
         {
             ola[i][j]=mxs(ola[i-][j],ola[i-][j+(<<(i-))]);
         }
     }
 }
 int lca(int a,int b)
 {
     if(p[a].oul>p[b].oul)
         swap(a,b);
     int lgg=lg[p[b].oul-p[a].oul+];
     return mxs(ola[lgg][p[a].oul],ola[lgg][p[b].oul-(<<lgg)+]);
 }
 void fdfs(int x,int f)
 {
     p[x].root=++T.ks;
     T.updte(rt[p[x].root],rt[p[f].root],p[x].a,p[x].b,,m,p[x].w);
     for(int i=p[x].hd;i;i=e[i].lst)
     {
         int to=e[i].twd;
         if(to-f)
         {
             fdfs(to,x);
         }
     }
 }
 int main()
 {
     freopen("c.in","r",stdin);
     freopen("c.out","w",stdout);
     scanf("%d%d%d",&n,&m,&q);
     for(int i=;i<=*n;i++)
     {
         lg[i]=lg[i/]+;
     }
     for(int i=;i<n;i++)
     {
         int x,y;
         scanf("%d%d",&x,&y);
         ade(x,y);
         ade(y,x);
     }
     for(int i=;i<=n;i++)
     {
         scanf("%d%d%d",&p[i].a,&p[i].b,&p[i].w);
         p[i].a=min(p[i].a,p[i].b);
         p[i].b=max(p[i].a,p[i].b);
     }
     dfs(,);
     kls();
     T.build(,m,rt[]);
     fdfs(,);
     for(int i=;i<=q;i++)
     {
         int u,v,l,r;
         scanf("%d%d%d%d",&u,&v,&l,&r);
         l=min(l,r);
         r=max(l,r);
         int la=lca(u,v);
         lnt ans=;
         int fl=p[la].fa;
         ans=T.sumls(rt[p[la].root],rt[p[u].root],l,r,,m);
         ans+=T.sumls(rt[p[fl].root],rt[p[v].root],l,r,,m);
         printf("%lld\n",ans);
     }
     return ;
 }