2020.1.15每日一题“快速幂求模”

@lebrongod

初来乍到

这是我的第一篇博客，如果写的有什么不对，欢迎提出意见，我会在日后的博客中改正。

快速幂

快速幂顾名思义就是快速地求一个幂指数计算式的值，比如：a^n对于这个式子我们第一时间想到的应该是循环n次累乘求解，显而易见他的时间复杂度为o(n)，显然这样的时间复杂度不是最优的。

这个时候我们就应该想办法如何降低他的循环次数，我们发现a^n =a^(n/2)*a(n/2)（前提是n为偶数，否则n/2会略去小数点导致等式不成立），而a^(n/2) 又可以变成a^ ((n/2)/2)*a^((n/2)/2),不断地像这样二分至n==1（即为递推的终止条件）。

不过在每一次分之前需要注意的是要先判断n的奇偶性，n为偶时，可以将其二分，累乘，并进入下一次二分，但当n为奇时，将（n-1）放入下一次二分，将之前累乘的结果乘入记录答案的变量（ans）。

核心代码如下：

快速幂取模

当我们在处理一个很大的由幂得到的数的模时，速度是非常慢的。

(ab)%c=(a%cb%c)%c

当数个a累乘起来的时候有

(a^n )%c=((a%c)^n)%c

这样一看我们要求最终的模就可以对每一次的a都做取余的步骤，
这样就避免了将一个很大的数进行求模，
而是将每次底数的规模缩小了，
这大大地提高了代码的效率。

核心代码如下：

以下为快速幂求模完整代码：

#include<stdio.h>
int kuaisumi(int n,int c)
{
    long long int a=2;
    long long int ans=1;
    a=a%c;
    while(n>0)
    {
        if(n%2==1)
            ans=(ans*a)%c;
        n=n/2;
        a=(a*a)%c;
    }
    return ans;
}
int  main()
{
    long long int s,i,f,n;
    while(~scanf("%lld",&n))
    {
        for(i=2,s=1; i<=n; i++)
        {
            f=kuaisumi(i-1,107);
            s+=f;
        }
        printf("%lld\n",s);
    }
}