Упрощение модульного возведения в степень C++

Я пытаюсь написать функцию расшифровки для системы шифрования RSA, все работает нормально для очень небольших чисел, однако иногда вывод просто неверен (я думаю, что причиной может быть ошибка с плавающей запятой или какое-то переполнение стека).

Процесс, который вызывает у меня проблемы, можно упростить до (11^23) mod 187, но я включу полный код на случай, если кто-то захочет его увидеть. Я знаю, что ответ должен быть 88, так как это пример, использованный в Приложении J "Книги кодов" доктора Саймона Сингха (я также проверил с помощью Wolfram Alpha). Тем не менее, я получаю результат 149. Однако с меньшими числами это согласуется с Wolfram Alpha.

Я думаю, что мне нужно упростить модульное возведение в степень, используя знание того, что:

a^b = a^c * a^d [где c + d = b]

Тем не менее, я до сих пор не уверен на 100%, если это проблема, это мое первое переполнение стека? (Я все еще не уверен на 100%, что это значит). Прежде чем кто-нибудь попытается на меня ответить, нет, это не какое-то домашнее задание, и мне жаль, если этот вопрос кажется тривиальным. Я открыт для использования gmp.h, если все думают, что это будет слишком сложно сделать, но я бы предпочел не делать этого, если буду полностью честен. Мой код приведен ниже (первая половина предназначена для вычисления закрытого ключа, который, как я считаю, не имеет отношения к моей проблеме, но я включил его на случай, если я ошибаюсь), я действительно надеюсь, что вы, ребята, можете помочь, спасибо вы очень заранее.

#include <iostream>
#include <math.h>

using namespace std;

unsigned int modinv(unsigned int u, unsigned int v)
{
    unsigned int inv, u1, u3, v1, v3, t1, t3, q;
    int iter;

    u1 = 1;
    u3 = u;
    v1 = 0;
    v3 = v;

    iter = 1;

    while (v3 != 0)
    {

        q = u3 / v3;
        t3 = u3 % v3;
        t1 = u1 + q * v1;

        u1 = v1; v1 = t1; u3 = v3; v3 = t3;
        iter = -iter;
    }

    if (u3 != 1)
        return 0;
    if (iter < 0)
        inv = v - u1;
    else
        inv = u1;
    return inv;
}

int main()
{ long unsigned int p = 17;
long unsigned int q = 11;
long unsigned int phi = (p-1)*(q-1);
long unsigned int e = 7;
long unsigned int c = 11;
long unsigned int n = p*q;
long unsigned int d = modinv (e,phi);
    {
         cout << fmod (pow (c, d), n);
    }
    return 0;
}