百度之星之路径交

Problem Description
给定一棵n个点的树，以及m条路径，每次询问第L条到第R条路径的交集部分的长度（如果一条边同时出现在2条路径上，那么它属于路径的交集）。
Input
第一行一个数n（n<=500,000）
接下来n-1行，每行三个数x，y，z，表示一条从x到y并且长度为z的边
第n+1行一个数m（m<=500,000）
接下来m行，每行两个数u，v，表示一条从u到v的路径
接下来一行一个数Q，表示询问次数（Q<=500,000）
接下来Q行，每行两个数L和R
Output
Q行，每行一个数表示答案。
Sample Input
4
1 2 5
2 3 2
1 4 3
2
1 2
3 4
1
1 2
Sample Output
5

话说我题我很早以前就想bac了。。
然而上段时间写过一次，还没提交电脑就关机学校的电脑啊
于是就一直没有时间。。
今天终于回来了。。
其实问题就是怎么快速地求两条路径的交。。
至于多段的话，弄一个线段树就好了

设两条路径分别为a1-b1,a2-b2，它们的lca是c1,c2，求出a1,b1-a2,b2两两之间的lca，设分别为d1,d2,d3,d4，假设dep[d1]<=dep[d2]<=dep[d3]<=dep[d4]，dep[c1]<=dep[c2]。那么如果两条路径有交，当且仅当dep[d1]>=dep[c1]且dep[d4]>=dep[d3]>=dep[c2]（自己画画就知道了）。并且如果有交，那么交一定是d3-d4。 from CQzhangyu

其实用我的话来讲，就是对于两条路径，每条路径有两个点，一共就是四个点
两两之前求一个lca。。
要是有交的话，深度最深的两个lca就是交。
这个手画一下大概可以理解一下
然后判有没有交的地方我也有点迷。。
CQzhangyu那里说画一下。。
然而，呵呵

update8.26
话说这几天在军训。。站军姿的时候也算是想通了
大概就是分了几种情况讨论一下。。
自己yy一下就好了

就当做我路径交的一个板子吧
code：

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
using namespace std;
const int N=500010;
struct qq
{
    int x,y,z,last;
}e[N*2];int num,last[N];
struct path
{
    int a,b,c;
    path(){}
    path(int A,int B,int C) {a=A,b=B,c=C;}
}p[N];//表示路径
struct qt
{
    int l,r;
    int s1,s2;
    path c;
}s[N*2];
int n;
void init (int x,int y,int z)
{
    num++;
    e[num].x=x;e[num].y=y;e[num].z=z;
    e[num].last=last[x];
    last[x]=num;
}
int dep[N];//节点的深度
int len[N];//这个点到根的路径长度 
int f[N][21];
void dfs (int x,int fa)
{
    f[x][0]=fa;
    for (int u=1;u<=20;u++)
        if (f[x][u-1]!=0)
            f[x][u]=f[f[x][u-1]][u-1];
    for (int u=last[x];u!=-1;u=e[u].last)
    {
        int y=e[u].y;
        if (y==fa) continue;
        len[y]=len[x]+e[u].z;dep[y]=dep[x]+1;
        dfs(y,x);
    }
}
int m;
int lca (int x,int y)
{
    if (dep[x]>dep[y]) swap(x,y);
    for (int u=20;u>=0;u--)
        if (dep[f[y][u]]>=dep[x])
            y=f[y][u];
    if (x==y) return x;
    for (int u=20;u>=0;u--)
        if (f[x][u]!=f[y][u])
            x=f[x][u],y=f[y][u];
    return f[x][0];
}
int Num=0;
int cs[10];
bool cmp (int a,int b){return dep[a]<dep[b];}
path mix (path x,path y)
{
    if (x.c==0||y.c==0) return path(0,0,0);
    cs[1]=lca(x.a,y.a);cs[2]=lca(x.a,y.b);
    cs[3]=lca(x.b,y.a);cs[4]=lca(x.b,y.b);
    sort(cs+1,cs+1+4,cmp);
    int md=max(dep[x.c],dep[y.c]),nd=min(dep[x.c],dep[y.c]);
    if (dep[cs[1]]<nd||dep[cs[3]]<md) return path(0,0,0);
    else return path(cs[3],cs[4],lca(cs[3],cs[4]));
}
void bt (int l,int r)
{
    int a=++Num;
    s[a].l=l;s[a].r=r;
    if (l==r)
    {
        s[a].c=p[l];
        return ;
    }
    int mid=(l+r)>>1;
    s[a].s1=Num+1;bt(l,mid);
    s[a].s2=Num+1;bt(mid+1,r);
    s[a].c=mix(s[s[a].s1].c,s[s[a].s2].c);
}
path get (int now,int l,int r)
{
    if (s[now].l==l&&s[now].r==r)
        return s[now].c;
    int s1=s[now].s1,s2=s[now].s2;
    int mid=(s[now].l+s[now].r)>>1;
    if (r<=mid) return get(s1,l,r);
    else if (l>mid) return get(s2,l,r);
    else return mix(get(s1,l,mid),get(s2,mid+1,r));
}
int main()
{
    num=0;memset(last,-1,sizeof(last));
    scanf("%d",&n);
    for (int u=1;u<n;u++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        init(x,y,z);init(y,x,z);
    }
    len[0]=0;dep[1]=1;dfs(1,0);
/*  for (int u=1;u<=n;u++)
        printf("%d %d %d\n",dep[u],len[u],f[u][0]);*/
    scanf("%d",&m);
    for (int u=1;u<=m;u++)
    {
        scanf("%d%d",&p[u].a,&p[u].b);
        p[u].c=lca(p[u].a,p[u].b);
    }
/*  for (int u=1;u<=m;u++)
        printf("%d %d %d\n",p[u].a,p[u].b,p[u].c);*/
    bt(1,m);
    int q;
    scanf("%d",&q);
    for (int u=1;u<=q;u++)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        path ans=get(1,a,b);
        printf("%d\n",len[ans.a]+len[ans.b]-2*len[ans.c]);
    }
    return 0;
}