uoj 198: [CTSC2016]时空旅行

       一道比较套路的数据结构题。而且写起来也挺顺的。

       首先可以发现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