快速幂与64位整数乘法

众所周知，求幂函数我们常用math.h里的pow函数。但是，大家有没有遇见这样的情况，我们用pow求5的3次方，得到的结果却是124。这是为什么呢？
其实，是因为pow返回值是double类型，我们转换为int类型时，存在精度损失，然后小数点后再被截断，结果就少了1。

那么，我们该怎么处理这个问题呢？
运用，int类型的四舍五入，就可以解决这个问题了。
int a=(int)(pow()+0.5);

但是，由于c语言数据类型最高只有64位(long long)，题目一般都会要求我们对结果取模，但如果是一个64位的数据的64位的数据次幂，那么pow是肯定用不了的。

此时，就需要用我接下来讲的 快速幂 了。

快速幂

快速幂，如其名，它很快(时间复杂度为O(logn))。 所谓的快速幂，实际上是快速幂取模的缩写，简单的说，就是快速的求一个幂式的模(余)。在程序设计过程中，经常要去求一些大数对于某个数的余数，为了得到更快、计算范围更大的算法，产生了快速幂取模算法。
(以下算法long long简写为ll)
算法1.首先直接地来设计这个算法：

ll ans = 1;
for(ll i = 1;i<=n;i++)
{
   
    ans = ans * a;
}
ans = ans % c;

这个算法的时间复杂度体现在for循环中，为O（n）.这个算法存在着明显的问题，如果a和n过大，很容易就会溢出。
那么，我们先来看看第一个改进方案：在讲这个方案之前，要先有这样一个公式：
.这个公式大家在离散数学或者数论当中应该学过，不过这里为了方便大家的阅读，还是给出证明：
引理1：

上面公式为下面公式的引理，即积的取余等于取余的积的取余。

证明了以上的公式以后，我们可以先让a关于c取余，这样可以大大减少a的大小，
于是不用思考的进行了改进：
算法2：

ll ans = 1;
a = a % c; //加上这一句
for(ll i = 1;i<=b;i++)
{
   
    ans = ans * a;
}
ans = ans % c;

聪明的读者应该可以想到，既然某个因子取余之后相乘再取余保持余数不变，那么新算得的ans也可以进行取余，所以得到比较良好的改进版本。

算法3：

ll ans = 1;
a = a % c; //加上这一句
for(ll i = 1;i<=b;i++)
{
   
    ans =