纪念

数据结构一百题50题了呢，该过半周年啦~~~~

LYC和WGY半年的努力让这个几乎玩笑一般的系列到了现在。

今后也请多多关照啦。

祝愿dp100p早日过半

51.CF1000F One Occurrence

给定一个长度为\(n\)序列，\(m\)个询问，每次询问给定一个区间\([l,r]\)，如果这个区间里存在只出现一次的数，输出这个数（如果有多个就输出任意一个），没有就输出0。\(n,m\le 5e5\)

题意简述

给定一个序列，每次询问输出一个区间 \([l,r]\) 中只出现一次的数。

题解

这道题其实比较有意思。他让你输出只出现一次的数，讲道理你让我维护个数都还好办，这个就比较好玩了。

我这里给出一种不用吸氧的卡常莫队做法。

其实莫队需要维护的东西很简单，就是一个数的出现次数，我们用一个数组 cnt 来记录。

每次我们算加贡献的时候，cnt 的计算方法很显然。我们同时维护一个栈，算加贡献的时候如果这个数的出现次数为1，我们就把他放到栈顶上去。我们顺手维护一个数组 pos 表示这个数在栈里的位置。

算减贡献的时候 cnt 的计算方法依旧显然。我们把栈顶元素放到删除的数的位置上即可。

每次询问的答案就是栈顶元素，由于不合法的情况输出零，所以没必要特判。

光这样是过不了的，会T飞。

我们需要一个针对于莫队的优化，叫做奇偶性排序优化（瞎编的一个名字），具体来说就是在对询问排序的时候分块的奇偶性，具体实现看代码。这种排序的方式理论来说有 \(\Theta(\frac{1}{2})\) 的常数。

这样还是会T。容易发现输出答案的循环可以用unroll循环展开，我设的unroll的参是10，已经能过了。

  1. #include <cstdio>
  2. #include <algorithm>
  4. using namespace std;
  6. const int Maxn = 5e5 + 5, Each = 720;
  7. int n, m, top, isa[ Maxn ], cnt[ Maxn ], pos[ Maxn ], htl[ Maxn ], ans[ Maxn ];
  8. struct Query_Node
  9. {
  10. 	int l, r, id, pos;
  11. 	Query_Node( int L = 0, int R = 0, int ID = 0 ) { l = L, r = R, id = ID; }
  12. } Q[ Maxn ];
  14. inline int read( )
  15. {
  16. 	int a = 0, minus = 0;
  17. 	char ch = getchar( );
  18. 	while( !isdigit( ch ) )
  19. 	{
  20. 		if( ch == '-' )	 minus = 1;
  21. 		ch = getchar( );
  22. 	}
  23. 	while( isdigit( ch ) )
  24. 	{
  25. 		a = a * 10 + ch - '0';
  26. 		ch = getchar( );
  27. 	}
  28. 	if( minus )	 return -a;
  29. 	else	return a;
  30. }
  32. inline void write( int x )
  33. {
  34. 	if( x < 0 )
  35. 	{
  36. 		x = -x;
  37. 		putchar( '-' );
  38. 	}
  39. 	if( x > 9 ) 	write( x / 10 );
  40. 	putchar( x % 10 + '0' );
  41. }
  43. inline bool cmp( Query_Node rhs, Query_Node shr )
  44. {
  45. 	if( rhs.pos != shr.pos )	return rhs.l < shr.l;
  46. 	else if( rhs.pos & 1 )	return rhs.r < shr.r;
  47. 	else 	return rhs.r > shr.r;
  48. }
  50. inline void Make_Cont( int x, int t )
  51. {
  52. 	x = isa[ x ];
  53. 	if( t == 1 ) ++ cnt[ x ];
  54. 	else	--cnt[ x ];
  55. 	if( cnt[ x ] == 1 )
  56. 	{
  57. 		htl[ ++ top ] = x;
  58. 		pos[ x ] = top;
  59. 	}
  60. 	else if( ( t == 1 ) ? ( cnt[ x ] == 2 ) : ( cnt[ x ] == 0 ) )
  61. 	{
  62. 		htl[ pos[ x ] ] = htl[ top ];
  63. 		pos[ htl[ top ] ] = pos[ x ];
  64. 		pos[ x ] = htl[ top -- ] = 0;
  65. 	}
  66. }
  68. inline void Contribute( )
  69. {
  70. 	int l = 1, r = 0;
  71. 	for( int i = 0; i < m; ++ i )
  72. 	{
  73. 		while( l < Q[ i ].l )	Make_Cont( l ++, 0 );
  74. 		while( l > Q[ i ].l )	Make_Cont( -- l, 1 );
  75. 		while( r < Q[ i ].r )	Make_Cont( ++ r, 1 );
  76. 		while( r > Q[ i ].r )	Make_Cont( r --, 0 );
  77. 		ans[ Q[ i ].id ] = htl[ top ];
  78. 	}
  79. }
  81. signed main( )
  82. {
  83. 	n = read( );
  84. 	for( int i = 1; i <= n; ++ i ) 		isa[ i ] = read( );
  85. 	m = read( );
  86. 	for( int i = 0; i < m; ++ i )
  87. 	{
  88. 		Q[ i ].l = read( );
  89. 		Q[ i ].r = read( );
  90. 		Q[ i ].id = i;
  91. 		Q[ i ].pos = Q[ i ].l / Each;
  92. 	}
  93. 	sort( Q, Q + m, cmp );
  94. 	Contribute( );
  95. 	if( m <= 10 )
  96. 	{
  97. 		for( int i = 0; i < m; ++ i )		write( ans[ i ] ), putchar( '\n' );
  98. 		return 0;
  99. 	}
  100. 	#pragma unroll 10
  101. 	for( int i = 0; i < m; ++ i )		write( ans[ i ] ), putchar( '\n' );
  102. 	return 0;
  103. }

52.P5046 [Ynoi2019模拟赛]Yuno loves sqrt technology I

给你一个长为 \(n\) 的排列，\(m\) 次询问，每次查询一个区间的逆序对数，强制在线。

题意简述

无修改区间求逆序对。

题解

首先有一个显然的 \(\Theta(N\sqrt{N}\log_{2}N)\) 做法，由于过不了所以我就不废话。

其实有了 \(\Theta(N\sqrt{N}\log_{2}N)\) 的过不去做法，我们就可以根据这个思路然后预处理解决问题。

我们需要处理的信息有：

  1. 散块的逆序对数量

  2. 以块为单位的区间逆序对数量

那么我们需要处理的数组就有以下几个：

  1. previous[i] 表示 \(i\) 到 该块开头的逆序对数量。

  2. suffix[i] 同理。

  3. block[i][j] 表示前 \(i\) 个块中 \(\leq j\) 元素个数。

  4. intervals[i][j] 表示以块为单位的区间 \([i,j]\) 中的逆序对数量。

讲讲预处理方法。

  1. previous[i] 和 suffix[i] 的处理方法都很显然，可以一直扫着然后FWT扫就行。

  2. block[i][j] 可以递推，递推式为 block[i][j]=block[i+1][j]+block[i][j-1]-block[i+1][j-1]+cont(i,j)。其中 cont(i,j) 表示计算对应块的逆序对数。

  3. intervals[i][j] 每次循环到块的开头继承上一个块的贡献即可。

计算贡献的方法很简单，归并即可。mrsrz讲得也挺清楚的，我这里就不再赘述，主要讲讲怎么卡常。

首先我们可以把主函数里的所有循环全部展开，经过实践参数传8的时候跑得比较快。

然后八聚氧先加上，luogu O2也开着。

再其次快读fread快输fwrite，这些都是卡常的标配。

然后就把能拿出来的结构体拿出来，实在不能就不管了。

然后去STL，pair vector能去就去。

然后long long开在正确的地方，不要无脑replace。

函数inline，循环register。虽然可能作用不大但是可以先加上。

然后调块长，经过无数次实践发现取150~170较为优秀。

然后加了过后发现就算rp再好也只有60pts。

然后谷歌搜索硫酸的化学式H₂SO₄，给评测机喂硫酸（idea来自SyadouHayami）。

然后本来交了5页都过不了，这下再交两次就过了。

  1. // 省略八聚氧
  2. #include <cstdio>
  3. #include <iostream>
  4. #include <algorithm>
  5. #include <cstring>
  7. using namespace std;
  9. const int Maxn = 1e5 + 5;
  10. const int Maxm = 650;
  11. const int each = 160;
  12. int n, m, blocks, Lp[Maxn], Rp[Maxn], isa[Maxn], head[Maxn], tail[Maxn], sorted[Maxn], belong[Maxn], previous[Maxn], suffix[Maxn], block[Maxm][Maxn];
  13. long long intervals[Maxm][Maxn];
  14. struct Holy_Pair
  15. {
  16. 	int first, second;
  18. 	bool operator < (const Holy_Pair& rhs) const
  19. 	{
  20. 		return first < rhs.first;
  21. 	}
  22. } current[Maxn];
  23. struct Fenwick_Tree
  24. {
  25. 	int fwt[Maxn];
  27. 	inline void Modify(int x, int v)
  28. 	{
  29. 		for (; x + 5 <= Maxn; x += x & -x)	fwt[x] += v;
  30. 	}
  32. 	inline int Query(int x)
  33. 	{
  34. 		int ret = 0;
  35. 		for (; x; x ^= x & -x)	ret += fwt[x];
  36. 		return ret;
  37. 	}
  38. } FWT;
  40. #define io_e '\0'
  41. #define io_s ' '
  42. #define io_l '\n'
  43. namespace Fast_IO
  44. {
  45. ... // 省略快读
  46. }  // namespace Fast_IO
  48. using Fast_IO::read;
  49. using Fast_IO::write;
  51. inline Holy_Pair make_pair(int first, int second)
  52. {
  53. 	Holy_Pair ret;
  54. 	ret.first = first;
  55. 	ret.second = second;
  56. 	return ret;
  57. }
  59. inline int Merge_Vct(int rhs[], int shr[], int szl, int szr)
  60. {
  61. 	int itl = 1, itr = 1;
  62. 	int ret = 0, ctr = 1;
  63. 	while (itl <= szl && itr <= szr)
  64. 	{
  65. 		if (rhs[itl] < shr[itr]) 	++itl, ++ctr;
  66. 		else
  67. 		{
  68. 			ret += szl - ctr + 1;
  69. 			++itr;
  70. 		}
  71. 	}
  72. 	return ret + szr - szr;
  73. }
  75. inline int Merge_Idx(int st1, int st2, int sz1, int sz2)
  76. {
  77. 	int ret = 0, id1 = st1 + 1, id2 = st2 + 1;
  78. 	sz1 += st1, sz2 += st2;
  79. 	while (id1 <= sz1 && id2 <= sz2)
  80. 	{
  81. 		if (sorted[id1] < sorted[id2])		++id1;
  82. 		else
  83. 		{
  84. 			ret += sz1 - id1 + 1;
  85. 			++id2;
  86. 		}
  87. 	}
  88. 	return ret;
  89. }
  91. inline void Behavior(int l, int r, long long &ans)
  92. {
  93. 	int itl = 0, itr = 0;
  94. 	if (belong[l] == belong[r])
  95. 	{
  96. 		for (int i = head[belong[l]]; i <= tail[belong[r]]; ++i)
  97. 		{
  98. 			if (current[i].second >= l && current[i].second <= r)	Rp[++itr] = sorted[i];
  99. 			else if (current[i].second < l)		Lp[++itl] = sorted[i];
  100. 		}
  101. 		if (l == head[belong[l]])	ans = previous[r] - Merge_Vct(Lp, Rp, itl, itr);
  102. 		else 	ans = previous[r] - previous[l - 1] - Merge_Vct(Lp, Rp, itl, itr);
  103. 	}
  104. 	else
  105. 	{
  106. 		ans = intervals[belong[l] + 1][belong[r] - 1] + previous[r] + suffix[l];
  107. 		for (int i = head[belong[l]]; i <= tail[belong[l]]; ++i)
  108. 		{
  109. 			if (current[i].second >= l)
  110. 			{
  111. 				Lp[++itl] = sorted[i];
  112. 				ans += block[belong[r] - 1][1] - block[belong[r] - 1][sorted[i]] - block[belong[l]][1] + block[belong[l]][sorted[i]];
  113. 			}
  114. 		}
  115. 		for (int i = head[belong[r]]; i <= tail[belong[r]]; ++i)
  116. 		{
  117. 			if (current[i].second <= r)
  118. 			{
  119. 				Rp[++itr] = sorted[i];
  120. 				ans += block[belong[r] - 1][sorted[i] + 1] - block[belong[l]][sorted[i] + 1];
  121. 			}
  122. 		}
  123. 		ans += Merge_Vct(Lp, Rp, itl, itr);
  124. 	}
  125. 	write(io_l, ans);
  126. }
  128. signed main()
  129. {
  130. 	read(n, m), blocks = (n - 1) / each + 1;
  131. 	if (n <= 8)
  132. 	{
  133. 		for (int i = 1; i <= n; ++i)
  134. 		{
  135. 			read(isa[i]);
  136. 			current[i] = make_pair(isa[i], i);
  137. 		}
  138. 	}
  139. 	else
  140. 	{
  141. 		#pragma unroll 8
  142. 		for (int i = 1; i <= n; ++i)
  143. 		{
  144. 			read(isa[i]);
  145. 			current[i] = make_pair(isa[i], i);
  146. 		}
  147. 	}
  148. 	if (blocks <= 8)
  149. 	{
  150. 		for (int i = 1; i <= blocks; ++i)
  151. 		{
  152. 			head[i] = tail[i - 1] + 1;
  153. 			tail[i] = tail[i - 1] + each;
  154. 			if (i == blocks) 	tail[i] = n;
  155. 		}
  156. 	}
  157. 	else
  158. 	{
  159. 		#pragma unroll 8
  160. 		for (int i = 1; i <= blocks; ++i)
  161. 		{
  162. 			head[i] = tail[i - 1] + 1;
  163. 			tail[i] = tail[i - 1] + each;
  164. 			if (i == blocks) 	tail[i] = n;
  165. 		}
  166. 	}
  167. 	if (blocks <= 8)
  168. 	{
  169. 		for (int i = 1; i <= blocks; ++i)
  170. 		{
  171. 			memcpy(block[i], block[i - 1], sizeof(block[0]));
  172. 			sort(current + head[i], current + 1 + tail[i]);
  173. 			for (int j = head[i]; j <= tail[i]; ++j)
  174. 			{
  175. 				++block[i][isa[j]];
  176. 				belong[j] = i;
  177. 				sorted[j] = current[j].first;
  178. 			}
  179. 			int satisfy = 0;
  180. 			for (int j = head[i]; j <= tail[i]; ++j)
  181. 			{
  182. 				FWT.Modify(isa[j], 1);
  183. 				satisfy += FWT.Query(n) - FWT.Query(isa[j]);
  184. 				previous[j] = satisfy;
  185. 			}
  186. 			intervals[i][i] = satisfy;
  187. 			for (int j = head[i]; j <= tail[i]; ++j)
  188. 			{
  189. 				suffix[j] = satisfy;
  190. 				FWT.Modify(isa[j], -1);
  191. 				satisfy -= FWT.Query(isa[j] - 1);
  192. 			}
  193. 		}
  194. 	}
  195. 	else
  196. 	{
  197. 		#pragma unroll 8
  198. 		for (int i = 1; i <= blocks; ++i)
  199. 		{
  200. 			memcpy(block[i], block[i - 1], sizeof(block[0]));
  201. 			sort(current + head[i], current + 1 + tail[i]);
  202. 			for (int j = head[i]; j <= tail[i]; ++j)
  203. 			{
  204. 				++block[i][isa[j]];
  205. 				belong[j] = i;
  206. 				sorted[j] = current[j].first;
  207. 			}
  208. 			int satisfy = 0;
  209. 			for (int j = head[i]; j <= tail[i]; ++j)
  210. 			{
  211. 				FWT.Modify(isa[j], 1);
  212. 				satisfy += FWT.Query(n) - FWT.Query(isa[j]);
  213. 				previous[j] = satisfy;
  214. 			}
  215. 			intervals[i][i] = satisfy;
  216. 			for (int j = head[i]; j <= tail[i]; ++j)
  217. 			{
  218. 				suffix[j] = satisfy;
  219. 				FWT.Modify(isa[j], -1);
  220. 				satisfy -= FWT.Query(isa[j] - 1);
  221. 			}
  222. 		}
  223. 	}
  224. 	if (blocks <= 8)
  225. 	{
  226. 		for (int dis = 1; dis <= blocks; ++dis)
  227. 		{
  228. 			for (int i = n - 1; i; --i)		block[dis][i] += block[dis][i + 1];
  229. 			for (int l = 1, r = dis + 1; r <= blocks + 1; ++l, ++r)
  230. 				intervals[l][r] = intervals[l + 1][r] + intervals[l][r - 1] - intervals[l + 1][r - 1] +
  231. 									Merge_Idx(head[l] - 1, head[r] - 1, tail[l] - head[l] + 1, tail[r] - head[r] + 1);
  232. 		}
  233. 	}
  234. 	else
  235. 	{
  236. 		#pragma unroll 8
  237. 		for (int dis = 1; dis <= blocks; ++dis)
  238. 		{
  239. 			for (int i = n - 1; i; --i)		block[dis][i] += block[dis][i + 1];
  240. 			for (int l = 1, r = dis + 1; r <= blocks + 1; ++l, ++r)
  241. 				intervals[l][r] = intervals[l + 1][r] + intervals[l][r - 1] - intervals[l + 1][r - 1] +
  242. 									Merge_Idx(head[l] - 1, head[r] - 1, tail[l] - head[l] + 1, tail[r] - head[r] + 1);
  243. 		}
  244. 	}
  246. 	if (m <= 8)
  247. 	{
  248. 		long long lastans = 0;
  249. 		for (int i = 0; i < m; ++i)
  250. 		{
  251. 			long long l, r;
  252. 			read(l, r);
  253. 			l ^= lastans;
  254. 			r ^= lastans;
  255. 			Behavior(l, r, lastans);
  256. 		}
  257. 	}
  258. 	else
  259. 	{
  260. 		long long lastans = 0;
  261. 		#pragma unroll 8
  262. 		for (int i = 0; i < m; ++i)
  263. 		{
  264. 			long long l, r;
  265. 			read(l, r);
  266. 			l ^= lastans;
  267. 			r ^= lastans;
  268. 			Behavior(l, r, lastans);
  269. 		}
  270. 	}
  271. 	return 0;
  272. }

53.SP20644 ZQUERY - Zero Query

长度为\(n\)的序列，序列中的值为\(1\)或\(-1\)

有\(m\)个询问，询问在\([L,R]\)中区间和为\(0\)的区间的最大长度

题意简述

求区间中区间和为0的最长区间。

题解

首先前缀和基本操作，然后问题转化为：

给定询问区间 \([l,r]\)，求区间中相同的数的最远间隔距离，即P5906。

考虑莫队来整。

可以维护两个东西，一个是 pre[x] 表示当前的询问区间 \([l,r]\) 中 \(x\) 最早出现的位置。

还有一个是 suf[x] 表示当前的询问区间 \([l,r]\) 中 \(x\) 最后出现的位置。

在莫队加贡献的时候这些很好整，分别像这样维护：

如果加贡献时出现了数 \(x\) 那么 pre[x]=p，\(p\) 表示加贡献时计算的位置。

suf[x] 同理。

那么询问的答案就是 \(\max_{i\in [l,r]}\{suf_{a_{x}}-pre_{a_{x}}\}\)。

中间有个特殊情况就是需要判断还没被更新过的时候，自己打代码的时候注意点。

然后发现一个神奇的事情：减贡献我们做不来了！

既然我们不会做减法，那就不做吧。

正常的莫队的询问区间左右指针 \(l,r\) 一般初值为 \(l=1,r=0\)。

我们这里换一种方式。

首先我们把询问分块，然后排序，这是莫队的常规操作。

然后再枚举每个块，我们首先处理询问的左端点在当前块中的询问，思考一下询问分块的排序方式就可以知道这种算法的正确性。

然后我们把指针的初值设为 \(l=R_{i}+1,r=R_{i}-1\)，其中数组 \(R\) 表示块的左端点。

然后如果询问区间在一个块中就暴力。

由于询问区间在一个块中的询问我们已经暴力处理过了，所以剩下的询问右端点一定在其他块。

然后 \(r\) 指针正常右移，这一步只需要加贡献。

所以 \(l\) 指针只需要左移，我们就可以合并两次的贡献，就是这次询问的答案。

由于排序的规则，剩下的询问右端点一定在当前询问右端点的右边，所以我们的 \(r\) 指针可以不用管，就放到当前询问的位置。

但是 \(l\) 是无序的，所以我们需要撤销计算 \(l\) 产生的贡献，这很简单，不赘述。

这种算法叫回滚莫队。

但是这种trick我在之前一道莫队题中想到过，不过觉得太麻烦就没打。

后来发现有这个东西的板子，又重学了一遍。

还是挺好玩的吧，数据结构。

  1. #include <cstdio>
  2. #include <iostream>
  3. #include <algorithm>
  4. #include <cstring>
  5. #include <queue>
  6. #include <cmath>
  8. using namespace std;
  10. const int Maxn = 2e5 + 5;
  11. int n, m, each, cube, isa[Maxn], pre[Maxn], nxt[Maxn], vio[Maxn], ans[Maxn], bel[Maxn], head[Maxn], tail[Maxn];
  12. struct Query_Node
  13. {
  14. 	int l, r, id;
  15. } Q[Maxn];
  17. bool cmp(Query_Node rhs, Query_Node shr)
  18. {
  19. 	if (bel[rhs.l] == bel[shr.l])	return rhs.r < shr.r;
  20. 	else 	return bel[rhs.l] < bel[shr.l];
  21. }
  23. void Contribute(int las = 0)
  24. {
  25. 	for (int k = 1; k <= cube; ++k)
  26. 	{
  27. 		int l = tail[k] + 1, r = tail[k], ret = 0;
  28. 		vector < int > vctr;
  29. 		for (int c = las + 1; bel[Q[c].l] == k; ++c)
  30. 		{
  31. 			int ql = Q[c].l, qr = Q[c].r, now = 0;
  32. 			if (bel[ql] == bel[qr])
  33. 			{
  34. 				for (int i = ql; i <= qr; ++i)	vio[isa[i]] = 0;
  35. 				for (int i = ql; i <= qr; ++i)
  36. 				{
  37. 					if (!vio[isa[i]])	vio[isa[i]] = i;
  38. 					else	now = max(now, i - vio[isa[i]]);
  39. 				}
  40. 				ans[Q[c].id] = now;
  41. 				las = c;
  42. 			}
  43. 			else
  44. 			{
  45. 				while (r < qr)
  46. 				{
  47. 					++r;
  48. 					nxt[isa[r]] = r;
  49. 					if (pre[isa[r]] == 0)
  50. 					{
  51. 						pre[isa[r]] = r;
  52. 						vctr.push_back(isa[r]);
  53. 					}
  54. 					ret = max(ret, nxt[isa[r]] - pre[isa[r]]);
  55. 				}
  56. 				now = ret;
  57. 				while (l > ql)
  58. 				{
  59. 					--l;
  60. 					if (nxt[isa[l]] == 0)	nxt[isa[l]] = l;
  61. 					else	ret = max(ret, nxt[isa[l]] - l);
  62. 				}
  63. 				while (l < tail[k] + 1)
  64. 				{
  65. 					if (nxt[isa[l]] == l)	nxt[isa[l]] = 0;
  66. 					++l;
  67. 				}
  68. 				ans[Q[c].id] = ret;
  69. 				ret = now;
  70. 				las = c;
  71. 			}
  72. 		}
  73. 		for (unsigned i = 0; i < vctr.size(); ++i)	pre[vctr[i]] = nxt[vctr[i]] = 0;
  74. 	}
  75. }
  77. signed main()
  78. {
  79. 	scanf("%d %d", &n, &m);
  80.     ++n;
  81. 	each = pow(n, 0.5);
  82. 	cube = (n - 1) / each + 1;
  83. 	for (int i = 2; i <= n; ++i)
  84. 	{
  85. 		scanf("%d", &isa[i]);
  86.         isa[i] += isa[i - 1];
  87. 	}
  88.     for (int i = 1; i <= n; ++i)    isa[i] += 50005;
  89. 	for (int i = 1; i <= cube; ++i)
  90. 	{
  91. 		head[i] = tail[i - 1] + 1;
  92. 		tail[i] = tail[i - 1] + each;
  93. 		for (int j = head[i]; j <= tail[i]; ++j)	bel[j] = i;
  94. 		if (i == cube)	tail[i] = n;
  95. 	}
  96. 	for (int i = 1; i <= m; ++i)
  97. 	{
  98. 		scanf("%d %d", &Q[i].l, &Q[i].r);
  99.         Q[i].r += 1;
  100. 		Q[i].id = i;
  101. 	}
  102. 	sort(Q + 1, Q + 1 + m, cmp);
  103. 	Contribute();
  104. 	for (int i = 1; i <= m; ++i)	printf("%d\n", ans[i]);
  105. 	return 0;
  106. }

54.P5268 [SNOI2017]一个简单的询问

给你一个长度为 \(N\) 的序列 \(a_i\)，\(1\leq i\leq N\)，和 \(q\) 组询问，每组询问读入 \(l_1,r_1,l_2,r_2\)，需输出

\[\sum\limits_{x=0}^\infty \text{get}(l_1,r_1,x)\times \text{get}(l_2,r_2,x)
\]

$ \text{get}(l,r,x)$ 表示计算区间 \([l,r]\) 中，数字 \(x\) 出现了多少次。

题意简述

求式子

题解

纯算式子nt行为，由于是出现次数，满足区间加减性，所以我们可以这样表达 \(\mathrm{get}(l,r,x)\)（省略 \(x\)）:

\[\mathrm{get}(l,r)=\mathrm{get}(1,r)-\mathrm{get}(1,l-1)
\]

那么我们代进原式，化一波式子：

\[\sum_{i=1}^{\infty}\mathrm{get}(l_{1},r_{1},x)\times\mathrm{get}(l_{2},r_{2},x)
\]
\[\sum_{i=1}^{\infty}(\mathrm{get}(1,r_{1})-\mathrm{get}(1,l_{1}-1))\times(\mathrm{get}(1,r_{2})-\mathrm{get}(1,l_{2}-1))
\]
\[\sum_{i=1}^{\infty}\mathrm{get}(r_{1})\times\mathrm{get}(r_{2})-\mathrm{get}(r_{1})\times\mathrm{get}(l_{2}-1)-\mathrm{get}(l_{1}-1)\times\mathrm{get}(r_{2})+\mathrm{get}(l_{1}-1))\times\mathrm{get}(l_{2}-1)
\]
\[\mathrm{let}\ F(x,y)=\sum_{d}\mathrm{get}(1,l,d)\times\mathrm{get}(1,r,d)
\]

则答案为:

\[F(r_{1},r_{2})-F(r_{1},l_{2}-1)-F(l_{1}-1,r_{2})+F(l_{1}-1,l_{2}-1)
\]

考虑怎么更新，比如从 \(l\) 更新到 \(l+1\)，则：

\[\mathrm{get(1,l)}\times\mathrm{get}(1,r)
\]
\[\mathrm{get(1,l+1)}\times\mathrm{get}(1,r)
\]
\[\mathrm{get(1,l)}\times\mathrm{get}(1,r)+\mathrm{cont}(a_{l})
\]

其中 \(\mathrm{cont}(a_{l})\) 表示 \(a_{l}\) 的出现次数。

则我们就知道怎么更新了，由于我们维护和的是前缀信息，所以姿势和普通莫队有点不一样。

维护两个数组 cntl[x] 和 cntr[y] 表示答案式子

\[F(r_{1},r_{2})-F(r_{1},l_{2}-1)-F(l_{1}-1,r_{2})+F(l_{1}-1,l_{2}-1)
\]

中 \(F(x,y)\) 的值。

式子和表达可能有些玄学，自己意会一下吧。

  1. #include <cstdio>
  2. #include <algorithm>
  4. using namespace std;
  6. const int Maxn = 50000 + 5;
  7. int n, m, q, each, cube, isa[Maxn], bel[Maxn], cntl[Maxn], cntr[Maxn];
  8. long long now, ans[Maxn];
  9. struct Query_Node
  10. {
  11. 	int l, r, id, typ;
  12. 	Query_Node(int tl = 0, int tr = 0, int tid = 0, int tpyt = 0)
  13. 	{
  14. 		l = tl;
  15. 		r = tr;
  16. 		id = tid;
  17. 		typ = tpyt;
  18. 	}
  19. } Q[Maxn << 2];
  21. bool cmp(Query_Node rhs, Query_Node shr)
  22. {
  23. 	if (bel[rhs.l] == bel[shr.l])	return rhs.r < shr.r;
  24. 	else 	return rhs.l < shr.l;
  25. }
  27. void Contribute()
  28. {
  29. 	int l = 0, r = 0;
  30. 	for (int i = 1; i <= m; ++i)
  31. 	{
  32. 		while (l < Q[i].l)
  33. 		{
  34. 			++cntl[isa[++l]];
  35. 			now += cntr[isa[l]];
  36. 		}
  37. 		while (l > Q[i].l)
  38. 		{
  39. 			--cntl[isa[l]];
  40. 			now -= cntr[isa[l--]];
  41. 		}
  42. 		while (r < Q[i].r)
  43. 		{
  44. 			++cntr[isa[++r]];
  45. 			now += cntl[isa[r]];
  46. 		}
  47. 		while (r > Q[i].r)
  48. 		{
  49. 			--cntr[isa[r]];
  50. 			now -= cntl[isa[r--]];
  51. 		}
  52. 		ans[Q[i].id] += Q[i].typ * now;
  53. 	}
  54. }
  56. signed main()
  57. {
  58. 	scanf("%d", &n);
  59. 	each = 230;
  60. 	for (int i = 1; i <= n; ++i)
  61. 	{
  62. 		scanf("%d", &isa[i]);
  63. 		bel[i] = (i - 1) / each + 1;
  64. 	}
  65. 	scanf("%d", &q);
  66. 	for (int i = 1; i <= q; ++i)
  67. 	{
  68. 		int l1, r1, l2, r2;
  69. 		scanf("%d %d %d %d", &l1, &r1, &l2, &r2);
  70. 		Q[++m] = Query_Node(r1, r2, i, 1);
  71. 		Q[++m] = Query_Node(r1, l2 - 1, i, -1);
  72. 		Q[++m] = Query_Node(l1 - 1, r2, i, -1);
  73. 		Q[++m] = Query_Node(l1 - 1, l2 - 1, i, 1);
  74. 	}
  75. 	for (int i = 1; i <= m; ++i)
  76. 	{
  77. 		if (Q[i].l < Q[i].r)	swap(Q[i].l, Q[i].r);
  78. 	}
  79. 	sort(Q + 1, Q + 1 + m, cmp);
  80. 	Contribute();
  81. 	for (int i = 1; i <= q; ++i)	printf("%lld\n", ans[i]);
  82. 	return 0;
  83. }

55.P3466 [POI2008]KLO-Building blocks

有 \(n\) 柱砖,每柱砖有一个高度，我们现在希望有连续 \(k\) 柱的高度是一样的。

你可以选择以下两个动作：

1:从某柱砖的顶端拿一块砖出来,丢掉不要了。

2:从仓库中拿出一块砖,放到另一柱，仓库是无限大的。

现在希望用最小次数的动作完成任务，除此之外你还要求输出结束状态时,每柱砖的高度。

题意简述

给定一个长度为 \(n\) 的序列，支持单点减一，单点加法的操作，问把一个长度为 \(k\) 的区间里的值变成同样的最小操作数。

题解

挺好的一道平衡树练手+猜结论题。（雾

化简一下题意，我们需要找到一个值 \(m\) 且 \(m\) 使得 \(\sum_{i=1}^{r}\mid a_{i}-x\mid\) 最小。

其中 \(r-l+1=k\)。

大约可以猜到这里需要取中位数。

感性理解一下，把 \(a_{l},a_{l+1},\cdots,a_{r}\) 排序后如果每个数要取最小只能取中位数。

可以当成一个结论式的东西积累起来。

那么我们就需要一个数据结构，支持插入，删除，查 \(k_{th}\)。

这里我选择的是Zip-Tree，也就是国内常说的FHQ-Treap。

具体的做法是，枚举每个长度为 \(k\) 的区间 \([l,r]\)，首先我们把 \(a_{1},a_{2},\cdots,a_{k}\) 插入平衡树，然后根据 \(k\) 的奇偶性用 \(k_{th}\) 操作找到中位数（序列的中位数根据序列的长度奇偶性不同有区别）。

然后计算贡献即可，顺手记录一下区间的位置，方便输出答案。

贡献的计算可以分成大于中位数，小于中位数，等于中位数的情况。

  1. #include <cstdio>
  2. #include <cstdlib>
  3. #include <ctime>
  4. #include <climits>
  5. #define int long long
  7. const int Maxn = 1e5 + 5;
  8. int n, k, rt, tot, isa[Maxn];
  9. struct Treap
  10. {
  11. 	int l, r, siz, sum, key, val;
  12. } t[Maxn];
  14. int newnode(int val)
  15. {
  16. 	t[++tot].val = val;
  17. 	t[tot].sum = val;
  18. 	t[tot].key = rand();
  19. 	t[tot].siz = 1;
  20. 	return tot;
  21. }
  23. void maintain(int x)
  24. {
  25. 	t[x].siz = t[t[x].l].siz + t[t[x].r].siz + 1;
  26. 	t[x].sum = t[t[x].l].sum + t[t[x].r].sum + t[x].val;
  27. }
  29. void Split(int now, int val, int &x, int &y)
  30. {
  31. 	if (now == 0)	x = y = 0;
  32. 	else
  33. 	{
  34. 		if (val >= t[now].val)
  35. 		{
  36. 			x = now;
  37. 			Split(t[now].r, val, t[now].r, y);
  38. 		}
  39. 		else
  40. 		{
  41. 			 y = now;
  42. 			 Split(t[now].l, val, x, t[now].l);
  43. 		}
  44. 		maintain(now);
  45. 	}
  46. }
  48. int Merge(int x, int y)
  49. {
  50. 	if (x == 0 || y == 0)	return x + y;
  51. 	else
  52. 	{
  53. 		if (t[x].key > t[y].key)
  54. 		{
  55. 			t[x].r = Merge(t[x].r, y);
  56. 			maintain(x);
  57. 			return x;
  58. 		}
  59. 		else
  60. 		{
  61. 			t[y].l = Merge(x, t[y].l);
  62. 			maintain(y);
  63. 			return y;
  64. 		}
  65. 	}
  66. }
  68. int rt1, rt2, rt3;
  69. void Insert(int val)
  70. {
  71. 	Split(rt, val, rt1, rt2);
  72. 	rt = Merge(rt1, Merge(newnode(val), rt2));
  73. }
  75. void Remove(int val)
  76. {
  77. 	Split(rt, val, rt1, rt3);
  78. 	Split(rt1, val - 1, rt1, rt2);
  79. 	rt2 = Merge(t[rt2].l, t[rt2].r);
  80. 	rt = Merge(rt1, Merge(rt2, rt3));
  81. }
  83. int Find_Kth(int rnk)
  84. {
  85. 	int now = rt;
  86. 	while (now)
  87. 	{
  88. 		if (t[t[now].l].siz + 1 == rnk) 	break;
  89. 		else if (t[t[now].l].siz >= rnk) 	now = t[now].l;
  90. 		else
  91. 		{
  92. 			rnk -= t[t[now].l].siz + 1;
  93. 			now = t[now].r;
  94. 		}
  95. 	}
  96. 	return t[now].val;
  97. }
  99. int Query(int mid)
  100. {
  101. 	Split(rt, mid, rt1, rt3);
  102. 	Split(rt1, mid - 1, rt1, rt2);
  103. 	int lth = rt1, mth = rt3;
  104. 	int cont_l = t[lth].siz * mid - t[lth].sum;
  105. 	int cont_r = t[mth].sum - t[mth].siz * mid;
  106. 	rt = Merge(rt1, Merge(rt2, rt3));
  107. 	return cont_l + cont_r;
  108. }
  110. signed main()
  111. {
  112. 	srand((unsigned)(time(NULL)));
  113. 	int MID, LEFT, RIGHT, ANS = LONG_LONG_MAX;
  114. 	scanf("%lld %lld", &n, &k);
  115. 	for (int i = 1; i <= n; ++i)	scanf("%lld", &isa[i]);
  116. 	for (int i = 1; i <= k; ++i) 	Insert(isa[i]);
  117. 	int l = 1, r = k;
  118. 	for (int i = k; i <= n; ++i)
  119. 	{
  120. 		int mid;
  121. 		if (k & 1)	mid = Find_Kth((k >> 1) + 1);
  122. 		else 	mid = (Find_Kth(k >> 1) + Find_Kth((k >> 1) + 1)) >> 1;
  123. 		int xxx = Query(mid);
  124. 		if (ANS > xxx)
  125. 		{
  126. 			ANS = xxx;
  127. 			MID = mid;
  128. 			LEFT = l - 1;
  129. 			RIGHT = r + 1;
  130. 		}
  131. 		l++, r++;
  132. 		Remove(isa[l - 1]);
  133. 		Insert(isa[r]);
  134. 	}
  135. 	printf("%lld\n", ANS);
  136. 	for (int i = 1; i <= LEFT; ++i) 	printf("%lld\n", isa[i]);
  137. 	for (int i = LEFT + 1; i <= RIGHT - 1; ++i) 	printf("%lld\n", MID);
  138. 	for (int i = RIGHT; i <= n; ++i) 	printf("%lld\n", isa[i]);
  139. 	return 0;
  140. }

56.P5610 [Ynoi2013]大学

一个长为 \(n\) 的非负整数序列 \(a\)，支持以下两个操作：

  • 1 l r x：把区间 \([l,r]\) 中所有 \(x\) 的倍数除以 \(x\)。
  • 2 l r：查询区间 \([l,r]\) 的和。

本题强制在线，每次的 \(l,r,x\) 需要 xor 上上次答案，如果之前没有询问，则上次答案为 \(0\)。

题意简述

区间查 \(x\) 的倍数并除掉，区间查和。

题解

平衡树。

首先有个基本的想法就是按 \(a_{i}\) 开平衡树，即对于每个 \(a_{i}\) 都开一棵平衡树，共5e5棵，第 \(i\) 棵平衡树维护的值是所有 \(a_{i}\) 的倍数在原数组中的下标，用处后面讲。

注意到一个小性质，一个正整数 \(A\) 最多被整除 \(\log_{2}A\) 次，这个很好想，每次都至少减少一半。可以当成一个 trick 记下来。

整个区间的数最多被除 \(\sum_{i=1}^{n}\log_{2}a_{i}\) 次，区间和的操作可以用树状数组操作实现，则整体的操作复杂度为 \(\Theta(\sum_{i=1}^{n}\log_{2}a_{i}+\log_{2}a_{i})\)。

开头提到了对于每个 \(a_{i}\) 我们都开一棵平衡树，作用就体现在这里。因为如果要保证正确的时间复杂度，我们需要快速的找到需要执行操作的数。

这里我采用的是 FHQ-Treap。

我们可以用两次 split 操作在 \(x\) 的平衡树中提取出当前的询问区间，由于我们是以下标为平衡树维护的权值，所以我们用按值分裂即可提取出区间。

然后我们就在提取出的子树中 DFS 遍历，然后暴力操作，把操作后依然是 \(x\) 的倍数的数记录下来，操作完后用这个数组再建一棵临时树然后和之前 split 出来的子树一起合并回去。

【记得配图，解释一下】

操作之前记得预处理每个数的所有约数，这个简单，直接用 vector 即可。

  1. #include <cstdio>
  2. #include <cstdlib>
  3. #include <algorithm>
  4. #include <ctime>
  5. #include <queue>
  7. using namespace std;
  8. typedef long long t_t;
  10. const int Maxn = 1e5 + 5;
  11. const int Maxa = 5e5 + 5;
  12. int n, m, xxx, zip, tot, isa[Maxn], bin[Maxa], root[Maxa];
  13. struct Treap
  14. {
  15. 	int l, r, key, val;
  16. } t[Maxa * 230];
  17. struct Fenwick_Tree
  18. {
  19. 	t_t bit[Maxn];
  21. 	void ins(int x, t_t v)
  22. 	{
  23. 		for (; x <= n; x += x & -x) 	bit[x] += v;
  24. 	}
  26. 	t_t sum(int x)
  27. 	{
  28. 		t_t ret = 0;
  29. 		for (; x; x -= x & -x)	ret += bit[x];
  30. 		return ret;
  31. 	}
  32. } fwt;
  33. vector < int > vec[Maxa];
  35. int newnode(int val)
  36. {
  37. 	t[++tot].val = val;
  38. 	t[tot].key = rand();
  39. 	return tot;
  40. }
  42. void split(int now, int val, int &x, int &y)
  43. {
  44. 	if (now == 0)	x = y = 0;
  45. 	else
  46. 	{
  47. 		if (t[now].val <= val)
  48. 		{
  49. 			x = now;
  50. 			split(t[now].r, val, t[now].r, y);
  51. 		}
  52. 		else
  53. 		{
  54. 			y = now;
  55. 			split(t[now].l, val, x, t[now].l);
  56. 		}
  57. 	}
  58. }
  60. int merge(int x, int y)
  61. {
  62. 	if (x == 0 || y == 0)	return x | y;
  63. 	if (t[x].key > t[y].key)
  64. 	{
  65. 		t[x].r = merge(t[x].r, y);
  66. 		return x;
  67. 	}
  68. 	else
  69. 	{
  70. 		t[y].l = merge(x, t[y].l);
  71. 		return y;
  72. 	}
  73. }
  75. int build(int l, int r)
  76. {
  77. 	if (l > r)	return 0;
  78. 	int mid = (l + r) >> 1;
  79. 	int now = newnode(bin[mid]);
  80. 	t[now].l = build(l, mid - 1);
  81. 	t[now].r = build(mid + 1, r);
  82. 	return now;
  83. }
  85. void dfs(int now, int val)
  86. {
  87. 	if (now == 0)	return ;
  88. 	dfs(t[now].l, val);
  89. 	if (isa[t[now].val] % val == 0)
  90. 	{
  91. 		fwt.ins(t[now].val, -isa[t[now].val]);
  92. 		isa[t[now].val] /= val;
  93. 		fwt.ins(t[now].val, isa[t[now].val]);
  94. 	}
  95. 	if (isa[t[now].val] % val == 0)	bin[++zip] = t[now].val;
  96. 	dfs(t[now].r, val);
  97. }
  99. int tx, ty, tp;
  100. void change(int l, int r, int x)
  101. {
  102. 	if (x == 1) 	return ;
  103. 	split(root[x], r, tx, tp);
  104. 	split(tx, l - 1, tx, ty);
  105. 	zip = 0;
  106. 	dfs(ty, x);
  107. 	root[x] = merge(tx, merge(build(1, zip), tp));
  108. }
  110. signed main()
  111. {
  112. 	srand((114514 - 1) / 3 - 4);
  113. 	scanf("%d %d", &n, &m);
  114. 	for (int i = 1; i <= n; ++i)
  115. 	{
  116. 		scanf("%d", &isa[i]);
  117. 		fwt.ins(i, isa[i]);
  118. 		xxx = max(xxx, isa[i]);
  119. 	}
  120. 	for (int i = 1; i <= n; ++i)
  121. 	{
  122. 		for (int j = 1; j * j <= isa[i]; ++j)
  123. 		{
  124. 			if (isa[i] % j == 0)
  125. 			{
  126. 				vec[j].push_back(i);
  127. 				if (j * j != isa[i])	vec[isa[i] / j].push_back(i);
  128. 			}
  129. 		}
  130. 	}
  131. 	for (int i = 1; i <= xxx; ++i)
  132. 	{
  133. 		zip = 0;
  134. 		for (unsigned j = 0; j < vec[i].size(); ++j)	bin[++zip] = vec[i][j];
  135. 		root[i] = build(1, zip);
  136. 	}
  137. 	for (int i = 0; i < m; ++i)
  138. 	{
  139. 		int t, l, r, x;
  140. 		scanf("%d %d %d", &t, &l, &r);
  141. 		if (t == 1)
  142. 		{
  143. 			scanf("%d", &x);
  144. 			change(l, r, x);
  145. 		}
  146. 		else	printf("%lld\n", fwt.sum(r) - fwt.sum(l - 1));
  147. 	}
  148. 	return 0;
  149. }

57.CF803G Periodic RMQ Problem

给你一个序列\(a\) 让你支持

\(1\) \(l\) \(r\) \(x\) 区间赋值

\(2\) \(l\) \(r\) 询问区间最小值

我们觉得这个问题太水了,所以我们不会给你序列\(a\)

而是给你序列一个长度为\(n\) 的序列\(b\) ,把\(b\) 复制粘贴\(k\) 次就可以得到\(a\)

\(n\le10^5,k\le10^4,q\le10^5,b_i\le10^9\)

\(1\le l\le r\le n\times k\)

题目解释：

有一个由给定序列不断首尾拼接而成的序列（很长，存不下）。

要求实现两个操作：区间赋值和查询区间最小值。

这个序列很长不好操作，所以我们要进行离散化，缩小一下操作的范围。

把每个操作的左端点减一和右端点加入离散化数组\(pri\)（这里下标从1开始）。

离散化后第i个块的左端点是\(pri[i-1]+1\),右端点是\(pri[i]\)。（离散化还要插入\(0\)和最后一个位置的下标，不然有些位置没有被覆盖到）

但是离散化之后我们不知道每个离散化出来的块内的最小值，这样就无法初始化。

我们可以利用之前给出用来首尾拼接的序列求出来。

对每个块分类讨论：（设给定用来拼接的序列长为\(n\)）

1).块长大于等于\(n\)：

此时这个块已经覆盖了整个给定序列，那它的最小值直接就是给定序列的最小值。

2).块长小于\(n\)：

又要分两种情况:

$\ \ \ \ $A).左端点投射到序列中的位置在右端点的后面（由于是由同一序列首尾拼接而成，我们可以把它们投射到一个序列上）。

$\ \ \ \ $就是说：这个块覆盖的给定序列的前面的一部分和后面的一部分。

$\ \ \ \ $就像这样：

$\ \ \ \ $投射到给定序列上就像这样：

$\ \ \ \ \(那么我们直接查询（涂色的）这两段的最小值再取个\)min$就好了。

$\ \ \ \ \(其实就是从他们中间的序列尾（就是更上面的那张图L,R中间那条黑线）分开来求两边的最小值再取\)min$，就得到了他们之间的最小值。

$\ \ \ \ $B).左端点投射到序列中的位置在右端点的前面。

$\ \ \ \ $仿照上个情况，投射到序列上查询他们之间的最小值即可。

由于是以操作的左右端点做的离散化，所以一个块内要么一起被询问，要么一起被修改。

所以我们用上面的方法处理出这个离散化后的序列每个块的初始值之后，就可以用线段树把它当板题做了！

代码：

  1. #include<cstdio>
  2. #include<vector>
  3. #include<algorithm>
  4. using namespace std;
  5. vector<int> pri;
  6. const int INF=1e9+1;
  7. int n,q,k,a[100010],ori[400010],op[100010],opl[100010],opr[100010],opx[100010],nodes[1000010],tag[2000010],s;
  8. void buildori(int l,int r,int x)//对于给定序列的线段树建树。
  9. {
  10. 	if(l^r)
  11. 	{
  12. 		int mid=(l+r)>>1;
  13. 		buildori(l,mid,x<<1);
  14. 		buildori(mid+1,r,x<<1|1);
  15. 		ori[x]=min(ori[x<<1],ori[x<<1|1]);
  16. 	}
  17. 	else	ori[x]=a[l];
  18. }
  19. int findori(int l,int r,int x,int fr,int ba)//查询给定序列里的区间最小值。
  20. {
  21. 	if(l>ba||r<fr)	return INF;
  22. 	if(l>=fr&&r<=ba)	return ori[x];
  23. 	else
  24. 	{
  25. 		int mid=(l+r)>>1;
  26. 		return min(findori(l,mid,x<<1,fr,ba),findori(mid+1,r,x<<1|1,fr,ba));
  27. 	}
  28. }
  29. void build(int l,int r,int x)//对离散化后的序列建树。
  30. {
  31. 	if(l^r)
  32. 	{
  33. 		int mid=(l+r)>>1;
  34. 		build(l,mid,x<<1);
  35. 		build(mid+1,r,x<<1|1);
  36. 		nodes[x]=min(nodes[x<<1],nodes[x<<1|1]);
  37. 	}
  38. 	else
  39. 	{
  40. 		if(pri[l]-pri[l-1]>=n)	nodes[x]=ori[1];//情况1).
  41. 		else
  42. 		{
  43. 			int one=(pri[l-1]+1)%n;
  44. 			if(one==0)	one=n;
  45. 			int two=pri[l]%n;
  46. 			if(two==0)	two=n;
  47. 			if(one>two)	nodes[x]=min(findori(1,n,1,one,n),findori(1,n,1,1,two));//情况2).A).
  48. 			else	nodes[x]=findori(1,n,1,one,two);//情况2).B).
  49. 		}
  50. 	}
  51. }
  52. //后面的就是模板操作。
  53. void pushdown(int x)
  54. {
  55. 	if(tag[x])
  56. 	{
  57. 		nodes[x]=tag[x];
  58. 		tag[x<<1]=tag[x<<1|1]=tag[x];
  59. 		tag[x]=0;
  60. 	}
  61. }
  62. void update(int l,int r,int x,int fr,int ba,int val)
  63. {
  64. 	if(l>ba||r<fr)	return;
  65. 	if(l>=fr&&r<=ba)	tag[x]=val;
  66. 	else
  67. 	{
  68. 		int mid=(l+r)>>1;
  69. 		pushdown(x);
  70. 		update(l,mid,x<<1,fr,ba,val);
  71. 		update(mid+1,r,x<<1|1,fr,ba,val);
  72. 		pushdown(x<<1);
  73. 		pushdown(x<<1|1);
  74. 		nodes[x]=min(nodes[x<<1],nodes[x<<1|1]);
  75. 	}
  76. }
  77. int find(int l,int r,int x,int fr,int ba)
  78. {
  79. 	if(l>ba||r<fr)	return INF;
  80. 	pushdown(x);
  81. 	if(l>=fr&&r<=ba)	return nodes[x];
  82. 	else
  83. 	{
  84. 		int mid=(l+r)>>1;
  85. 		return min(find(l,mid,x<<1,fr,ba),find(mid+1,r,x<<1|1,fr,ba));
  86. 	}
  87. }
  88. int main()
  89. {
  90. 	scanf("%d %d",&n,&k);
  91. 	for(int i=1;i<=n;++i)	scanf("%d",&a[i]);
  92. 	buildori(1,n,1);
  93. 	scanf("%d",&q);
  94. 	pri.push_back(0);
  95. 	pri.push_back(n*k);//插入0和最后一个位置的下标。
  96. 	for(int i=1;i<=q;++i)
  97. 	{
  98. 		scanf("%d %d %d",&op[i],&opl[i],&opr[i]);
  99. 		if(op[i]==1)	scanf("%d",&opx[i]);
  100. 		pri.push_back(opl[i]-1);
  101. 		pri.push_back(opr[i]);//把每个操作的左端点减一和右端点加入离散化数组pri。
  102. 	}
  103. 	sort(pri.begin(),pri.end());
  104. 	pri.erase(unique(pri.begin(),pri.end()),pri.end());
  105. 	s=pri.size();
  106. 	build(1,s,1);//0下标的位置没有前一个，所以从1开始建树。
  107. 	for(int i=1;i<=q;++i)
  108. 	{
  109. 		//lower_bound返回大于等于给定值的下标，插入操作时是插入的opl[i]-1和opr[i]，所以查询大于等于opl[i]的下标可以查询到opl[i]所在的块。
  110. 		if(op[i]==1)	update(1,s,1,lower_bound(pri.begin(),pri.end(),opl[i])-pri.begin(),lower_bound(pri.begin(),pri.end(),opr[i])-pri.begin(),opx[i]);
  111. 		else	printf("%d\n",find(1,s,1,lower_bound(pri.begin(),pri.end(),opl[i])-pri.begin(),lower_bound(pri.begin(),pri.end(),opr[i])-pri.begin()));
  112. 	}
  113. 	return 0;
  114. }

58.P5305 [GXOI/GZOI2019]旧词

> 浮生有梦三千场

> 穷尽千里诗酒荒

> 徒把理想倾倒

> 不如早还乡

>

> 温一壶风尘的酒

> 独饮往事迢迢

> 举杯轻思量

> 泪如潮青丝留他方

>

> ——乌糟兽/愚青《旧词》

你已经解决了五个问题，不妨在这大树之下，吟唱旧词一首抒怀。最后的问题就是关于这棵树的，它的描述很简单。

给定一棵 \(n\) 个点的有根树，节点标号 \(1 \sim n\)，\(1\) 号节点为根。

给定常数 \(k\)。

给定 \(Q\) 个询问，每次询问给定 \(x,y\)。

求：

\[\sum\limits_{i \le x} \text{depth}(\text{lca}(i,y))^k
\]

\(\text{lca}(x,y)\) 表示节点 \(x\) 与节点 \(y\) 在有根树上的最近公共祖先。

\(\text{depth}(x)\) 表示节点 \(x\) 的深度，根节点的深度为 \(1\)。

由于答案可能很大，你只需要输出答案模 \(998244353\) 的结果。

这道题目非常简洁：要求 \(\sum_{i\le x}depth(lca(i,y))^k\) （都写在题面里了）

要直接解决这个问题是有点困难的，那么——

我们先看它的弱化版：[LNOI2014]LCA

要求 \(\sum_{i=l}^rdepth(lca(i,z))\)

少了个k次方呢！

首先转化一下 \(\sum_{i=l}^rdepth(lca(i,z))=\sum_{i=1}^rdepth(lca(i,z))-\sum_{i=1}^{l-1}depth(lca(i,z))\)

可以用前缀和来解决这个问题。

把每个询问拆成 \((1,l-1)\) 和 \((1,r)\) 分别解决。

那么我们可以按顺序把 \(1\) 到 \(n\) 的点加进来，同时计算加到每个点时的询问答案。

那么现在考虑加到第 \(i\) 个点时如何计算贡献。

先上张图举个例子~

其中涂成红色的点就是已经加入的点。

我们先计算一下每个节点的子树中有多少已经加入的点，用 \(siz[x]\) 表示。没有加入的节点对当前询问肯定没有贡献。所以我们直接统计加入的节点就好了

那么 \(siz[1]=4\) , \(siz[2]=2\), \(siz[3]=1\), \(siz[4]=1\) , \(siz[5]=1\) , \(siz[6]=0\) .

现在考虑计算当前加入的节点到4的LCA深度和。

\(ans=siz[4]\times dep[4]+(siz[2]-siz[4])\times dep[2]+(siz[1]-siz[2])\times dep[1]\)

即加入的节点中在 \(4\) 的子树中的到 \(4\) 的 \(LCA\) 肯定为 \(4\)。在 \(2\) 的子树但不在 \(4\) 的子树中的到 \(4\) 的 \(LCA\) 肯定为 \(2\)，在 \(1\) 的子树但不在 \(2\) 的子树中的到 \(4\) 的 \(LCA\) 肯定为 \(1\)。

把式子拆开化简得到:

\(ans=siz[4]\times (dep[4]-dep[2])+siz[2]\times (dep[2]-dep[1])+siz[1]\)

\(ans=siz[4]+siz[2]+siz[1]\)

于是我们发现其实答案就是该节点到根路径的节点子树大小之和。

普遍形式：设查询节点到根的路径为 \(\{v_1,v_2,v_3\dots v_k\}\)。

\(ans=siz[v_1]\times dep[v_1]+(siz[v_2]-siz[v_1])\times dep[v_2]+(siz[v_3]-siz[v_2])\times dep[v_3]+\dots +(siz[v_k]-siz[v_{k-1}])\times dep[v_k]\)

由于是到根节点的路径，所以 \(dep[v_i]=dep[v_{i+1}]+1\)

于是化简得

\(ans=siz[v_1]\times (dep[v_1]-dep[v_2])+siz[v_2]\times (dep[v_2]-dep[v_3])+siz[v_3]\times (dep[v_3]-dep[v_4])+\dots +siz[v_k]\times dep[v_k]\)

\(ans=siz[v_1]+siz[v_2]+siz[v_3]+\dots +siz[v_k]\)

就此得出结论：答案就是该节点到根路径的节点子树大小之和。

查询直接查询要查询的节点到根路径上的节点的子树大小之和。

加入节点时只有它到根的路径上的节点的子树大小加了一。

使用树链剖分加线段树维护即可。

回到本题

本题没有 \(l\) 的限制，所以我们不用前缀和拆询问处理了。

那么我们抬出之前的式子改成这道题的样子。

\(ans=siz[v1]\times dep[v1]^k+(siz[v2]-siz[v1])\times dep[v2]^k+(siz[v3]-siz[v2])\times dep[v3]^k+\dots+(siz[vk]-siz[vk-1])\times dep[vk]^k\)

化简一下，得到：

\(ans=siz[v_1]\times (dep[v_1]^k-dep[v_2]^k)+siz[v_2]\times (dep[v_2]^k-dep[v_3]^k)+siz[v_3]\times (dep[v_3]^k-dep[v4]^k)+\dots siz[v_k]\times dep[v_k]^k\)

可以发现，\(dep[v_2]\) 总等于 \(dep[v_1]-1\)（因为是从查询节点到根节点的路径嘛）。

所以对于\((dep[v_i]^k-dep[v_{i+1}]^k)\)这种东西，我们可以对它进行预处理。

设对于点i的预处理答案为 \(val[i]\)。

则答案为:

\(ans=siz[v_1]\times val[v_1]+siz[v_2]\times val[v_2]+siz[v_3]\times val[v_3]+...siz[v_k]\times val[v_k]\)

就是在线段树维护时多加一个权值的问题了。

在线段树的每个节点上附加一个权值，表示这个区间里所有点的 \(val\) 之和。

这样就可以区间修改和下传标记了。

跟上面那道弱化版差得不多吧~

代码：

  1. //直接用LCA代码改的，可能有点迷惑（？
  2. #include<cstdio>
  3. #include<vector>
  4. #include<algorithm>
  5. using namespace std;
  6. vector<long long> e[50010];
  7. struct node
  8. {
  9. 	long long bas,sum,tag;//bas为附加的权值
  10. }nodes[400010];
  11. struct ask
  12. {
  13. 	long long pos,val,ID,nag;
  14. }q[100010];
  15. bool cmp(ask one,ask ano)
  16. {
  17. 	return one.pos<ano.pos;
  18. }
  19. long long n,m,s,opl,opr,opz,f,cnt,tot,waste,ans[50010],kkk,dep[50010],fa[50010],son[50010],siz[50010],hb[50010],dfn[50010],ton[50010],power[50010],bruh;
  20. const long long mod=998244353;
  21. long long qpow(long long bas,long long tim)//快速幂用来处理val
  22. {
  23. 	long long res=1,fur=bas;
  24. 	while(tim)
  25. 	{
  26. 		if(tim&1)	res=(res*fur)%mod;
  27. 		fur=(fur*fur)%mod;
  28. 		tim>>=1;
  29. 	}
  30. 	return res;
  31. }
  32. void dfs(long long x,long long las)
  33. {
  34. 	fa[x]=las;
  35. 	dep[x]=dep[las]+1;
  36. 	siz[x]=1;
  37. 	long long b=0,s=0;
  38. 	for(long long i=0;i<e[x].size();++i)
  39. 	{
  40. 		long long y=e[x][i];
  41. 		dfs(y,x);
  42. 		siz[x]+=siz[y];
  43. 		if(siz[y]>b)
  44. 		{
  45. 			b=siz[y];
  46. 			s=y;
  47. 		}
  48. 	}
  49. 	son[x]=s;
  50. }
  51. void dfs2(long long x,long long las,long long heavy)//树剖dfs
  52. {
  53. 	if(heavy)	hb[x]=hb[las];
  54. 	else	hb[x]=x;
  55. 	dfn[x]=++cnt;
  56. 	ton[cnt]=x;
  57. 	if(son[x])	dfs2(son[x],x,1);
  58. 	for(long long i=0;i<e[x].size();++i)
  59. 	{
  60. 		long long y=e[x][i];
  61. 		if(y^son[x])	dfs2(y,x,0);
  62. 	}
  63. }
  64. void pushdown(long long x)
  65. {
  66. 	if(nodes[x].tag)
  67. 	{
  68. 		//(siz[x]+a)*val[x]+(siz[y]+a)*val[y]+...+(siz[z]+a)*val[z]
  69. 		//=siz[x]*val[x]+siz[y]*val[y]+...+siz[z]*val[z]+a*(val[x]+val[y]+...+val[z])
  70. 		nodes[x].sum+=(nodes[x].bas*nodes[x].tag);
  71. 		nodes[x].sum%=mod;
  72. 		nodes[x<<1].tag+=nodes[x].tag;
  73. 		nodes[x<<1].tag%=mod;
  74. 		nodes[x<<1|1].tag+=nodes[x].tag;
  75. 		nodes[x<<1|1].tag%=mod;
  76. 		nodes[x].tag=0;
  77. 	}
  78. }
  79. void build(long long l,long long r,long long x)//预处理val区间和
  80. {
  81. 	if(l^r)
  82. 	{
  83. 		long long mid=(l+r)>>1;
  84. 		build(l,mid,x<<1);
  85. 		build(mid+1,r,x<<1|1);
  86. 		nodes[x].bas=(nodes[x<<1].bas+nodes[x<<1|1].bas)%mod;
  87. 	}
  88. 	else	nodes[x].bas=power[dep[ton[l]]]-power[dep[ton[l]]-1];
  89. //	printf(" %lld %lld %lld\n",l,r,nodes[x].bas);
  90. }
  91. void update(long long l,long long r,long long x,long long fr,long long ba)
  92. {
  93. 	if(l>ba||r<fr)	return;
  94. 	if(l>=fr&&r<=ba)	nodes[x].tag=(nodes[x].tag+1)%mod;
  95. 	else
  96. 	{
  97. 		pushdown(x);
  98. 		long long mid=(l+r)>>1;
  99. 		update(l,mid,x<<1,fr,ba);
  100. 		update(mid+1,r,x<<1|1,fr,ba);
  101. 		pushdown(x<<1);
  102. 		pushdown(x<<1|1);
  103. 		nodes[x].sum=(nodes[x<<1].sum+nodes[x<<1|1].sum)%mod;
  104. 	}
  105. }
  106. long long find(long long l,long long r,long long x,long long fr,long long ba)
  107. {
  108. 	if(l>ba||r<fr)	return 0;
  109. 	pushdown(x);
  110. 	if(l>=fr&&r<=ba)	return nodes[x].sum;
  111. 	else
  112. 	{
  113. 		long long mid=(l+r)>>1;
  114. 		return (find(l,mid,x<<1,fr,ba)+find(mid+1,r,x<<1|1,fr,ba))%mod;
  115. 	}
  116. }
  117. void output(long long l,long long r,long long x)
  118. {
  119. 	pushdown(x);
  120. 	printf(" %lld %lld %lld\n",l,r,nodes[x].sum);
  121. 	if(l^r)
  122. 	{
  123. 		long long mid=(l+r)>>1;
  124. 		output(l,mid,x<<1);
  125. 		output(mid+1,r,x<<1|1);
  126. 	}
  127. }
  128. void update_LCA(long long x,long long y)
  129. {
  130. 	long long fx=hb[x],fy=hb[y];
  131. 	while(fx^fy)
  132. 	{
  133. 		if(dep[fx]<dep[fy])
  134. 		{
  135. 			swap(fx,fy);
  136. 			swap(x,y);
  137. 		}
  138. 		update(1,s,1,dfn[fx],dfn[x]);
  139. 		x=fa[fx];
  140. 		fx=hb[x];
  141. 	}
  142. 	update(1,s,1,min(dfn[x],dfn[y]),max(dfn[x],dfn[y]));
  143. }
  144. long long find_LCA(long long x,long long y)
  145. {
  146. 	long long res=0;
  147. 	long long fx=hb[x],fy=hb[y];
  148. 	while(fx^fy)
  149. 	{
  150. 		if(dep[fx]<dep[fy])
  151. 		{
  152. 			swap(fx,fy);
  153. 			swap(x,y);
  154. 		}
  155. 		res+=find(1,s,1,dfn[fx],dfn[x]);
  156. 		res%=mod;
  157. 		x=fa[fx];
  158. 		fx=hb[x];
  159. 	}
  160. 	res+=find(1,s,1,min(dfn[x],dfn[y]),max(dfn[x],dfn[y]));
  161. 	res%=mod;
  162. 	return res;
  163. }
  164. int main()
  165. {
  166. 	scanf("%lld %lld %lld",&n,&m,&bruh);
  167. 	s=n;
  168. 	for(long long i=1;i<=n;++i)	power[i]=qpow(i,bruh);
  169. 	for(long long i=2;i<=n;++i)
  170. 	{
  171. 		scanf("%lld",&f);
  172. 		e[f].push_back(i);
  173. 	}
  174. 	dfs(1,1);
  175. 	dfs2(1,1,0);
  176. 	build(1,s,1);
  177. 	for(long long i=1;i<=m;++i)
  178. 	{
  179. 		scanf("%lld %lld",&opr,&opz);
  180. 		q[++tot].ID=i;
  181. 		q[tot].nag=1;
  182. 		q[tot].pos=opr;
  183. 		q[tot].val=opz;
  184. 	}
  185. 	sort(q+1,q+1+tot,cmp);
  186. 	while(q[kkk].pos<=0&&kkk<=tot)
  187. 	{
  188. 		ans[q[kkk].ID]+=(q[kkk].nag*0);
  189. 		++kkk;
  190. 	}
  191. 	for(long long i=1;i<=n;++i)
  192. 	{
  193. 		update_LCA(1,i);
  194. 		while(q[kkk].pos<=i&&kkk<=tot)
  195. 		{
  196. 			ans[q[kkk].ID]+=(q[kkk].nag*find_LCA(1,q[kkk].val));
  197. 			ans[q[kkk].ID]%=mod;
  198. 			++kkk;
  199. 		}
  200. //		puts("");
  201. //		output(1,s,1);
  202. //		puts("\n");
  203. 	}
  204. 	for(long long i=1;i<=m;++i)	printf("%lld\n",((ans[i]%mod)+mod)%mod);
  205. 	return 0;
  206. }

59.CF643G Choosing Ads

  • 给定一个长度为 \(n\) 的序列和一个整数 \(p\)。
  • 有 \(m\) 个操作，操作要么是区间赋值，要么是询问区间内出现次数至少占 \(p\%\) 的数。
  • 输出询问的答案时，可以包含错的数，也可以重复输出，但对的数一定要在答案中，且输出的数的个数不超过 \(\lfloor \frac{100}{p} \rfloor\)。

题意简述

够简单了。

题解

我觉得很有意思的一道题目。

第一个操作，嗯，不讲。

我们着重来看看第二个操作。

这个查询操作要我们输出所有出现次数的频率 \(\ge p\%\) 的数字。

先考虑一种特殊情况，即当 \(p>50\) 的情况，这时候这个问题的本质就成了求区间的绝对众数。

区间绝对众数有一个经典的 \(\Theta(1)\) 空间的做法，摩尔投票。

这里简单介绍一下摩尔投票法。

举个例子，序列 \(\texttt{7 5 7 6 7 7 4 7}\)。

我们定义两个变量 major,cnt。

我们遍历一遍序列：

起初 major=7,cnt=1。

然后 \(5\neq 7\)，cnt--。

此时 cnt==0，我们令 major=5,cnt=1。

后面以此类推。

那么拓展到这道题目，我们的 p 最小是 20。

和普通的摩尔投票法类似，我们把 major,cnt 替换为 major[5],cnt[5]，用数组来求出分别的众数。

不过这样做有可能会出现出现频率不符合答案规范的数存进来，不过反正良心出题人允许我们输出无意义内容，所以就不用管了。

然后考虑用一个数据结构来加速这个过程。

区间推平、查询首选线段树。

那么要维护的信息就很明显了。

每个结点维护一个摩尔投票的结构体，然后再维护一个区间推平的 tag。

最后 pushup 合并左右儿子贡献时，直接把右区间合并到左边来即可。

下传标记很简单，不过要注意下传时因为时区间推平，所以摩尔投票的结构体里的 cnt 部分要设为 r-l+1。

  1. #include <cstdio>
  2. #include <algorithm>
  3. #include <cstring>
  4. #include <assert.h>
  6. using namespace std;
  8. const int Maxn = 15e4 + 5;
  9. int n, m, p, isa[Maxn];
  10. struct Major_Voting // 摩尔投票的结构体
  11. {
  12. 	int major[5], cnt[5];
  13. 	Major_Voting()
  14. 	{
  15. 		memset(major, 0, sizeof(major));
  16. 		memset(cnt, 0, sizeof(cnt));
  17. 	}
  18. };
  19. struct Segment_Tree
  20. {
  21. 	Major_Voting val;
  22. 	int tag;
  23. } t[Maxn << 2];
  25. Major_Voting make_vote(int major) // 以 major 为关键字创建一个新的摩尔投票
  26. {
  27. 	Major_Voting ret;
  28. 	memset(ret.major, 0, sizeof(ret.major));
  29. 	memset(ret.cnt, 0, sizeof(ret.cnt));
  30. 	ret.major[0] = major;
  31. 	ret.cnt[0] = 1;
  32. 	return ret;
  33. }
  35. Major_Voting merge(Major_Voting one, Major_Voting ano) // 合并区间，这一步的细节较多
  36. {
  37. 	Major_Voting ret = one;
  38. 	for (int i = 0; i < p; ++i)
  39. 	{
  40. 		if (ano.cnt[i])
  41. 		{
  42. 			bool flag = 0;
  43. 			for (int j = 0; j < p; ++j)
  44. 			{
  45. 				if (ret.major[j] == ano.major[i])
  46. 				{
  47. 					ret.cnt[j] += ano.cnt[i];
  48. 					flag = 1;
  49. 					break;
  50. 				}
  51. 			}
  52. 			if (flag == 0)
  53. 			{
  54. 				for (int j = 0; j < p; ++j)
  55. 				{
  56. 					if (ret.cnt[j] == 0)
  57. 					{
  58. 						ret.major[j] = ano.major[i];
  59. 						ret.cnt[j] = ano.cnt[i];
  60. 						flag = 1;
  61. 						break;
  62. 					}
  63. 				}
  64. 			} // 把右区间的贡献算到左边来
  65. 			if (flag)	continue;
  66. 			int line = 2e9;
  67. 			for (int j = 0; j < p; ++j)	line = min(line, ret.cnt[j]);
  68. 			if (line >= ano.cnt[i])
  69. 			{
  70. 				for (int j = 0; j < p; ++j) 	ret.cnt[j] -= ano.cnt[i];
  71. 				continue;
  72. 			}
  73. 			bool book = 0;
  74. 			for (int j = 0; j < p; ++j)
  75. 			{
  76. 				if (ret.cnt[j] == line && book == 0)
  77. 				{
  78. 					ret.major[j] = ano.major[i];
  79. 					ret.cnt[j] = ano.cnt[i] - line;
  80. 					book = 1;
  81. 				}
  82. 				else	ret.cnt[j] -= line;
  83. 			} // 分类讨论最小值的处理
  84. 		}
  85. 	}
  86. 	return ret;
  87. }
  89. void maintain(int p)
  90. {
  91. 	t[p].val = merge(t[p << 1].val, t[p << 1 | 1].val);
  92. }
  94. void spread(int p, int l, int r)
  95. {
  96. 	if (t[p].tag)
  97. 	{
  98. 		int mid = (l + r) >> 1;
  99. 		t[p << 1].tag = t[p].tag;
  100. 		t[p << 1 | 1].tag = t[p].tag;
  101. 		t[p << 1].val = make_vote(t[p].tag);
  102. 		t[p << 1].val.cnt[0] = mid - l + 1;
  103. 		t[p << 1 | 1].val = make_vote(t[p].tag);
  104. 		t[p << 1 | 1].val.cnt[0] = r - mid;
  105. 		t[p].tag = 0;
  106. 	}
  107. }
  109. void build(int p, int l, int r)
  110. {
  111. 	if (l == r)
  112. 	{
  113. 		t[p].val = make_vote(isa[l]);
  114. 		return ;
  115. 	}
  116. 	int mid = (l + r) >> 1;
  117. 	build(p << 1, l, mid);
  118. 	build(p << 1 | 1, mid + 1, r);
  119. 	maintain(p);
  120. }
  122. void change(int p, int l, int r, int x, int y, int v)
  123. {
  124. 	if (l > y || r < x) 	return ;
  125. 	if (l >= x && r <= y)
  126. 	{
  127. 		t[p].tag = v;
  128. 		t[p].val = make_vote(v);
  129. 		t[p].val.cnt[0] = r - l + 1;
  130. 		return ;
  131. 	}
  132. 	int mid = (l + r) >> 1;
  133. 	spread(p, l, r);
  134. 	if (mid >= x)	change(p << 1, l, mid, x, y, v);
  135. 	if (mid < y)	change(p << 1 | 1, mid + 1, r, x, y, v);
  136. 	maintain(p);
  137. }
  139. Major_Voting query(int p, int l, int r, int x, int y)
  140. {
  141. 	if (l > y || r < x) 	return make_vote(0);
  142. 	if (l >= x && r <= y)	return t[p].val;
  143. 	int mid = (l + r) >> 1;
  144. 	Major_Voting ret;
  145. 	spread(p, l, r);
  146. 	if (mid >= x)	ret = merge(ret, query(p << 1, l, mid, x, y));
  147. 	if (mid < y)	ret = merge(ret, query(p << 1 | 1, mid + 1, r, x, y));
  148. 	return ret;
  149. } // 以上板子
  151. signed main()
  152. {
  153. 	scanf("%d %d %d", &n, &m, &p);
  154. 	assert(p >= 20 && p <= 100);
  155. 	p = 100 / p;
  156. 	for (int i = 1; i <= n; ++i)	scanf("%d", &isa[i]);
  157. 	build(1, 1, n);
  158. 	for (int i = 0, t, l, r, x; i < m; ++i)
  159. 	{
  160. 		scanf("%d %d %d", &t, &l, &r);
  161. 		if (t == 1)
  162. 		{
  163. 			scanf("%d", &x);
  164. 			change(1, 1, n, l, r, x);
  165. 		}
  166. 		else
  167. 		{
  168. 			Major_Voting ans = query(1, 1, n, l, r);
  169. 			printf("%d", p);
  170. 			for (int i = 0; i < p; ++i) 	printf(" %d", ans.major[i]); // 因为题目的特殊性，我们可以随便乱输出
  171. 			putchar('\n');
  172. 		}
  173. 	}
  174. 	return 0;
  175. }

60.P5071 [Ynoi2015]此时此刻的光辉

珂朵莉给你了一个长为 \(n\) 的序列，有 \(m\) 次查询，每次查询一段区间的乘积的约数个数 \(\bmod 19260817\) 的值。

我们知道一个结论：

若 \(n=\Pi_{i=1}^kp_i^{a_i}\)。

则 \(n\) 的约数个数为 \(\Pi_{i=1}^k(a_i+1)\)。

知道这个结论之后，我们就有一个 自 然 而 然 的想法：

用莫队维护每个质因数的出现次数和答案。

每加入或删除一个质因数，就先乘上它之前出现的次数加一的逆元，然后把出现次数加上或减去一之后再乘。

一算时间复杂度，每个数最多会有 \(9\) 个质因数，那么序列中会有 \(9\times 10^5\) 个数。莫队跑不过去。

但我们发现：小的质因数很多。我们可以抓住这一点进行优化。

由于只用知道每个质因数出现的次数就可以了。

我们可以把这些小的质因数都筛出来，做个前缀和，每次询问都单独统计乘上去。

这里我是把 \(\sqrt{10^9}\) 以内的质数都筛出来了（用来筛其他的大的质因数），然后把前 \(200\) 个的出现次数拿来做前缀和。

之前人傻了写线段树卡了半天常

代码：

  1. #include<cstdio>
  2. #include<vector>
  3. #include<algorithm>
  4. using namespace std;
  5. #define lowbit(x) ((x)&(-(x)))
  6. vector<int> p;
  7. int prime[100010],cntot,tag[100010];
  8. void search(int x)//线性筛质数
  9. {
  10. 	tag[1]=1;
  11. 	for(int i=2;i<=x;++i)
  12. 	{
  13. 		if(!tag[i])	prime[++cntot]=i;
  14. 		for(int j=1;j<=cntot&&prime[j]*i<=x;++j)
  15. 		{
  16. 			tag[i*prime[j]]=1;
  17. 			if(i%prime[j]==0)	break;
  18. 		}
  19. 	}
  20. }
  21. struct ask
  22. {
  23. 	int l,r,ID,lyc,wgy;//lyc和wgy表示这个询问在原序列上的左右端点，用来之后查询小的质因数的出现个数
  24. }as[100010];
  25. const int each=700,mod=19260817;
  26. int n,m,a[100010],cnt[100010][201],nodes[100010][201],ton[201],pri[200010],belong[200010],ans[100010],in[100010],out[100010],inv[200010],tot,opl,opr,nowl,nowr,app[200010],cur=1;
  27. bool cmp(ask one,ask ano)
  28. {
  29. 	if(belong[one.l]^belong[ano.l])	return belong[one.l]<belong[ano.l];
  30. 	if(belong[one.l]&1)	return one.r<ano.r;//奇偶排序
  31. 	else	return one.r>ano.r;
  32. }
  33. void read(int &hhh)
  34. {
  35. 	int x=0;
  36. 	char c=getchar();
  37. 	while((c<'0')|(c>'9')&&c^'-')	c=getchar();
  38. 	if(c^'-')	x=c^'0';
  39. 	char d=getchar();
  40. 	while((d<='9')&(d>='0'))
  41. 	{
  42. 		x=(x<<1)+(x<<3)+(d^'0');
  43. 		d=getchar();
  44. 	}
  45. 	if(c^'-')	hhh=x;
  46. 	else	hhh=-x;
  47. }
  48. void writing(int x)
  49. {
  50. 	if(!x)	return;
  51. 	writing(x/10);
  52. 	putchar((x%10)+'0');
  53. }
  54. void write(int x)
  55. {
  56. 	if(x<0)
  57. 	{
  58. 		putchar('-');
  59. 		x=-x;
  60. 	}
  61. 	else if(!x)
  62. 	{
  63. 		putchar('0');
  64. 		putchar('\n');
  65. 		return;
  66. 	}
  67. 	writing(x);
  68. 	putchar('\n');
  69. }
  70. int main()
  71. {
  72. 	read(n);
  73. 	read(m);
  74. 	search(33000);
  75. 	for(int i=1;i^201;++i)	ton[i]=prime[i];
  76. 	for(int i=1;i<=n;++i)
  77. 	{
  78. 		read(a[i]);
  79. 		in[i]=tot+1;//in[],out[]记录每个原序列的元素质因数在新的序列中的开头结尾
  80. 		for(int j=1;j^201;++j)
  81. 		{
  82. 			while(a[i]%ton[j]==0)
  83. 			{
  84. 				a[i]/=ton[j];
  85. 				++cnt[i][j];
  86. 			}
  87. 		}
  88. 		if(a[i]>1)
  89. 		{
  90. 			for(int j=201;prime[j]*prime[j]<=a[i];++j)//把大的质因数拿出来作为一个新的序列
  91. 			{
  92. 				while(a[i]%prime[j]==0)
  93. 				{
  94. 					a[i]/=prime[j];
  95. 					pri[++tot]=prime[j];
  96. 				}
  97. 			}
  98. 			if(a[i]>1)	pri[++tot]=a[i];
  99. 		}
  100. 		out[i]=tot;
  101. 	}
  102. 	inv[1]=1;
  103. 	for(int i=2;i<=tot;++i)	inv[i]=(((-(long long)(mod/i)*inv[mod%i])%mod)+mod)%mod;//线性处理逆元
  104. 	for(int i=1;i<=tot;++i)
  105. 	{
  106. 		p.push_back(pri[i]);
  107. 		belong[i]=(i-1)/each+1;
  108. 	}
  109. 	sort(p.begin(),p.end());
  110. 	p.erase(unique(p.begin(),p.end()),p.end());
  111. 	for(int i=1;i<=tot;++i)	pri[i]=lower_bound(p.begin(),p.end(),pri[i])-p.begin()+1;
  112. 	for(int i=1;i^201;++i)//对小的质因数做前缀和
  113. 	{
  114. 		for(int j=1;j<=n;++j)	nodes[j][i]=nodes[j-1][i]+cnt[j][i];
  115. 	}
  116. 	for(int i=1;i<=m;++i)
  117. 	{
  118. 		read(opl);
  119. 		read(opr);
  120. 		as[i].l=in[opl];
  121. 		as[i].r=out[opr];//莫队的询问左右端点，把原序列中查询区间的质因数全部覆盖掉。
  122. 		as[i].ID=i;
  123. 		as[i].lyc=opl;
  124. 		as[i].wgy=opr;
  125. 	}
  126. 	sort(as+1,as+1+m,cmp);
  127. 	nowl=1;
  128. 	for(int i=1;i<=m;++i)
  129. 	{
  130. 		while(nowl>as[i].l)
  131. 		{
  132. 			--nowl;
  133. 			cur=((long long)cur*inv[app[pri[nowl]]+1])%mod;//先乘上以前的逆元，消除之前的答案贡献
  134. 			++app[pri[nowl]];
  135. 			cur=((long long)cur*(app[pri[nowl]]+1))%mod;//乘上现在的答案贡献
  136. 		}
  137. 		while(nowr<as[i].r)
  138. 		{
  139. 			++nowr;
  140. 			cur=((long long)cur*inv[app[pri[nowr]]+1])%mod;
  141. 			++app[pri[nowr]];
  142. 			cur=((long long)cur*(app[pri[nowr]]+1))%mod;
  143. 		}
  144. 		while(nowl<as[i].l)
  145. 		{
  146. 			cur=((long long)cur*inv[app[pri[nowl]]+1])%mod;
  147. 			--app[pri[nowl]];
  148. 			cur=((long long)cur*(app[pri[nowl]]+1))%mod;
  149. 			++nowl;
  150. 		}
  151. 		while(nowr>as[i].r)
  152. 		{
  153. 			cur=((long long)cur*inv[app[pri[nowr]]+1])%mod;
  154. 			--app[pri[nowr]];
  155. 			cur=((long long)cur*(app[pri[nowr]]+1))%mod;
  156. 			--nowr;
  157. 		}
  158. 		int tmp=cur;
  159. 		for(int j=1;j^201;++j)	tmp=((long long)tmp*(nodes[as[i].wgy][j]-nodes[as[i].lyc-1][j]+1))%mod;//临时乘上小的质因数对答案的贡献
  160. 		ans[as[i].ID]=tmp;
  161. 	}
  162. 	for(int i=1;i<=m;++i)	write(ans[i]);
  163. 	return 0;
  164. }

Ds100p -「数据结构百题」51~60的更多相关文章

  1. [luogu_P1251][LOJ#6008]「网络流 24 题」餐巾计划

    [luogu_P1251][LOJ#6008]「网络流 24 题」餐巾计划 试题描述 一个餐厅在相继的 \(N\) 天里,第 \(i\) 天需要 \(R_i\) 块餐巾 \((i=l,2,-,N)\) ...

  2. LOJ6003 - 「网络流 24 题」魔术球

    原题链接 Description 假设有根柱子,现要按下述规则在这根柱子中依次放入编号为的球. 每次只能在某根柱子的最上面放球. 在同一根柱子中,任何2个相邻球的编号之和为完全平方数. 试设计一个算法 ...

  3. LOJ6002 - 「网络流 24 题」最小路径覆盖

    原题链接 Description 求一个DAG的最小路径覆盖,并输出一种方案. Solution 模板题啦~ Code //「网络流 24 题」最小路径覆盖 #include <cstdio&g ...

  4. LOJ6001 - 「网络流 24 题」太空飞行计划

    原题链接 Description 有个实验和个仪器,做实验有报酬买仪器有花费.每个实验都需要一些仪器,求最大净收益(实验报酬仪器花费),并输出一组方案. Solution 实验向所需仪器连边,实验的点 ...

  5. LOJ6000 - 「网络流 24 题」搭配飞行员

    原题链接 题意简述 求二分图的最大匹配. 题解 这里写的是匈牙利算法. 表示节点的当前匹配. 为真表示在这一轮匹配中,无法给节点一个新的匹配.所以如果为真就不用再dfs它了,直接continue就好. ...

  6. LibreOJ #6014. 「网络流 24 题」最长 k 可重区间集

    #6014. 「网络流 24 题」最长 k 可重区间集 内存限制:256 MiB时间限制:1000 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: 匿名 提交提交记录统计讨论测试数据   ...

  7. Libre 6013 「网络流 24 题」负载平衡 （网络流，最小费用最大流）

    Libre 6013 「网络流 24 题」负载平衡 (网络流,最小费用最大流) Description G 公司有n 个沿铁路运输线环形排列的仓库,每个仓库存储的货物数量不等.如何用最少搬运量可以使n ...

  8. Libre 6012 「网络流 24 题」分配问题 （网络流，费用流）

    Libre 6012 「网络流 24 题」分配问题 (网络流,费用流) Description 有n件工作要分配给n个人做.第i个人做第j件工作产生的效益为\(c_{ij}\).试设计一个将n件工作分 ...

  9. Libre 6011 「网络流 24 题」运输问题 （网络流，最小费用最大流）

    Libre 6011 「网络流 24 题」运输问题 (网络流,最小费用最大流) Description W 公司有m个仓库和n个零售商店.第i个仓库有\(a_i\)个单位的货物:第j个零售商店需要\( ...

  10. LibreOJ #6013. 「网络流 24 题」负载平衡 最小费用最大流 供应平衡问题

    #6013. 「网络流 24 题」负载平衡 内存限制:256 MiB时间限制:1000 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: 匿名 提交提交记录统计讨论测试数据   题目描述 ...

随机推荐

  1. Dapr在Java中的实践 之 状态管理

    状态管理 状态管理(State Management)使用键值对作为存储机制,可以轻松的使长时运行.高可用的有状态服务和无状态服务共同运行在我们的服务中. 我们的服务可以利用Dapr的状态管理API在 ...

  2. shell编程-发送消息

    需求:利用 Linux 自带的 mesg 和 write 工具,编写一个向用户快速发送消息的脚本,输入用户名作为第一个参数,消息内容为第二个参数.脚本需要检测用户是否登录,是否打开消息功能,以及当前发 ...

  3. docker镜像的原理

    docker镜像的原理 docker镜像是由特殊的文件系统叠加而成 最低端是bootfs,并使用宿主机的bootfs 第二层是root文件系统rootfs,称之为base image 再往上是可叠加的 ...

  4. Python潮流周刊#6：Python 3.12 有我贡献的代码！

    你好,我是猫哥.这里记录每周值得分享的 Python 及通用技术内容,部分为英文,已在小标题注明.(标题取自其中一则分享,不代表全部内容都是该主题,特此声明.) 首发于我的博客,https://pyt ...

  5. @Target元注解的使用

    @Target注解标记另外的注解用于限制此注解可以应用哪种Java元素类型.先看Java SE 8中@Target是如何声明的: package java.lang.annotation; publi ...

  6. XTTS系列之三：中转空间的选择和优化

    通常选择XTTS做迁移的数据库都不会太小的,至少都是几T.几十T这样的规模,这种级别的数据量原有空间不够用,所以在迁移过程临时用作存放迁移数据库备份文件的空间也是需要提前考虑规划的问题. 最近就有客户 ...

  7. 提升性能的利器：深入解析SectionReader

    一. 简介 本文将介绍 Go 语言中的 SectionReader,包括 SectionReader的基本使用方法.实现原理.使用注意事项.从而能够在合适的场景下,更好得使用SectionReader ...

  8. JDBC的增删改-结果集的元数据-Class反射-JDBC查询封装

    一.使用JDBC批量添加 ​ 知识点复习: ​1.JDBC的六大步骤 (导入jar包, 加载驱动类,获取连接对象, 获取sql执行器.执行sql与并返回结果, 关闭数据库连接) 2.​封装了一个DBU ...

  9. 2023-07-17：给定一个数组arr，长度为n， 再给定一个数字k，表示一定要将arr划分成k个集合， 每个数字只能进一个集合。 返回每个集合内部的平均值都累加起来最小的值。 平均值向下取整。 1

    2023-07-17:给定一个数组arr,长度为n, 再给定一个数字k,表示一定要将arr划分成k个集合, 每个数字只能进一个集合. 返回每个集合内部的平均值都累加起来最小的值. 平均值向下取整. 1 ...

  10. 投个 3D 冰壶，上班玩一玩

    ​本篇文章将介绍如何使用物理引擎和图扑 3D 可视化技术来呈现冰壶运动的模拟. Oimo.js 物理引擎 Oimo.js 是一个轻量级的物理引擎,它使用 JavaScript 语言编写,并且基于 Oi ...