Lib.rs

› 算法 | 科学 › 数学
#fraction #precision #arithmetic #version #big-int #python #big-integer

算法

任意精度算术

所有权归 Azat Ibrakov。

35个稳定版本 (14个主要版本)

14.1.0 2023年10月17日
14.0.0 2023年7月24日
13.2.0 2023年7月5日
12.0.0 2023年6月8日
0.0.0 2021年7月14日

#274 in 算法

Download history 1/week @ 2024-03-12 2/week @ 2024-03-19 2/week @ 2024-03-26 11/week @ 2024-04-02

每月607次下载
在 rene 中使用

MIT 许可证

575KB
16K SLoC

算法

在下文中 python 是 python3.8 或 pypy3.8 或任何更高版本 (python3.9，pypy3.9 等等) 的别名。

安装

安装最新的 pip 和 setuptools 包版本

python -m pip install --upgrade pip setuptools

用户

从 PyPI 仓库下载并安装最新稳定版本

python -m pip install --upgrade rithm

开发者

从 GitHub 仓库下载最新版本

git clone https://github.com/lycantropos/rithm.git
cd rithm

安装依赖

python -m pip install -r requirements.txt

安装

python setup.py install

使用

Python

任意精度整数

使用设置

>>> from rithm.integer import Int

我们可以

  • 构建
    >>> Int()
Int(0)
>>> Int(9)
Int(9)
>>> Int('9')
Int(9)
>>> Int('0b1001', 2)
Int(9)
>>> Int('0o11', 8)
Int(9)
>>> Int('0x9', 16)
Int(9)
>>> Int('1001', 2)
Int(9)
>>> Int('0o11', 8)
Int(9)
>>> Int('9', 16)
Int(9)
>>> Int(9.99)
Int(9)
  • 比较
    >>> Int(9) == Int(9)
True
>>> Int(9) >= Int(9)
True
>>> Int(9) > Int(8)
True
>>> Int(9) <= Int(9)
True
>>> Int(9) < Int(10)
True
  • 计算
    >>> abs(Int(-9))
Int(9)
>>> Int(4) + Int(5)
Int(9)
>>> Int(9) & Int(11)
Int(9)
>>> Int(19) // Int(2)
Int(9)
>>> ~Int(-10)
Int(9)
>>> Int(19) % Int(10)
Int(9)
>>> Int(3) * Int(3)
Int(9)
>>> -Int(-9)
Int(9)
>>> Int(1) | Int(8)
Int(9)
>>> Int(3) ** Int(2)
Int(9)
>>> Int(5) << Int(1)
Int(10)
>>> Int(5) >> Int(1)
Int(2)
>>> Int(25) - Int(16)
Int(9)
>>> Int(18) / Int(2)
Fraction(Int(9), Int(1))
>>> Int(2) ^ Int(11)
Int(9)

精确分数

使用设置

>>> from rithm.fraction import Fraction

我们可以

  • 构建
    >>> Fraction()
Fraction(Int(0), Int(1))
>>> Fraction(1)
Fraction(Int(1), Int(1))
>>> Fraction(1, 2)
Fraction(Int(1), Int(2))
>>> Fraction(50, 100)
Fraction(Int(1), Int(2))
>>> Fraction(0.5)
Fraction(Int(1), Int(2))
  • 比较
    >>> Fraction(1, 2) == Fraction(1, 2)
True
>>> Fraction(1, 2) >= Fraction(1, 2)
True
>>> Fraction(1, 2) > Fraction(1, 3)
True
>>> Fraction(1, 2) < Fraction(2, 3)
True
>>> Fraction(1, 2) != Fraction(1, 3)
True
  • 计算
    >>> abs(Fraction(-1, 2))
Fraction(Int(1), Int(2))
>>> Fraction(1, 3) + Fraction(1, 6)
Fraction(Int(1), Int(2))
>>> Fraction(3, 2) // Fraction(1)
Int(1)
>>> Fraction(3, 2) % Fraction(1)
Fraction(Int(1), Int(2))
>>> Fraction(1, 3) * Fraction(3, 2)
Fraction(Int(1), Int(2))
>>> -Fraction(-1, 2)
Fraction(Int(1), Int(2))
>>> Fraction(1, 2) ** 2
Fraction(Int(1), Int(4))
>>> Fraction(3, 2) - Fraction(1)
Fraction(Int(1), Int(2))
>>> Fraction(1, 3) / Fraction(2, 3)
Fraction(Int(1), Int(2))

Rust

任意精度整数

/// With setup
use std::convert::TryFrom;
use traiter::numbers::{
    Abs, DivEuclid, FromStrRadix, Pow, RemEuclid, Zero
};
use rithm::big_int;

#[cfg(target_arch = "x86")]
type Digit = u16;
#[cfg(not(target_arch = "x86"))]
type Digit = u32;
const DIGIT_BITNESS: usize = (Digit::BITS - 1) as usize;
const _: () = assert!(big_int::is_valid_digit_bitness::<Digit, DIGIT_BITNESS>());
type BigInt = big_int::BigInt<Digit, DIGIT_BITNESS>;
/// we can:
/// - construct
assert_eq!(BigInt::zero(), 0);
assert_eq!(BigInt::from(9), 9);
assert_eq!(BigInt::try_from("9").unwrap(), 9);
assert_eq!(BigInt::try_from("0b1001").unwrap(), 9);
assert_eq!(BigInt::try_from("0o11").unwrap(), 9);
assert_eq!(BigInt::try_from("0x9").unwrap(), 9);
assert_eq!(BigInt::from_str_radix("1001", 2).unwrap(), 9);
assert_eq!(BigInt::from_str_radix("11", 8).unwrap(), 9);
assert_eq!(BigInt::from_str_radix("9", 16).unwrap(), 9);
assert_eq!(BigInt::try_from(9.99).unwrap(), 9);
/// - compare
assert!(BigInt::from(9) == BigInt::from(9));
assert!(BigInt::from(9) >= BigInt::from(9));
assert!(BigInt::from(9) > BigInt::from(8));
assert!(BigInt::from(9) <= BigInt::from(9));
assert!(BigInt::from(9) < BigInt::from(10));
/// - calculate
assert_eq!(BigInt::from(-9).abs(), 9);
assert_eq!(BigInt::from(4) + BigInt::from(5), 9);
assert_eq!(BigInt::from(9) & BigInt::from(11), 9);
assert_eq!(BigInt::from(1) | BigInt::from(8), 9);
assert_eq!(BigInt::from(2) ^ BigInt::from(11), 9);
assert_eq!(BigInt::from(19) / BigInt::from(2), 9);
assert_eq!(BigInt::from(19).div_euclid(BigInt::from(2)), 9);
assert_eq!(BigInt::from(3) * BigInt::from(3), 9);
assert_eq!(-BigInt::from(-9), 9);
assert_eq!(!BigInt::from(-10), 9);
assert_eq!(BigInt::from(3).pow(BigInt::from(2)), 9);
assert_eq!(BigInt::from(19) % BigInt::from(10), 9);
assert_eq!(BigInt::from(19).rem_euclid(BigInt::from(10)), 9);
assert_eq!(BigInt::from(5) << 1, 10);
assert_eq!(BigInt::from(5) >> 1, 2);
assert_eq!(BigInt::from(25) - BigInt::from(16), 9);

精确分数

/// With setup
use std::convert::TryFrom;
use traiter::numbers::{Abs, DivEuclid, One, Pow, RemEuclid, Zero};
use rithm::fraction;

type Fraction = fraction::Fraction<i8>;
/// we can:
/// - construct
assert_eq!(Fraction::zero(), 0);
assert_eq!(Fraction::one(), 1);
assert_eq!(Fraction::new(1, 2), Some(Fraction::from(1) / 2));
assert_eq!(Fraction::new(50, 100), Fraction::new(1, 2));
assert_eq!(Fraction::try_from(0.5).unwrap(), Fraction::new(1, 2).unwrap());
/// - compare
assert!(Fraction::new(1, 2).unwrap() == Fraction::new(1, 2).unwrap());
assert!(Fraction::new(1, 2).unwrap() >= Fraction::new(1, 2).unwrap());
assert!(Fraction::new(1, 2).unwrap() > Fraction::new(1, 3).unwrap());
assert!(Fraction::new(1, 2).unwrap() <= Fraction::new(1, 2).unwrap());
assert!(Fraction::new(1, 2).unwrap() < Fraction::new(2, 3).unwrap());
assert!(Fraction::new(1, 2).unwrap() != Fraction::new(1, 3).unwrap());
/// - calculate
assert_eq!(Fraction::new(-1, 2).unwrap().abs(), Fraction::new(1, 2).unwrap());
assert_eq!(Fraction::new(1, 3).unwrap() + Fraction::new(1, 6).unwrap(),
           Fraction::new(1, 2).unwrap());
assert_eq!(Fraction::new(1, 3).unwrap() / Fraction::new(2, 3).unwrap(),
           Fraction::new(1, 2).unwrap());
assert_eq!(Fraction::new(3, 2).unwrap().div_euclid(Fraction::from(1)), 1);
assert_eq!(Fraction::new(1, 3).unwrap() * Fraction::new(3, 2).unwrap(),
           Fraction::new(1, 2).unwrap());
assert_eq!(-Fraction::new(-1, 2).unwrap(), Fraction::new(1, 2).unwrap());
assert_eq!(Fraction::new(1, 2).unwrap().pow(2), Fraction::new(1, 4).unwrap());
assert_eq!(Fraction::new(3, 2).unwrap() % Fraction::from(1),
           Fraction::new(1, 2).unwrap());
assert_eq!(Fraction::new(3, 2).unwrap().rem_euclid(Fraction::from(1)),
           Fraction::new(1, 2).unwrap());
assert_eq!(Fraction::new(3, 2).unwrap() - Fraction::from(1),
           Fraction::new(1, 2).unwrap());

开发

版本升级

准备

安装 bump2version。

预发布

根据 semver 规范 选择要升级的版本号类别。

测试版本升级

bump2version --dry-run --verbose $CATEGORY

其中 $CATEGORY 是目标版本号类别名称，可能的值是 patch/minor/major。

升级版本

bump2version --verbose $CATEGORY

这将版本设置为 major.minor.patch-alpha。

发布

测试版本升级

bump2version --dry-run --verbose release

升级版本

bump2version --verbose release

这将设置版本为 major.minor.patch。

运行测试

安装依赖

python -m pip install -r requirements-tests.txt

普通

pytest

在 Docker 容器内部

  • 使用 CPython
    docker-compose --file docker-compose.cpython.yml up
  • 使用 PyPy
    docker-compose --file docker-compose.pypy.yml up

Bash 脚本

  • 使用 CPython

    ./run-tests.sh

    或

    ./run-tests.sh cpython

  • 使用 PyPy

    ./run-tests.sh pypy

PowerShell 脚本

  • 使用 CPython
    .\run-tests.ps1
    或
    .\run-tests.ps1 cpython
  • 使用 PyPy
    .\run-tests.ps1 pypy

依赖项

~3.5–8.5MB
~71K SLoC