整形数快速幂

整形数快速幂：

计算整形数的幂的普通算法时间复杂度为O（N）,快速幂的时间复杂度为O(log N)。

基本思想：

以x^10 为例：可以看做(x^5)^2,还可以继续拆分为:(x*(x^2)^2)^5;

1. 令m=n/2。 假设已经知道如何计算x^m，再来求xⁿ。

2. 分两种情况：

如果n为偶数，xⁿ=(x^m)²

如果n为奇数，xⁿ=x(x^m)²

根据这个思想，可以得出递归的算法：

伪代码：

我们推导出迭代的算法。

推导过程：

求实数x的n次方法：

1. 令n表示为二进制d_k，d_k-1，...d₀。n=2^k+2^k-1+...+2⁰

2. 令y=1.

3. 从左到右扫描二进制数字，如果当前数字为0，y赋值y²，如果当前数字为1，y赋值xy²。

伪代码：

具体实现代码：

#include<iostream>
using namespace std;
int Pow(int a,int m){
	int temp;
	temp = 1;
	while(m){
		if(m&1){
			temp*=a;
		}
		a*=a;
		m/=2;
	}
	return temp;
} 
int main(){
	int n,m;
	cin>>n>>m;
	cout<<Pow(n,m);
	return 0;
} 

幂数较大，最终结果超出整数范围：

#include<iostream>
using namespace std;
int a[10000];
int main(){
	int n,m;
	scanf("%d%d",&n,&m);
	a[0] = 1;
	int len=1,jin=0;
	for(int i = 1;i <= m;i++){
		for(int j = 0;j < len;j++){
			int temp = a[j]*n+jin;
			a[j] = temp%10;
			jin = temp/10;
		}
		while(jin){
			a[len++] = jin%10;
			jin /= 10;
		}
	}
	for(int i = len-1;i>=0;i--){
		printf("%d",a[i]);
	}
	return 0;
}