Большая арифметика от MinGW провалилась

Question

Моя общая задача - вычислить По словам Вольфрама,

Я ожидал, что следующее вернет то же самое:

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

unsigned long long f(const std::vector<unsigned> &coefficients)
{
    unsigned short exponent = coefficients.size() - 1;
    unsigned long long sum;
    for (unsigned i : coefficients)
    {
        sum += i * pow(36, exponent);
        --exponent;
    }
    return sum;
}

int main()
{
    std::cout << f({9,13,19,6,7,8,2});
}

но вместо этого он возвращает 20416905041, По словам Алекса Б, минимальная емкость длинного целого числа без знака - от 0 до 18446744073709551615, поэтому проблема не в емкости.

---

Характеристики:

• Компилятор: x86_64-w64-mingw32 из набора компиляторов Windows TDM-GCC
• Составлено через g++ mwe.cpp -std=c++11 -omwe
• ОС: Windows 10 Home

Answer 1

Из комментария М.М. к вопросу:

pow является функцией с плавающей запятой и поэтому может быть подвержен ошибкам округления. Вместо этого вы должны использовать целочисленные переменные (сделать повторное умножение на 36).