凸包+可持久化栈 Lost My Music

这玩意看上去好恶心。。。实际上并没有。

原来的式子是(c[top]-c[i])/(dep[i]-dep[top])—-> -(c[top]-c[i])/(dep[top]-dep[i]).这样就变成了斜率的形式了。要原来的式子最小，同理就是要现在的式子最大，并且这个斜率满足单调性。如果把斜率从小到大排序，然后放在一起，发现就成了一个下凸壳。考虑维护一个栈，dfs，每到一个点就退栈直到第一个满足的位置（也就是满足斜率单调的下一个位置），推入栈中，并用前一个点的算出ans。dfs完自己的子节点后把退出去的栈补回来。

明显很慢QAQ考虑优化。退栈过程可以从一个一个退—>二分。而且并没必要真的把栈全退掉，只要把相应的一位改成当前节点，dfs完子节点后再补回来。维护栈长度的话，就没有必要把后面的点真的退掉了。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <vector>
#define N 500005
using namespace std;
int n,f[N],a[N],dep[N],zhan[N];
double ans[N];
vector<int> son[N];
int get(int x,int h)
{
    int l=2,r=h,mid;double k1,k2;
    while(l<=r)
    {
        mid=l+r>>1;
        k1=-1.0*(a[zhan[mid]]-a[zhan[mid-1]])/(dep[zhan[mid]]-dep[zhan[mid-1]]);
        k2=-1.0*(a[x]-a[zhan[mid]])/(dep[x]-dep[zhan[mid]]);
        if(k1<k2)r=mid-1;
        else l=mid+1;
    }
    return l-1;
}
void dfs(int x,int h)
{
    dep[x]=dep[f[x]]+1;
    int k=get(x,h)+1,id=zhan[k],y=zhan[k-1];
    if(x==1)k=1;
    ans[x]=-1.0*(a[y]-a[x])/(dep[y]-dep[x]);
    zhan[k]=x;
    for(int i=0;i<son[x].size();i++)dfs(son[x][i],k);
    zhan[k]=id;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%d",&a[i]);
    for(int i=2;i<=n;i++)scanf("%d",&f[i]),son[f[i]].push_back(i);
    dfs(1,0);
    for(int i=2;i<=n;i++)printf("%.9lf\n",ans[i]);
}