C ++:实现Modular Exponentiation

时间:2017-03-08 04:57:03

标签: c++

我正在使用这个新的和改进的代码我更正了以解决这个问题。

我正在使用模块化的Exponentiation来使用公式[a ^ k mod n]来获得我必须做的任务的答案,我需要在两个步骤中对其进行编码。

首先必须将 int k转换为二进制 表示K由0和1的列表组成。其次,必须执行Modular Exponentiation 使用anK[]作为参数..

早些时候我的代码不正确并且能够纠正它。

我现在面临的问题是,当我使用5^3 % 13的模块化指数的在线计算器时,它应== 8

我从代码中得到的结果是5。 我试图了解代码中是否有一些小问题,或者我的数学错了?谢谢

#include <iostream>
#include <vector>
using namespace std;

vector <int> BinaryK(int k);
int ModularExpo(int a, vector <int> & k, int n);

int main()
{
    int a = 0;
    int k = 0;
    int n = 0;

    cout << "a^k % n" << endl;
    cout << "a = ";
    cin >> a;

    cout << "k = ";
    cin >> k;

    cout << "n = ";
    cin >> n;

    vector<int> B = BinaryK(k);
    int result = ModularExpo(a, B, n);
    cout << "a ^ k mod n == " << result << endl;

    return 0;
}

// c == b^e % m

vector<int> BinaryK(int k)
{
    vector<int> K; //hint: make K a vector
    int tmp = k;
    while (tmp > 0)
    {
        K.push_back(tmp % 2); //hint: use pushback
        tmp = tmp / 2;
    }

    return K;
}

int ModularExpo(int a, vector<int> & K, int n)
{
    if (n == 1)
        return 0;

    int b = 1;
    if (K.size() == 0)
        return b;
    int A = a;
    if (K[0] == 1)
        b = a;

    for (int i = 1; i < K.size() - 1; i++)
    {
        A = A * A % n;
        if (K[i] == 1)
            b = A*b % n;
    }
    return (b);
}

1 个答案:

答案 0 :(得分:1)

改变这一行:

    for (int i = 1; i < K.size(); i++)   // K.size() not K.size()-1