BZOJ5077: [Ctsc2016]时空旅行（线段树+凸包）

传送门

题解：
首先答案为$\min {(x-x_i)^2+c_i}$，移项得$x^2 -2x_ix+x_i^2+c$
其实是给一堆直线和给定的$x$，求与直线最下方的交点。

那么显然是维护凸壳了，因为一个星球只对部分时空产生影响，那么我们在线段树上把这些直线丢进去，维护线段树上每个节点的凸壳。

注意先将所有直线排序然后依次加入，否则会多出排序的时间。

查询类似，先将查询排序然后每个凸壳维护一个指针，单调挪动就好了，时间复杂度$O(n\log n)$。

不过我用的vector，在bzoj上是卡着时间过去的，最好自己手写一个链表。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef long double LD;
typedef pair<int,int> pii;
inline LL rd() {
    char ch=getchar(); LL i=0,f=1;
    while(!isdigit(ch)) {if(ch=='-')f=-1; ch=getchar();}
    while(isdigit(ch)) {i=(i<<1)+(i<<3)+ch-'0'; ch=getchar();}
    return i*f;
}
const int N=5e5+50;
const LL INF=0x3f3f3f3f3f3f3f3f;
const LD eps=1e-8;
inline int sgn(LD x) {return (x>eps)-(x<-eps);}
int n,m,tot,a[N],p[N],vis[N],dfn[N],sze[N],ind,p2[N],nowpos[N*4];
LL c[N],ans[N];
pair < int,pair<int,int> > st[N];
set <pii> seg[N];
struct data {int op,id;}op[N];
vector <int> edge[N];
vector <int> cv[N*4];
typedef set <pii>::iterator it;
inline bool cmp(const int &x,const int &y) {
    return a[x]>a[y] || (a[x]==a[y] && c[x]<c[y]);
}
inline int dec(set <pii> &s,int l,int r) {
    it p=--s.upper_bound(make_pair(l,r));
    pii t=*p; s.erase(p);
    if(t.first!=l) s.insert(make_pair(t.first,l-1));
    if(r!=t.second) s.insert(make_pair(r+1,t.second)); 
}
inline int inc(set <pii> &s,int l,int r) {s.insert(make_pair(l,r));}
inline void dfs(int x,int f) {
    dfn[x]=++ind; sze[x]=1;
    for(int e=edge[x].size()-1;e>=0;e--) {
        int v=edge[x][e]; if(v==f) continue;
        dfs(v,x); sze[x]+=sze[v];
    }
} 
inline bool judge(int i,int j,int k) {
    LD x=(LD)(c[k]-c[i])/(a[i]-a[k]);
    LD yc1=a[k]*x+c[k],yc2=a[j]*x+c[j];
    return sgn(yc1-yc2)<=0;
}
inline void build_cv(vector<int> &cx,int l) {
    if(cx.size() && a[l]==a[cx[cx.size()-1]]) return;
    while(cx.size()>=2 && judge(cx[cx.size()-2],cx[cx.size()-1],l)) cx.pop_back();
    cx.push_back(l);
}
inline void inc(int k,int l,int r,int L,int R,int id) {
    if(L<=l&&r<=R) {build_cv(cv[k],id); return;}
    int mid=(l+r)>>1;
    if(R<=mid) inc(k<<1,l,mid,L,R,id);
    else if(L>mid) inc(k<<1|1,mid+1,r,L,R,id);
    else inc(k<<1,l,mid,L,R,id),inc(k<<1|1,mid+1,r,L,R,id);
}
inline LL calc(int p,int id) {return (LL)a[id]*p+c[id];}
inline void solve(int &pos,const vector <int> &cx,int p,int id) {
    if(!cx.size()) return;
    while(pos<cx.size()-1 && calc(p,cx[pos])>=calc(p,cx[pos+1])) ++pos;
    ans[id]=min(ans[id],calc(p,cx[pos]));
}
inline void inc2(int k,int l,int r,int p,int pos,int id) {
    solve(nowpos[k],cv[k],pos,id);
    if(l==r) return;
    int mid=(l+r)>>1;
    if(p<=mid) inc2(k<<1,l,mid,p,pos,id);
    else inc2(k<<1|1,mid+1,r,p,pos,id);
}
int main() {
    n=rd(), m=rd();
    memset(ans,0x3f,sizeof(ans));
    a[1]=0; c[1]=rd(); vis[1]=1;
    inc(seg[1],1,n);
    for(int i=2;i<=n;i++) {
        op[i].op=rd(); int f=rd()+1;
        edge[f].push_back(i);
        edge[i].push_back(f);
        op[i].id=rd()+1;
        if(!op[i].op) {
            vis[op[i].id]=1;
            int x=rd(); rd(); rd();
            a[op[i].id]=-2*x;
            c[op[i].id]=rd()+(LL)x*x;
        }
    }
    dfs(1,0);
    for(int i=2;i<=n;i++) {
        if(!op[i].op) inc(seg[op[i].id],dfn[i],dfn[i]+sze[i]-1);
        else dec(seg[op[i].id],dfn[i],dfn[i]+sze[i]-1);
    }
    for(int i=1;i<=n;i++) if(vis[i]) p[++tot]=i;
    sort(p+1,p+tot+1,cmp);
    for(int i=1;i<=tot;i++) {
        int x=p[i];
        for(it q=seg[x].begin();q!=seg[x].end();++q) 
            inc(1,1,n,q->first,q->second,x);
    }
    for(int i=1;i<=m;i++) 
        st[i].second.first=rd()+1, st[i].first=rd(), p2[i]=st[i].first, st[i].second.second=i;
    sort(st+1,st+m+1);
    for(int i=1;i<=m;i++)
        inc2(1,1,n,dfn[st[i].second.first],st[i].first,st[i].second.second);
    for(int i=1;i<=m;i++) 
        printf("%lld\n",ans[i]+(LL)p2[i]*p2[i]);
}