1995-11-06 02:34:19 -08:00
|
|
|
(***********************************************************************)
|
|
|
|
(* *)
|
1996-04-30 07:53:58 -07:00
|
|
|
(* Objective Caml *)
|
1995-11-06 02:34:19 -08:00
|
|
|
(* *)
|
|
|
|
(* Valerie Menissier-Morain, projet Cristal, INRIA Rocquencourt *)
|
|
|
|
(* *)
|
1996-04-30 07:53:58 -07:00
|
|
|
(* Copyright 1996 Institut National de Recherche en Informatique et *)
|
1999-11-17 10:59:06 -08:00
|
|
|
(* en Automatique. All rights reserved. This file is distributed *)
|
2001-12-07 05:41:02 -08:00
|
|
|
(* under the terms of the GNU Library General Public License, with *)
|
|
|
|
(* the special exception on linking described in file ../../LICENSE. *)
|
1995-11-06 02:34:19 -08:00
|
|
|
(* *)
|
|
|
|
(***********************************************************************)
|
|
|
|
|
|
|
|
(* $Id$ *)
|
|
|
|
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Operations on arbitrary-precision integers.
|
1995-11-06 02:34:19 -08:00
|
|
|
|
2001-12-05 09:04:56 -08:00
|
|
|
Big integers (type [big_int]) are signed integers of arbitrary size.
|
|
|
|
*)
|
1995-11-06 02:34:19 -08:00
|
|
|
|
|
|
|
open Nat
|
|
|
|
|
|
|
|
type big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** The type of big integers. *)
|
1995-11-06 02:34:19 -08:00
|
|
|
|
|
|
|
val zero_big_int : big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** The big integer [0]. *)
|
1995-11-06 02:34:19 -08:00
|
|
|
val unit_big_int : big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** The big integer [1]. *)
|
2001-07-24 01:45:31 -07:00
|
|
|
|
2001-12-28 15:15:23 -08:00
|
|
|
(** {6 Arithmetic operations} *)
|
2001-07-24 01:45:31 -07:00
|
|
|
|
1995-11-06 02:34:19 -08:00
|
|
|
val minus_big_int : big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Unary negation. *)
|
1995-11-06 02:34:19 -08:00
|
|
|
val abs_big_int : big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Absolute value. *)
|
2001-07-24 01:45:31 -07:00
|
|
|
val add_big_int : big_int -> big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Addition. *)
|
2001-07-24 01:45:31 -07:00
|
|
|
val succ_big_int : big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Successor (add 1). *)
|
2001-07-24 01:45:31 -07:00
|
|
|
val add_int_big_int : int -> big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Addition of a small integer to a big integer. *)
|
2001-07-24 01:45:31 -07:00
|
|
|
val sub_big_int : big_int -> big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Subtraction. *)
|
2001-07-24 01:45:31 -07:00
|
|
|
val pred_big_int : big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Predecessor (subtract 1). *)
|
2001-07-24 01:45:31 -07:00
|
|
|
val mult_big_int : big_int -> big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Multiplication of two big integers. *)
|
2001-07-24 01:45:31 -07:00
|
|
|
val mult_int_big_int : int -> big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Multiplication of a big integer by a small integer *)
|
2001-07-24 01:45:31 -07:00
|
|
|
val square_big_int: big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Return the square of the given big integer *)
|
2001-07-24 01:45:31 -07:00
|
|
|
val sqrt_big_int: big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** [sqrt_big_int a] returns the integer square root of [a],
|
2001-07-24 01:45:31 -07:00
|
|
|
that is, the largest big integer [r] such that [r * r <= a].
|
|
|
|
Raise [Invalid_argument] if [a] is negative. *)
|
|
|
|
val quomod_big_int : big_int -> big_int -> big_int * big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Euclidean division of two big integers.
|
2001-07-24 01:45:31 -07:00
|
|
|
The first part of the result is the quotient,
|
|
|
|
the second part is the remainder.
|
|
|
|
Writing [(q,r) = quomod_big_int a b], we have
|
|
|
|
[a = q * b + r] and [0 <= r < |b|].
|
|
|
|
Raise [Division_by_zero] if the divisor is zero. *)
|
|
|
|
val div_big_int : big_int -> big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Euclidean quotient of two big integers.
|
2001-07-24 01:45:31 -07:00
|
|
|
This is the first result [q] of [quomod_big_int] (see above). *)
|
|
|
|
val mod_big_int : big_int -> big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Euclidean modulus of two big integers.
|
2001-07-24 01:45:31 -07:00
|
|
|
This is the second result [r] of [quomod_big_int] (see above). *)
|
|
|
|
val gcd_big_int : big_int -> big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Greatest common divisor of two big integers. *)
|
2001-07-24 01:45:31 -07:00
|
|
|
val power_int_positive_int: int -> int -> big_int
|
|
|
|
val power_big_int_positive_int: big_int -> int -> big_int
|
|
|
|
val power_int_positive_big_int: int -> big_int -> big_int
|
|
|
|
val power_big_int_positive_big_int: big_int -> big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Exponentiation functions. Return the big integer
|
2002-03-14 12:12:54 -08:00
|
|
|
representing the first argument [a] raised to the power [b]
|
|
|
|
(the second argument). Depending
|
|
|
|
on the function, [a] and [b] can be either small integers
|
|
|
|
or big integers. Raise [Invalid_argument] if [b] is negative. *)
|
2001-07-24 01:45:31 -07:00
|
|
|
|
2001-12-28 15:15:23 -08:00
|
|
|
(** {6 Comparisons and tests} *)
|
2001-07-24 01:45:31 -07:00
|
|
|
|
|
|
|
val sign_big_int : big_int -> int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Return [0] if the given big integer is zero,
|
2001-07-24 01:45:31 -07:00
|
|
|
[1] if it is positive, and [-1] if it is negative. *)
|
1995-11-06 02:34:19 -08:00
|
|
|
val compare_big_int : big_int -> big_int -> int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** [compare_big_int a b] returns [0] if [a] and [b] are equal,
|
2001-07-24 01:45:31 -07:00
|
|
|
[1] if [a] is greater than [b], and [-1] if [a] is smaller
|
|
|
|
than [b]. *)
|
1995-11-06 02:34:19 -08:00
|
|
|
val eq_big_int : big_int -> big_int -> bool
|
|
|
|
val le_big_int : big_int -> big_int -> bool
|
|
|
|
val ge_big_int : big_int -> big_int -> bool
|
|
|
|
val lt_big_int : big_int -> big_int -> bool
|
|
|
|
val gt_big_int : big_int -> big_int -> bool
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Usual boolean comparisons between two big integers. *)
|
1995-11-06 02:34:19 -08:00
|
|
|
val max_big_int : big_int -> big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Return the greater of its two arguments. *)
|
1995-11-06 02:34:19 -08:00
|
|
|
val min_big_int : big_int -> big_int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Return the smaller of its two arguments. *)
|
2001-07-24 01:45:31 -07:00
|
|
|
val num_digits_big_int : big_int -> int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Return the number of machine words used to store the
|
2001-07-24 01:45:31 -07:00
|
|
|
given big integer. *)
|
|
|
|
|
2001-12-28 15:15:23 -08:00
|
|
|
(** {6 Conversions to and from strings} *)
|
2001-07-24 01:45:31 -07:00
|
|
|
|
|
|
|
val string_of_big_int : big_int -> string
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Return the string representation of the given big integer,
|
2001-07-24 01:45:31 -07:00
|
|
|
in decimal (base 10). *)
|
|
|
|
val big_int_of_string : string -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Convert a string to a big integer, in decimal.
|
2001-07-24 01:45:31 -07:00
|
|
|
The string consists of an optional [-] or [+] sign,
|
|
|
|
followed by one or several decimal digits. *)
|
|
|
|
|
2001-12-28 15:15:23 -08:00
|
|
|
(** {6 Conversions to and from other numerical types} *)
|
2001-07-24 01:45:31 -07:00
|
|
|
|
1995-11-06 02:34:19 -08:00
|
|
|
val big_int_of_int : int -> big_int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Convert a small integer to a big integer. *)
|
1995-11-06 02:34:19 -08:00
|
|
|
val is_int_big_int : big_int -> bool
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Test whether the given big integer is small enough to
|
2001-07-24 01:45:31 -07:00
|
|
|
be representable as a small integer (type [int])
|
|
|
|
without loss of precision. On a 32-bit platform,
|
|
|
|
[is_int_big_int a] returns [true] if and only if
|
2001-12-05 09:04:56 -08:00
|
|
|
[a] is between 2{^30} and 2{^30}-1. On a 64-bit platform,
|
2001-07-24 01:45:31 -07:00
|
|
|
[is_int_big_int a] returns [true] if and only if
|
2001-12-05 09:04:56 -08:00
|
|
|
[a] is between -2{^62} and 2{^62}-1. *)
|
2001-07-24 01:45:31 -07:00
|
|
|
val int_of_big_int : big_int -> int
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Convert a big integer to a small integer (type [int]).
|
2001-07-24 01:45:31 -07:00
|
|
|
Raises [Failure "int_of_big_int"] if the big integer
|
|
|
|
is not representable as a small integer. *)
|
|
|
|
val float_of_big_int : big_int -> float
|
2001-12-05 09:04:56 -08:00
|
|
|
(** Returns a floating-point number approximating the
|
2001-07-24 01:45:31 -07:00
|
|
|
given big integer. *)
|
|
|
|
|
2001-12-05 09:04:56 -08:00
|
|
|
(**/**)
|
2001-07-24 01:45:31 -07:00
|
|
|
|
2001-12-28 15:15:23 -08:00
|
|
|
(** {6 For internal use} *)
|
1995-11-06 02:34:19 -08:00
|
|
|
val nat_of_big_int : big_int -> nat
|
|
|
|
val big_int_of_nat : nat -> big_int
|
|
|
|
val base_power_big_int: int -> int -> big_int -> big_int
|
|
|
|
val sys_big_int_of_string: string -> int -> int -> big_int
|
|
|
|
val round_futur_last_digit : string -> int -> int -> bool
|
|
|
|
val approx_big_int: int -> big_int -> string
|