【USACO】 Balanced Lineup

【题目链接】

点击打开链接

【算法】

这是一道经典的最值查询（RMQ)问题。

我们首先想到线段树。但有没有更快的方法呢？对于这类问题，我们可以用ST表（稀疏表）算法求解。

稀疏表算法。其实也是一种动态规划的算法。是先做一遍预处理，然后O（1）求出答案。
        
        设计状态 : f[i][j] 表示从第i个数开始连续2^j次方个数（包括第i个数），中的（最大或最小值）

那么状态转移方程是什么呢？

我们知道 2^j可分解为两个2^(j-1)，所以f[i][j] = max或min(f[i][j-1],f[i+2^(j-1)][j-1])

做完预处理，我们又该如何查询呢？

首先我们证明 ： 2^log(x) > x / 2

先将x表示为2^i+k的形式，则log(x)=i

得到： 2^i > 2^(i-1) + k / 2
        2^(i-1) > k / 2
        2^i > k

我们知道2^i一定大于k，所以成立。

所以在查询（最大或最小值）时，我们只需找两个长度为2^log(x)的区间求最大最小值就可以了

此题堪称ST表裸题

【代码】

#include<bits/stdc++.h>

using namespace std;

int i,j,N,M,k,x,y;
int a[],f_max[][],f_min[][];

int main() {

        scanf("%d%d",&N,&M);
        for (i = ; i <= N; i++) {
                scanf("%d",&a[i]);
                f_max[i][] = f_min[i][] = a[i];
        }

        for (i = N; i >= ; i--) {
                for (j = ; i + (<<j) -  <= N; j++) {
                        f_max[i][j] = max(f_max[i][j-],f_max[i+(<<(j-))][j-]);
                        f_min[i][j] = min(f_min[i][j-],f_min[i+(<<(j-))][j-]);
                }
        }

        for (i = ; i <= M; i++) {
                scanf("%d%d",&x,&y);
                k = (int)(log(y - x + ) / log(2.0));
                printf("%d\n",max(f_max[x][k],f_max[y-(<<k)+][k]) - min(f_min[x][k],f_min[y-(<<k)+][k]));
        }

        return ;

}