洛谷 P3865 【模板】ST表(倍增, ST表)

倍增, ST表

ST算法， 在 O(n*log(n)) 的时间预处理后，以 O(1) 的时间回答：
数列 A 的区间 [l, r] 之间的最大值。
本题要点：
1、f[i][j] 表示从i开始，2^j 个数范围内，区间最大值
2、把长度是 2^j 的区间，分为左右两半， 各长 2^(j-1), 2^j 的区间的最大值，就是两个区间的较大值。
f[i][j] = max(f[i][j - 1], f[i + (1 << (j - 1))][j - 1])
3、查询的时候，寻找到第一个 k， 使得 2^k <= r - l + 1 < 2^(k+1), 然后
f[l][r] = max(f[l][k], f[r - (1 << k) + 1][k])

#include <cstdio>
#include <cstring>
#include <iostream>
#include <cmath>
using namespace std;
const int MaxN = 100010;
int a[MaxN];
int f[MaxN][20];	// f[i][j] 表示从i开始，2^j 个数范围内，区间最大值
int n, m;

inline int read()
{
	int x = 0,f = 1; char ch = getchar();
	while (!isdigit(ch))
	{
		if (ch == '-') 
			f = -1; ch = getchar();
	}
	while (isdigit(ch))
	{
		x = x * 10 + ch - 48; ch = getchar();
	}
	return x * f;
}

void ST_prework()
{
	for(int i = 1; i <= n; ++i)
	{
		f[i][0] = a[i];
	}
	int t = log(n) / log(2) + 1;
	for(int j = 1; j < t; ++j)
	{
		for(int i = 1; i <= n - (1 << j) + 1; ++i)
		{
			f[i][j] = max(f[i][j - 1], f[i + (1 << (j - 1))][j - 1]);
		}
	}
}

int ST_query(int l, int r)
{
	int k = log(r - l + 1) / log(2);
	return max(f[l][k], f[r - (1 << k) + 1][k]);
}

int main()
{
	int x, y;
	n = read(), m = read();
	for(int i = 1; i <= n; ++i)
		a[i] = read();
	ST_prework();
	for(int i = 0; i < m; ++i)
	{
		x = read(), y = read();
		printf("%d\n", ST_query(x, y));
	}
	return 0;
}

/*
8 8
9 3 1 7 5 6 0 8
1 6
1 5
2 7
2 6
1 8
4 8
3 7
1 8
*/

/*
9
9
7
7
9
8
7
9
*/