[NOI2014] 购票 题解

合集 - 题解(82)

1.[JSOI2008]火星人 题解2024-01-192.[BZOJ3786] 星系探索 题解2024-01-193.[ZJOI2015]幻想乡战略游戏 题解2024-02-034.[SDOI2015] 寻宝游戏2024-02-205.[HNOI2011] 数学作业 题解2024-03-116.[HDU5396] Expression 题解2024-03-237.[NOIP2018] 旅行 题解2024-04-148.[POI2012] Rendezvous 题解2024-04-149.[BZOJ3037] 创世纪 题解2024-04-2010.[HEOI2014]大工程 题解2024-04-2011.[ZJOI2019] 语言 题解2024-04-2112.[国家集训队] 矩阵乘法 题解2024-05-0413.[THUPC2017] 天天爱射击 题解2024-05-0414.[SCOI2007] 蜥蜴 题解2024-05-0515.[SDOI2015] 星际战争 题解2024-05-0516.[国家集训队] happiness 题解2024-05-0917.[SDOI2016] 数字配对 题解2024-05-1918.[SDOI2009] 晨跑 题解2024-05-1919.[SCOI2007] 修车2024-05-1920.acwing329 围栏障碍训练场 题解2024-06-0421.无限之环 题解2024-06-0822.[BZOJ4350] 括号序列再战猪猪侠 题解2024-07-0823.[TJOI2015] 弦论 题解2024-07-1024.[SDOI2008] Sandy的卡片 题解2024-07-1025.[AHOI2013] 差异 题解2024-07-1026.[Ynoi2016] 镜中的昆虫 题解2024-08-1527.[JOISC 2023 Day3] Tourism 题解2024-08-2128.[PA2021] Od deski do deski 题解2024-10-1429.[TJOI2019] 甲苯先生的字符串 题解2024-10-1430.[POI2014] HOT-Hotels 加强版题解2024-10-1631.[Ynoi2015] 盼君勿忘 题解2024-10-2532.[luogu2123] 皇后游戏2024-11-1033.[JXOI2017] 加法 题解2024-11-1034.[luogu1248] 加工生产调度 题解2024-11-1035.[NOIP2018] 赛道修建2024-11-1036.[Ynoi2015] 我回来了 题解2024-11-2937.[HNOI2009] 图的同构计数2024-12-1438.[BZOJ3811] 玛里苟斯 题解2024-12-1439.[BZOJ3569] DZY Loves Chinese II 题解2024-12-1540.[SCOI2016] 幸运数字 题解2024-12-1641.[HDU5603] the soldier of love 题解2024-12-1842.[QOJ8672][PKUSC2024] 排队2024-12-1843.[BZOJ3489] A simple rmq problem2024-12-1944.[HAOI2010] 软件安装 题解2024-12-1945.[IOI2020] 连接擎天树 题解2024-12-2046.[LOJ6669] Nauuo and Binary Tree 题解2024-12-2047.[JOISC2019] 聚会 题解2024-12-2048.[CERC2014] Parades 题解2024-12-2049.[SHOI2017] 摧毁“树状图”2024-12-21

50.[NOI2014] 购票 题解2024-12-23

51.[WC2018] 通道 题解2024-12-2452.[BZOJ2741][FOTILE模拟赛] L 题解2024-12-2453.[THUSC2015] 异或运算 题解2024-12-2454.[BZOJ4771] 七彩树 题解2024-12-2555.[BZOJ4605] 崂山白花蛇草水 题解2024-12-2956.[BZOJ3600] 没有人的算术 题解2024-12-3057.[WC2014] 紫荆花之恋 题解01-0558.[SDOI2008] 洞穴勘测 题解01-0659.[POJ3237] 树的维护 题解01-0660.[国家集训队] Tree2 题解01-0861.[WC2006] 水管局长 题解01-0862.[BZOJ3514] [Codechef MARCH14] GERALD07加强版 题解01-0863.[luogu4114] Qtree1 题解01-0864.[THUWC2017] 在美妙的数学王国中畅游 题解（内附求导小技巧）01-0965.[BZOJ3159] 决战 题解01-0966.[BZOJ2194] 快速傅立叶之二 题解01-1867.[ZJOI2014] 力 题解01-1868.[BZOJ3451] Normal 题解01-1869.[BZOJ3771] Triple 题解01-1970.[BZOJ3160] 万径人踪灭 题解01-2071.[SDOI2015] 序列统计 题解01-2172.[联合省选 2020A] 组合数问题 题解01-2373.[HDU4625] JZPTREE+[国家集训队] Crash 的文明世界 题解01-2374.[TJOI/HEOI2016] 求和 题解01-2375.[BZOJ5093] 图的价值 题解01-2376.[FJOI2016] 建筑师 题解01-2377.[BZOJ4665] 小w的喜糖 题解01-2378.[BZOJ3622] 已经没有什么好害怕的了 题解01-2379.[BZOJ4671] 异或图 题解01-2380.[BZOJ4833] 最小公倍佩尔数 题解01-2381.[BalticOI 2022] Uplifting Excursion (Day1) 题解02-2182.[JLOI2016] 成绩比较 题解02-25

收起

首先发现 \(p_x\times dis(x,y)+q_x\) 异常像是能斜率优化的样子，那先把求 \(f_x\) 的式子写出来（下设 \(d_x\) 表示 \(x\) 到根的距离）：

\[f_x=\min_{lca(x,y)=y,y\ne x}(p_x\times(d_x-d_y)+q_x+f_y) \]

提出公共部分 \(p_x\times d_x+q_x\)，得：

\[f_x=p_x\times d_x+q_x+\min_{lca(x,y)=y,y\ne x}(f_y-p_x\times d_y) \]

那当 \(y\) 优于 \(z\) 时（设 \(d_y<d_z\)），一定满足：

\[f_y-p_x\times d_y<f_z-p_x\times d_z \]

\[f_y-f_z<p_x\times(d_y-d_z) \]

\[\dfrac{f_y-f_z}{d_y-d_z}>p_x \]

那这里就形成了一个斜率优化 \(dp\) 的式子，得用单调栈维护凸包。

但是这里就出现了一个很现实的问题：实际上你要维护的是多个滑动窗口内的凸包，也就意味着你原来淘汰的点可能在后面的决策中再次出现！

问题转化为维护后缀凸包，考虑树状数组维护单调栈，这样就能在 \(O(n\log^2n)\) 的时间中解决问题。

最后发现由于是在树上，所以要支持可撤销，但是明显的，单调栈的时间复杂度是均摊的，单次修改仍然有可能达到 \(O(n)\)，所以单调队列插入时，应当使用二分找到插入点，保证单次时间复杂度 \(O(\log n)\)。

时间复杂度 \(O(n\log^2n)\)。注意单调栈要用 \(vector\)。

#include<bits/stdc++.h>
#define int long long
#define db double
using namespace std;
const int N=2e5+5,M=4e6+5;
struct del{
	int id,tp,ps;
}sk[M];vector<int>g[N];
int n,m,idx[N],ln[N],ft;
int d[N],p[N],q[N],f[N];
struct mstack{
	int tp=-1;vector<int>st;
	db sp(int x,int y){
		return 1.0*(f[x]-f[y])/(d[x]-d[y]);
	}int cmp(int x,int y,int z){
		return sp(x,y)>=sp(y,z);
	}void add(int x,int id){
		int l=1,r=tp,ans=tp+1;
		while(l<=r){
			int mid=(l+r)/2;
			if(cmp(st[mid-1],st[mid],x))
				ans=mid,r=mid-1;
			else l=mid+1;
		}if(ans>tp) st.push_back(0);
		sk[++ft]={id,tp,st[ans]},st[tp=ans]=x; 
	}int ans(int x){
		if(tp<0) return (int)9e18;
		int l=0,r=tp-1,ans=st[tp];
		while(l<=r){
			int mid=(l+r)/2;
			if(sp(st[mid],st[mid+1])>x)
				ans=st[mid],r=mid-1;
			else l=mid+1;
		}return f[ans]-d[ans]*x;
	}
};namespace BIT{
	mstack c[N];
	void clear(int lft){
		while(ft>lft){
			del x=sk[ft--];
			c[x.id].st[c[x.id].tp]=x.ps;
			c[x.id].tp=x.tp;
		}
	}void add(int x,int v){
		x=n-x+1;
		for(;x<=n;x+=x&-x)
			c[x].add(v,x);
	}int ans(int x,int v){
		int mn=9e18;x=n-x+1;
		for(;x;x-=x&-x)
			mn=min(mn,c[x].ans(v));
		return mn;
	}
}using namespace BIT;
void dfs(int x,int fa){
	idx[++m]=(d[x]+=d[fa]);
	int cc=lower_bound(idx+1,idx+m+1,d[x]-ln[x])-idx,lft=ft;
	if(x>1) f[x]=d[x]*p[x]+q[x]+ans(cc,p[x]);add(m,x);
	for(auto y:g[x]) dfs(y,x);clear(lft),m--;
}signed main(){
	ios::sync_with_stdio(0);
	cin.tie(0),cout.tie(0);
	cin>>n>>ft,ft=0;
	for(int i=2;i<=n;i++){
		int fa;cin>>fa>>d[i]>>p[i]>>q[i];
		cin>>ln[i],g[fa].push_back(i);
	}dfs(1,0);
	for(int i=2;i<=n;i++)
		cout<<f[i]<<"\n";
	return 0;
}