倍增思想和ST算法

这是我从书上抄下来的一段话，括号内是自己的理解。

当我们递推的时候状态空间很大，线性递推无法满足时间与空间复杂度，那么我们可以通过成倍增长的方式只递推出状态空间中在2的整数次幂位置上的值作为代表 （这个代表能够代表之前倍增达到的代表值或者代表这一段到上一个2的整数次幂区域的代表值，根据题目而异） ， （如果这个位置不能够代表或者需要其他位置作为代表） ，那么可以通过“任意整数都可以表示成若干个2的次幂项的和”这一性质，使用之前的代表值拼成所需要的值，所以倍增算法要求我们递推的问题状态空间关于2的次幂具有可划分性 （像是值递增数组） 。

问题：
给定数列A[n]（A[i]>0 ），S是A的前缀和，每次给定一个T，求最大的k，使得S[k] <= T。
1.最朴素的算法就是从前到后枚举
2.二分
3.如果k靠前那么二分就会耗费很大的时间，还不如枚举，所以可以用倍增的思想

大致代码就是这样的

#include<iostream>
using namespace std;
int main()
{
	int sum=0,p=1,k=0,T;
	int a[100]={0},s[100]={0};
	for(int i=1;i<=10;i++){
		a[i] = i;
		s[i] += s[i-1]+a[i];
	}
	cin>>T;
	while(1){
		if(!p) break;
		sum = s[k+p];
		if(sum > T) p /= 2,sum = s[k];
		else k += p,p <<= 1;
	}
	cout<<sum;
	
	return 0;
 } 

接下来讲ST算法

ST算法是一个能在O(NlogN)的时间预处理以后，以O(1)的时间复杂度在线回答“区间最值”的算法。
一个序列的子区间有N²个，根据倍增思想，我们首先在这个规模为O(N²)的状态空间里选择一些2的整数次幂的位置作为代表值。

我们可以设F[i,j]表示在序列A中在下标在子区间[i,i+2^(j-1)]里的数的最大值(通俗来讲这个i就是从A的那个下标开始算，也就是L，j就是表示从i开始的 2^j个数)，递推边界就是F[i,0]=A[i]，每次将等倍数的范围求预处理，再翻倍求2倍于前面的倍数的范围，这样就得到递推公式
F[i,j] = max(F[i,j-1],F[i+2^(j-1),j-1]) 2²
max 里面的“j-1”意思就是小于F[i,j]一倍倍数范围的范围，i和i+2^(j-1)就是F[i,j]中的左半边子区间和右半边子区间，因为区间大小都是以2的整数次幂来计算的，所以就有了此公式。

//预处理区间 
void ST_pre(){
	for(int i=1;i<=n;i++) F[i][0] = A[i];
	//while(!temp) temp<<1,t++;
	t = log(n) / log(2) + 1;
	for(int i=1;i<=t;i++){
		for(int j=1;j<=n-(1<<j)+1;j++){
			F[i][j] = max(F[i][j-1],F[i+(1<<(j-1))][j-1]);
		}
	} 
} 

查询：给出区间 [i,j] 首先先求s=(r-l+1)，s是这个区间的范围，所以只需要求出一个k满足2^k <= s <=2^(k+1) 因为这个k+1的范围是k范围的两倍，所以k的前一段和k的后一段就是k+1的那一段，而且s<=2^(k+1) 所以s中的 l ~ l+2^k 和 r-2^k+1 ~ r 一定包含了 l ~ r 。

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]);
}