今天,栋栋大神又在做题了,宝宝看他一直没做出来,就来凑凑热闹,。。
果然,,我也做了好久。
快仨小时了都;;
先上题目:
题目描述 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";
}
}
噫。~~