2 pmod - integral power of an integer modulo a specified integer
15 pmod(x, n, md) returns the integer value of the canonical residue
16 of x^n modulo md, where the set of canonical residues is determined
17 by md and bits 0, 2, and 4 of config("mod") (other bits are ignored).
19 If md is zero, the value is simply x^n.
21 For nonzero md, the canonical residues v modulo md are as follows:
23 config("mod") md > 0 md < 0
25 0 0 < v < md md < v < 0
26 1 -md < v < 0 0 < v < -md
27 4 0 < v < md 0 < v < -md
28 5 -md < v < 0 md < v < 0
29 16 -md/2 < v <= md/2 md/2 <= v < -md/2
30 17 -md/2 <= v < md/2 md/2 < v <= -md/2
31 20 -md/2 < v <= md/2 md/2 < v <= -md/2
32 21 -md/2 <= v < md/2 md/2 <= v < -md/2
36 > print pmod(2,3,10), pmod(2,5,10), pmod(2,3,-10), pod(2,5,-10)
39 > c = config("mod",16)
40 > print pmod(2,3,10), pmod(2,5,10), pmod(2,3,-10), pmod(2,5,-10)
47 NUMBER *qpowermod(NUMBER *x, NUMBER *n, NUMBER *md)