数学-快速幂与矩阵快速幂

快速幂

二进制取幂的想法是，我们将取幂的任务按照指数的 二进制表示 来分割成更小的任务

快速幂的递归版本

long long fastpower(long long a, long long b) {
  if (!b) return 1;
  long long res = fastpower(a, b / 2);//递归调用；
  if (b % 2)
    return res * res * a;
  else
    return res * res;
}

while循环版本 相对于原来的递归版本更快（解决你的TLE）

long long fastpower(long long base,long long power)
{
	long long result =1;
	while(power>0)
	{
		if(power&1)//判断指数是否是奇数
			result=result*base;//将指数分出1
			power>>=1;//指数>>1
			base=(base*base);//基底平方
	}
	return result;
}

中间取模

long long fastpower(long long base,long long power)
{
	long long result =1;
	while(power>0)
	{
		if(power&1)//判断指数是否是奇数
			result=（result*base）%mod;//将指数分出1
			power>>=1;//指数>>1
			base=(base*base)%mod;//基底平方
	}
	return result;
}

练习题：HDU2035传送门
AC代码

#include <cstdio>
#include <iostream> 
const int mod=1000; 
using namespace std; 
long long qp(long long a,int b) {
    long long ans=1;
    while(b){
        if(b&1) ans=(ans*a)%mod;
        a=(a*a)%mod;
        b>>=1;
    }
    return ans; } int main() {

    int a,b;
    while(cin>>a>>b,a&&b){
        cout<<qp(a,b)<<endl;
    }
    return 0; }

练习题：
UVA374
UVA1230
UVA11029 //试试科学计数法

UVA11029AC
方法分析

#include <cstdio>
#include <iostream>
#include<cmath>
 using namespace std; 
 int mod; long long qp(long long a,int b) {
    long long ans=1;
    while(b){
        if(b&1) ans=(ans*a)%mod;
        a=(a*a)%mod;
        b>>=1;
    }
    return ans; } int main() {

    int n;
    cin>>n;					
    int a,b,m;
    mod=1000;
    while(n--){
        cin>>a>>b;
        printf("%d...%03d\n",(int)pow(10,2+fmod(b*log10(a),1)),qp(a,b));
    }
    return 0;
} 

总结1：取模运算
(a+b)%n=((a mod n)+(b mod n)) mod n;
(a-b)%n=((a mod n)-(b mod n)) mod n;
(ab)%n=((a mod n)(b mod n)) mod n;
取对数的利用

矩阵快速幂

矩阵是一个按照长方阵列排列的复数或实数集合

矩阵的基本运算

其中矩阵乘法应该遵循以下的条件

前一个矩阵的列数等于后一个矩阵的行数

计算规则是，第一个矩阵第一行的每个数字（1,2,3），各自乘以第二个矩阵第一列对应位置的数字（7,9,10），然后将乘积相加（1x7+2x9+3x10），得到结果矩阵左上角的那个值58

快速明白矩阵乘法

利用结构体和操作符重载实现矩阵快速幂

洛谷P3390

#include <bits/stdc++.h>
using namespace std;
const int MOD=1e9+7;
struct matrix{
	long long C[101][101];//申请二维数组 模拟矩阵
}A,I;
long long  n,k;
matrix operator *(const matrix &x,const matrix &y)
{
	matrix a;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)//矩阵初始化
			a.C[i][j]=0;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)//三重循环模拟矩阵乘法
			for(int k=1;k<=n;k++)
			{
				a.C[i][j]+=x.C[i][k]*y.C[k][j]%MOD;//矩阵乘法公式 
				a.C[i][j]%=MOD;	
			}
	return a;
}
void fp(matrix a,long long b)//fastpower部分
{
	while(b>0)//非递归版本快速幂
	{
		if(b&1)
		I=I*a;
		b>>=1;
		a=a*a;
	}
}
int main()
{
	cin>>n>>k;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			cin>>A.C[i][j];//读入矩阵
	for(int i=1;i<=n;i++)
		I.C[i][i]=1; //初始化答案矩阵
	fp(A,k);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++)
			cout<<I.C[i][j]<<" ";//输出答案
			cout<<endl;
	}
	return 0;
}