codevs 3305 水果姐逛水果街Ⅱ 题解报告

今天，栋栋大神又在做题了，宝宝看他一直没做出来，就来凑凑热闹，。。
果然，，我也做了好久。
快仨小时了都；；

先上题目：
题目描述 Description
水果姐第二天心情也很不错，又来逛水果街。

突然，cgh又出现了。cgh施展了魔法，水果街变成了树结构（店与店之间只有一条唯一的路径）。

同样还是n家水果店，编号为1~n，每家店能买水果也能卖水果，并且同一家店卖与买的价格一样。

cgh给出m个问题，每个问题要求水果姐从第x家店出发到第y家店，途中只能选一家店买一个水果，然后选一家店（可以是同一家店，但不能往回走）卖出去。求最多可以赚多少钱。

水果姐向学过oi的你求助。

输入描述 Input Description
第一行n，表示有n家店

下来n个正整数，表示每家店一个苹果的价格。

下来n-1行，每行两个整数x，y，表示第x家店和第y家店有一条边。

下来一个整数m，表示下来有m个询问。

下来有m行，每行两个整数x和y，表示从第x家店出发到第y家店。

输出描述 Output Description
有m行。

每行对应一个询问，一个整数，表示面对cgh的每次询问，水果姐最多可以赚到多少钱。

样例输入 Sample Input
10
16 5 1 15 15 1 8 9 9 15
1 2
1 3
2 4
2 5
2 6
6 7
4 8
1 9
1 10
6
9 1
5 1
1 7
3 3
1 1
3 6

样例输出 Sample Output
7
11
7
0
0
15

N<=200000 m<=10000

和T1要求的内容一样，不过不是线性的图，而是一颗树；
考虑LCA 倍增。
在求两点LCA的同时处理出最大利益；
对于这棵树，维护j[x][y] 为从点x向上跳2^y 距离的点是谁；
ma[N][20],mi[N][20],lr[N][20],rl[N][20]；
分别为 区间最大 区间最小 从左到右最大利润 从右到左最大利润
递推式为：
lr[x][i]=max(ma[j[x][i-1]][i-1]-mi[x][i-1],max(lr[x][i-1],lr[j[x][i-1]][i-1]));
rl[x][i]=max(ma[x][i-1]-mi[j[x][i-1]][i-1],max(rl[x][i-1],rl[j[x][i-1]][i-1]));
然后dfs

void dfs(int x,int s)
{
    deep[x]=s;
    for(int i=1;(1<<i)<s;i++)
    {
        j[x][i]=j[j[x][i-1]][i-1];
        ma[x][i]=max(ma[x][i-1],ma[j[x][i-1]][i-1]);
        mi[x][i]=min(mi[x][i-1],mi[j[x][i-1]][i-1]);
        lr[x][i]=max(ma[j[x][i-1]][i-1]-mi[x][i-1],max(lr[x][i-1],lr[j[x][i-1]][i-1]));
        rl[x][i]=max(ma[x][i-1]-mi[j[x][i-1]][i-1],max(rl[x][i-1],rl[j[x][i-1]][i-1]));
    }
    int k=h[x];
    while(k)
    {
        if(vv[v[k]])
        {
            k=nex[k];
            continue;
        }
            vv[v[k]]=1;
            j[v[k]][0]=x;
            ma[v[k]][0]=max(a[x],a[v[k]]);
            mi[v[k]][0]=min(a[x],a[v[k]]);
            lr[v[k]][0]=a[x]-a[v[k]];
            rl[v[k]][0]=a[v[k]]-a[x];
            dfs(v[k],s+1);
            k=nex[k];
    }
}

dfs后进行倍增求LCA
找一条从X到Y的道路，维护X-LCA(X,Y)这一条链上从左到右最大利益和LCA(X,Y)-Y这一条链上从右到左最大利益；
最后再考虑 点LCA(X,Y)和两侧最大最小值的差是否比求得两侧最大利益大。
四者取max
即为答案；

代码：：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<stack>
#include<cstdlib>
#include<string>
#include<bitset>
#include<iomanip>
#include<deque>
#define INF 1000000000
#define fi first
#define se second
#define N 1000005
#define P 1000000007
#define debug(x) cerr<<#x<<"="<<x<<endl
#define MP(x,y) make_pair(x,y)
using namespace std;
int n,m;
int a[200010],p=0,nex[N*2],h[N*2],v[N*2],j[N][20],ma[N][20],mi[N][20],lr[N][20],rl[N][20],deep[N];
void add(int aa,int b)
{
    p++;
    v[p]=b;
    nex[p]=h[aa];
    h[aa]=p;
}
bool vv[N];
inline int get_num()
{
int num = 0;
char c;
bool flag = false;
while ((c = getchar()) == ' ' || c == '\n' || c == '\r');
if (c == '-') flag = true;
else num = c - '0';
while (isdigit(c = getchar()))
num = num * 10 + c - '0';
return (flag ? -1 : 1) * num;
}
void dfs(int x,int s)
{
    deep[x]=s;
    for(int i=1;(1<<i)<s;i++)
    {
        j[x][i]=j[j[x][i-1]][i-1];
        ma[x][i]=max(ma[x][i-1],ma[j[x][i-1]][i-1]);
        mi[x][i]=min(mi[x][i-1],mi[j[x][i-1]][i-1]);
        lr[x][i]=max(ma[j[x][i-1]][i-1]-mi[x][i-1],max(lr[x][i-1],lr[j[x][i-1]][i-1]));
        rl[x][i]=max(ma[x][i-1]-mi[j[x][i-1]][i-1],max(rl[x][i-1],rl[j[x][i-1]][i-1]));
    }
    int k=h[x];
    while(k)
    {
        if(vv[v[k]])
        {
            k=nex[k];
            continue;
        }
            vv[v[k]]=1;
            j[v[k]][0]=x;
            ma[v[k]][0]=max(a[x],a[v[k]]);
            mi[v[k]][0]=min(a[x],a[v[k]]);
            lr[v[k]][0]=a[x]-a[v[k]];
            rl[v[k]][0]=a[v[k]]-a[x];
            dfs(v[k],s+1);
            k=nex[k];
    }
}
int just_doit(int x,int y)
{
    int ans=0,xmi=a[x],yma=a[y];
    if(x==y)
    {
        return 0;
    }
    if(deep[x]>deep[y])
    {
        int mii=a[x],maa=0;
        for(int i=15;i>=0;i--)
        {
        if(deep[x]-(1<<i)>=deep[y])
        {
            ans=max(ans,max(lr[x][i],ma[x][i]-xmi));
            xmi=min(xmi,mi[x][i]);
            x=j[x][i];
        }
        }
    }else
    {
        if(deep[x]<deep[y])
        {
        int mii=2147483647,maa=a[y];
        for(int i=15;i>=0;i--)
        {
        if(deep[y]-(1<<i)>=deep[x])
        {
            ans=max(ans,max(rl[y][i],yma-mi[y][i]));
            yma=max(yma,ma[y][i]);
            y=j[y][i];
        }
        }
        }
        else
        {
            xmi=a[x];
            yma=a[y];
        }
    }
    if(x==y)
    {
        return ans;
    }
    int may=0,mix=a[x];
    for(int i=15;i>=0;i--)
    {
        if(deep[x]-(1<<i)>=1&&j[x][i]!=j[y][i])
        {
            ans=max(ans,max(lr[x][i],ma[x][i]-xmi));
            ans=max(ans,max(rl[y][i],yma-mi[y][i]));
            xmi=min(xmi,mi[x][i]);
            yma=max(yma,ma[y][i]);
            x=j[x][i];
            y=j[y][i];
        }
    }
    ans=max(ans,max(yma-xmi,max(a[j[y][0]]-xmi,yma-a[j[x][0]])));
    //刚开始就卡在这里一个多小时，注意！！
    return ans;
}
int main()
{
    n=get_num();
    for(int i=1;i<=n;i++)
    {
        a[i]=get_num();
    }
    for(int i=1;i<n;i++)
    {
        int q,w;
        q=get_num();w=get_num();
        add(q,w);
        add(w,q);
    }
    vv[1]=1;
    dfs(1,1);
    //cout<<rl[6][1];
    cin>>m;
    for(int i=1;i<=m;i++)
    {
        int q,w;
        q=get_num();w=get_num();
        cout<<just_doit(q,w)<<"\n";
    }

}

噫。~~