multiprec-test/main.cpp at master · jamse101/multiprec-test · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
#define NDEBUG

#include <iostream>
#include <chrono>
#include <cmath>

#include <boost/multiprecision/cpp_int.hpp>
#include <boost/multiprecision/cpp_bin_float.hpp>
#include <boost/multiprecision/float128.hpp>
#include <boost/multiprecision/mpfr.hpp>

namespace mp = boost::multiprecision;

void taxi_cab(mp::uint128_t tc)
{
    mp::uint128_t raja = (mp::cbrt(tc.convert_to<mp::float128>()) + 0.5).convert_to<mp::uint128_t>();

    std::cout << raja << "\n";

    mp::uint128_t temp = 0;

    mp::uint128_t i = 0;
    mp::uint128_t j = 0;

    for ( i = 1 ; i < raja ; i++ )
    {
        mp::uint128_t iii = i*i*i;

        //mp::float128 f = (tc-iii).convert_to<mp::float128>();

        //j = (mp::round(mp::cbrt(f))).convert_to<mp::uint128_t>();

        double f = (tc-iii).convert_to<double>();

        j = mp::uint128_t(std::round(std::cbrt(f)));

        if (iii + j*j*j == tc)
        {
            if (i == temp)
                break;

            std::cout << i << " " << j << "\n";

            temp = j;
        }
    }
}

int main()
{
    auto t0 = std::chrono::steady_clock::now();

    taxi_cab(mp::uint128_t("24153319581254312065344"));

    auto t1 = std::chrono::steady_clock::now();

    std::chrono::duration<double> dur = t1 - t0;

    std::cout << dur.count() << "\n";

    return 0;
}