【题解】最小值

题目描述

有一个长度为 (n) 的序列,初始时序列中的数全为 (2^{31}-1)

(m) 次操作,第 (i) 次操作为将序列中第 (a_i) 个数修改为 (b_i)

记第 (i) 次操作后序列中的最小值为 (s_i),你需要输出 (sumlimits_{i=1}^m s_i imes 10099^i)

(a_i)(b_i) 用以下方法确定:

输入整数 (x_0)(x_1)(a)(b)(c),令 (x_i=(ax_{i-2}+bx_{i-1}+c) mod 2^{32}quad (ige 2)),则 (a_i=leftlfloordfrac{x_{2i-1}}{4} ight floor mod n)(b_i=leftlfloordfrac{x_{2i}}{4} ight floor)

输入格式

一行七个整数 (n)(m)(x_0)(x_1)(a)(b)(c)

输出格式

一行一个整数表示答案。

数据范围

测试时间限制 (1000 mathrm{ms}),空间限制 (256 mathrm{MiB})

  • 对于 (10\%) 的数据,(1le n,mle 1000)
  • 对于 (50\%) 的数据,(1le n,mle 10^5)
  • 对于 (100\%) 的数据,(1le n,mle 10^7)(x_0)(x_1)(a)(b)(c)(left[0,2^{31}-1 ight)) 中均匀随机。

分析

注意到这题要求的操作很少,所以我们要考虑一些比较大胆的做法。

(mathtt{10 pts})

很简单,记录一个数组,每一次修改以后的暴力扫一遍最小值,复杂度 (Theta(1)-Theta(n)),稳稳超时。

(mathtt{50 pts})

也很简单,随便套一个 (mathcal{O}(log n)) 的数据结构乱搞就行了。

(mathtt{100 pts})

这个就有点难想到了。

根据数据范围,我们八成要设计一种查询 (Theta(1)) 的算法。也就是单纯的修改 (Theta(1)),查询 (Theta(1))

在脑中搜刮一下后,好像支持修改、查询,修改是 (Theta(1)) 的就只有暴力了。

但在一般的题目中,毒瘤出题人都会设计数据卡暴力,使之在查询时跑得奇慢无比。

但是……我们发现,这题的操作都是伪随机生成的。难道……

考虑这样的算法:

维护数组 (A),记录第 (i) 位的值 (A_i)

修改:直接修改 (A_i)

查询:

  • 如果不是最小值所在的位置,则更新最小值和其所在位置。

  • 如果是的话,则暴力更新最小值。

利用伪随机的性质,修改到最小值的概率为 (dfrac{m}{n}),则复杂度的数学期望

[egin{aligned}E(x)&=dfrac{m}{n} imes mathcal{O}(n)+left(m-dfrac{m}{n} ight) imesmathcal{O}(1)\&=mathcal{O}(m)+mathcal{O}(m)\&=mathcal{O}(m)end{aligned} ]

可以通过本题。

Code

顺带一提,这道题有点卡空间……开 long long 就炸了。

#include <cstdio>
#include <climits>
using namespace std;

typedef unsigned int ui;
const int max_n = 10000000;
const ui mul = 10099;

ui nums[max_n];

int main()
{
	int n, m, min_pos = -1;
	ui x0, x1, a, b, c, ans = 0, tk = 1, xt, ai, bi, min_val = INT_MAX;
	
	scanf("%d%d%u%u%u%u%u", &n, &m, &x0, &x1, &a, &b, &c);
	
	for (int i = 0; i < n; i++)
		nums[i] = INT_MAX;
	
	for (int i = 0; i < m; i++)
	{
		ai = (x1 / 4) % n;
		
		xt = a * x0 + b * x1 + c;
		bi = xt / 4, x0 = x1, x1 = xt;
		
		xt = a * x0 + b * x1 + c;
		x0 = x1, x1 = xt;
		
		nums[ai] = bi;
		if (ai == min_pos && bi > min_val)
		{
			min_val = INT_MAX, min_pos = -1;
			for (int i = 0; i < n; i++)
				if (nums[i] < min_val)
				{
					min_val = nums[i];
					min_pos = i;
				}
		}
		else if (bi < min_val)
		{
			min_val = bi;
			min_pos = ai;
		}
		
		tk *= mul;
		ans += tk * min_val;
	}
	
	printf("%u
", ans);
	
	return 0;
}

后记

这道题告诉了我们一个很有意思的道理——复杂度越大,适用范围越广。

同样地,复杂度小的,适用范围会小一点。

解题时,一般要选择适合的算法,即题目在适用范围之内,而且复杂度较优。这题就是一个很好的例子。

原文地址:https://www.cnblogs.com/5ab-juruo/p/solution-20200229-min.html