随机算法 之模幂函数

在随机素数测试算法中要用到模幂运算，在O(lgn)的时间内产生模幂结果是非常有用的。
在诸如RSA等算法中都要用到求a^n mod p的运算，例如费马小定理（a^(n-1) mod n = 1,p是a的非素数因子)
及rsa算法用到的费马定理的推广(a^(y(n))mod n = 1，y(n)为n的欧拉函数)等等都需要用到模幂运算，那么
怎么能快速的到模幂运算结果其实原理很简单，这是我用english写的，希望没写错，呵呵
Basing on the two simple theories:
 1,(x*y)%z=((x%z)*(y%z))%z
 2,(x^y)%z=((x%z)^y)%z
To solve the y=(g^x)%p problem left shifting x until x become 0,compute g^x by computing
k=(g^(x/2)%p) and g^x=k*k or g^x=k*k*g 。so this algorithm is O(lgn) which is also linear。
This is a simple and recursive algorithm.
这是我写的参考程序 （仅供参考）
class MMath{
 
 public:
  
 static int power(int g,int p,int n){

  int s,u;

  if(n==0)return 1;
  else {
   s = power(g,p,n>>1);
   u = s*s;
   if(n%2==0)
    return u%p;
   else
    return u*g%p;
  }
 }

 static bool prime_judge(int n){
  if(power(2,n,n-1)!=1)return false;
  else return true;
 }
}
另外附一个根据费马定理写的素数测试程序（会产生错误，伪素数如1387等）
MillerRabin测试还没有想好怎么写程序呢。