概述
倍增算法是一种优化算法,经常用于指数幂的问题的优化,快速幂就是它的直接体现。不仅仅是快速幂,在很多算法题中,都会需要倍增的思想,而且往往不是想快速幂那样直接,而是隐性的,需要挖掘题意。下面给出一道例题,体会倍增思想。
例题
问题描述
给定A,B,C,Q,进行以下操作 Q 次:
将A变为
⌊
A
B
⌋
\lfloor \frac{A}{B} \rfloor
⌊BA⌋ + C;
请你输出 Q 次操作后 A 的权值。
输入格式
第一行包含 4 个正整数 A,B,C,Q。
输出格式
输出共1行,包含1个整数,表示最终答案。
样例输入
5 2 1 1
样例输出
3
评测数据规模
对于所有测评数据,1 ≤ A, B, C < 10^5,1 ≤ Q ≤ 10^9
运行限制
c++ 时间:2s 空间:256mb
解题思路
由于Q的最大取值为 10^9 ,所以直接模拟肯定超时。观察题意我们发现,A无论变多少次,都不会超出[1,2*10^5]的这个区间,也就是它的值域是一定的且并不大,那我们可以考虑倍增思想,我们用f[a][b]来表示A执行 2^b 次后的值,通过f[a][b] = f[f[a][b-1][b-1]我们就可以求出值域A的任意次操作之后的值了。
点击查看代码#include<bits/stdc++.h>
using namespace std;
const int N=1e6,M=31;
int a,b,c,d,F[N][M];
int main(){
cin>>a>>b>>c>>d;
if(b==1){
cout<<a+1ll*c*d<<endl;
return 0;
}
for(int i=0;i<=200000;i++) F[i][0]=(i/b)+c;
for(int x=1;x<=30;x++)
for(int i=0;i<=200000;i++) F[i][x]=F[F[i][x-1]][x-1];
for(int i=0;i<=30;i++)
if((d&(1<<i) a=F[a][i];
cout<<a<<endl;
return 0;
}