poj 3728 The merchant

题意：

求树上x->y的路径最大差值，且大的权在小的权之后。

题解：

挺好的一道经典陈题。
首先好像有种差分+树剖的做法，但太恶心。
扫一眼题解，发现要维护一大堆东西，于是YY了下，好像可以。
大概是这样的：

struct trnode{    int dis,fa[20],max[20],min[20],up[20],down[20];}tr[50010];

分别是，x的2i的祖先，这一段上最大值，最小值，从x->2ifa的答案，从2ifa->x的答案。
写考虑暴力做法，枚举路径，用当前值-已有的最小值更新答案。
那么现在可以倍增求这些东西，分两段，x->lca与lca->y。
我也不知道在口胡什么，具体代码
还有种离线的做法，然而我并不会，丢链接：ozy

总结：

想做这种题最好别想太多，变一变就好。
有时就差一点灵感。
剩下99%都是码代码。
但灵感更重要。
code：

#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>using namespace std;const int inf=(1<<28);struct node{    int y,next;}a[100010];int len=0,last[50010];int n,w[50010];struct trnode{    int dis,fa[20],max[20],min[20],up[20],down[20];}tr[50010];void ins(int x,int y){    a[++len].y=y;    a[len].next=last[x];last[x]=len;}void cmax(int &x,int y){x=max(x,y);}void cmin(int &x,int y){x=min(x,y);}void dfs(int x,int fa){    tr[x].dis=tr[fa].dis+1;    tr[x].fa[0]=fa;    if(fa!=0) cmax(tr[x].up[0],w[fa]-w[x]);    if(fa!=0) cmax(tr[x].down[0],w[x]-w[fa]);    tr[x].max[0]=tr[x].min[0]=w[x];    if(fa!=0) cmax(tr[x].max[0],w[fa]);    if(fa!=0) cmin(tr[x].min[0],w[fa]);    for(int i=1;(1<<i)<=tr[x].dis;i++)    {        int y=tr[x].fa[i-1];        tr[x].fa[i]=tr[y].fa[i-1];        cmax(tr[x].up[i],max(tr[x].up[i-1],tr[y].up[i-1]));        cmax(tr[x].up[i],tr[y].max[i-1]-tr[x].min[i-1]);        cmax(tr[x].down[i],max(tr[x].down[i-1],tr[y].down[i-1]));        cmax(tr[x].down[i],tr[x].max[i-1]-tr[y].min[i-1]);        tr[x].max[i]=max(tr[x].max[i-1],tr[y].max[i-1]);        tr[x].min[i]=min(tr[x].min[i-1],tr[y].min[i-1]);    }    for(int i=last[x];i;i=a[i].next)    {        int y=a[i].y;        if(y!=fa) dfs(y,x);    }}int getlca(int x,int y){    if(tr[x].dis<tr[y].dis) swap(x,y);    for(int i=18;i>=0;i--)        if((1<<i)<=tr[x].dis-tr[y].dis) x=tr[x].fa[i];    if(x==y) return x;    for(int i=18;i>=0;i--)        if((1<<i)<=tr[x].dis&&tr[x].fa[i]!=tr[y].fa[i]) x=tr[x].fa[i],y=tr[y].fa[i];    return tr[x].fa[0];}int solve(int x,int y){    int lca=getlca(x,y);    int ans=0;    int mx=w[y],mn=w[x];    for(int i=18;i>=0;i--)        if((1<<i)<=tr[x].dis-tr[lca].dis)        {            cmax(ans,tr[x].up[i]);            cmax(ans,tr[x].max[i]-mn);            cmin(mn,tr[x].min[i]);            x=tr[x].fa[i];        }    for(int i=18;i>=0;i--)        if((1<<i)<=tr[y].dis-tr[lca].dis)        {            cmax(ans,tr[y].down[i]);            cmax(ans,mx-tr[y].min[i]);            cmax(mx,tr[y].max[i]);            y=tr[y].fa[i];        }    cmax(ans,mx-mn);    return ans;}int main(){    scanf("%d",&n);    for(int i=1;i<=n;i++) scanf("%d",&w[i]);    for(int i=1;i<n;i++)    {        int x,y;scanf("%d %d",&x,&y);        ins(x,y);ins(y,x);    }    tr[0].dis=-1;dfs(1,0);    int q;scanf("%d",&q);    while(q--)    {        int x,y;scanf("%d %d",&x,&y);        printf("%d\n",solve(x,y));    }}