信息学奥赛一本通 提高篇 序列第k个数 及 快速幂

Posted 2021-11-30 lztzs

我是传送门

这个题首先是先判断是等差还是等比数列

等差的话非常简单:

前后两个数是等差的,举个栗子:

3 6 9 12

这几个数,(我感觉 1 2 3 4并说明不了什么)

每次都加3嘛,很容易看出,第一个数是3 * 1,第二个是3 * 2....以此类推

第k个数 = (第2个数 - 第1个数) * k ;

(z - y) * k % 200907

% 200907 的原因是题目要求

 

但是这样并不能过

hack一下

4 7 10 13

用原先的公式:
(7 - 4) * 4 % 200907 = 12;

明显不对啊

但是我们将序列每个都减一

就变成了: 3 6 9 12

熟悉滴感觉,但我们并不需要真的都减一,只需要加上差就可以了

x + (z - y) * (k - 1) % 200907

这个和我说的有些出入,但仔细想想就会发现是一样的

 

等比的话稍微难一点:

前后两个数是等比的,举个栗子:

3 6 18 27

第一个数是3^1 

第二个数是3^2

第三个数是3^3

emmm,好像也很简单

这数据 不能乘10^9次吧

100000000000000% 超时

 

快速幂不错

快速幂怎么写呢

我在这里说一种非递归写法(其实本蒟蒻根本不明白递归的实现)

把幂次拆成二进制数,只需要乘log n 次,非常优秀的算法

比如 2 ^ 10 的二进制是 2的1010次幂

ans 就是  2 ^ 2 * 2 ^ 8 ;

当然数学上来讲  2 ^ 2 * 2 ^ 8 = 2 ^ (2 + 10)

long long ksm(long long a,long long b,long long n){//a是底数,b是幂数,n是mod数
    long long ans = 1;
    while(b){//等价于while(b != 0) ,dalao勿喷,讲给我这样的蒟蒻听
        if ( b & 1){//b & 1 等价于 b % 2 == 1
            ans *= a;//是的话乘起来
            ans %= 200907;//随时取膜性质
            /* code */
        }
        a *= a;//下一位
        a %= 200907;//膜
        b >>= 1;//右移一位  等价于 b /= 2;
    }
    return ans % 200907;
}