Make.Z
Integers.
This modules provides arbitrary-precision integers. Small integers internally use a regular OCaml int
. When numbers grow too large, we switch transparently to GMP numbers (mpn
numbers fully allocated on the OCaml heap).
This interface is rather similar to that of Int32
and Int64
, with some additional functions provided natively by GMP (GCD, square root, pop-count, etc.).
This file is part of the Zarith library http://forge.ocamlcore.org/projects/zarith . It is distributed under LGPL 2 licensing, with static linking exception. See the LICENSE file included in the distribution.
Copyright (c) 2010-2011 Antoine Miné, Abstraction project. Abstraction is part of the LIENS (Laboratoire d'Informatique de l'ENS), a joint laboratory by: CNRS (Centre national de la recherche scientifique, France), ENS (École normale supérieure, Paris, France), INRIA Rocquencourt (Institut national de recherche en informatique, France).
For an optimal experience with the ocaml
interactive toplevel, the magic commands are:
#load "zarith.cma";;
#install_printer Z.pp_print;;
Alternatively, using the new Zarith_top
toplevel module, simply:
#require "zarith.top";;
Raised by conversion functions when the value cannot be represented in the destination type.
val zero : t
The number 0.
val one : t
The number 1.
val minus_one : t
The number -1.
val of_int : int -> t
Converts from a base integer.
val of_int32 : int32 -> t
Converts from a 32-bit (signed) integer.
val of_int64 : int64 -> t
Converts from a 64-bit (signed) integer.
val of_string : string -> t
Converts a string to an integer. An optional -
prefix indicates a negative number, while a +
prefix is ignored. An optional prefix 0x
, 0o
, or 0b
(following the optional -
or +
prefix) indicates that the number is, represented, in hexadecimal, octal, or binary, respectively. Otherwise, base 10 is assumed. (Unlike C, a lone 0
prefix does not denote octal.) Raises an Invalid_argument
exception if the string is not a syntactically correct representation of an integer.
val of_substring : string -> pos:int -> len:int -> t
of_substring s ~pos ~len
is the same as of_string (String.sub s
pos len)
val of_string_base : int -> string -> t
Parses a number represented as a string in the specified base, with optional -
or +
prefix. The base must be between 2 and 16.
val of_substring_base : int -> string -> pos:int -> len:int -> t
of_substring_base base s ~pos ~len
is the same as of_string_base
base (String.sub s pos len)
Integer division. The result is truncated towards zero and obeys the rule of signs. Raises Division_by_zero
if the divisor (second argument) is 0.
Integer remainder. Can raise a Division_by_zero
. The result of rem a b
has the sign of a
, and its absolute value is strictly smaller than the absolute value of b
. The result satisfies the equality a = b * div a b + rem a b
.
Computes both the integer quotient and the remainder. div_rem a b
is equal to (div a b, rem a b)
. Raises Division_by_zero
if b = 0
.
Integer division with rounding towards +oo (ceiling). Can raise a Division_by_zero
.
Integer division with rounding towards -oo (floor). Can raise a Division_by_zero
.
Euclidean division and remainder. ediv_rem a b
returns a pair (q, r)
such that a = b * q + r
and 0 <= r < |b|
. Raises Division_by_zero
if b = 0
.
Euclidean division. ediv a b
is equal to fst (ediv_rem a b)
. The result satisfies 0 <= a - b * ediv a b < |b|
. Raises Division_by_zero
if b = 0
.
Euclidean remainder. erem a b
is equal to snd (ediv_rem a b)
. The result satisfies 0 <= erem a b < |b|
and a = b * ediv a b + erem a b
. Raises Division_by_zero
if b = 0
.
divexact a b
divides a
by b
, only producing correct result when the division is exact, i.e., when b
evenly divides a
. It should be faster than general division. Can raise a Division_by_zero
.
divisible a b
returns true
if a
is exactly divisible by b
. Unlike the other division functions, b = 0
is accepted (only 0 is considered divisible by 0).
congruent a b c
returns true
if a
is congruent to b
modulo c
. Unlike the other division functions, c = 0
is accepted (only equal numbers are considered equal congruent 0).
For all bit-level operations, negative numbers are considered in 2's complement representation, starting with a virtual infinite number of 1s.
Shifts to the left. Equivalent to a multiplication by a power of 2. The second argument must be nonnegative.
Shifts to the right. This is an arithmetic shift, equivalent to a division by a power of 2 with rounding towards -oo. The second argument must be nonnegative.
Shifts to the right, rounding towards 0. This is equivalent to a division by a power of 2, with truncation. The second argument must be nonnegative.
val numbits : t -> int
Returns the number of significant bits in the given number. If x
is zero, numbits x
returns 0. Otherwise, numbits x
returns a positive integer n
such that 2^{n-1} <= |x| < 2^n
. Note that numbits
is defined for negative arguments, and that numbits (-x) = numbits x
.
val trailing_zeros : t -> int
Returns the number of trailing 0 bits in the given number. If x
is zero, trailing_zeros x
returns max_int
. Otherwise, trailing_zeros x
returns a nonnegative integer n
which is the largest n
such that 2^n
divides x
evenly. Note that trailing_zeros
is defined for negative arguments, and that trailing_zeros (-x) = trailing_zeros x
.
val testbit : t -> int -> bool
testbit x n
return the value of bit number n
in x
: true
if the bit is 1, false
if the bit is 0. Bits are numbered from 0. Raise Invalid_argument
if n
is negative.
val popcount : t -> int
Counts the number of bits set. Raises Overflow
for negative arguments, as those have an infinite number of bits set.
Counts the number of different bits. Raises Overflow
if the arguments have different signs (in which case the distance is infinite).
Note that, when converting to an integer type that cannot represent the converted value, an Overflow
exception is raised.
val to_int : t -> int
Converts to a signed OCaml int
. Raises an Overflow
if the value does not fit in a signed OCaml int
.
val to_int32 : t -> int32
Converts to a signed 32-bit integer int32
. Raises an Overflow
if the value does not fit in a signed int32
.
val to_int64 : t -> int64
Converts to a signed 64-bit integer int64
. Raises an Overflow
if the value does not fit in a signed int64
.
val to_string : t -> string
Gives a human-readable, decimal string representation of the argument.
val format : string -> t -> string
Gives a string representation of the argument in the specified printf-like format. The general specification has the following form:
% [flags] [width] type
Where the type actually indicates the base:
i
, d
, u
: decimalb
: binaryo
: octalx
: lowercase hexadecimalX
: uppercase hexadecimalSupported flags are:
+
: prefix positive numbers with a +
sign-
: left-justify (default is right justification)0
: pad with zeroes (instead of spaces)#
: alternate formatting (actually, simply output a literal-like prefix: 0x
, 0b
, 0o
)Unlike the classic printf
, all numbers are signed (even hexadecimal ones), there is no precision field, and characters that are not part of the format are simply ignored (and not copied in the output).
val fits_int : t -> bool
Whether the argument fits in an OCaml signed int
.
val fits_int32 : t -> bool
Whether the argument fits in a signed int32
.
val fits_int64 : t -> bool
Whether the argument fits in a signed int64
.
val pp_print : Format.formatter -> t -> unit
Prints the argument on the specified formatter. Can be used as %a
format printer in Format.printf
and as argument to #install_printer
in the top-level.
Comparison. compare x y
returns 0 if x
equals y
, -1 if x
is smaller than y
, and 1 if x
is greater than y
.
Note that Pervasive.compare can be used to compare reliably two integers only on OCaml 3.12.1 and later versions.
val sign : t -> int
Returns -1, 0, or 1 when the argument is respectively negative, null, or positive.
val is_even : t -> bool
Returns true if the argument is even (divisible by 2), false if odd.
val is_odd : t -> bool
Returns true if the argument is odd, false if even.
pow base exp
raises base
to the exp
power. exp
must be nonnegative. Note that only exponents fitting in a machine integer are supported, as larger exponents would surely make the result's size overflow the address space.
Returns the square root. The result is truncated (rounded down to an integer). Raises an Invalid_argument
on negative arguments.
Returns the square root truncated, and the remainder. Raises an Invalid_argument
on negative arguments.
root x n
computes the n
-th root of x
. n
must be positive and, if n
is even, then x
must be nonnegative. Otherwise, an Invalid_argument
is raised.
rootrem x n
computes the n
-th root of x
and the remainder x-root**n
. n
must be positive and, if n
is even, then x
must be nonnegative. Otherwise, an Invalid_argument
is raised.
val perfect_power : t -> bool
True if the argument has the form a^b
, with b>1
val perfect_square : t -> bool
True if the argument has the form a^2
.
val log2 : t -> int
Returns the base-2 logarithm of its argument, rounded down to an integer. If x
is positive, log2 x
returns the largest n
such that 2^n <= x
. If x
is negative or zero, log2 x
raise the Invalid_argument
exception.
val log2up : t -> int
Returns the base-2 logarithm of its argument, rounded up to an integer. If x
is positive, log2up x
returns the smallest n
such that x <= 2^n
. If x
is negative or zero, log2up x
raise the Invalid_argument
exception.
val size : t -> int
Returns the number of machine words used to represent the number.
extract a off len
returns a nonnegative number corresponding to bits off
to off
+len
-1 of b
. Negative a
are considered in infinite-length 2's complement representation.
signed_extract a off len
extracts bits off
to off
+len
-1 of b
, as extract
does, then sign-extends bit len-1
of the result (that is, bit off + len - 1
of a
). The result is between - 2{^[len]-1}
(included) and 2{^[len]-1}
(excluded), and equal to extract a off len
modulo 2{^len}
.
val to_bits : t -> string
Returns a binary representation of the argument. The string result should be interpreted as a sequence of bytes, corresponding to the binary representation of the absolute value of the argument in little endian ordering. The sign is not stored in the string.
val of_bits : string -> t
Constructs a number from a binary string representation. The string is interpreted as a sequence of bytes in little endian order, and the result is always positive. We have the identity: of_bits (to_bits x) = abs x
. However, we can have to_bits (of_bits s) <> s
due to the presence of trailing zeros in s.