calx

コード例

Int — 多倍長整数

#include <math/core/mp/Int.hpp>
#include <iostream>
using namespace calx;

int main() {
    Int a("123456789012345678901234567890");
    Int b("987654321098765432109876543210");

    std::cout << "a + b = " << a + b << std::endl;
    std::cout << "a * b = " << a * b << std::endl;
    std::cout << "gcd   = " << gcd(a, b) << std::endl;

    // 素数判定 (Miller-Rabin)
    Int p("170141183460469231731687303715884105727");  // 2^127 - 1
    std::cout << p << " is prime: " << std::boolalpha
              << isProbablePrime(p) << std::endl;

    // NaN/Infinity 安全伝播
    Int zero;
    Int nan = Int(1) / zero;
    std::cout << "1/0 = " << nan << std::endl;  // NaN
}

Float — 多倍長浮動小数点

#include <math/core/mp/Float.hpp>
#include <iostream>
using namespace calx;

int main() {
    // 数学定数 (スレッドローカルキャッシュ、2 回目以降は即座に返る)
    Float pi = Float::pi(1000);    // π を 1000 桁で計算 (Chudnovsky)
    Float e  = Float::e(1000);     // e を 1000 桁で計算

    // 超越関数 (精度は引数から自動決定)
    Float y = exp(pi);             // e^π
    Float z = log(y);              // log(e^π) ≈ π

    std::cout << "exp(pi) = " << y << std::endl;
    std::cout << "log(exp(pi)) = " << z << std::endl;

    // 三角関数
    Float s = sin(pi / 6);         // sin(π/6) = 0.5
    std::cout << "sin(pi/6) = " << s << std::endl;
    std::cout << "pi(1000 digits) = " << pi << std::endl;
}