一、倍增算法
要了解倍增之前,强烈建议大家先看一下这位大佬对倍增的解释:【白话系列】倍增算法
看完以后相信你已经对倍增有了大致初步的了解,下面给出倍增的定义
倍增 从字面的上意思看就是成倍的增长 ,这是指我们在进行递推时,如果状态空间很大,通常的线性递推无法满足时间和空间复杂度的要求 ,那么我们就可以通过成倍的增长,只递推状态空间中在 2 的整数次幂位置上的值作为代表 。当需要其他位置上的值时,我们只需要通过" 任意整数可以表示成若干个2的次幂项的和 " 这一性质( 13 = 2 3 + 2 2 + 2 0 13 = 2^3 + 2^2 +2^0 13=23+22+20 ), 使用之前求出的代表值拼成所需的值。
给定一个长度为 N 的数列 A ,然后进行若干次查询 , 每一次给定一个整数 T , 求出最大的 k , 满足 ∑ 1 = 1 k A [ i ] < = T \sum_{1 = 1 }^{k} A[i] <=T ∑1=1kA[i]<=T . 算法必须是在线的(每给一次询问,就给出结果) ;
我们当然可以先 预处理前缀和 , ( A [ i ] > 0 A[i] > 0 A[i]>0 )然后用 二分找到 这个k , 这个复杂度最坏是 O ( n ) O(n) O(n) , 因为如果我们的 询问的T 太小 ,就还不如从前往后枚举。
我们可以设计这样一个 倍增算法 :
1 . 令 p = 1 , k = 0 , sum = 0 ;
2 . 比较" A 数组中 k 之后的 p 个数的 和" 与 T 的关系 , 也就是说 ,如果 s u m + s [ k + p ] − s [ k ] < = T sum + s[k+p] - s[k] <=T sum+s[k+p]−s[k]<=T , 则令 s u m + = s [ k + p ] − s [ k ] , k + = p , p ∗ = 2 ; sum += s[k+p] -s[k] , k+=p , p*=2 ; sum+=s[k+p]−s[k],k+=p,p∗=2; 也就是累加上这 p 个数的和 ,然后把 p 的跨度增长了一倍。如果 s u m + s [ k + p ] − s [ k ] > T sum+ s[k+p] - s[k] > T sum+s[k+p]−s[k]>T , 则令 p / = 2 p /=2 p/=2 .
3 重复上一步,直到 p 的值变为0 , 此时 k 就是 答案 .
注意这里我加了一个特判是因为如果p经过神奇的运算会小于k,那么代码就会直接死循环,不信你把那个特判删了然后试试我编的这组数据:
输入:
5
1 2 3 4 5
10
输出:
4
然后就是代码了(注意这里只是简单的应用帮助理解倍增,并不一定是本题的正确解法)
#include<iostream>
#include<string.h>
#include<cstdio>
#include<math.h>
#include<map>
#define ls (p<<1)
#define rs (p<<1|1)
#define over(i,s,t) for(register long long i=s;i<=t;++i)
#define lver