#include <mp.h>
#include <stdio.h>
typedef struct mint { int len; short *val; } MINT;
madd(a, b, c)
msub(a, b, c)
mult(a, b, c)
mdiv(a, b, q, r)
pow(a, b, m, c)
gcd(a, b, c)
invert(a, b, c)
rpow(a, n, c)
msqrt(a, b, r)
mcmp(a, b)
move(a, b)
min(a)
omin(a)
fmin(a, f)
m_in(a, n, f)
mout(a)
omout(a)
fmout(a, f)
m_out(a, n, f)
MINT *a, *b, *c, *m, *q, *r;
FILE *f;
int n;
sdiv(a, n, q, r)
MINT *a, *q;
short n;
short *r;
MINT *itom(n)
madd, msub and mult assign to c the sum, difference and product, respectively, of a and b. mdiv assigns to q and r the quotient and remainder obtained from dividing a by b. sdiv is like mdiv except that the divisor is a short integer n and the remainder is placed in a short whose address is given as r. msqrt produces the integer square root of a in b and places the remainder in r. rpow calculates in c the value of a raised to the (``regular'' integral) power n, while pow calculates this with a full multiple precision exponent b and the result is reduced modulo m. gcd returns the greatest common denominator of a and b in c, and invert computes c such that a*c mod b = 1, for a and b relatively prime. mcmp returns a negative, zero or positive integer value when a is less than, equal to or greater than b, respectively. move copies a to b. min and mout do decimal input and output while omin and omout do octal input and output. More generally, fmin and fmout do decimal input and output using file f, and m_in and m_out do I/O with arbitrary radix n. On input, records should have the form of strings of digits terminated by a newline; output records have a similar form.
Programs which use the multiple-precision arithmetic library must be loaded using the loader flag -lmp.
dc(1) and bc(1) don't use this library.
The input and output routines are a crock.
pow is also the name of a standard math library routine.