标签：龟速 光速 int res sqrt base long 数学知识 mod

1. 快速幂

考虑求 $a^b \operatorname{mod} p$ ，$p$ 是质数

用乘法累乘实在是太慢了，所以我们要找出更优秀的算法

不妨将 $b$ 分解为二进制，比如 $(11)_{10}$ 分解成 $(1011)_2$

那么 $11=8+2+1$ ，也就是 $a^{11}=a^{8+2+1}=a^8a^2a^1$

又发现 $a^8=(a^4)^2$ ，$a^4=(a^2)^2 \cdots$

那么自然我们可以用 $a\times a$ 算出 $a^2$ ，再让 $a^2 \times a^2$ 得到 $a^4 \cdots$

假设底数是 $a$ ，指数是 $b$ ，模数是 $mod$

那么有以下代码

　

#define int long long
int quick_pow(int a,int b,int mod){
    int base=a,res=1;
    while(b){
        if(b&1) res=res*base%mod;
        base=base*base%mod;
        b>>=1;
    }
    return res;
}

 

（感觉好没有意义

2. 龟速乘

比起计算机自带的乘法，龟速乘的的运行速度还要慢上一些。

但是，它可以有效地保证你的 $\mathbf{long}$ $\mathbf{long}$ 不会炸掉并送给你一个神奇的数字。

#define int long long
int slow_time(int a,int b,int mod){
    int base=a,res=1;
    while(b){
        if(b&1) res=(res+base)%mod;
        base=(base+base)%mod;
        b>>=1;
    }
    return res;
}

我们可以发现它的原理和快速幂是基本相同的，只不过是把乘换成了加，仅此而已

原理是显然的，在每次自加的时候可以进行一个取模，这样就可以保证不会爆掉

 

3. 光速幂

还是设底数是 $a$ ，指数是 $b$ ，模数是 $mod$

主要思想就是分块，预处理出 $a^1 , a^2 , a^3 , \cdots ,a^{\sqrt{n}}$ 和 $a^{\sqrt{n}} , a^{2\sqrt{n}} , a^{3\sqrt{n}} , \cdots ,a^n$ （这里的 $n$ 比最大的 $b$ 要大）

然后在询问的时候取模直接查询即可

用的范围比较窄 $\cdots$

$\Theta (\sqrt{n})$ 预处理， $\Theta(1)$ 查询

struct Lightspeed_Pow{
    int base1[WR],basesqrt[WR];
    int BL;//block_len不是别的什么东西
    int maxn=1e12;//最大的可能模数，建议小于等于1e12否则可能炸空间
    void init(int x){
        BL=sqrt(maxn)+1;
        base1[0]=1;
        for(int i=1;i<=BL;i++){
            base1[i]=base1[i-1]*x%mod;
        }//处理常数次幂
        basesqrt[0]=1;
        for(int i=1;i<=BL;i++){
            basesqrt[i]=basesqrt[i-1]*base1[BL]%mod;
        }//处理根号n的倍数次幂
    }
    int calc(int x){
        return basesqrt[x/BL]*base1[x%BL]%mod;
    }
}pw;

 

标签：龟速,光速,int,res,sqrt,base,long,数学知识,mod
来源： https://www.cnblogs.com/WintersRain/p/16505292.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。