LOJ6285 数列分块入门9(分块毒瘤)

题目描述:

给出一个长为 的数列,以及 个操作,操作涉及询问区间的最小众数。

输入格式

第一行输入一个数字 。

第二行输入 个数字,第 个数字为 ,以空格隔开。

接下来输入 行询问,每行输入两个数字 、,以空格隔开。

表示查询位于 的数字的众数。

输出格式

对于每次询问,输出一行一个数字表示答案。

样例

样例输入

4
1 2 2 4
1 2
1 4
2 4
3 4

样例输出

1
2
2
2

题意:

给出区间内的最小众数

分块,离散化每个数,开 v e c t o r vector vector 记录每个数 p p p 出现的位置,这样就能二分出 L , R L,R LR 以内 p p p 的个数了。先预处理出任意两块的众数 f [ i ] [ j ] f[i][j] f[i][j],这样众数就是 f [ i ] [ j ] f[i][j] f[i][j] 和旁边两块数中的其中一个了,直接遍历这些数即可。但是这道题不像普通的分块开 n \sqrt{n} n 就行,这个得把块再缩小,把时间复杂度优化。

AC代码:

const int N = 1e5 + 10;
struct Block
{
	int l, r;
} b[N];
int a[N], belong[N];
int f[10000][10000];
int num[N];
int n, block;
vector<int> vv;
vector<int> pos[N];
void init()
{
	rep(i, 1, belong[n])
	{
		rep(j, 0, vv.size())
			num[j] = 0;
		int mode = INF, NUM = 0;
		rep(j, b[i].l, n)
		{
			num[a[j]]++;
			if (num[a[j]] > NUM || (num[a[j]] == NUM && a[j] < mode))
			{
				mode = a[j];
				NUM = num[a[j]];
			}
			f[i][belong[j]] = mode;
		}
	}
}
int getNum(int l, int r, int v)
{
	int t = upper_bound(pos[v].begin(), pos[v].end(), r) - lower_bound(pos[v].begin(), pos[v].end(), l);
	return t;
}
int query(int l, int r)
{
	int bl = belong[l], br = belong[r];
	int ans = INF, NUM = 0;
	if (bl == br)
	{
		rep(i, l, r)
		{
			int tot = getNum(l, r, a[i]);
			if (tot > NUM || (tot == NUM && a[i] < ans))
			{
				ans = a[i];
				NUM = tot;
			}
		}
	}
	else
	{
		rep(i, l, b[bl].r)
		{
			int tot = getNum(l, r, a[i]);
			if (tot > NUM || (tot == NUM && a[i] < ans))
			{
				ans = a[i];
				NUM = tot;
			}
		}
		if (bl + 1 <= br - 1)
		{
			int v = f[bl + 1][br - 1];
			int tot = getNum(l, r, v);
			if (tot > NUM || (tot == NUM && v < ans))
			{
				ans = v;
				NUM = tot;
			}
		}
		rep(i, b[br].l, r)
		{
			int tot = getNum(l, r, a[i]);
			if (tot > NUM || (tot == NUM && a[i] < ans))
			{
				ans = a[i];
				NUM = tot;
			}
		}
	}
	return vv[ans - 1];
}

int main()
{
	sd(n);
	vv.clear();
	rep(i, 1, n)
	{
		sd(a[i]);
		vv.pb(a[i]);
	}
	sort(vv.begin(), vv.end());
	vv.erase(unique(vv.begin(), vv.end()), vv.end());
	rep(i, 1, n)
	{
		a[i] = lower_bound(vv.begin(), vv.end(), a[i]) - vv.begin() + 1;
	}
	rep(i, 0, vv.size())
	{
		pos[i].clear();
	}
	block = 30;
	rep(i, 1, n)
	{
		belong[i] = (i - 1) / block + 1;
		pos[a[i]].push_back(i);
	}
	rep(i,1,belong[n])
	{
		b[i].l = (i - 1) * block + 1;
		b[i].r = b[i].l + block - 1;
	}
	b[belong[n]].r = n;
	init();
	rep(i, 1, n)
	{
		int l, r;
		sdd(l,r);
		pd(query(l, r));
	}
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值