倍增搜索

二分搜索适合既有下界也有上界,对于只知道一边的界情况,可以用倍增搜索法

int search(vector<int> A,  int x) {
    for (int i = 0; A[i] <= x;) {
        if (A[i] == x) return i;
        if (i + 1 == A.size() || A[i + 1] > x) return -1;
        for (int k = 1; i + k < A.size() && A[i + k] <= x;  k += k)
            i += k;
    }
    return -1;
}

类似问题有用除法的定义(只用加法)实现除法。

循环不变式:如果数组中存在x, 它必然在大于等于位置i 的地方

算法框架:

1)如果起点处i的值A[i] > x,否则说明这个范围内不存在x, 返回-1

2)如果起点A[i]==x,找到,返回位置i

      否则就是A[i] < x:

          这时候如果 A[i + 1] > x, 也说明不存在,返回-1

  否则进行倍增查找,当i跳某个步长k后的值A[i+k] 大于x ,本轮倍增结束,一直保持着不变式 A[i] <= x ,主循环下一轮相当于是解决一个相同的问题,只不过i已经推进了很多,离目标近了一些。用递归看更清晰:

int search(vector<int> &A,  int i, int x) {
    if (A[i] > x) return -1;
    if (A[i] == x) return i;
    if (i + 1 == A.size() || A[i + 1] > x) return -1;
    for (int k = 1; i + k < A.size() && A[i + k] <= x; i += k, k += k);
    return search(A, i, x);
}


另一种更通用的模板:


int search(vector<int> A,  int x) {
    int i = -1; //i is the farest position currently searched
    for (;;) {
        int k = 1;
        for (; i + k < A.size() && A[i + k] <= x;  k += k) {
            i += k;
        }
        if (i >= 0 && A[i] == x) return i;
        if (k == 1) return -1;
    }
}


bool isPowerOfThree(int n) {
        int x = 1; //closed position that not exceed
        for (;;) {
            int step = 3;
            for (; n / step >= x; step *= step) {
                x *= step;
            }//exit: could not add the current step
            if (x == n) return true;
            if (step == 3) return false; //could not go one step further, or will exceed
        }
    }





倍增LCA(Least Common Ancestor,最近公共祖先)是求解树上最近公共祖先问题的一种经典算法,它是朴素算法的改进算法,通过预处理和巧妙的倍增思想,大幅提升了求解效率[^1]。 ### 原理 - **倍增思想**:倍增法是一种每次将情况翻倍从而将线性处理转化为对数级处理,进而极大优化时间复杂度的方法。对于树上的节点,通过预处理出每个节点的第 $2^j$ 个祖先,在查找最近公共祖先时,利用这些预处理信息进行快速跳转,避免了朴素算法中一步一步向上查找的低效过程。 - **预处理**:构建一个二维数组 $f[i][j]$,其中 $f[i][j]$ 表示节点 $i$ 的第 $2^j$ 个祖先。这个数组可以通过深度优先搜索(DFS)进行预处理。 ### 代码示例 以下是一个使用C++实现的倍增LCA算法的基本模板: ```cpp #include<bits/stdc++.h> using namespace std; const int N = 5e4 + 10; int a[N];//原始输入数组 //st表,表示需要查询的数组的从下标i到下标i+2^j-1的最值 int mx[N][30];//最大值 int mn[N][30];//最小值 //预处理 void init(int n) { for(int i = 1; i <= n; i++){ mx[i][0] = a[i]; mn[i][0] = a[i]; } for(int j = 1; (1 << j) <= n; j++) { for(int i = 1; i + (1 << j) <= n + 1; i++){ mn[i][j] = min(mn[i][j - 1], mn[i + (1 << (j - 1))][j - 1]); mx[i][j] = max(mx[i][j - 1], mx[i + (1 << (j - 1))][j - 1]); } } } //查询函数 int search(int l, int r)//l和r范围为1~n { int k = log2(r - l + 1); int t1 = min(mn[l][k], mn[r - (1 << k) + 1][k]);//区间最小值 int t2 = max(mx[l][k], mx[r - (1 << k) + 1][k]);//区间最大值 return t2 - t1; } int main(){ // freopen("in.txt","r",stdin); int n, q; scanf("%d%d", &n, &q); for(int i = 1; i <= n; i++) scanf("%d", &a[i]); init(n); while(q--){ int l, r; scanf("%d%d", &l, &r); printf("%d\n", search(l, r)); } return 0; } ``` ### 复杂度分析 - **时间复杂度**:预处理的时间复杂度为 $O(n \log n)$,每次查询的时间复杂度为 $O(\log n)$,其中 $n$ 是树中节点的数量。 - **空间复杂度**:主要用于存储预处理的 $f[i][j]$ 数组,空间复杂度为 $O(n \log n)$。 ### 应用场景 在解决树上的路径查询、距离计算等问题时,经常需要用到LCA。例如,在处理树节点的第K个祖先、在传球游戏中最大化函数值、林卡树等问题中,都可以使用树上倍增LCA的方法[^2]。
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值