一道比较套路的数据结构题。而且写起来也挺顺的。
首先可以发现y和z就是来卖萌的,无视即可。那么考虑一个点(星球)位置为x0,费用为c0,那么询问x到它的总花费为(x-x0)^2+c0=-2x0*x+(c0+x0^2) +x^2,那么我们可以将一个点看做一条直线y=-2x0*x+(c0+x0^2),那么查询就相当于是求在当前空间中所有直线在x点的值的最小值。
直接斜率优化,相当于维护一个凸壳。求一遍dfs序之后可以发现一个星球存在的地方一定是一些在dfs序上连续的段,并且段的总数量是O(N)。那么考虑用线段树维护,每一个维护对应区间的一个凸壳。凸壳只需要按照直线的斜率排序在去求就好了。
查询的时候,按照x排序之后,在线段树上对应的区间像斜率优化的单调队列那样跑就好了。
思路是很清晰的,写起来也就不恶心辣≧.≦
AC代码如下:
#include<bits/stdc++.h>
#define ll long long
#define N 500005
#define M 1200005
#define cant(x,y,z) ((b[y]-b[x])*(a[z]-a[y])<=(b[z]-b[y])*(a[y]-a[x]))
using namespace std;
int n,m,cnt,tot,dfsclk,q[N],tk[N],tx[N],pnt[M],nxt[M],lf[N],rg[N],id[N];
int h[50000005],pos[N],head[M],tail[M],sz[M]; ll a[N],b[N],ans[N];
ll read(){
ll x=0; char cr=getchar(); bool flag=0;
while (cr<'0' || cr>'9'){ if (cr=='-') flag=1; cr=getchar(); }
while (cr>='0' && cr<='9'){ x=x*10+cr-'0'; cr=getchar(); }
return (flag)?-x:x;
}
struct graph{
int fst[N];
void add(int x,int y){
pnt[++tot]=y; nxt[tot]=fst[x]; fst[x]=tot;
}
}g0,g1;
void dfs(int x){
int i; lf[x]=++dfsclk;
for (i=g0.fst[x]; i; i=nxt[i]) dfs(pnt[i]);
rg[x]=dfsclk;
}
void mdy(int k,int l,int r,int x,int y){
if (x<=l && r<=y){ sz[k]++; return; }
int mid=l+r>>1;
if (x<=mid) mdy(k<<1,l,mid,x,y); if (y>mid) mdy(k<<1|1,mid+1,r,x,y);
}
void ins(int k,int l,int r,int x,int y,int z){
if (x<=l && r<=y){
int &i=tail[k];
for (; head[k]<i && cant(h[i-1],h[i],z); i--); h[++i]=z;
return;
}
int mid=l+r>>1;
if (x<=mid) ins(k<<1,l,mid,x,y,z); if (y>mid) ins(k<<1|1,mid+1,r,x,y,z);
}
ll qry(int k,int x){
int i; ll ans=1000000000000000000ll;
for (i=pos[k]; i; i>>=1){
int &j=head[i];
for (; j<tail[i] && (a[h[j+1]]-a[h[j]])*x<=b[h[j+1]]-b[h[j]]; j++);
if (j<=tail[i]) ans=min(ans,a[h[j]]*x-b[h[j]]);
}
return ans;
}
void build(int k,int l,int r){
head[k]=cnt+1; tail[k]=cnt; cnt+=sz[k];
if (l==r){ pos[l]=k; return; }
int mid=l+r>>1;
build(k<<1,l,mid); build(k<<1|1,mid+1,r);
}
bool cmp0(int x,int y){ return lf[x]<lf[y]; }
bool cmp1(int x,int y){ return a[x]<a[y] || a[x]==a[y] && b[x]<b[y]; }
bool cmp2(int x,int y){ return tx[x]>tx[y]; }
int main(){
int tt=clock();
n=read(); m=read();
b[1]=-read(); id[1]=1;
int i,j,k,l,op,x,y;
for (i=2; i<=n; i++){
op=read(); x=read()+1; g0.add(x,i);
y=read()+1;
if (op) g1.add(y,i); else{
id[y]=i;
a[y]=read(); b[y]=read(); b[y]=read(); b[y]=-read();
b[y]-=a[y]*a[y]; a[y]=-(a[y]<<1);
}
}
dfs(1);
for (i=1; i<=n; i++) if (id[i]){
for (j=g1.fst[i],cnt=0; j; j=nxt[j]) q[++cnt]=pnt[j];
sort(q+1,q+cnt+1,cmp0);
for (j=g1.fst[i],k=lf[id[i]],cnt=0; j; j=nxt[j],k=rg[x]+1){
x=pnt[j]=q[++cnt];
if (k<lf[x]) mdy(1,1,n,k,lf[x]-1);
}
if (k<=rg[id[i]]) mdy(1,1,n,k,rg[id[i]]);
}
cnt=0; build(1,1,n);
for (i=1; i<=n; i++) q[i]=i;
sort(q+1,q+n+1,cmp1);
for (l=1; l<=n; l++) if (id[i=q[l]]){
for (j=g1.fst[i],k=lf[id[i]]; j; j=nxt[j],k=rg[x]+1){
x=pnt[j];
if (k<lf[x]) ins(1,1,n,k,lf[x]-1,i);
}
if (k<=rg[id[i]]) ins(1,1,n,k,rg[id[i]],i);
}
for (i=1; i<=m; i++){
q[i]=i;
tk[i]=lf[read()+1]; tx[i]=read();
}
sort(q+1,q+m+1,cmp2);
for (l=1; l<=m; l++){
i=q[l];
ans[i]=qry(tk[i],tx[i])+(ll)tx[i]*tx[i];
}
for (i=1; i<=m; i++) printf("%lld\n",ans[i]);
return 0;
}
by lych
2016.12.8