倍增法和ST算法 个人学习笔记&代码

倍增法

P4155 国旗计划

cpp 复制代码
#include<bits/stdc++.h>
#define endl '\n'
#define forr(i,l,r) for(int i=l;i<=r;i++)
#define reforr(i,l,r) for(int i=r;i>=l;i--)
#define int long long
using namespace std;
const int N=4e5+10,M=20;
struct wa
{
    int id,l,r;
}a[N*2];
int n,m,nt;
int go[N][M],res[N];
void init(){
    int nxt=1;
    forr(i,1,nt){
        while (nxt<=nt&&a[nxt].l<=a[i].r)//找下一个最优区间
        {
            nxt++;
        }
        go[i][0]=nxt-1;//2^0=1 记录跳一步的最优位置
    }
    //递推
    for(int i=1;(1<<i)<=n;i++){//倍增
        forr(s,1,nt){
            go[s][i]=go[go[s][i-1]][i-1];//跳到2^(i-1)的跳板 再跳2^(i-1)
        }
    }
}
void fans(int x){
    int len=a[x].l+m;//绕一圈
    int ans=0,cur=x;
    reforr(i,0,(int)log2(N)){//从最大数开始跳
        int pos=go[cur][i];
        if(pos&&a[pos].l<len){
            ans+=(1<<i);
            cur=pos;//跳到新位置
        }
    }
    res[a[x].id]=ans+1;//排过序 不是以前的顺序了 输出答案要还原顺序
}
void solve()
{
	cin>>n>>m;
    forr(i,1,n){
        a[i].id=i;
        cin>>a[i].l>>a[i].r;
        if(a[i].r<a[i].l)a[i].r+=m;
    }
    sort(a+1,a+n+1,[](wa x,wa y){
        return x.l<y.l;
    });
    nt=n;
    forr(i,1,n){
        nt++,a[nt].l=a[i].l+m,a[nt].r=a[i].r+m;//拆环成链
    }
    init();
    forr(i,1,n) fans(i);
    forr(i,1,n)cout<<res[i]<<' ';

}
signed main()
{
	ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0);

	int t=1;
    // cin>>t;
	while(t--) solve();
    
	return 0;
}

ST算法

用于求解静态数组RMQ(区间最值查询)问题

P2880 Balanced Lineup G

cpp 复制代码
const int N=5e4+10,M=20;
int n,m;
int a[N],l2[N];
int dpmn[N][M],dpmx[N][M];
void init(){
    l2[0]=-1;//log2递推计算
    forr(i,1,N)l2[i]=l2[i>>1]+1;
    forr(i,1,n){//初始化
        dpmx[i][0]=a[i];
        dpmn[i][0]=a[i];
    }
    int p=log2(n);
    //递推
    forr(i,1,p){//区间长度 小区间凑出大区间
        for(int s=1;s+(1<<i)<=n+1;s++){//左端点位置
            dpmx[s][i]=max(dpmx[s][i-1],dpmx[s+(1<<(i-1))][i-1]);
            dpmn[s][i]=min(dpmn[s][i-1],dpmn[s+(1<<(i-1))][i-1]);
        }
    }
}
int ans(int l,int r){
    int k=l2[r-l+1];//2^k<=len 2^(k+1)>=len 
    int mx=max(dpmx[l][k],dpmx[r-(1<<k)+1][k]);
    int mn=min(dpmn[l][k],dpmn[r-(1<<k)+1][k]);
    return mx-mn;
}
void solve()
{
	cin>>n>>m;
    forr(i,1,n){
        cin>>a[i];
    }
    init();
    forr(i,1,m){
        int l,r;
        cin>>l>>r;
        cout<<ans(l,r)<<endl;
    }
}
相关推荐
秋说1 小时前
【PTA数据结构 | C语言版】一元多项式求导
c语言·数据结构·算法
Maybyy1 小时前
力扣61.旋转链表
算法·leetcode·链表
兴趣使然_3 小时前
【笔记】使用 html 创建网址快捷方式
笔记·html·js
卡卡卡卡罗特3 小时前
每日mysql
数据结构·算法
chao_7894 小时前
二分查找篇——搜索旋转排序数组【LeetCode】一次二分查找
数据结构·python·算法·leetcode·二分查找
aramae4 小时前
C++ -- STL -- vector
开发语言·c++·笔记·后端·visual studio
lifallen4 小时前
Paimon 原子提交实现
java·大数据·数据结构·数据库·后端·算法
lixzest4 小时前
C++ Lambda 表达式详解
服务器·开发语言·c++·算法
EndingCoder5 小时前
搜索算法在前端的实践
前端·算法·性能优化·状态模式·搜索算法
丶小鱼丶5 小时前
链表算法之【合并两个有序链表】
java·算法·链表