UOJ #7 【NOI2014】 购票

题目链接：购票

　　这道题我调了好久啊……主要还是因为这种用(CDQ)分治来搞斜率优化的题已经很久没写过了……上一次要追溯到去年暑假去了……

　　看下面这些东西之前你需要先自己推出斜率优化的式子……

　　这道题我们先来考虑一下如果不是树，在序列上怎么做。用(CDQ)分治的思想，先递归处理左半区间，然后用左半区间的值来更新右半区间，最后递归处理右半区间。这样的话，每次递归处理到一个区间，我们只需要考虑这个区间左边的点对右边的点的影响，对左边的点维护一个凸包，每次在凸包上二分最优解即可。

　　考虑如何把这个拓展到树上。我们树分治的时候每次是找了一个重心出来，然后递归处理剩下的部分。因此，我们在做这道题的时候，只需先递归处理剩下的子树中重心父亲所在的块，然后更新其他点的答案，最后递归其他联通块即可。注意重心需要单独考虑一下。

　　还有一个细节：先把所有祖先拿出来构凸包，然后每次在上面二分是不对的。因为这样可能会出现这种情况：对于一个点(x)，它的最优决策是(u)，在(fa_u)加入了凸包后点(u)会被从凸包中踢掉。但是如果(x)只能到(u)而到不了(fa_u)，就挂掉了……因为这里我调了好久……所以我们需要把点抠出来按能够到达的最小深度排序，然后把祖先一个一个往凸包里面加……

　　不知道为什么我的程序好慢……

　　下面贴代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)
#define maxn 200010
#define INF (1ll<<60)

using namespace std;
typedef long long llg;

int n,fa[maxn],p[maxn],T;
int head[maxn],next[maxn<<1],to[maxn<<1],tt;
llg d[maxn],f[maxn],q[maxn],L[maxn];

struct point{
	llg x,y;
	bool operator < (const point &h)const{return x>h.x;}
	point(llg a=0,llg b=0):x(a),y(b){}
}s[maxn];

struct Convex{//维护下凸包
	point s[maxn];int ls;
	double R(point a,point b){return (double)(b.y-a.y)/(b.x-a.x);}
	void push(point x){
		while(ls>1 && R(s[ls-1],s[ls])<=R(s[ls],x)) ls--;
		s[++ls]=x;
	}
	llg find(int x){//二分斜率
		if(!ls) return INF;
		int l=1,r=ls,mid;
		while(l<r){
			mid=(l+r)>>1;
			if(R(s[mid],s[mid+1])<x) r=mid;
			else l=mid+1;
		}
		return s[l].y-1ll*x*s[l].x;
	}
}tu;

void link(int x,int y){
	to[++tt]=y;next[tt]=head[x];head[x]=tt;
	to[++tt]=x;next[tt]=head[y];head[y]=tt;
}

int a[maxn],la,val[maxn],siz[maxn],dax[maxn];
bool vis[maxn];

void dfs(int u){
	vis[u]=1; a[++la]=u; siz[u]=1; dax[u]=0;
	for(int i=head[u],v;v=to[i],i;i=next[i])
		if(!vis[v]) dfs(v),siz[u]+=siz[v],dax[u]=max(dax[u],siz[v]);
	vis[u]=0;
}

void solve(int ff,int fr){
	la=0; dfs(ff); int k=a[1],ls=0;
	for(int i=1,u;u=a[i],i<=la;i++){
		val[u]=max(dax[u],siz[ff]-siz[u]);
		if(val[u]<val[k]) k=u;
	}
	
	vis[k]=1; if(!vis[fa[k]]) solve(fa[k],fr);//先处理父亲

	a[la=1]=k;
	for(int i=head[k],v;v=to[i],i;i=next[i])
		if(!vis[v] && v!=fa[k]) dfs(v);
	for(int i=1,u;u=a[i],i<=la;i++) s[i]=point(d[u]-L[u],u);
	ls=la; sort(s+1,s+ls+1); la=0;//把所有需要被更新的点抠出来排序
	
	for(int u=fa[k];!vis[u];u=fa[u]) a[++la]=u;
	if(fr) a[++la]=fr; tu.ls=0;//上一层的重心没有在上一层考虑，这一层要加进来

	for(int i=1,j=1,u;u=s[i].y,i<=ls;i++){//更新
		while(j<=la && d[a[j]]>=s[i].x) tu.push(point(d[a[j]],f[a[j]])),j++;
		f[u]=min(f[u],tu.find(p[u])+1ll*p[u]*d[u]+q[u]);
	}
	
	for(int i=head[k],v;v=to[i],i;i=next[i])
		if(!vis[v] && v!=fa[k]) solve(v,k);
	vis[k]=0;
}

int main(){
	File("a");
	scanf("%d %d",&n,&T);
	for(int i=2;i<=n;i++){
		scanf("%d %lld",&fa[i],&d[i]);
		scanf("%d %lld %lld",&p[i],&q[i],&L[i]);
		d[i]+=d[fa[i]]; link(i,fa[i]); f[i]=INF;
	}
	vis[0]=1; solve(1,0);
	for(int i=2;i<=n;i++) printf("%lld
",f[i]);
	return 0;
}

 　　BZOJ提交链接：BZOJ 3672 购票