求a的j次幂模p

2019-04-14 16:07发布

// exponent.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"
#include
using namespace std;
int ModularExponent(int a,int j,int p);

int _tmain(int argc, _TCHAR* argv[])
{
    //int k=ModularExponent(2,3,3);
    //cout<     int a,j,p;
    do
    {
        cout<<"Input a,j,p: ";   
        cin>>a;
        if(a<0)
            break;
        cin>>j>>p;
        cout<     } while (a>0);
    return 0;
}
int ModularExponent(int a, int j, int p)
{
    //求方幂模s=aj mod p, 注意先求aj可能会溢出
    int s=1;
    while (j>0)
    {
        if(j%2)
            s=(s*a%p);//if (j is odd) s ← s·a mod p;
        a=(a*a)%p;//a ← a2 mod p;
        j=j/2;//j ← j div 2;
    }
    return s;
}