求教 a的b次方(1<=A,B<=10000)

[xxx0624]

C语言
根据二分的思想求解，常规的算法会超时
例如：A*A*A==((A*A)%C)*A

[jcadre]

常规想法 for循环啊 不过应该会越界啦。

[xxx0624]

会啊。

[Gloveing]

大数运算

[pp25210]

#include<iostream>
　　using namespace std;
　　int f(int a,int b,int n);
　　int main()
　　{
　　cout<<f(19999,2999896,17)<<endl;
　　return 0;
　　}
　　int f(int a,int b,int n)
　　{
　　if(a<0||b<0||n<0)return -1;
　　int binary[32];
　　int index=0;
　　while(b)
　　{
　　if(b&1)
　　{
　　binary[index]=1;
　　}
　　else
　　{
　　binary[index]=0;
　　}
　　index++;
　　b>>=1;
　　}
　　index--;
　　int result=1;
　　while(index>=0)
　　{
　　result=(result*result)%n;
　　if(binary[index])
　　{
　　result=(result*a)%n;
　　}
　　index--;
　　}
　　return result;
　　}

这是我刚刚找到的一个算法，依据《算法导论》（第二版中文版，机械工业出版社出版）第539页的算法来实现，速度非常快。即使是对非常大的参数也是如此

[elated]

#b为正整数
def power(a, b):
	if b == 1:
		return a

	remainder = b % 2
	half = (b - remainder) / 2
	value = power(a, half)
	value *= value
	
	if remainder == 1:
		value *= a
	
	return value