RMQ 区间最值问题

rmq：指区间最值问题。就是给你一个数列，然后不断地询问在区间[a,b]中的最大值/最小值是什么。

 

可以看例题：

tyvj 的 1297题

 

数据规模：

40%的数据：1<=n,p<=1,000

100%的数据：1<=n<=50,000    1<=p<=20,000

n是数列大小，p是询问数

 

对于这类问题：

1.朴素算法：遍历区间[a,b]找最值 O(np)

 

2.线段树：不会

 

3.简单分段法：

 

例如上面的例题，预处理以100个元素为1段，最后剩下不足100的为1段。

设pos(n) = n/100既n所在的段数。那么对于[a,b],pa = pos(a)+1,pb = pos[b]-1

pa*100到pb*100可以一段一段地找，a到pa*100与pb*100到b则遍历。

注意如果a,b很近，没有完整的段，就要直接遍历了。

 

这个方法容易想到，效率也有较不错的改进。例题90分。最后一个点超了

 

#include <iostream> #include <stdio.h> using namespace std; #define POS(n) (n/100) const int _DEF_MAXNUM = 0x7fffffff; int _fmax( int a,int b ) { return a > b ? a : b; } int _fmin( int a,int b ) { return a < b ? a : b; } struct _cinfo { int max,min; }; int _ghigh[60000],_gn,_gp; _cinfo _gset[600]; int main() { cin >> _gn >> _gp; int max = 0,min = _DEF_MAXNUM,cnt = 0; for( int i = 1;i <= _gn;++i ) { if( i%100 == 0 ) //new { _gset[cnt].min = min,_gset[cnt].max = max; ++cnt; max = 0,min = _DEF_MAXNUM; } scanf("%d",&_ghigh[i]); min = _fmin( min,_ghigh[i] ); max = _fmax( max,_ghigh[i] ); } _gset[cnt].min = min,_gset[cnt].max = max; ++cnt; int a,b; for( int i = 0;i < _gp;++i ) { scanf("%d %d",&a,&b); int min = _DEF_MAXNUM,max = 0,pa = POS(a)+1,pb = POS(b)-1; if( b - a < 300 ) { while( a <= b ) { min = _fmin( _ghigh[a],min ); max = _fmax( _ghigh[a],max ); ++a; } } else { while( pa <= pb ) { min = _fmin( _gset[pa].min,min ); max = _fmax( _gset[pa].max,max ); ++pa; } for( int beg = a;beg <= pa*100;++beg ) { min = _fmin( min,_ghigh[beg] ); max = _fmax( max,_ghigh[beg] ); } for( int beg = pb*100;beg <= b;++beg ) { min = _fmin( min,_ghigh[beg] ); max = _fmax( max,_ghigh[beg] ); } } printf("%d %d/n",max,min); } return 0; }  

 

4.st算法：

 

其实是动态规划。设d[i,j]表示[i,i+2^j-1]这个区间内的最大值，那么在询问到[a,b]区间的最大值时答案就是max(d[a,k], d[b-2^k+1,k])，其中k是满足2^k<=b-a的最大的k,即k=[log(b-a+1)/log(2)]

d的求法可以用动态规划，d[i,j]=max(d[i,j-1],d[i+2^(j-1),j-1])

 

边界:d[i,0] 表示 i 到 i.所以d[i,0] = h[i]

 

效率O(nlogn).100分

 

#include <iostream> #include <stdio.h> #include <math.h> using namespace std; int _fmax( int a,int b ) { return a > b ? a : b; } int _fmin( int a,int b ) { return a < b ? a : b; } int _gstmax[60000][30],_gstmin[60000][30],_gn,_gp; int main() { scanf("%d %d",&_gn,&_gp); for( int i = 1;i <= _gn;++i ) { scanf("%d",&_gstmax[i][0]); _gstmin[i][0] = _gstmax[i][0]; } int lev = (int)( log(_gn) / log(2) ); //初始化 for( int k = 1;k <= lev;++k ) { for( int i = 1;i <= _gn-(1<<k)+1;++i ) { _gstmax[i][k] = _fmax( _gstmax[i][k-1],_gstmax[i+(1<<(k-1))][k-1] ); _gstmin[i][k] = _fmin( _gstmin[i][k-1],_gstmin[i+(1<<(k-1))][k-1] ); } } int a,b,max,min; for( int i = 0;i < _gp;++i ) { scanf("%d %d",&a,&b); int t = (int)( log(b-a+1)/log(2) ); max = _fmax( _gstmax[a][t],_gstmax[b-(1<<t)+1][t] ); min = _fmin( _gstmin[a][t],_gstmin[b-(1<<t)+1][t] ); printf("%d %d/n",max,min); } return 0; }