4515: [Sdoi2016]游戏

最新推荐文章于 2021-10-09 13:32:09 发布
最新推荐文章于 2021-10-09 13:32:09 发布
阅读量1.1k 收藏
点赞数 1
分类专栏： 李超树 线段树 树链剖分
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/CRZbulabula/article/details/52244666
版权

4515: [Sdoi2016]游戏

Time Limit: 40 Sec   Memory Limit: 256 MB
Submit: 396   Solved: 178
[ Submit][ Status][ Discuss]

Description

Alice 和 Bob 在玩一个游戏。
游戏在一棵有 n 个点的树上进行。最初，每个点上都只有一个数字，那个数字是 123456789123456789。
有时，Alice 会选择一条从 s 到 t 的路径，在这条路径上的每一个点上都添加一个数字。对于路径上的一个点 r，
若 r 与 s 的距离是 dis，那么 Alice 在点 r 上添加的数字是 a×dis+b。有时，Bob 会选择一条从 s 到 t 的路径。
他需要先从这条路径上选择一个点，再从那个点上选择一个数字。
Bob 选择的数字越小越好，但大量的数字让 Bob 眼花缭乱。Bob 需要你帮他找出他能够选择的最小的数字。

Input

第一行两个数字 n、m，表示树的点数和进行的操作数。
接下来 n−1 行，每行三个数字 u、v、w，表示树上有一条连接 u、v 的边，长度是 w。
接下来 m 行。每行第一个数字是 1 或 2。
若第一个数是 1，表示 Alice 进行操作，接下来四个数字 s、t、a、b。
若第一个数是 2，表示 Bob 进行操作，接下来四个数字 s、t。

Output

每当 Bob 进行操作，输出一行一个数，表示他能够选择的最小的数字

Sample Input

3 5
1 2 10
2 3 20
2 1 3
1 2 3 5 6
2 2 3
1 2 3 -5 -6
2 2 3

Sample Output

123456789123456789
6
-106

HINT

 n≤100000，m≤100000，∣a∣≤10000,0<=w，|b|<=10^9

Source

鸣谢Menci上传

[ Submit][ Status][ Discuss]

今天算是栽在这上面了，，日，，细节真多
总体思路显然：树链剖分+李超树
从s到Lca一段，从Lca到t一段
可分别用不同的斜率和截距表示，便于处理
即将原数据稍加修改得(详见代码)
然后是，，维护一个区间最低点
更新时刻有:
本区间被新线段覆盖
更新了子区间
记得更新区间最低点时不能舍去以前的点

查询答案的时候，从上往下一边找最低点一边用覆盖区间的线段尝试更新


日啊，，，调了半天，，真的是半天(from 12:00 to 20:30) 
#include<iostream>
#include<cstdio>
#include<queue>
#include<vector>
#include<bitset>
#include<algorithm>
#include<cstring>
#include<map>
#include<stack>
#include<set>
#include<cmath>
#include<ext/pb_ds/priority_queue.hpp>
using namespace std;

const int maxn = 1E5 + 10;
const int T = 20;
typedef long long LL;
const LL Inf = 123456789123456789;
typedef double DB;

struct E{
	int to,w;
	E (int _to = 0,int _w = 0) {to = _to; w = _w;}
};

int n,m,cnt,Root,Lca,I,fa[maxn][20],len[maxn],Nex[maxn],L[maxn]
	,First[maxn],pos[maxn],siz[maxn],root[maxn],lc[maxn*T],rc[maxn*T];
LL ans,Min[maxn*T],Light[maxn],k[maxn*T],b[maxn*T],dis2[maxn],dis[maxn];
bool Huge[maxn];

vector <E> v[maxn];
vector <int> Road[maxn];

int New() {++cnt; k[cnt] = 0; b[cnt] = Min[cnt] = Inf; return cnt;}
int Getlc(int o) {return lc[o]?lc[o]:lc[o] = New();}
int Getrc(int o) {return rc[o]?rc[o]:rc[o] = New();}

void dfs(int x,int from)
{
	Light[x] = Inf;
	siz[x] = 1; int Max = 0,po;
	for (int i = 1; i < 20; i++) fa[x][i] = fa[fa[x][i-1]][i-1];
	for (int i = 0; i < v[x].size(); i++) {
		int to = v[x][i].to;
		if (to == from) continue;
		fa[to][0] = x; dis[to] = dis[x] + 1LL*v[x][i].w;
		L[to] = L[x] + 1; dfs(to,x);
		siz[x] += siz[to];
		if (siz[to] > Max) Max = siz[to],po = to;
	}
	if (siz[x] > 1) Huge[po] = 1,Nex[x] = po;
}

void dfs2(int x,int from,int Len)
{
	if (siz[x] == 1) {
		if (Huge[x]) {
			len[x] = Len + 1; First[x] = x;
			root[x] = New(); pos[x] = 1;
			k[cnt] = 0; b[cnt] = Inf;
			dis2[x] = 0;
			Road[x].push_back(0);
			Road[x].push_back(x);
		}
		return;
	}
	Len = Huge[x]?++Len:0;
	for (int i = 0; i < v[x].size(); i++) {
		int to = v[x][i].to;
		if (to == from) continue;
		dfs2(to,x,Len);
	}
	if (Huge[x]) {
		root[x] = root[Nex[x]]; len[x] = len[Nex[x]]; First[x] = First[Nex[x]];
		pos[x] = pos[Nex[x]] + 1; dis2[x] = dis2[Nex[x]] + dis[Nex[x]] - dis[x];
		Road[First[x]].push_back(x);
	}
}

void pushdown(int o,int belong,int l,int r,LL K,LL B)
{
	int numl = Road[belong][l];
	int numr = Road[belong][r];
	LL A1 = dis2[numl]*K + B,A2 = dis2[numl]*k[o] + b[o];
	LL B1 = dis2[numr]*K + B,B2 = dis2[numr]*k[o] + b[o];
	Min[o] = min(Min[o],min(A1,B1));
	if (A1 < A2 && B1 < B2) {
		Min[o] = min(Min[o],min(A1,B1));
		k[o] = K; b[o] = B;
		return;
	}
	if (A1 >= A2 && B1 >= B2) return;
	DB x = (DB)(b[o] - B)/(DB)(K - k[o]);
	int Mid = Road[belong][(l+r)>>1]; LL mid = dis2[Mid]; Mid = pos[Mid];
	if (x <= mid) {
		if (A1 > A2) pushdown(Getlc(o),belong,l,Mid,k[o],b[o]),k[o] = K,b[o] = B;
		else pushdown(Getlc(o),belong,l,Mid,K,B);
	}
	else {
		if (A1 > A2) pushdown(Getrc(o),belong,Mid+1,r,K,B);
		else pushdown(Getrc(o),belong,Mid+1,r,k[o],b[o]),k[o] = K,b[o] = B;
	}
	Min[o] = min(Min[o],min(Min[Getlc(o)],Min[Getrc(o)]));
}

void modify(int o,int belong,int l,int r,int ml,int mr,LL K,LL B)
{
 	if (ml <= l && r <= mr) {pushdown(o,belong,l,r,K,B); return;}
	int mid = (l + r) >> 1;
	if (ml <= mid) modify(Getlc(o),belong,l,mid,ml,mr,K,B);
	if (mr > mid) modify(Getrc(o),belong,mid+1,r,ml,mr,K,B);
	Min[o] = min(Min[o],min(Min[Getlc(o)],Min[Getrc(o)]));
}

int Quickfa(int x,int y)
{
	for (int num = 0; y; y >>= 1,++num) 
		if (y & 1) x = fa[x][num]; 
	return x;
}

void Modify1(int x,LL K,LL B)
{
	LL Dis = 0;
	for (; x != Lca; Dis += (dis[x] - dis[fa[x][0]]),x = fa[x][0]) {
		if (!Huge[x]) Light[x] = min(Light[x],K*Dis + B);
		else {
			if (pos[x] == 1) {
				if (L[x] - (len[x] - pos[x]) > L[Lca]) {
					modify(root[x],First[x],1,len[x],pos[x],len[x],K,B+Dis*K);
					Dis += dis[x];
					x = Quickfa(x,len[x] - pos[x]);
					Dis -= dis[x];
				}
				else {
					modify(root[x],First[x],1,len[x],pos[x],pos[Lca]-1,K,B+Dis*K); 
					return;
				}
			}
			else {
				if (L[x] - (len[x] - pos[x]) > L[Lca]) {
					modify(root[x],First[x],1,len[x],pos[x],len[x],K,B+(Dis-dis2[x])*K);
					Dis += dis[x];
					x = Quickfa(x,len[x] - pos[x]);
					Dis -= dis[x];
				}
				else {
					modify(root[x],First[x],1,len[x],pos[x],pos[Lca]-1,K,B+(Dis-dis2[x])*K); 
					return;
				}
			}
		}
	}
}

void Modify2(int x,LL K,LL B,LL Dis)
{
	Dis += (dis[x] - dis[Lca]);
	for (; x != Lca; Dis -= (dis[x] - dis[fa[x][0]]),x = fa[x][0]) {
		if (!Huge[x]) Light[x] = min(Light[x],K*Dis + B);
		else {
			if (pos[x] == 1) {
				if (L[x] - (len[x] - pos[x]) > L[Lca]) {
					modify(root[x],First[x],1,len[x],pos[x],len[x],-K,B+Dis*K);
					Dis -= dis[x];
					x = Quickfa(x,len[x] - pos[x]);
					Dis += dis[x];
				}
				else {
					modify(root[x],First[x],1,len[x],pos[x],pos[Lca]-1,-K,B+Dis*K); 
					return;
				}
			}
			else {
				if (L[x] - (len[x] - pos[x]) > L[Lca]) {
					modify(root[x],First[x],1,len[x],pos[x],len[x],-K,B+(Dis+dis2[x])*K);
					Dis -= dis[x];
					x = Quickfa(x,len[x] - pos[x]);
					Dis += dis[x];
				}
				else {
					modify(root[x],First[x],1,len[x],pos[x],pos[Lca]-1,-K,B+(Dis+dis2[x])*K); 
					return;
				}
			}
		}	
	}
}

LL query(int o,int belong,int l,int r,int ql,int qr)
{
	LL QL = Road[belong][max(ql,l)]; QL = dis2[QL]; 
	LL QR = Road[belong][min(qr,r)]; QR = dis2[QR]; 
	LL A1 = QL*k[o] + b[o];
	LL B1 = QR*k[o] + b[o];
	if (ql <= l && r <= qr) return min(Min[o],min(A1,B1));
	int mid = (l + r) >> 1;
	LL ret = Inf;
	if (ql <= mid) ret = min(ret,query(Getlc(o),belong,l,mid,ql,qr));
	if (qr > mid) ret = min(ret,query(Getrc(o),belong,mid+1,r,ql,qr));
	return min(ret,min(A1,B1));
}

void Query(int x)
{
	for (; x != Lca; x = fa[x][0]) {
		if (!Huge[x]) ans = min(ans,Light[x]);
		else {
			if (L[x] - (len[x] - pos[x]) > L[Lca]) {
				ans = min(ans,query(root[x],First[x],1,len[x],pos[x],len[x]));
				x = Quickfa(x,len[x] - pos[x]);
			}
			else {
				ans = min(ans,query(root[x],First[x],1,len[x],pos[x],pos[Lca]-1)); 
				return;
			}
		}
	}
}

int getint()
{
	char ch = getchar();
	int ret = 0,multi = 1;
	while (ch < '0' || '9' < ch) {
		if (ch == '-') multi = -1;
		ch = getchar();
	}
	while ('0' <= ch && ch <= '9') ret = ret*10 + ch - '0',ch = getchar();
	return ret*multi;
}

int LCA(int p,int q)
{
	if (L[p] < L[q]) swap(p,q);
	int Log; for (Log = 0; L[p] - (1<<Log) >= 1; Log++); --Log;
	for (int j = Log; j >= 0; j--)
		if (L[p] - (1<<j) >= L[q])
			p = fa[p][j];
	if (p == q) return p;
	for (int j = Log; j >= 0; j--)
		if (fa[p][j] != fa[q][j])
			p = fa[p][j],q = fa[q][j];
	return fa[p][0];
}

int main()
{
	#ifdef DMC
		//freopen("DMC.txt","r",stdin);
		freopen("game20.in","r",stdin);
		freopen("test.txt","w",stdout);
	#endif
	
	n = getint(); m = getint();
	for (int i = 1; i < n; i++) {
		int x = getint(),y = getint(),z = getint();
		v[x].push_back(E(y,z));
		v[y].push_back(E(x,z));
	}
	Root = n/2; L[Root] = 1;
	dfs(Root,0); dfs2(Root,0,0);
	for (I = 1; I <= m; I++) {
		int typ = getint(),s = getint(),t = getint();
		Lca = LCA(s,t);
		if (typ == 1) {
			LL K = getint(),B = getint();
			if (Lca != s) Modify1(s,K,B); 
			if (!Huge[Lca]) Light[Lca] = min(Light[Lca],(dis[s] - dis[Lca])*K + B);
			else {
				DB NewB = (dis[s] - dis[Lca] - dis2[Lca])*K+B;
				modify(root[Lca],First[Lca],1,len[Lca],pos[Lca],pos[Lca],K,NewB);
			}
			if (Lca != t) Modify2(t,K,B,dis[s] - dis[Lca]);
		}
		else {
			ans = Inf;
			if (s != Lca) Query(s); 
			if (t != Lca) Query(t);
			if (!Huge[Lca]) ans = min(ans,Light[Lca]);
			else ans = min(ans,query(root[Lca],First[Lca],1,len[Lca],pos[Lca],pos[Lca]));
			printf("%lld\n",ans);
		}
	}
	return 0; 
}

CRZbulabula
关注
专栏目录
[SDOI2016] 游戏
C202044zxy的博客
11-29 218
一、题目 点此看题 二、解法
SDOI2016 总结
qq_42835841的博客
04-22 455
4644【SDOI2016】储能表 考虑二进制下数位dpdpdp，令>k>k>k表示满足条件，计算满足条件的数个数及和 记录f[N][0/1][0/1][0/1]f[N][0/1][0/1][0/1]f[N][0/1][0/1][0/1]表示f[step][顶着n][顶着m][>k]f[step][顶着n][顶着m][>k]f[step][顶着n][顶着m][>...
BZOJ 4515|SDOI 2016|游戏|树链剖分
とある黄鱼のOIブログ
04-28 577
未完待续#include <cstdio> #include <algorithm> #include <cstring> typedef long long ll; using namespace std; const int N = 100005, M = N * 2, S = M * 2; int h[N], p[M], v[M], sz[N], fa[N], cnt, ts; int dep
bzoj4515: [Sdoi2016]游戏
ACoder
02-19 592
树链剖分+标记永久化线段树
[SDOI2016]游戏
qq_52383696的博客
10-09 694
[SDOI2016]游戏 思路： 树剖+李超线段树 李超线段树模板题，把对一条链的操作变成李超线段树上一段区间的操作。维护李超就像普通的维护直线一样，但是每个区间加个标记，代表在自己下面的区间直线所出现的最小值，这样就保证了复杂度。 把一条链拆成从sss到lcalcalca，从lcalcalca到ttt的两个线段。 普普通通的AC吧； #include<bits/stdc++.h> using namespace std; #define ll long long int n,m; struc
bzoj 4515: [Sdoi2016]游戏 树链剖分
lych的博客
04-15 3196
       首先可以把修改转换一下，因为那个dis非常不爽。显然s~t的路径有s~lca和lca~t组成。令d[x]表示x的深度，对于s~lca上面的点，修改的值相当于a*(d[s]-d[x])+b=-a*d[x]+(b+a*d[s])，lca~t上面的点的值相当于a*(d[s]+d[x]-2*d[lca])+b=a*d[x]+(b+a*(d[s]-d[lca]*2))，这样就可以得到新的a'和...
SDOI 2016 排列计数
WedsonLin的博客
12-04 165
SDOI 2016 排列计数 题目链接 设满足题目条件且长度为n恰好有m个数稳定的序列数为ans(n,m)ans(n,m)ans(n,m) 那么易得题目答案ANS可表示为： ANS=ans(n,m)=C(n,m)∗ans(n−m,0)ANS=ans(n,m)=C(n,m)*ans(n-m,0)ANS=ans(n,m)=C(n,m)∗ans(n−m,0) 因为题目规定n,m≤1e6n,m\le1e6n,m≤1e6，所以C(n,m)可以通过预处理阶乘+乘法逆元快速求得。 至于ans(n-m,0)，又被称作为全错
[SDOI2016] 游戏（树链剖分+李超线段树）
hipamp的博客
10-20 245
题意 其中，n,m≤100000,∣a∣≤10000n,m\leq 100000,|a|\leq 10000n,m≤100000,∣a∣≤10000。 分析 设 iii 到根的距离为 disidis_idisi​。 首先化简一下式子，从 sss 到 lca(s,t)lca(s,t)lca(s,t) 路径： add=a×(diss−disi)+b=−a×disi+(b+a×diss)add=a\times (dis_s-dis_i)+b=-a\times dis_i+(b+a\times dis_s)add
[SDOI2016]储能表
C202044zxy的博客
03-08 263
一、题目 点此看题 二、解法 不要把这道题从图形的角度考虑，我们考虑它的本质就是求： ∑i=0n−1∑j=0m−1max⁡((i xor j)−k,0)\sum_{i=0}^{n-1}\sum_{j=0}^{m-1}\max((i\space xor\space j)-k,0)i=0∑n−1​j=0∑m−1​max((i xor j)−k,0)可以考虑数位d...
【BZOJ4515】【SDOI2016】游戏
cz_xuyixuan的博客
05-19 318
【题目链接】点击打开链接【思路要点】先树链剖分，然后用李超线段树维护序列，支持一段区间对某一次函数取Min，以及查询区间最小值。由于需要区间定位，单次1号操作的时间复杂度为\(O(Log^3N)\)，单次2号操作的时间复杂度为\(O(Log^2N)\)。总时间复杂度\(O(MLog^3N)\)，但数据没有将该做法的复杂度卡满，实际上做到这一点也很难，因此该做法可以通过本题，甚至跑得很快。【代码】#...
【LuoguP4069】 [SDOI2016]游戏
Yu Gi Oh!
02-10 247
题目链接 题意 树上路径插入以到起点距离为 xxx 的一次函数,询问路径 min 。边有边权 Sol 李超线段树+树剖 简单讲李超线段树: 支持区间加入一次函数,询问区间的一次函数最值。 修改复杂度为 O(log2n)O(log^2n)O(log2n) , 如果是全局加入就是O(log&nbsp;n)O(log\ n)O(log&nbsp;n) 查询复杂度为 O(log&nbsp;n)O(log\...
bzoj 4515: [Sdoi2016]游戏（树链剖分+线段树）
clover_hxy的博客
10-18 593
4515: [Sdoi2016]游戏 Time Limit: 40 Sec  Memory Limit: 256 MB Submit: 417  Solved: 185 [Submit][Status][Discuss] Description Alice 和 Bob 在玩一个游戏。 游戏在一棵有 n 个点的树上进行。最初，每个点上都只有一个数字，那个数字是 12345678912
Bzoj 4515 线段树 SDOI2016 游戏 Game
zzzzzfy的博客
04-20 1202
树链剖分 区间维护 直线  相交的话判下怎么传即可。。 注意相交上下取整问题 以及一系列long long 我求线段端点 写成了 b*x+b真是愚蠢至极 #include #include #include #include #define LL long long using namespace std; const int N=100010; const LL INF=123456789
BZOJ4515: [Sdoi2016]游戏
L_0_Forever_LF的专栏
12-01 455
传送门（这篇写得比我好qwq s到t的链可以拆成两条纵链 令dis[x]表示x到根的距离 对于一条纵链上的一个点x，添加的数就可以看成 (dis[x]-dis[t])* k+b =dis[x]*k-dis[t]*k+b =k*dis[x]+B 很像直线的表达式y=kx+b 树剖，问题变成线段树维护多条直线在区间内的最小值 对于一条链，因为这个区间内的dis是递增的即x坐标递增，这个
[SDOI 2016]游戏
HT008的Blog
03-09 226
题目描述： 雾. 题目分析： 首先可以把修改转换一下，因为那个dis非常不爽。显然s~t的路径有s~lca和lca~t组成。令d[x]表示x的深度，对于s~lca上面的点，修改的值相当于a*(d[s]-d[x])+b=-a*d[x]+(b-a*d[s]) lca~t上面的点的值相当于a*(d[s]+d[x]-2*d[lca])+b=a*d[x]+(b+a*(d[s]-d[lca]*2)...
BZOJ 4515: [Sdoi2016]游戏
月下沙茶树
05-18 594
树链剖分 将每次更新拆成若干条链上关于depth的一次函数 于是就变成了线段树维护半平面交中x轴上一段区间的y最小值 考虑每个点代表一条覆盖整个区间的直线，以此来进行标记永久化 若两个标记，即两条直线覆盖了同一区间 Case 1:其中一个标记处处比另一个优，那么显然取优的那个 Case 2:在该区间中有交点，于是取更优区间较大的那个为此区间的标记，另一个递归到相应的区间去继续更新 由
SDOI2016 游戏
g21wcr的博客
10-18 143
传送门 仔细一看：每个节点动态维护半平面交？？ 一搜，原来是李超线段树。 把原树剖分一下，就转化为二维平面内如下操作： ①加上一条线段。 ②询问一段区间内的最小值。 该线段树利用了标记永久化的思想。 每个区间[l,r][l,r][l,r]只保留一条线段：k∗mid+bk*mid+bk∗mid+b最小的那条。也就是与直线x=midx=midx=mid相交的最下面的一条。插入线段的时候根据交点维护这个...
洛谷·[SDOI2016]游戏
樱狸✨愿携清风归来处
08-22 175
洛谷P4069 [SDOI2016]游戏 李超线段树
[BZOJ4515][SDOI2016]游戏（线段树+链剖+差分）
薇小薇
04-23 270
题目： 我是超链接 题解： 我们先随意分析一下题目，这个处理路径上的最小值应该是树链剖分。 还有这个距离，一开始想的是lca，但我发现这个s作为一条路径的开始不能白费，而且对于每个节点和s求lca显然太慢了，我们考虑差分，即： 上升的路径：A(dis[s]−dis[r])+B=−Adis[r]+Adis[s]+BA(dis[s]−dis[r])+B=−Adis[r]+Adis[s]+...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值