/ Usenet 1994 January / usenetsourcesnewsgroupsinfomagicjanuary1994.iso / sources / unix / volume26 / calc / part16

< prev

next >

MacOS 8.1

Win98

DOS

browse contents

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 4e 65 77 73 67 72 6f 75 | 70 73 3a 20 63 6f 6d 70 |Newsgrou|ps: comp|
|00000010| 2e 73 6f 75 72 63 65 73 | 2e 75 6e 69 78 0a 46 72 |.sources|.unix.Fr|
|00000020| 6f 6d 3a 20 64 62 65 6c | 6c 40 70 64 61 63 74 2e |om: dbel|l@pdact.|
|00000030| 70 64 2e 6e 65 63 69 73 | 61 2e 6f 7a 2e 61 75 20 |pd.necis|a.oz.au |
|00000040| 28 44 61 76 69 64 20 49 | 2e 20 42 65 6c 6c 29 0a |(David I|. Bell).|
|00000050| 53 75 62 6a 65 63 74 3a | 20 76 32 36 69 30 34 32 |Subject:| v26i042|
|00000060| 3a 20 43 41 4c 43 20 2d | 20 41 6e 20 61 72 62 69 |: CALC -| An arbi|
|00000070| 74 72 61 72 79 20 70 72 | 65 63 69 73 69 6f 6e 20 |trary pr|ecision |
|00000080| 43 2d 6c 69 6b 65 20 63 | 61 6c 63 75 6c 61 74 6f |C-like c|alculato|
|00000090| 72 2c 20 50 61 72 74 31 | 36 2f 32 31 0a 53 65 6e |r, Part1|6/21.Sen|
|000000a0| 64 65 72 3a 20 75 6e 69 | 78 2d 73 6f 75 72 63 65 |der: uni|x-source|
|000000b0| 73 2d 6d 6f 64 65 72 61 | 74 6f 72 40 70 61 2e 64 |s-modera|tor@pa.d|
|000000c0| 65 63 2e 63 6f 6d 0a 41 | 70 70 72 6f 76 65 64 3a |ec.com.A|pproved:|
|000000d0| 20 76 69 78 69 65 40 70 | 61 2e 64 65 63 2e 63 6f | vixie@p|a.dec.co|
|000000e0| 6d 0a 0a 53 75 62 6d 69 | 74 74 65 64 2d 42 79 3a |m..Submi|tted-By:|
|000000f0| 20 64 62 65 6c 6c 40 70 | 64 61 63 74 2e 70 64 2e | dbell@p|dact.pd.|
|00000100| 6e 65 63 69 73 61 2e 6f | 7a 2e 61 75 20 28 44 61 |necisa.o|z.au (Da|
|00000110| 76 69 64 20 49 2e 20 42 | 65 6c 6c 29 0a 50 6f 73 |vid I. B|ell).Pos|
|00000120| 74 69 6e 67 2d 4e 75 6d | 62 65 72 3a 20 56 6f 6c |ting-Num|ber: Vol|
|00000130| 75 6d 65 20 32 36 2c 20 | 49 73 73 75 65 20 34 32 |ume 26, |Issue 42|
|00000140| 0a 41 72 63 68 69 76 65 | 2d 4e 61 6d 65 3a 20 63 |.Archive|-Name: c|
|00000150| 61 6c 63 2f 70 61 72 74 | 31 36 0a 0a 23 21 20 2f |alc/part|16..#! /|
|00000160| 62 69 6e 2f 73 68 0a 23 | 20 54 68 69 73 20 69 73 |bin/sh.#| This is|
|00000170| 20 61 20 73 68 65 6c 6c | 20 61 72 63 68 69 76 65 | a shell| archive|
|00000180| 2e 20 20 52 65 6d 6f 76 | 65 20 61 6e 79 74 68 69 |.  Remov|e anythi|
|00000190| 6e 67 20 62 65 66 6f 72 | 65 20 74 68 69 73 20 6c |ng befor|e this l|
|000001a0| 69 6e 65 2c 20 74 68 65 | 6e 20 75 6e 70 61 63 6b |ine, the|n unpack|
|000001b0| 0a 23 20 69 74 20 62 79 | 20 73 61 76 69 6e 67 20 |.# it by| saving |
|000001c0| 69 74 20 69 6e 74 6f 20 | 61 20 66 69 6c 65 20 61 |it into |a file a|
|000001d0| 6e 64 20 74 79 70 69 6e | 67 20 22 73 68 20 66 69 |nd typin|g "sh fi|
|000001e0| 6c 65 22 2e 20 20 54 6f | 20 6f 76 65 72 77 72 69 |le".  To| overwri|
|000001f0| 74 65 20 65 78 69 73 74 | 69 6e 67 0a 23 20 66 69 |te exist|ing.# fi|
|00000200| 6c 65 73 2c 20 74 79 70 | 65 20 22 73 68 20 66 69 |les, typ|e "sh fi|
|00000210| 6c 65 20 2d 63 22 2e 20 | 20 59 6f 75 20 63 61 6e |le -c". | You can|
|00000220| 20 61 6c 73 6f 20 66 65 | 65 64 20 74 68 69 73 20 | also fe|ed this |
|00000230| 61 73 20 73 74 61 6e 64 | 61 72 64 20 69 6e 70 75 |as stand|ard inpu|
|00000240| 74 20 76 69 61 0a 23 20 | 75 6e 73 68 61 72 2c 20 |t via.# |unshar, |
|00000250| 6f 72 20 62 79 20 74 79 | 70 69 6e 67 20 22 73 68 |or by ty|ping "sh|
|00000260| 20 3c 66 69 6c 65 22 2c | 20 65 2e 67 2e 2e 20 20 | <file",| e.g..  |
|00000270| 49 66 20 74 68 69 73 20 | 61 72 63 68 69 76 65 20 |If this |archive |
|00000280| 69 73 20 63 6f 6d 70 6c | 65 74 65 2c 20 79 6f 75 |is compl|ete, you|
|00000290| 0a 23 20 77 69 6c 6c 20 | 73 65 65 20 74 68 65 20 |.# will |see the |
|000002a0| 66 6f 6c 6c 6f 77 69 6e | 67 20 6d 65 73 73 61 67 |followin|g messag|
|000002b0| 65 20 61 74 20 74 68 65 | 20 65 6e 64 3a 0a 23 09 |e at the| end:.#.|
|000002c0| 09 22 45 6e 64 20 6f 66 | 20 61 72 63 68 69 76 65 |."End of| archive|
|000002d0| 20 31 36 20 28 6f 66 20 | 32 31 29 2e 22 0a 23 20 | 16 (of |21).".# |
|000002e0| 43 6f 6e 74 65 6e 74 73 | 3a 20 20 7a 6d 6f 64 2e |Contents|:  zmod.|
|000002f0| 63 0a 23 20 57 72 61 70 | 70 65 64 20 62 79 20 64 |c.# Wrap|ped by d|
|00000300| 62 65 6c 6c 40 65 6c 6d | 20 6f 6e 20 54 75 65 20 |bell@elm| on Tue |
|00000310| 46 65 62 20 32 35 20 31 | 35 3a 32 31 3a 31 34 20 |Feb 25 1|5:21:14 |
|00000320| 31 39 39 32 0a 50 41 54 | 48 3d 2f 62 69 6e 3a 2f |1992.PAT|H=/bin:/|
|00000330| 75 73 72 2f 62 69 6e 3a | 2f 75 73 72 2f 75 63 62 |usr/bin:|/usr/ucb|
|00000340| 20 3b 20 65 78 70 6f 72 | 74 20 50 41 54 48 0a 69 | ; expor|t PATH.i|
|00000350| 66 20 74 65 73 74 20 2d | 66 20 27 7a 6d 6f 64 2e |f test -|f 'zmod.|
|00000360| 63 27 20 2d 61 20 22 24 | 7b 31 7d 22 20 21 3d 20 |c' -a "$|{1}" != |
|00000370| 22 2d 63 22 20 3b 20 74 | 68 65 6e 20 0a 20 20 65 |"-c" ; t|hen .  e|
|00000380| 63 68 6f 20 73 68 61 72 | 3a 20 57 69 6c 6c 20 6e |cho shar|: Will n|
|00000390| 6f 74 20 63 6c 6f 62 62 | 65 72 20 65 78 69 73 74 |ot clobb|er exist|
|000003a0| 69 6e 67 20 66 69 6c 65 | 20 5c 22 27 7a 6d 6f 64 |ing file| \"'zmod|
|000003b0| 2e 63 27 5c 22 0a 65 6c | 73 65 0a 65 63 68 6f 20 |.c'\".el|se.echo |
|000003c0| 73 68 61 72 3a 20 45 78 | 74 72 61 63 74 69 6e 67 |shar: Ex|tracting|
|000003d0| 20 5c 22 27 7a 6d 6f 64 | 2e 63 27 5c 22 20 5c 28 | \"'zmod|.c'\" \(|
|000003e0| 33 32 34 36 38 20 63 68 | 61 72 61 63 74 65 72 73 |32468 ch|aracters|
|000003f0| 5c 29 0a 73 65 64 20 22 | 73 2f 5e 58 2f 2f 22 20 |\).sed "|s/^X//" |
|00000400| 3e 27 7a 6d 6f 64 2e 63 | 27 20 3c 3c 27 45 4e 44 |>'zmod.c|' <<'END|
|00000410| 5f 4f 46 5f 46 49 4c 45 | 27 0a 58 2f 2a 0a 58 20 |_OF_FILE|'.X/*.X |
|00000420| 2a 20 43 6f 70 79 72 69 | 67 68 74 20 28 63 29 20 |* Copyri|ght (c) |
|00000430| 31 39 39 32 20 44 61 76 | 69 64 20 49 2e 20 42 65 |1992 Dav|id I. Be|
|00000440| 6c 6c 0a 58 20 2a 20 50 | 65 72 6d 69 73 73 69 6f |ll.X * P|ermissio|
|00000450| 6e 20 69 73 20 67 72 61 | 6e 74 65 64 20 74 6f 20 |n is gra|nted to |
|00000460| 75 73 65 2c 20 64 69 73 | 74 72 69 62 75 74 65 2c |use, dis|tribute,|
|00000470| 20 6f 72 20 6d 6f 64 69 | 66 79 20 74 68 69 73 20 | or modi|fy this |
|00000480| 73 6f 75 72 63 65 2c 0a | 58 20 2a 20 70 72 6f 76 |source,.|X * prov|
|00000490| 69 64 65 64 20 74 68 61 | 74 20 74 68 69 73 20 63 |ided tha|t this c|
|000004a0| 6f 70 79 72 69 67 68 74 | 20 6e 6f 74 69 63 65 20 |opyright| notice |
|000004b0| 72 65 6d 61 69 6e 73 20 | 69 6e 74 61 63 74 2e 0a |remains |intact..|
|000004c0| 58 20 2a 0a 58 20 2a 20 | 52 6f 75 74 69 6e 65 73 |X *.X * |Routines|
|000004d0| 20 74 6f 20 64 6f 20 6d | 6f 64 75 6c 6f 20 61 72 | to do m|odulo ar|
|000004e0| 69 74 68 6d 65 74 69 63 | 20 62 6f 74 68 20 6e 6f |ithmetic| both no|
|000004f0| 72 6d 61 6c 6c 79 20 61 | 6e 64 20 61 6c 73 6f 20 |rmally a|nd also |
|00000500| 75 73 69 6e 67 20 74 68 | 65 20 52 45 44 43 0a 58 |using th|e REDC.X|
|00000510| 20 2a 20 61 6c 67 6f 72 | 69 74 68 6d 20 67 69 76 | * algor|ithm giv|
|00000520| 65 6e 20 62 79 20 50 65 | 74 65 72 20 4c 2e 20 4d |en by Pe|ter L. M|
|00000530| 6f 6e 74 67 6f 6d 65 72 | 79 20 69 6e 20 4d 61 74 |ontgomer|y in Mat|
|00000540| 68 65 6d 61 74 69 63 73 | 20 6f 66 20 43 6f 6d 70 |hematics| of Comp|
|00000550| 75 74 61 74 69 6f 6e 2c | 0a 58 20 2a 20 76 6f 6c |utation,|.X * vol|
|00000560| 75 6d 65 20 34 34 2c 20 | 6e 75 6d 62 65 72 20 31 |ume 44, |number 1|
|00000570| 37 30 20 28 41 70 72 69 | 6c 2c 20 31 39 38 35 29 |70 (Apri|l, 1985)|
|00000580| 2e 20 20 46 6f 72 20 6d | 75 6c 74 69 70 6c 65 20 |.  For m|ultiple |
|00000590| 6d 75 6c 74 69 70 6c 69 | 65 73 20 75 73 69 6e 67 |multipli|es using|
|000005a0| 0a 58 20 2a 20 74 68 65 | 20 73 61 6d 65 20 6c 61 |.X * the| same la|
|000005b0| 72 67 65 20 6d 6f 64 75 | 6c 75 73 2c 20 74 68 65 |rge modu|lus, the|
|000005c0| 20 52 45 44 43 20 61 6c | 67 6f 72 69 74 68 6d 20 | REDC al|gorithm |
|000005d0| 61 76 6f 69 64 73 20 74 | 68 65 20 75 73 75 61 6c |avoids t|he usual|
|000005e0| 20 64 69 76 69 73 69 6f | 6e 0a 58 20 2a 20 62 79 | divisio|n.X * by|
|000005f0| 20 74 68 65 20 6d 6f 64 | 75 6c 75 73 2c 20 69 6e | the mod|ulus, in|
|00000600| 73 74 65 61 64 20 72 65 | 70 6c 61 63 69 6e 67 20 |stead re|placing |
|00000610| 69 74 20 77 69 74 68 20 | 74 77 6f 20 6d 75 6c 74 |it with |two mult|
|00000620| 69 70 6c 69 65 73 20 6f | 72 20 65 6c 73 65 20 61 |iplies o|r else a|
|00000630| 0a 58 20 2a 20 73 70 65 | 63 69 61 6c 20 61 6c 67 |.X * spe|cial alg|
|00000640| 6f 72 69 74 68 6d 2e 20 | 20 57 68 65 6e 20 74 68 |orithm. | When th|
|00000650| 65 73 65 20 74 77 6f 20 | 6d 75 6c 74 69 70 6c 69 |ese two |multipli|
|00000660| 65 73 20 6f 72 20 74 68 | 65 20 73 70 65 63 69 61 |es or th|e specia|
|00000670| 6c 20 61 6c 67 6f 72 69 | 74 68 6d 0a 58 20 2a 20 |l algori|thm.X * |
|00000680| 61 72 65 20 66 61 73 74 | 65 72 20 74 68 65 6e 20 |are fast|er then |
|00000690| 74 68 65 20 64 69 76 69 | 73 69 6f 6e 2c 20 74 68 |the divi|sion, th|
|000006a0| 65 6e 20 74 68 65 20 52 | 45 44 43 20 61 6c 67 6f |en the R|EDC algo|
|000006b0| 72 69 74 68 6d 20 69 73 | 20 62 65 74 74 65 72 2e |rithm is| better.|
|000006c0| 0a 58 20 2a 2f 0a 58 0a | 58 23 69 6e 63 6c 75 64 |.X */.X.|X#includ|
|000006d0| 65 20 22 6d 61 74 68 2e | 68 22 0a 58 0a 58 0a 58 |e "math.|h".X.X.X|
|000006e0| 23 64 65 66 69 6e 65 09 | 50 4f 57 42 49 54 53 09 |#define.|POWBITS.|
|000006f0| 34 09 09 2f 2a 20 62 69 | 74 73 20 66 6f 72 20 70 |4../* bi|ts for p|
|00000700| 6f 77 65 72 20 63 68 75 | 6e 6b 73 20 28 6d 75 73 |ower chu|nks (mus|
|00000710| 74 20 64 69 76 69 64 65 | 20 42 41 53 45 42 29 20 |t divide| BASEB) |
|00000720| 2a 2f 0a 58 23 64 65 66 | 69 6e 65 09 50 4f 57 4e |*/.X#def|ine.POWN|
|00000730| 55 4d 53 09 28 31 3c 3c | 50 4f 57 42 49 54 53 29 |UMS.(1<<|POWBITS)|
|00000740| 09 2f 2a 20 6e 75 6d 62 | 65 72 20 6f 66 20 70 6f |./* numb|er of po|
|00000750| 77 65 72 73 20 6e 65 65 | 64 65 64 20 69 6e 20 74 |wers nee|ded in t|
|00000760| 61 62 6c 65 20 2a 2f 0a | 58 0a 58 0a 58 4c 45 4e |able */.|X.X.XLEN|
|00000770| 20 5f 70 6f 77 32 5f 20 | 3d 20 50 4f 57 5f 41 4c | _pow2_ |= POW_AL|
|00000780| 47 32 3b 09 09 2f 2a 20 | 6d 6f 64 75 6c 6f 20 73 |G2;../* |modulo s|
|00000790| 69 7a 65 20 74 6f 20 75 | 73 65 20 52 45 44 43 20 |ize to u|se REDC |
|000007a0| 66 6f 72 20 70 6f 77 65 | 72 73 20 2a 2f 0a 58 4c |for powe|rs */.XL|
|000007b0| 45 4e 20 5f 72 65 64 63 | 32 5f 20 3d 20 52 45 44 |EN _redc|2_ = RED|
|000007c0| 43 5f 41 4c 47 32 3b 09 | 2f 2a 20 6d 6f 64 75 6c |C_ALG2;.|/* modul|
|000007d0| 6f 20 73 69 7a 65 20 74 | 6f 20 75 73 65 20 73 65 |o size t|o use se|
|000007e0| 63 6f 6e 64 20 52 45 44 | 43 20 61 6c 67 6f 72 69 |cond RED|C algori|
|000007f0| 74 68 6d 20 2a 2f 0a 58 | 0a 58 73 74 61 74 69 63 |thm */.X|.Xstatic|
|00000800| 20 52 45 44 43 20 2a 70 | 6f 77 65 72 6d 6f 64 72 | REDC *p|owermodr|
|00000810| 65 64 63 20 3d 20 4e 55 | 4c 4c 3b 09 2f 2a 20 52 |edc = NU|LL;./* R|
|00000820| 45 44 43 20 69 6e 66 6f | 20 66 6f 72 20 72 61 69 |EDC info| for rai|
|00000830| 73 69 6e 67 20 74 6f 20 | 70 6f 77 65 72 20 2a 2f |sing to |power */|
|00000840| 0a 58 0a 58 23 69 66 20 | 30 0a 58 65 78 74 65 72 |.X.X#if |0.Xexter|
|00000850| 6e 20 76 6f 69 64 20 7a | 61 64 64 6d 6f 64 20 70 |n void z|addmod p|
|00000860| 72 6f 74 6f 28 28 5a 56 | 41 4c 55 45 20 7a 31 2c |roto((ZV|ALUE z1,|
|00000870| 20 5a 56 41 4c 55 45 20 | 7a 32 2c 20 5a 56 41 4c | ZVALUE |z2, ZVAL|
|00000880| 55 45 20 7a 33 2c 20 5a | 56 41 4c 55 45 20 2a 72 |UE z3, Z|VALUE *r|
|00000890| 65 73 29 29 3b 0a 58 65 | 78 74 65 72 6e 20 76 6f |es));.Xe|xtern vo|
|000008a0| 69 64 20 7a 6e 65 67 6d | 6f 64 20 70 72 6f 74 6f |id znegm|od proto|
|000008b0| 28 28 5a 56 41 4c 55 45 | 20 7a 31 2c 20 5a 56 41 |((ZVALUE| z1, ZVA|
|000008c0| 4c 55 45 20 7a 32 2c 20 | 5a 56 41 4c 55 45 20 2a |LUE z2, |ZVALUE *|
|000008d0| 72 65 73 29 29 3b 0a 58 | 0a 58 2f 2a 0a 58 20 2a |res));.X|.X/*.X *|
|000008e0| 20 4d 75 6c 74 69 70 6c | 79 20 74 77 6f 20 6e 75 | Multipl|y two nu|
|000008f0| 6d 62 65 72 73 20 74 6f | 67 65 74 68 65 72 20 61 |mbers to|gether a|
|00000900| 6e 64 20 74 68 65 6e 20 | 6d 6f 64 20 74 68 65 20 |nd then |mod the |
|00000910| 72 65 73 75 6c 74 20 77 | 69 74 68 20 61 20 74 68 |result w|ith a th|
|00000920| 69 72 64 20 6e 75 6d 62 | 65 72 2e 0a 58 20 2a 20 |ird numb|er..X * |
|00000930| 54 68 65 20 74 77 6f 20 | 6e 75 6d 62 65 72 73 20 |The two |numbers |
|00000940| 74 6f 20 62 65 20 6d 75 | 6c 74 69 70 6c 69 65 64 |to be mu|ltiplied|
|00000950| 20 63 61 6e 20 62 65 20 | 6e 65 67 61 74 69 76 65 | can be |negative|
|00000960| 20 6f 72 20 6f 75 74 20 | 6f 66 20 6d 6f 64 75 6c | or out |of modul|
|00000970| 6f 20 72 61 6e 67 65 2e | 0a 58 20 2a 20 54 68 65 |o range.|.X * The|
|00000980| 20 72 65 73 75 6c 74 20 | 77 69 6c 6c 20 62 65 20 | result |will be |
|00000990| 69 6e 20 74 68 65 20 72 | 61 6e 67 65 20 30 20 74 |in the r|ange 0 t|
|000009a0| 6f 20 74 68 65 20 6d 6f | 64 75 6c 75 73 20 2d 20 |o the mo|dulus - |
|000009b0| 31 2e 0a 58 20 2a 2f 0a | 58 76 6f 69 64 0a 58 7a |1..X */.|Xvoid.Xz|
|000009c0| 6d 75 6c 6d 6f 64 28 7a | 31 2c 20 7a 32 2c 20 7a |mulmod(z|1, z2, z|
|000009d0| 33 2c 20 72 65 73 29 0a | 58 09 5a 56 41 4c 55 45 |3, res).|X.ZVALUE|
|000009e0| 20 7a 31 3b 09 09 2f 2a | 20 66 69 72 73 74 20 6e | z1;../*| first n|
|000009f0| 75 6d 62 65 72 20 74 6f | 20 62 65 20 6d 75 6c 74 |umber to| be mult|
|00000a00| 69 70 6c 69 65 64 20 2a | 2f 0a 58 09 5a 56 41 4c |iplied *|/.X.ZVAL|
|00000a10| 55 45 20 7a 32 3b 09 09 | 2f 2a 20 73 65 63 6f 6e |UE z2;..|/* secon|
|00000a20| 64 20 6e 75 6d 62 65 72 | 20 74 6f 20 62 65 20 6d |d number| to be m|
|00000a30| 75 6c 74 69 70 6c 69 65 | 64 20 2a 2f 0a 58 09 5a |ultiplie|d */.X.Z|
|00000a40| 56 41 4c 55 45 20 7a 33 | 3b 09 09 2f 2a 20 6e 75 |VALUE z3|;../* nu|
|00000a50| 6d 62 65 72 20 74 6f 20 | 74 61 6b 65 20 6d 6f 64 |mber to |take mod|
|00000a60| 20 77 69 74 68 20 2a 2f | 0a 58 09 5a 56 41 4c 55 | with */|.X.ZVALU|
|00000a70| 45 20 2a 72 65 73 3b 09 | 09 2f 2a 20 72 65 73 75 |E *res;.|./* resu|
|00000a80| 6c 74 20 2a 2f 0a 58 7b | 0a 58 09 5a 56 41 4c 55 |lt */.X{|.X.ZVALU|
|00000a90| 45 20 74 6d 70 3b 0a 58 | 09 46 55 4c 4c 20 70 72 |E tmp;.X|.FULL pr|
|00000aa0| 6f 64 3b 0a 58 09 46 55 | 4c 4c 20 64 69 67 69 74 |od;.X.FU|LL digit|
|00000ab0| 3b 0a 58 09 42 4f 4f 4c | 20 6e 65 67 3b 0a 58 0a |;.X.BOOL| neg;.X.|
|00000ac0| 58 09 69 66 20 28 69 73 | 7a 65 72 6f 28 7a 33 29 |X.if (is|zero(z3)|
|00000ad0| 20 7c 7c 20 69 73 6e 65 | 67 28 7a 33 29 29 0a 58 | || isne|g(z3)).X|
|00000ae0| 09 09 65 72 72 6f 72 28 | 22 4d 6f 64 20 6f 66 20 |..error(|"Mod of |
|00000af0| 6e 6f 6e 2d 70 6f 73 69 | 74 69 76 65 20 69 6e 74 |non-posi|tive int|
|00000b00| 65 67 65 72 22 29 3b 0a | 58 09 69 66 20 28 69 73 |eger");.|X.if (is|
|00000b10| 7a 65 72 6f 28 7a 31 29 | 20 7c 7c 20 69 73 7a 65 |zero(z1)| || isze|
|00000b20| 72 6f 28 7a 32 29 20 7c | 7c 20 69 73 75 6e 69 74 |ro(z2) ||| isunit|
|00000b30| 28 7a 33 29 29 20 7b 0a | 58 09 09 2a 72 65 73 20 |(z3)) {.|X..*res |
|00000b40| 3d 20 5f 7a 65 72 6f 5f | 3b 0a 58 09 09 72 65 74 |= _zero_|;.X..ret|
|00000b50| 75 72 6e 3b 0a 58 09 7d | 0a 58 0a 58 09 2f 2a 0a |urn;.X.}|.X.X./*.|
|00000b60| 58 09 20 2a 20 49 66 20 | 74 68 65 20 6d 6f 64 75 |X. * If |the modu|
|00000b70| 6c 75 73 20 69 73 20 61 | 20 73 69 6e 67 6c 65 20 |lus is a| single |
|00000b80| 64 69 67 69 74 20 6e 75 | 6d 62 65 72 2c 20 74 68 |digit nu|mber, th|
|00000b90| 65 6e 20 64 6f 20 74 68 | 65 20 72 65 73 75 6c 74 |en do th|e result|
|00000ba0| 0a 58 09 20 2a 20 63 68 | 65 61 70 6c 79 2e 20 20 |.X. * ch|eaply.  |
|00000bb0| 43 68 65 63 6b 20 65 73 | 70 65 63 69 61 6c 6c 79 |Check es|pecially|
|00000bc0| 20 66 6f 72 20 61 20 73 | 6d 61 6c 6c 20 70 6f 77 | for a s|mall pow|
|00000bd0| 65 72 20 6f 66 20 74 77 | 6f 2e 0a 58 09 20 2a 2f |er of tw|o..X. */|
|00000be0| 0a 58 09 69 66 20 28 69 | 73 74 69 6e 79 28 7a 33 |.X.if (i|stiny(z3|
|00000bf0| 29 29 20 7b 0a 58 09 09 | 6e 65 67 20 3d 20 28 7a |)) {.X..|neg = (z|
|00000c00| 31 2e 73 69 67 6e 20 21 | 3d 20 7a 32 2e 73 69 67 |1.sign !|= z2.sig|
|00000c10| 6e 29 3b 0a 58 09 09 64 | 69 67 69 74 20 3d 20 7a |n);.X..d|igit = z|
|00000c20| 33 2e 76 5b 30 5d 3b 0a | 58 09 09 69 66 20 28 28 |3.v[0];.|X..if ((|
|00000c30| 64 69 67 69 74 20 26 20 | 2d 64 69 67 69 74 29 20 |digit & |-digit) |
|00000c40| 3d 3d 20 64 69 67 69 74 | 29 20 7b 09 2f 2a 20 4e |== digit|) {./* N|
|00000c50| 45 45 44 53 20 32 27 53 | 20 43 4f 4d 50 20 2a 2f |EEDS 2'S| COMP */|
|00000c60| 0a 58 09 09 09 70 72 6f | 64 20 3d 20 28 28 46 55 |.X...pro|d = ((FU|
|00000c70| 4c 4c 29 20 7a 31 2e 76 | 5b 30 5d 29 20 2a 20 28 |LL) z1.v|[0]) * (|
|00000c80| 28 46 55 4c 4c 29 20 7a | 32 2e 76 5b 30 5d 29 3b |(FULL) z|2.v[0]);|
|00000c90| 0a 58 09 09 09 70 72 6f | 64 20 26 3d 20 28 64 69 |.X...pro|d &= (di|
|00000ca0| 67 69 74 20 2d 20 31 29 | 3b 0a 58 09 09 7d 20 65 |git - 1)|;.X..} e|
|00000cb0| 6c 73 65 20 7b 0a 58 09 | 09 09 7a 31 2e 73 69 67 |lse {.X.|..z1.sig|
|00000cc0| 6e 20 3d 20 30 3b 0a 58 | 09 09 09 7a 32 2e 73 69 |n = 0;.X|...z2.si|
|00000cd0| 67 6e 20 3d 20 30 3b 0a | 58 09 09 09 70 72 6f 64 |gn = 0;.|X...prod|
|00000ce0| 20 3d 20 28 46 55 4c 4c | 29 20 7a 6d 6f 64 69 28 | = (FULL|) zmodi(|
|00000cf0| 7a 31 2c 20 28 6c 6f 6e | 67 29 20 64 69 67 69 74 |z1, (lon|g) digit|
|00000d00| 29 3b 0a 58 09 09 09 70 | 72 6f 64 20 2a 3d 20 28 |);.X...p|rod *= (|
|00000d10| 46 55 4c 4c 29 20 7a 6d | 6f 64 69 28 7a 32 2c 20 |FULL) zm|odi(z2, |
|00000d20| 28 6c 6f 6e 67 29 20 64 | 69 67 69 74 29 3b 0a 58 |(long) d|igit);.X|
|00000d30| 09 09 09 70 72 6f 64 20 | 25 3d 20 64 69 67 69 74 |...prod |%= digit|
|00000d40| 3b 0a 58 09 09 7d 0a 58 | 09 09 69 66 20 28 6e 65 |;.X..}.X|..if (ne|
|00000d50| 67 20 26 26 20 70 72 6f | 64 29 0a 58 09 09 09 70 |g && pro|d).X...p|
|00000d60| 72 6f 64 20 3d 20 64 69 | 67 69 74 20 2d 20 70 72 |rod = di|git - pr|
|00000d70| 6f 64 3b 0a 58 09 09 69 | 74 6f 7a 28 28 6c 6f 6e |od;.X..i|toz((lon|
|00000d80| 67 29 20 70 72 6f 64 2c | 20 72 65 73 29 3b 0a 58 |g) prod,| res);.X|
|00000d90| 09 09 72 65 74 75 72 6e | 3b 0a 58 09 7d 0a 58 0a |..return|;.X.}.X.|
|00000da0| 58 09 2f 2a 0a 58 09 20 | 2a 20 54 68 65 20 6d 6f |X./*.X. |* The mo|
|00000db0| 64 75 6c 75 73 20 69 73 | 20 6d 6f 72 65 20 74 68 |dulus is| more th|
|00000dc0| 61 6e 20 6f 6e 65 20 64 | 69 67 69 74 2e 0a 58 09 |an one d|igit..X.|
|00000dd0| 20 2a 20 41 63 74 75 61 | 6c 6c 79 20 64 6f 20 74 | * Actua|lly do t|
|00000de0| 68 65 20 6d 75 6c 74 69 | 70 6c 79 20 61 6e 64 20 |he multi|ply and |
|00000df0| 64 69 76 69 64 65 20 69 | 66 20 6e 65 63 65 73 73 |divide i|f necess|
|00000e00| 61 72 79 2e 0a 58 09 20 | 2a 2f 0a 58 09 7a 6d 75 |ary..X. |*/.X.zmu|
|00000e10| 6c 28 7a 31 2c 20 7a 32 | 2c 20 26 74 6d 70 29 3b |l(z1, z2|, &tmp);|
|00000e20| 0a 58 09 69 66 20 28 69 | 73 70 6f 73 28 74 6d 70 |.X.if (i|spos(tmp|
|00000e30| 29 20 26 26 20 28 28 74 | 6d 70 2e 6c 65 6e 20 3c |) && ((t|mp.len <|
|00000e40| 20 7a 33 2e 6c 65 6e 29 | 20 7c 7c 20 28 28 74 6d | z3.len)| || ((tm|
|00000e50| 70 2e 6c 65 6e 20 3d 3d | 20 7a 33 2e 6c 65 6e 29 |p.len ==| z3.len)|
|00000e60| 20 26 26 0a 58 09 09 28 | 74 6d 70 2e 76 5b 74 6d | &&.X..(|tmp.v[tm|
|00000e70| 70 2e 6c 65 6e 2d 31 5d | 20 3c 20 7a 32 2e 76 5b |p.len-1]| < z2.v[|
|00000e80| 7a 33 2e 6c 65 6e 2d 31 | 5d 29 29 29 29 0a 58 09 |z3.len-1|])))).X.|
|00000e90| 7b 0a 58 09 09 2a 72 65 | 73 20 3d 20 74 6d 70 3b |{.X..*re|s = tmp;|
|00000ea0| 0a 58 09 09 72 65 74 75 | 72 6e 3b 0a 58 09 7d 0a |.X..retu|rn;.X.}.|
|00000eb0| 58 09 7a 6d 6f 64 28 74 | 6d 70 2c 20 7a 33 2c 20 |X.zmod(t|mp, z3, |
|00000ec0| 72 65 73 29 3b 0a 58 09 | 66 72 65 65 68 28 74 6d |res);.X.|freeh(tm|
|00000ed0| 70 2e 76 29 3b 0a 58 7d | 0a 58 0a 58 0a 58 2f 2a |p.v);.X}|.X.X.X/*|
|00000ee0| 0a 58 20 2a 20 53 71 75 | 61 72 65 20 61 20 6e 75 |.X * Squ|are a nu|
|00000ef0| 6d 62 65 72 20 61 6e 64 | 20 74 68 65 6e 20 6d 6f |mber and| then mo|
|00000f00| 64 20 74 68 65 20 72 65 | 73 75 6c 74 20 77 69 74 |d the re|sult wit|
|00000f10| 68 20 61 20 73 65 63 6f | 6e 64 20 6e 75 6d 62 65 |h a seco|nd numbe|
|00000f20| 72 2e 0a 58 20 2a 20 54 | 68 65 20 6e 75 6d 62 65 |r..X * T|he numbe|
|00000f30| 72 20 74 6f 20 62 65 20 | 73 71 75 61 72 65 64 20 |r to be |squared |
|00000f40| 63 61 6e 20 62 65 20 6e | 65 67 61 74 69 76 65 20 |can be n|egative |
|00000f50| 6f 72 20 6f 75 74 20 6f | 66 20 6d 6f 64 75 6c 6f |or out o|f modulo|
|00000f60| 20 72 61 6e 67 65 2e 0a | 58 20 2a 20 54 68 65 20 | range..|X * The |
|00000f70| 72 65 73 75 6c 74 20 77 | 69 6c 6c 20 62 65 20 69 |result w|ill be i|
|00000f80| 6e 20 74 68 65 20 72 61 | 6e 67 65 20 30 20 74 6f |n the ra|nge 0 to|
|00000f90| 20 74 68 65 20 6d 6f 64 | 75 6c 75 73 20 2d 20 31 | the mod|ulus - 1|
|00000fa0| 2e 0a 58 20 2a 2f 0a 58 | 76 6f 69 64 0a 58 7a 73 |..X */.X|void.Xzs|
|00000fb0| 71 75 61 72 65 6d 6f 64 | 28 7a 31 2c 20 7a 32 2c |quaremod|(z1, z2,|
|00000fc0| 20 72 65 73 29 0a 58 09 | 5a 56 41 4c 55 45 20 7a | res).X.|ZVALUE z|
|00000fd0| 31 3b 09 09 2f 2a 20 6e | 75 6d 62 65 72 20 74 6f |1;../* n|umber to|
|00000fe0| 20 62 65 20 73 71 75 61 | 72 65 64 20 2a 2f 0a 58 | be squa|red */.X|
|00000ff0| 09 5a 56 41 4c 55 45 20 | 7a 32 3b 09 09 2f 2a 20 |.ZVALUE |z2;../* |
|00001000| 6e 75 6d 62 65 72 20 74 | 6f 20 74 61 6b 65 20 6d |number t|o take m|
|00001010| 6f 64 20 77 69 74 68 20 | 2a 2f 0a 58 09 5a 56 41 |od with |*/.X.ZVA|
|00001020| 4c 55 45 20 2a 72 65 73 | 3b 09 09 2f 2a 20 72 65 |LUE *res|;../* re|
|00001030| 73 75 6c 74 20 2a 2f 0a | 58 7b 0a 58 09 5a 56 41 |sult */.|X{.X.ZVA|
|00001040| 4c 55 45 20 74 6d 70 3b | 0a 58 09 46 55 4c 4c 20 |LUE tmp;|.X.FULL |
|00001050| 70 72 6f 64 3b 0a 58 09 | 46 55 4c 4c 20 64 69 67 |prod;.X.|FULL dig|
|00001060| 69 74 3b 0a 58 0a 58 09 | 69 66 20 28 69 73 7a 65 |it;.X.X.|if (isze|
|00001070| 72 6f 28 7a 32 29 20 7c | 7c 20 69 73 6e 65 67 28 |ro(z2) ||| isneg(|
|00001080| 7a 32 29 29 0a 58 09 09 | 65 72 72 6f 72 28 22 4d |z2)).X..|error("M|
|00001090| 6f 64 20 6f 66 20 6e 6f | 6e 2d 70 6f 73 69 74 69 |od of no|n-positi|
|000010a0| 76 65 20 69 6e 74 65 67 | 65 72 22 29 3b 0a 58 09 |ve integ|er");.X.|
|000010b0| 69 66 20 28 69 73 7a 65 | 72 6f 28 7a 31 29 20 7c |if (isze|ro(z1) ||
|000010c0| 7c 20 69 73 75 6e 69 74 | 28 7a 32 29 29 20 7b 0a || isunit|(z2)) {.|
|000010d0| 58 09 09 2a 72 65 73 20 | 3d 20 5f 7a 65 72 6f 5f |X..*res |= _zero_|
|000010e0| 3b 0a 58 09 09 72 65 74 | 75 72 6e 3b 0a 58 09 7d |;.X..ret|urn;.X.}|
|000010f0| 0a 58 0a 58 09 2f 2a 0a | 58 09 20 2a 20 49 66 20 |.X.X./*.|X. * If |
|00001100| 74 68 65 20 6d 6f 64 75 | 6c 75 73 20 69 73 20 61 |the modu|lus is a|
|00001110| 20 73 69 6e 67 6c 65 20 | 64 69 67 69 74 20 6e 75 | single |digit nu|
|00001120| 6d 62 65 72 2c 20 74 68 | 65 6e 20 64 6f 20 74 68 |mber, th|en do th|
|00001130| 65 20 72 65 73 75 6c 74 | 0a 58 09 20 2a 20 63 68 |e result|.X. * ch|
|00001140| 65 61 70 6c 79 2e 20 20 | 43 68 65 63 6b 20 65 73 |eaply.  |Check es|
|00001150| 70 65 63 69 61 6c 6c 79 | 20 66 6f 72 20 61 20 73 |pecially| for a s|
|00001160| 6d 61 6c 6c 20 70 6f 77 | 65 72 20 6f 66 20 74 77 |mall pow|er of tw|
|00001170| 6f 2e 0a 58 09 20 2a 2f | 0a 58 09 69 66 20 28 69 |o..X. */|.X.if (i|
|00001180| 73 74 69 6e 79 28 7a 32 | 29 29 20 7b 0a 58 09 09 |stiny(z2|)) {.X..|
|00001190| 64 69 67 69 74 20 3d 20 | 7a 32 2e 76 5b 30 5d 3b |digit = |z2.v[0];|
|000011a0| 0a 58 09 09 69 66 20 28 | 28 64 69 67 69 74 20 26 |.X..if (|(digit &|
|000011b0| 20 2d 64 69 67 69 74 29 | 20 3d 3d 20 64 69 67 69 | -digit)| == digi|
|000011c0| 74 29 20 7b 09 2f 2a 20 | 4e 45 45 44 53 20 32 27 |t) {./* |NEEDS 2'|
|000011d0| 53 20 43 4f 4d 50 20 2a | 2f 0a 58 09 09 09 70 72 |S COMP *|/.X...pr|
|000011e0| 6f 64 20 3d 20 28 46 55 | 4c 4c 29 20 7a 31 2e 76 |od = (FU|LL) z1.v|
|000011f0| 5b 30 5d 3b 0a 58 09 09 | 09 70 72 6f 64 20 3d 20 |[0];.X..|.prod = |
|00001200| 28 70 72 6f 64 20 2a 20 | 70 72 6f 64 29 20 26 20 |(prod * |prod) & |
|00001210| 28 64 69 67 69 74 20 2d | 20 31 29 3b 0a 58 09 09 |(digit -| 1);.X..|
|00001220| 7d 20 65 6c 73 65 20 7b | 0a 58 09 09 09 7a 31 2e |} else {|.X...z1.|
|00001230| 73 69 67 6e 20 3d 20 30 | 3b 0a 58 09 09 09 70 72 |sign = 0|;.X...pr|
|00001240| 6f 64 20 3d 20 28 46 55 | 4c 4c 29 20 7a 6d 6f 64 |od = (FU|LL) zmod|
|00001250| 69 28 7a 31 2c 20 28 6c | 6f 6e 67 29 20 64 69 67 |i(z1, (l|ong) dig|
|00001260| 69 74 29 3b 0a 58 09 09 | 09 70 72 6f 64 20 3d 20 |it);.X..|.prod = |
|00001270| 28 70 72 6f 64 20 2a 20 | 70 72 6f 64 29 20 25 20 |(prod * |prod) % |
|00001280| 64 69 67 69 74 3b 0a 58 | 09 09 7d 0a 58 09 09 69 |digit;.X|..}.X..i|
|00001290| 74 6f 7a 28 28 6c 6f 6e | 67 29 20 70 72 6f 64 2c |toz((lon|g) prod,|
|000012a0| 20 72 65 73 29 3b 0a 58 | 09 09 72 65 74 75 72 6e | res);.X|..return|
|000012b0| 3b 0a 58 09 7d 0a 58 0a | 58 09 2f 2a 0a 58 09 20 |;.X.}.X.|X./*.X. |
|000012c0| 2a 20 54 68 65 20 6d 6f | 64 75 6c 75 73 20 69 73 |* The mo|dulus is|
|000012d0| 20 6d 6f 72 65 20 74 68 | 61 6e 20 6f 6e 65 20 64 | more th|an one d|
|000012e0| 69 67 69 74 2e 0a 58 09 | 20 2a 20 41 63 74 75 61 |igit..X.| * Actua|
|000012f0| 6c 6c 79 20 64 6f 20 74 | 68 65 20 73 71 75 61 72 |lly do t|he squar|
|00001300| 65 20 61 6e 64 20 64 69 | 76 69 64 65 20 69 66 20 |e and di|vide if |
|00001310| 6e 65 63 65 73 73 61 72 | 79 2e 0a 58 09 20 2a 2f |necessar|y..X. */|
|00001320| 0a 58 09 7a 73 71 75 61 | 72 65 28 7a 31 2c 20 26 |.X.zsqua|re(z1, &|
|00001330| 74 6d 70 29 3b 0a 58 09 | 69 66 20 28 28 74 6d 70 |tmp);.X.|if ((tmp|
|00001340| 2e 6c 65 6e 20 3c 20 7a | 32 2e 6c 65 6e 29 20 7c |.len < z|2.len) ||
|00001350| 7c 0a 58 09 09 28 28 74 | 6d 70 2e 6c 65 6e 20 3d ||.X..((t|mp.len =|
|00001360| 3d 20 7a 32 2e 6c 65 6e | 29 20 26 26 20 28 74 6d |= z2.len|) && (tm|
|00001370| 70 2e 76 5b 74 6d 70 2e | 6c 65 6e 2d 31 5d 20 3c |p.v[tmp.|len-1] <|
|00001380| 20 7a 32 2e 76 5b 7a 32 | 2e 6c 65 6e 2d 31 5d 29 | z2.v[z2|.len-1])|
|00001390| 29 29 20 7b 0a 58 09 09 | 09 2a 72 65 73 20 3d 20 |)) {.X..|.*res = |
|000013a0| 74 6d 70 3b 0a 58 09 09 | 09 72 65 74 75 72 6e 3b |tmp;.X..|.return;|
|000013b0| 0a 58 09 7d 0a 58 09 7a | 6d 6f 64 28 74 6d 70 2c |.X.}.X.z|mod(tmp,|
|000013c0| 20 7a 32 2c 20 72 65 73 | 29 3b 0a 58 09 66 72 65 | z2, res|);.X.fre|
|000013d0| 65 68 28 74 6d 70 2e 76 | 29 3b 0a 58 7d 0a 58 0a |eh(tmp.v|);.X}.X.|
|000013e0| 58 0a 58 2f 2a 0a 58 20 | 2a 20 41 64 64 20 74 77 |X.X/*.X |* Add tw|
|000013f0| 6f 20 6e 75 6d 62 65 72 | 73 20 74 6f 67 65 74 68 |o number|s togeth|
|00001400| 65 72 20 61 6e 64 20 74 | 68 65 6e 20 6d 6f 64 20 |er and t|hen mod |
|00001410| 74 68 65 20 72 65 73 75 | 6c 74 20 77 69 74 68 20 |the resu|lt with |
|00001420| 61 20 74 68 69 72 64 20 | 6e 75 6d 62 65 72 2e 0a |a third |number..|
|00001430| 58 20 2a 20 54 68 65 20 | 74 77 6f 20 6e 75 6d 62 |X * The |two numb|
|00001440| 65 72 73 20 74 6f 20 62 | 65 20 61 64 64 65 64 20 |ers to b|e added |
|00001450| 63 61 6e 20 62 65 20 6e | 65 67 61 74 69 76 65 20 |can be n|egative |
|00001460| 6f 72 20 6f 75 74 20 6f | 66 20 6d 6f 64 75 6c 6f |or out o|f modulo|
|00001470| 20 72 61 6e 67 65 2e 0a | 58 20 2a 20 54 68 65 20 | range..|X * The |
|00001480| 72 65 73 75 6c 74 20 77 | 69 6c 6c 20 62 65 20 69 |result w|ill be i|
|00001490| 6e 20 74 68 65 20 72 61 | 6e 67 65 20 30 20 74 6f |n the ra|nge 0 to|
|000014a0| 20 74 68 65 20 6d 6f 64 | 75 6c 75 73 20 2d 20 31 | the mod|ulus - 1|
|000014b0| 2e 0a 58 20 2a 2f 0a 58 | 73 74 61 74 69 63 20 76 |..X */.X|static v|
|000014c0| 6f 69 64 0a 58 7a 61 64 | 64 6d 6f 64 28 7a 31 2c |oid.Xzad|dmod(z1,|
|000014d0| 20 7a 32 2c 20 7a 33 2c | 20 72 65 73 29 0a 58 09 | z2, z3,| res).X.|
|000014e0| 5a 56 41 4c 55 45 20 7a | 31 3b 09 09 2f 2a 20 66 |ZVALUE z|1;../* f|
|000014f0| 69 72 73 74 20 6e 75 6d | 62 65 72 20 74 6f 20 62 |irst num|ber to b|
|00001500| 65 20 61 64 64 65 64 20 | 2a 2f 0a 58 09 5a 56 41 |e added |*/.X.ZVA|
|00001510| 4c 55 45 20 7a 32 3b 09 | 09 2f 2a 20 73 65 63 6f |LUE z2;.|./* seco|
|00001520| 6e 64 20 6e 75 6d 62 65 | 72 20 74 6f 20 62 65 20 |nd numbe|r to be |
|00001530| 61 64 64 65 64 20 2a 2f | 0a 58 09 5a 56 41 4c 55 |added */|.X.ZVALU|
|00001540| 45 20 7a 33 3b 09 09 2f | 2a 20 6e 75 6d 62 65 72 |E z3;../|* number|
|00001550| 20 74 6f 20 74 61 6b 65 | 20 6d 6f 64 20 77 69 74 | to take| mod wit|
|00001560| 68 20 2a 2f 0a 58 09 5a | 56 41 4c 55 45 20 2a 72 |h */.X.Z|VALUE *r|
|00001570| 65 73 3b 09 09 2f 2a 20 | 72 65 73 75 6c 74 20 2a |es;../* |result *|
|00001580| 2f 0a 58 7b 0a 58 09 5a | 56 41 4c 55 45 20 74 6d |/.X{.X.Z|VALUE tm|
|00001590| 70 3b 0a 58 09 46 55 4c | 4c 20 73 75 6d 64 69 67 |p;.X.FUL|L sumdig|
|000015a0| 69 74 3b 0a 58 09 46 55 | 4c 4c 20 6d 6f 64 64 69 |it;.X.FU|LL moddi|
|000015b0| 67 69 74 3b 0a 58 0a 58 | 09 69 66 20 28 69 73 7a |git;.X.X|.if (isz|
|000015c0| 65 72 6f 28 7a 33 29 20 | 7c 7c 20 69 73 6e 65 67 |ero(z3) ||| isneg|
|000015d0| 28 7a 33 29 29 0a 58 09 | 09 65 72 72 6f 72 28 22 |(z3)).X.|.error("|
|000015e0| 4d 6f 64 20 6f 66 20 6e | 6f 6e 2d 70 6f 73 69 74 |Mod of n|on-posit|
|000015f0| 69 76 65 20 69 6e 74 65 | 67 65 72 22 29 3b 0a 58 |ive inte|ger");.X|
|00001600| 09 69 66 20 28 28 69 73 | 7a 65 72 6f 28 7a 31 29 |.if ((is|zero(z1)|
|00001610| 20 26 26 20 69 73 7a 65 | 72 6f 28 7a 32 29 29 20 | && isze|ro(z2)) |
|00001620| 7c 7c 20 69 73 75 6e 69 | 74 28 7a 33 29 29 20 7b ||| isuni|t(z3)) {|
|00001630| 0a 58 09 09 2a 72 65 73 | 20 3d 20 5f 7a 65 72 6f |.X..*res| = _zero|
|00001640| 5f 3b 0a 58 09 09 72 65 | 74 75 72 6e 3b 0a 58 09 |_;.X..re|turn;.X.|
|00001650| 7d 0a 58 09 69 66 20 28 | 69 73 74 77 6f 28 7a 32 |}.X.if (|istwo(z2|
|00001660| 29 29 20 7b 0a 58 09 09 | 69 66 20 28 28 7a 31 2e |)) {.X..|if ((z1.|
|00001670| 76 5b 30 5d 20 2b 20 7a | 32 2e 76 5b 30 5d 29 20 |v[0] + z|2.v[0]) |
|00001680| 26 20 30 78 31 29 0a 58 | 09 09 09 2a 72 65 73 20 |& 0x1).X|...*res |
|00001690| 3d 20 5f 6f 6e 65 5f 3b | 0a 58 09 09 65 6c 73 65 |= _one_;|.X..else|
|000016a0| 0a 58 09 09 09 2a 72 65 | 73 20 3d 20 5f 7a 65 72 |.X...*re|s = _zer|
|000016b0| 6f 5f 3b 0a 58 09 09 72 | 65 74 75 72 6e 3b 0a 58 |o_;.X..r|eturn;.X|
|000016c0| 09 7d 0a 58 09 7a 61 64 | 64 28 7a 31 2c 20 7a 32 |.}.X.zad|d(z1, z2|
|000016d0| 2c 20 26 74 6d 70 29 3b | 0a 58 09 69 66 20 28 69 |, &tmp);|.X.if (i|
|000016e0| 73 6e 65 67 28 74 6d 70 | 29 20 7c 7c 20 28 74 6d |sneg(tmp|) || (tm|
|000016f0| 70 2e 6c 65 6e 20 3e 20 | 7a 33 2e 6c 65 6e 29 29 |p.len > |z3.len))|
|00001700| 20 7b 0a 58 09 09 7a 6d | 6f 64 28 74 6d 70 2c 20 | {.X..zm|od(tmp, |
|00001710| 7a 33 2c 20 72 65 73 29 | 3b 0a 58 09 09 66 72 65 |z3, res)|;.X..fre|
|00001720| 65 68 28 74 6d 70 2e 76 | 29 3b 0a 58 09 09 72 65 |eh(tmp.v|);.X..re|
|00001730| 74 75 72 6e 3b 0a 58 09 | 7d 0a 58 09 73 75 6d 64 |turn;.X.|}.X.sumd|
|00001740| 69 67 69 74 20 3d 20 74 | 6d 70 2e 76 5b 74 6d 70 |igit = t|mp.v[tmp|
|00001750| 2e 6c 65 6e 20 2d 20 31 | 5d 3b 0a 58 09 6d 6f 64 |.len - 1|];.X.mod|
|00001760| 64 69 67 69 74 20 3d 20 | 7a 33 2e 76 5b 7a 33 2e |digit = |z3.v[z3.|
|00001770| 6c 65 6e 20 2d 20 31 5d | 3b 0a 58 09 69 66 20 28 |len - 1]|;.X.if (|
|00001780| 28 74 6d 70 2e 6c 65 6e | 20 3c 20 7a 33 2e 6c 65 |(tmp.len| < z3.le|
|00001790| 6e 29 20 7c 7c 20 28 73 | 75 6d 64 69 67 69 74 20 |n) || (s|umdigit |
|000017a0| 3c 20 6d 6f 64 64 69 67 | 69 74 29 29 20 7b 0a 58 |< moddig|it)) {.X|
|000017b0| 09 09 2a 72 65 73 20 3d | 20 74 6d 70 3b 0a 58 09 |..*res =| tmp;.X.|
|000017c0| 09 72 65 74 75 72 6e 3b | 0a 58 09 7d 0a 58 09 69 |.return;|.X.}.X.i|
|000017d0| 66 20 28 73 75 6d 64 69 | 67 69 74 20 3c 20 32 20 |f (sumdi|git < 2 |
|000017e0| 2a 20 6d 6f 64 64 69 67 | 69 74 29 20 7b 0a 58 09 |* moddig|it) {.X.|
|000017f0| 09 7a 73 75 62 28 74 6d | 70 2c 20 7a 33 2c 20 72 |.zsub(tm|p, z3, r|
|00001800| 65 73 29 3b 0a 58 09 09 | 66 72 65 65 68 28 74 6d |es);.X..|freeh(tm|
|00001810| 70 2e 76 29 3b 0a 58 09 | 09 72 65 74 75 72 6e 3b |p.v);.X.|.return;|
|00001820| 0a 58 09 7d 0a 58 09 7a | 6d 6f 64 28 74 6d 70 2c |.X.}.X.z|mod(tmp,|
|00001830| 20 7a 32 2c 20 72 65 73 | 29 3b 0a 58 09 66 72 65 | z2, res|);.X.fre|
|00001840| 65 68 28 74 6d 70 2e 76 | 29 3b 0a 58 7d 0a 58 0a |eh(tmp.v|);.X}.X.|
|00001850| 58 0a 58 2f 2a 0a 58 20 | 2a 20 53 75 62 74 72 61 |X.X/*.X |* Subtra|
|00001860| 63 74 20 74 77 6f 20 6e | 75 6d 62 65 72 73 20 74 |ct two n|umbers t|
|00001870| 6f 67 65 74 68 65 72 20 | 61 6e 64 20 74 68 65 6e |ogether |and then|
|00001880| 20 6d 6f 64 20 74 68 65 | 20 72 65 73 75 6c 74 20 | mod the| result |
|00001890| 77 69 74 68 20 61 20 74 | 68 69 72 64 20 6e 75 6d |with a t|hird num|
|000018a0| 62 65 72 2e 0a 58 20 2a | 20 54 68 65 20 74 77 6f |ber..X *| The two|
|000018b0| 20 6e 75 6d 62 65 72 73 | 20 74 6f 20 62 65 20 73 | numbers| to be s|
|000018c0| 75 62 74 72 61 63 74 20 | 63 61 6e 20 62 65 20 6e |ubtract |can be n|
|000018d0| 65 67 61 74 69 76 65 20 | 6f 72 20 6f 75 74 20 6f |egative |or out o|
|000018e0| 66 20 6d 6f 64 75 6c 6f | 20 72 61 6e 67 65 2e 0a |f modulo| range..|
|000018f0| 58 20 2a 20 54 68 65 20 | 72 65 73 75 6c 74 20 77 |X * The |result w|
|00001900| 69 6c 6c 20 62 65 20 69 | 6e 20 74 68 65 20 72 61 |ill be i|n the ra|
|00001910| 6e 67 65 20 30 20 74 6f | 20 74 68 65 20 6d 6f 64 |nge 0 to| the mod|
|00001920| 75 6c 75 73 20 2d 20 31 | 2e 0a 58 20 2a 2f 0a 58 |ulus - 1|..X */.X|
|00001930| 76 6f 69 64 0a 58 7a 73 | 75 62 6d 6f 64 28 7a 31 |void.Xzs|ubmod(z1|
|00001940| 2c 20 7a 32 2c 20 7a 33 | 2c 20 72 65 73 29 0a 58 |, z2, z3|, res).X|
|00001950| 09 5a 56 41 4c 55 45 20 | 7a 31 3b 09 09 2f 2a 20 |.ZVALUE |z1;../* |
|00001960| 6e 75 6d 62 65 72 20 74 | 6f 20 62 65 20 73 75 62 |number t|o be sub|
|00001970| 74 72 61 63 74 65 64 20 | 66 72 6f 6d 20 2a 2f 0a |tracted |from */.|
|00001980| 58 09 5a 56 41 4c 55 45 | 20 7a 32 3b 09 09 2f 2a |X.ZVALUE| z2;../*|
|00001990| 20 6e 75 6d 62 65 72 20 | 74 6f 20 62 65 20 73 75 | number |to be su|
|000019a0| 62 74 72 61 63 74 65 64 | 20 2a 2f 0a 58 09 5a 56 |btracted| */.X.ZV|
|000019b0| 41 4c 55 45 20 7a 33 3b | 09 09 2f 2a 20 6e 75 6d |ALUE z3;|../* num|
|000019c0| 62 65 72 20 74 6f 20 74 | 61 6b 65 20 6d 6f 64 20 |ber to t|ake mod |
|000019d0| 77 69 74 68 20 2a 2f 0a | 58 09 5a 56 41 4c 55 45 |with */.|X.ZVALUE|
|000019e0| 20 2a 72 65 73 3b 09 09 | 2f 2a 20 72 65 73 75 6c | *res;..|/* resul|
|000019f0| 74 20 2a 2f 0a 58 7b 0a | 58 09 69 66 20 28 69 73 |t */.X{.|X.if (is|
|00001a00| 7a 65 72 6f 28 7a 33 29 | 20 7c 7c 20 69 73 6e 65 |zero(z3)| || isne|
|00001a10| 67 28 7a 33 29 29 0a 58 | 09 09 65 72 72 6f 72 28 |g(z3)).X|..error(|
|00001a20| 22 4d 6f 64 20 6f 66 20 | 6e 6f 6e 2d 70 6f 73 69 |"Mod of |non-posi|
|00001a30| 74 69 76 65 20 69 6e 74 | 65 67 65 72 22 29 3b 0a |tive int|eger");.|
|00001a40| 58 09 69 66 20 28 69 73 | 7a 65 72 6f 28 7a 32 29 |X.if (is|zero(z2)|
|00001a50| 29 20 7b 0a 58 09 09 7a | 6d 6f 64 28 7a 31 2c 20 |) {.X..z|mod(z1, |
|00001a60| 7a 33 2c 20 72 65 73 29 | 3b 0a 58 09 09 72 65 74 |z3, res)|;.X..ret|
|00001a70| 75 72 6e 3b 0a 58 09 7d | 0a 58 09 69 66 20 28 69 |urn;.X.}|.X.if (i|
|00001a80| 73 7a 65 72 6f 28 7a 31 | 29 29 20 7b 0a 58 09 09 |szero(z1|)) {.X..|
|00001a90| 7a 6e 65 67 6d 6f 64 28 | 7a 32 2c 20 7a 33 2c 20 |znegmod(|z2, z3, |
|00001aa0| 72 65 73 29 3b 0a 58 09 | 09 72 65 74 75 72 6e 3b |res);.X.|.return;|
|00001ab0| 0a 58 09 7d 0a 58 09 69 | 66 20 28 28 7a 31 2e 73 |.X.}.X.i|f ((z1.s|
|00001ac0| 69 67 6e 20 3d 3d 20 7a | 32 2e 73 69 67 6e 29 20 |ign == z|2.sign) |
|00001ad0| 26 26 20 28 7a 31 2e 6c | 65 6e 20 3d 3d 20 7a 32 |&& (z1.l|en == z2|
|00001ae0| 2e 6c 65 6e 29 20 26 26 | 0a 58 09 09 28 7a 31 2e |.len) &&|.X..(z1.|
|00001af0| 76 5b 30 5d 20 3d 3d 20 | 7a 32 2e 76 5b 30 5d 29 |v[0] == |z2.v[0])|
|00001b00| 20 26 26 20 28 7a 63 6d | 70 28 7a 31 2c 20 7a 32 | && (zcm|p(z1, z2|
|00001b10| 29 20 3d 3d 20 30 29 29 | 20 7b 0a 58 09 09 09 2a |) == 0))| {.X...*|
|00001b20| 72 65 73 20 3d 20 5f 7a | 65 72 6f 5f 3b 0a 58 09 |res = _z|ero_;.X.|
|00001b30| 09 09 72 65 74 75 72 6e | 3b 0a 58 09 7d 0a 58 09 |..return|;.X.}.X.|
|00001b40| 7a 32 2e 73 69 67 6e 20 | 3d 20 21 7a 32 2e 73 69 |z2.sign |= !z2.si|
|00001b50| 67 6e 3b 0a 58 09 7a 61 | 64 64 6d 6f 64 28 7a 31 |gn;.X.za|ddmod(z1|
|00001b60| 2c 20 7a 32 2c 20 7a 33 | 2c 20 72 65 73 29 3b 0a |, z2, z3|, res);.|
|00001b70| 58 7d 0a 58 0a 58 0a 58 | 2f 2a 0a 58 20 2a 20 43 |X}.X.X.X|/*.X * C|
|00001b80| 61 6c 63 75 6c 61 74 65 | 20 74 68 65 20 6e 65 67 |alculate| the neg|
|00001b90| 61 74 69 76 65 20 6f 66 | 20 61 20 6e 75 6d 62 65 |ative of| a numbe|
|00001ba0| 72 20 6d 6f 64 75 6c 6f | 20 61 6e 6f 74 68 65 72 |r modulo| another|
|00001bb0| 20 6e 75 6d 62 65 72 2e | 0a 58 20 2a 20 54 68 65 | number.|.X * The|
|00001bc0| 20 6e 75 6d 62 65 72 20 | 74 6f 20 62 65 20 6e 65 | number |to be ne|
|00001bd0| 67 61 74 65 64 20 63 61 | 6e 20 62 65 20 6e 65 67 |gated ca|n be neg|
|00001be0| 61 74 69 76 65 20 6f 72 | 20 6f 75 74 20 6f 66 20 |ative or| out of |
|00001bf0| 6d 6f 64 75 6c 6f 20 72 | 61 6e 67 65 2e 0a 58 20 |modulo r|ange..X |
|00001c00| 2a 20 54 68 65 20 72 65 | 73 75 6c 74 20 77 69 6c |* The re|sult wil|
|00001c10| 6c 20 62 65 20 69 6e 20 | 74 68 65 20 72 61 6e 67 |l be in |the rang|
|00001c20| 65 20 30 20 74 6f 20 74 | 68 65 20 6d 6f 64 75 6c |e 0 to t|he modul|
|00001c30| 75 73 20 2d 20 31 2e 0a | 58 20 2a 2f 0a 58 73 74 |us - 1..|X */.Xst|
|00001c40| 61 74 69 63 20 76 6f 69 | 64 0a 58 7a 6e 65 67 6d |atic voi|d.Xznegm|
|00001c50| 6f 64 28 7a 31 2c 20 7a | 32 2c 20 72 65 73 29 0a |od(z1, z|2, res).|
|00001c60| 58 09 5a 56 41 4c 55 45 | 20 7a 31 3b 09 09 2f 2a |X.ZVALUE| z1;../*|
|00001c70| 20 6e 75 6d 62 65 72 20 | 74 6f 20 74 61 6b 65 20 | number |to take |
|00001c80| 6e 65 67 61 74 69 76 65 | 20 6f 66 20 2a 2f 0a 58 |negative| of */.X|
|00001c90| 09 5a 56 41 4c 55 45 20 | 7a 32 3b 09 09 2f 2a 20 |.ZVALUE |z2;../* |
|00001ca0| 6e 75 6d 62 65 72 20 74 | 6f 20 74 61 6b 65 20 6d |number t|o take m|
|00001cb0| 6f 64 20 77 69 74 68 20 | 2a 2f 0a 58 09 5a 56 41 |od with |*/.X.ZVA|
|00001cc0| 4c 55 45 20 2a 72 65 73 | 3b 09 09 2f 2a 20 72 65 |LUE *res|;../* re|
|00001cd0| 73 75 6c 74 20 2a 2f 0a | 58 7b 0a 58 09 69 6e 74 |sult */.|X{.X.int|
|00001ce0| 20 73 69 67 6e 3b 0a 58 | 09 69 6e 74 20 63 76 3b | sign;.X|.int cv;|
|00001cf0| 0a 58 0a 58 09 69 66 20 | 28 69 73 7a 65 72 6f 28 |.X.X.if |(iszero(|
|00001d00| 7a 32 29 20 7c 7c 20 69 | 73 6e 65 67 28 7a 32 29 |z2) || i|sneg(z2)|
|00001d10| 29 0a 58 09 09 65 72 72 | 6f 72 28 22 4d 6f 64 20 |).X..err|or("Mod |
|00001d20| 6f 66 20 6e 6f 6e 2d 70 | 6f 73 69 74 69 76 65 20 |of non-p|ositive |
|00001d30| 69 6e 74 65 67 65 72 22 | 29 3b 0a 58 09 69 66 20 |integer"|);.X.if |
|00001d40| 28 69 73 7a 65 72 6f 28 | 7a 31 29 20 7c 7c 20 69 |(iszero(|z1) || i|
|00001d50| 73 75 6e 69 74 28 7a 32 | 29 29 20 7b 0a 58 09 09 |sunit(z2|)) {.X..|
|00001d60| 2a 72 65 73 20 3d 20 5f | 7a 65 72 6f 5f 3b 0a 58 |*res = _|zero_;.X|
|00001d70| 09 09 72 65 74 75 72 6e | 3b 0a 58 09 7d 0a 58 09 |..return|;.X.}.X.|
|00001d80| 69 66 20 28 69 73 74 77 | 6f 28 7a 32 29 29 20 7b |if (istw|o(z2)) {|
|00001d90| 0a 58 09 09 69 66 20 28 | 7a 31 2e 76 5b 30 5d 20 |.X..if (|z1.v[0] |
|00001da0| 26 20 30 78 31 29 0a 58 | 09 09 09 2a 72 65 73 20 |& 0x1).X|...*res |
|00001db0| 3d 20 5f 6f 6e 65 5f 3b | 0a 58 09 09 65 6c 73 65 |= _one_;|.X..else|
|00001dc0| 0a 58 09 09 09 2a 72 65 | 73 20 3d 20 5f 7a 65 72 |.X...*re|s = _zer|
|00001dd0| 6f 5f 3b 0a 58 09 09 72 | 65 74 75 72 6e 3b 0a 58 |o_;.X..r|eturn;.X|
|00001de0| 09 7d 0a 58 0a 58 09 2f | 2a 0a 58 09 20 2a 20 49 |.}.X.X./|*.X. * I|
|00001df0| 66 20 74 68 65 20 61 62 | 73 6f 6c 75 74 65 20 76 |f the ab|solute v|
|00001e00| 61 6c 75 65 20 6f 66 20 | 74 68 65 20 6e 75 6d 62 |alue of |the numb|
|00001e10| 65 72 20 69 73 20 77 69 | 74 68 69 6e 20 74 68 65 |er is wi|thin the|
|00001e20| 20 6d 6f 64 75 6c 6f 20 | 72 61 6e 67 65 2c 0a 58 | modulo |range,.X|
|00001e30| 09 20 2a 20 74 68 65 6e | 20 74 68 65 20 72 65 73 |. * then| the res|
|00001e40| 75 6c 74 20 69 73 20 6a | 75 73 74 20 61 20 63 6f |ult is j|ust a co|
|00001e50| 70 79 20 6f 72 20 61 20 | 73 75 62 74 72 61 63 74 |py or a |subtract|
|00001e60| 69 6f 6e 2e 20 20 4f 74 | 68 65 72 77 69 73 65 20 |ion.  Ot|herwise |
|00001e70| 67 6f 0a 58 09 20 2a 20 | 61 68 65 61 64 20 61 6e |go.X. * |ahead an|
|00001e80| 64 20 6e 65 67 61 74 65 | 20 61 6e 64 20 72 65 64 |d negate| and red|
|00001e90| 75 63 65 20 74 68 65 20 | 72 65 73 75 6c 74 2e 0a |uce the |result..|
|00001ea0| 58 09 20 2a 2f 0a 58 09 | 73 69 67 6e 20 3d 20 7a |X. */.X.|sign = z|
|00001eb0| 31 2e 73 69 67 6e 3b 0a | 58 09 7a 31 2e 73 69 67 |1.sign;.|X.z1.sig|
|00001ec0| 6e 20 3d 20 30 3b 0a 58 | 09 63 76 20 3d 20 7a 72 |n = 0;.X|.cv = zr|
|00001ed0| 65 6c 28 7a 31 2c 20 7a | 32 29 3b 0a 58 09 69 66 |el(z1, z|2);.X.if|
|00001ee0| 20 28 63 76 20 3d 3d 20 | 30 29 20 7b 0a 58 09 09 | (cv == |0) {.X..|
|00001ef0| 2a 72 65 73 20 3d 20 5f | 7a 65 72 6f 5f 3b 0a 58 |*res = _|zero_;.X|
|00001f00| 09 09 72 65 74 75 72 6e | 3b 0a 58 09 7d 0a 58 09 |..return|;.X.}.X.|
|00001f10| 69 66 20 28 63 76 20 3c | 20 30 29 20 7b 0a 58 09 |if (cv <| 0) {.X.|
|00001f20| 09 69 66 20 28 73 69 67 | 6e 29 0a 58 09 09 09 7a |.if (sig|n).X...z|
|00001f30| 63 6f 70 79 28 7a 31 2c | 20 72 65 73 29 3b 0a 58 |copy(z1,| res);.X|
|00001f40| 09 09 65 6c 73 65 0a 58 | 09 09 09 7a 73 75 62 28 |..else.X|...zsub(|
|00001f50| 7a 32 2c 20 7a 31 2c 20 | 72 65 73 29 3b 0a 58 09 |z2, z1, |res);.X.|
|00001f60| 09 72 65 74 75 72 6e 3b | 0a 58 09 7d 0a 58 09 7a |.return;|.X.}.X.z|
|00001f70| 31 2e 73 69 67 6e 20 3d | 20 21 73 69 67 6e 3b 0a |1.sign =| !sign;.|
|00001f80| 58 09 7a 6d 6f 64 28 7a | 31 2c 20 7a 32 2c 20 72 |X.zmod(z|1, z2, r|
|00001f90| 65 73 29 3b 0a 58 7d 0a | 58 23 65 6e 64 69 66 0a |es);.X}.|X#endif.|
|00001fa0| 58 0a 58 0a 58 2f 2a 0a | 58 20 2a 20 43 61 6c 63 |X.X.X/*.|X * Calc|
|00001fb0| 75 6c 61 74 65 20 74 68 | 65 20 6e 75 6d 62 65 72 |ulate th|e number|
|00001fc0| 20 63 6f 6e 67 72 75 65 | 6e 74 20 74 6f 20 74 68 | congrue|nt to th|
|00001fd0| 65 20 67 69 76 65 6e 20 | 6e 75 6d 62 65 72 20 77 |e given |number w|
|00001fe0| 68 6f 73 65 20 61 62 73 | 6f 6c 75 74 65 0a 58 20 |hose abs|olute.X |
|00001ff0| 2a 20 76 61 6c 75 65 20 | 69 73 20 6d 69 6e 69 6d |* value |is minim|
|00002000| 61 6c 2e 20 20 54 68 65 | 20 6e 75 6d 62 65 72 20 |al.  The| number |
|00002010| 74 6f 20 62 65 20 72 65 | 64 75 63 65 64 20 63 61 |to be re|duced ca|
|00002020| 6e 20 62 65 20 6e 65 67 | 61 74 69 76 65 20 6f 72 |n be neg|ative or|
|00002030| 20 6f 75 74 20 6f 66 0a | 58 20 2a 20 6d 6f 64 75 | out of.|X * modu|
|00002040| 6c 6f 20 72 61 6e 67 65 | 2e 20 20 54 68 65 20 72 |lo range|.  The r|
|00002050| 65 73 75 6c 74 20 77 69 | 6c 6c 20 62 65 20 77 69 |esult wi|ll be wi|
|00002060| 74 68 69 6e 20 74 68 65 | 20 72 61 6e 67 65 20 2d |thin the| range -|
|00002070| 69 6e 74 28 28 6d 6f 64 | 75 6c 75 73 2d 31 29 2f |int((mod|ulus-1)/|
|00002080| 32 29 0a 58 20 2a 20 74 | 6f 20 69 6e 74 28 6d 6f |2).X * t|o int(mo|
|00002090| 64 75 6c 75 73 2f 32 29 | 20 69 6e 63 6c 75 73 69 |dulus/2)| inclusi|
|000020a0| 76 65 2e 20 20 46 6f 72 | 20 65 78 61 6d 70 6c 65 |ve.  For| example|
|000020b0| 2c 20 66 6f 72 20 6d 6f | 64 75 6c 75 73 20 37 2c |, for mo|dulus 7,|
|000020c0| 20 6e 75 6d 62 65 72 73 | 20 61 72 65 0a 58 20 2a | numbers| are.X *|
|000020d0| 20 72 65 64 75 63 65 64 | 20 74 6f 20 74 68 65 20 | reduced| to the |
|000020e0| 72 61 6e 67 65 20 5b 2d | 33 2c 20 33 5d 2c 20 61 |range [-|3, 3], a|
|000020f0| 6e 64 20 66 6f 72 20 6d | 6f 64 75 6c 75 73 20 38 |nd for m|odulus 8|
|00002100| 2c 20 6e 75 6d 62 65 72 | 73 20 61 72 65 20 72 65 |, number|s are re|
|00002110| 64 75 63 65 64 20 74 6f | 0a 58 20 2a 20 74 68 65 |duced to|.X * the|
|00002120| 20 72 61 6e 67 65 20 5b | 2d 33 2c 20 34 5d 2e 0a | range [|-3, 4]..|
|00002130| 58 20 2a 2f 0a 58 76 6f | 69 64 0a 58 7a 6d 69 6e |X */.Xvo|id.Xzmin|
|00002140| 6d 6f 64 28 7a 31 2c 20 | 7a 32 2c 20 72 65 73 29 |mod(z1, |z2, res)|
|00002150| 0a 58 09 5a 56 41 4c 55 | 45 20 7a 31 3b 09 09 2f |.X.ZVALU|E z1;../|
|00002160| 2a 20 6e 75 6d 62 65 72 | 20 74 6f 20 66 69 6e 64 |* number| to find|
|00002170| 20 6d 69 6e 69 6d 75 6d | 20 63 6f 6e 67 72 75 65 | minimum| congrue|
|00002180| 6e 63 65 20 6f 66 20 2a | 2f 0a 58 09 5a 56 41 4c |nce of *|/.X.ZVAL|
|00002190| 55 45 20 7a 32 3b 09 09 | 2f 2a 20 6e 75 6d 62 65 |UE z2;..|/* numbe|
|000021a0| 72 20 74 6f 20 74 61 6b | 65 20 6d 6f 64 20 77 69 |r to tak|e mod wi|
|000021b0| 74 68 20 2a 2f 0a 58 09 | 5a 56 41 4c 55 45 20 2a |th */.X.|ZVALUE *|
|000021c0| 72 65 73 3b 09 09 2f 2a | 20 72 65 73 75 6c 74 20 |res;../*| result |
|000021d0| 2a 2f 0a 58 7b 0a 58 09 | 5a 56 41 4c 55 45 20 74 |*/.X{.X.|ZVALUE t|
|000021e0| 6d 70 31 2c 20 74 6d 70 | 32 3b 0a 58 09 69 6e 74 |mp1, tmp|2;.X.int|
|000021f0| 20 73 69 67 6e 3b 0a 58 | 09 69 6e 74 20 63 76 3b | sign;.X|.int cv;|
|00002200| 0a 58 0a 58 09 69 66 20 | 28 69 73 7a 65 72 6f 28 |.X.X.if |(iszero(|
|00002210| 7a 32 29 20 7c 7c 20 69 | 73 6e 65 67 28 7a 32 29 |z2) || i|sneg(z2)|
|00002220| 29 0a 58 09 09 65 72 72 | 6f 72 28 22 4d 6f 64 20 |).X..err|or("Mod |
|00002230| 6f 66 20 6e 6f 6e 2d 70 | 6f 73 69 74 69 76 65 20 |of non-p|ositive |
|00002240| 69 6e 74 65 67 65 72 22 | 29 3b 0a 58 09 69 66 20 |integer"|);.X.if |
|00002250| 28 69 73 7a 65 72 6f 28 | 7a 31 29 20 7c 7c 20 69 |(iszero(|z1) || i|
|00002260| 73 75 6e 69 74 28 7a 32 | 29 29 20 7b 0a 58 09 09 |sunit(z2|)) {.X..|
|00002270| 2a 72 65 73 20 3d 20 5f | 7a 65 72 6f 5f 3b 0a 58 |*res = _|zero_;.X|
|00002280| 09 09 72 65 74 75 72 6e | 3b 0a 58 09 7d 0a 58 09 |..return|;.X.}.X.|
|00002290| 69 66 20 28 69 73 74 77 | 6f 28 7a 32 29 29 20 7b |if (istw|o(z2)) {|
|000022a0| 0a 58 09 09 69 66 20 28 | 69 73 6f 64 64 28 7a 31 |.X..if (|isodd(z1|
|000022b0| 29 29 0a 58 09 09 09 2a | 72 65 73 20 3d 20 5f 6f |)).X...*|res = _o|
|000022c0| 6e 65 5f 3b 0a 58 09 09 | 65 6c 73 65 0a 58 09 09 |ne_;.X..|else.X..|
|000022d0| 09 2a 72 65 73 20 3d 20 | 5f 7a 65 72 6f 5f 3b 0a |.*res = |_zero_;.|
|000022e0| 58 09 09 72 65 74 75 72 | 6e 3b 0a 58 09 7d 0a 58 |X..retur|n;.X.}.X|
|000022f0| 0a 58 09 2f 2a 0a 58 09 | 20 2a 20 44 6f 20 61 20 |.X./*.X.| * Do a |
|00002300| 71 75 69 63 6b 20 63 68 | 65 63 6b 20 74 6f 20 73 |quick ch|eck to s|
|00002310| 65 65 20 69 66 20 74 68 | 65 20 6e 75 6d 62 65 72 |ee if th|e number|
|00002320| 20 69 73 20 76 65 72 79 | 20 73 6d 61 6c 6c 20 63 | is very| small c|
|00002330| 6f 6d 70 61 72 65 64 0a | 58 09 20 2a 20 74 6f 20 |ompared.|X. * to |
|00002340| 74 68 65 20 6d 6f 64 75 | 6c 75 73 2e 20 20 49 66 |the modu|lus.  If|
|00002350| 20 73 6f 2c 20 74 68 65 | 6e 20 74 68 65 20 72 65 | so, the|n the re|
|00002360| 73 75 6c 74 20 69 73 20 | 6f 62 76 69 6f 75 73 2e |sult is |obvious.|
|00002370| 0a 58 09 20 2a 2f 0a 58 | 09 69 66 20 28 7a 31 2e |.X. */.X|.if (z1.|
|00002380| 6c 65 6e 20 3c 20 7a 32 | 2e 6c 65 6e 20 2d 20 31 |len < z2|.len - 1|
|00002390| 29 20 7b 0a 58 09 09 7a | 63 6f 70 79 28 7a 31 2c |) {.X..z|copy(z1,|
|000023a0| 20 72 65 73 29 3b 0a 58 | 09 09 72 65 74 75 72 6e | res);.X|..return|
|000023b0| 3b 0a 58 09 7d 0a 58 0a | 58 09 2f 2a 0a 58 09 20 |;.X.}.X.|X./*.X. |
|000023c0| 2a 20 4e 6f 77 20 6d 61 | 6b 65 20 73 75 72 65 20 |* Now ma|ke sure |
|000023d0| 74 68 65 20 69 6e 70 75 | 74 20 6e 75 6d 62 65 72 |the inpu|t number|
|000023e0| 20 69 73 20 77 69 74 68 | 69 6e 20 74 68 65 20 6d | is with|in the m|
|000023f0| 6f 64 75 6c 6f 20 72 61 | 6e 67 65 2e 0a 58 09 20 |odulo ra|nge..X. |
|00002400| 2a 20 49 66 20 6e 6f 74 | 2c 20 74 68 65 6e 20 72 |* If not|, then r|
|00002410| 65 64 75 63 65 20 69 74 | 20 74 6f 20 62 65 20 77 |educe it| to be w|
|00002420| 69 74 68 69 6e 20 72 61 | 6e 67 65 20 61 6e 64 20 |ithin ra|nge and |
|00002430| 6d 61 6b 65 20 74 68 65 | 0a 58 09 20 2a 20 71 75 |make the|.X. * qu|
|00002440| 69 63 6b 20 63 68 65 63 | 6b 20 61 67 61 69 6e 2e |ick chec|k again.|
|00002450| 0a 58 09 20 2a 2f 0a 58 | 09 73 69 67 6e 20 3d 20 |.X. */.X|.sign = |
|00002460| 7a 31 2e 73 69 67 6e 3b | 0a 58 09 7a 31 2e 73 69 |z1.sign;|.X.z1.si|
|00002470| 67 6e 20 3d 20 30 3b 0a | 58 09 63 76 20 3d 20 7a |gn = 0;.|X.cv = z|
|00002480| 72 65 6c 28 7a 31 2c 20 | 7a 32 29 3b 0a 58 09 69 |rel(z1, |z2);.X.i|
|00002490| 66 20 28 63 76 20 3d 3d | 20 30 29 20 7b 0a 58 09 |f (cv ==| 0) {.X.|
|000024a0| 09 2a 72 65 73 20 3d 20 | 5f 7a 65 72 6f 5f 3b 0a |.*res = |_zero_;.|
|000024b0| 58 09 09 72 65 74 75 72 | 6e 3b 0a 58 09 7d 0a 58 |X..retur|n;.X.}.X|
|000024c0| 09 74 6d 70 31 20 3d 20 | 7a 31 3b 0a 58 09 69 66 |.tmp1 = |z1;.X.if|
|000024d0| 20 28 63 76 20 3e 20 30 | 29 20 7b 0a 58 09 09 7a | (cv > 0|) {.X..z|
|000024e0| 31 2e 73 69 67 6e 20 3d | 20 28 42 4f 4f 4c 29 73 |1.sign =| (BOOL)s|
|000024f0| 69 67 6e 3b 0a 58 09 09 | 7a 6d 6f 64 28 7a 31 2c |ign;.X..|zmod(z1,|
|00002500| 20 7a 32 2c 20 26 74 6d | 70 31 29 3b 0a 58 09 09 | z2, &tm|p1);.X..|
|00002510| 69 66 20 28 74 6d 70 31 | 2e 6c 65 6e 20 3c 20 7a |if (tmp1|.len < z|
|00002520| 32 2e 6c 65 6e 20 2d 20 | 31 29 20 7b 0a 58 09 09 |2.len - |1) {.X..|
|00002530| 09 2a 72 65 73 20 3d 20 | 74 6d 70 31 3b 0a 58 09 |.*res = |tmp1;.X.|
|00002540| 09 09 72 65 74 75 72 6e | 3b 0a 58 09 09 7d 0a 58 |..return|;.X..}.X|
|00002550| 09 09 73 69 67 6e 20 3d | 20 30 3b 0a 58 09 7d 0a |..sign =| 0;.X.}.|
|00002560| 58 0a 58 09 2f 2a 0a 58 | 09 20 2a 20 4e 6f 77 20 |X.X./*.X|. * Now |
|00002570| 63 61 6c 63 75 6c 61 74 | 65 20 74 68 65 20 64 69 |calculat|e the di|
|00002580| 66 66 65 72 65 6e 63 65 | 20 6f 66 20 74 68 65 20 |fference| of the |
|00002590| 6d 6f 64 75 6c 75 73 20 | 61 6e 64 20 74 68 65 20 |modulus |and the |
|000025a0| 61 62 73 6f 6c 75 74 65 | 0a 58 09 20 2a 20 76 61 |absolute|.X. * va|
|000025b0| 6c 75 65 20 6f 66 20 74 | 68 65 20 6f 72 69 67 69 |lue of t|he origi|
|000025c0| 6e 61 6c 20 6e 75 6d 62 | 65 72 2e 20 20 43 6f 6d |nal numb|er.  Com|
|000025d0| 70 61 72 65 20 74 68 65 | 20 6f 72 69 67 69 6e 61 |pare the| origina|
|000025e0| 6c 20 6e 75 6d 62 65 72 | 20 77 69 74 68 0a 58 09 |l number| with.X.|
|000025f0| 20 2a 20 74 68 65 20 64 | 69 66 66 65 72 65 6e 63 | * the d|ifferenc|
|00002600| 65 2c 20 61 6e 64 20 72 | 65 74 75 72 6e 20 74 68 |e, and r|eturn th|
|00002610| 65 20 6f 6e 65 20 77 69 | 74 68 20 74 68 65 20 73 |e one wi|th the s|
|00002620| 6d 61 6c 6c 65 73 74 20 | 61 62 73 6f 6c 75 74 65 |mallest |absolute|
|00002630| 0a 58 09 20 2a 20 76 61 | 6c 75 65 2c 20 77 69 74 |.X. * va|lue, wit|
|00002640| 68 20 74 68 65 20 63 6f | 72 72 65 63 74 20 73 69 |h the co|rrect si|
|00002650| 67 6e 2e 20 20 49 66 20 | 74 68 65 20 74 77 6f 20 |gn.  If |the two |
|00002660| 76 61 6c 75 65 73 20 61 | 72 65 20 65 71 75 61 6c |values a|re equal|
|00002670| 2c 20 74 68 65 6e 0a 58 | 09 20 2a 20 72 65 74 75 |, then.X|. * retu|
|00002680| 72 6e 20 74 68 65 20 70 | 6f 73 69 74 69 76 65 20 |rn the p|ositive |
|00002690| 72 65 73 75 6c 74 2e 0a | 58 09 20 2a 2f 0a 58 09 |result..|X. */.X.|
|000026a0| 7a 73 75 62 28 7a 32 2c | 20 74 6d 70 31 2c 20 26 |zsub(z2,| tmp1, &|
|000026b0| 74 6d 70 32 29 3b 0a 58 | 09 63 76 20 3d 20 7a 72 |tmp2);.X|.cv = zr|
|000026c0| 65 6c 28 74 6d 70 31 2c | 20 74 6d 70 32 29 3b 0a |el(tmp1,| tmp2);.|
|000026d0| 58 09 69 66 20 28 63 76 | 20 3c 20 30 29 20 7b 0a |X.if (cv| < 0) {.|
|000026e0| 58 09 09 66 72 65 65 68 | 28 74 6d 70 32 2e 76 29 |X..freeh|(tmp2.v)|
|000026f0| 3b 0a 58 09 09 74 6d 70 | 31 2e 73 69 67 6e 20 3d |;.X..tmp|1.sign =|
|00002700| 20 28 42 4f 4f 4c 29 73 | 69 67 6e 3b 0a 58 09 09 | (BOOL)s|ign;.X..|
|00002710| 69 66 20 28 74 6d 70 31 | 2e 76 20 3d 3d 20 7a 31 |if (tmp1|.v == z1|
|00002720| 2e 76 29 0a 58 09 09 09 | 7a 63 6f 70 79 28 74 6d |.v).X...|zcopy(tm|
|00002730| 70 31 2c 20 72 65 73 29 | 3b 0a 58 09 09 65 6c 73 |p1, res)|;.X..els|
|00002740| 65 0a 58 09 09 09 2a 72 | 65 73 20 3d 20 74 6d 70 |e.X...*r|es = tmp|
|00002750| 31 3b 0a 58 09 7d 20 65 | 6c 73 65 20 7b 0a 58 09 |1;.X.} e|lse {.X.|
|00002760| 09 69 66 20 28 63 76 29 | 0a 58 09 09 09 74 6d 70 |.if (cv)|.X...tmp|
|00002770| 32 2e 73 69 67 6e 20 3d | 20 21 73 69 67 6e 3b 0a |2.sign =| !sign;.|
|00002780| 58 09 09 69 66 20 28 74 | 6d 70 31 2e 76 20 21 3d |X..if (t|mp1.v !=|
|00002790| 20 7a 31 2e 76 29 0a 58 | 09 09 09 66 72 65 65 68 | z1.v).X|...freeh|
|000027a0| 28 74 6d 70 31 2e 76 29 | 3b 0a 58 09 09 2a 72 65 |(tmp1.v)|;.X..*re|
|000027b0| 73 20 3d 20 74 6d 70 32 | 3b 0a 58 09 7d 0a 58 7d |s = tmp2|;.X.}.X}|
|000027c0| 0a 58 0a 58 0a 58 2f 2a | 0a 58 20 2a 20 43 6f 6d |.X.X.X/*|.X * Com|
|000027d0| 70 61 72 65 20 74 77 6f | 20 6e 75 6d 62 65 72 73 |pare two| numbers|
|000027e0| 20 66 6f 72 20 65 71 75 | 61 6c 69 74 79 20 6d 6f | for equ|ality mo|
|000027f0| 64 75 6c 6f 20 61 20 74 | 68 69 72 64 20 6e 75 6d |dulo a t|hird num|
|00002800| 62 65 72 2e 0a 58 20 2a | 20 54 68 65 20 74 77 6f |ber..X *| The two|
|00002810| 20 6e 75 6d 62 65 72 73 | 20 74 6f 20 62 65 20 63 | numbers| to be c|
|00002820| 6f 6d 70 61 72 65 64 20 | 63 61 6e 20 62 65 20 6e |ompared |can be n|
|00002830| 65 67 61 74 69 76 65 20 | 6f 72 20 6f 75 74 20 6f |egative |or out o|
|00002840| 66 20 6d 6f 64 75 6c 6f | 20 72 61 6e 67 65 2e 0a |f modulo| range..|
|00002850| 58 20 2a 20 52 65 74 75 | 72 6e 73 20 54 52 55 45 |X * Retu|rns TRUE|
|00002860| 20 69 66 20 74 68 65 20 | 6e 75 6d 62 65 72 73 20 | if the |numbers |
|00002870| 61 72 65 20 6e 6f 74 20 | 63 6f 6e 67 72 75 65 6e |are not |congruen|
|00002880| 74 2c 20 61 6e 64 20 46 | 41 4c 53 45 20 69 66 20 |t, and F|ALSE if |
|00002890| 74 68 65 79 20 61 72 65 | 0a 58 20 2a 20 63 6f 6e |they are|.X * con|
|000028a0| 67 72 75 65 6e 74 2e 0a | 58 20 2a 2f 0a 58 42 4f |gruent..|X */.XBO|
|000028b0| 4f 4c 0a 58 7a 63 6d 70 | 6d 6f 64 28 7a 31 2c 20 |OL.Xzcmp|mod(z1, |
|000028c0| 7a 32 2c 20 7a 33 29 0a | 58 09 5a 56 41 4c 55 45 |z2, z3).|X.ZVALUE|
|000028d0| 20 7a 31 3b 09 09 2f 2a | 20 66 69 72 73 74 20 6e | z1;../*| first n|
|000028e0| 75 6d 62 65 72 20 74 6f | 20 62 65 20 63 6f 6d 70 |umber to| be comp|
|000028f0| 61 72 65 64 20 2a 2f 0a | 58 09 5a 56 41 4c 55 45 |ared */.|X.ZVALUE|
|00002900| 20 7a 32 3b 09 09 2f 2a | 20 73 65 63 6f 6e 64 20 | z2;../*| second |
|00002910| 6e 75 6d 62 65 72 20 74 | 6f 20 62 65 20 63 6f 6d |number t|o be com|
|00002920| 70 61 72 65 64 20 2a 2f | 0a 58 09 5a 56 41 4c 55 |pared */|.X.ZVALU|
|00002930| 45 20 7a 33 3b 09 09 2f | 2a 20 6d 6f 64 75 6c 75 |E z3;../|* modulu|
|00002940| 73 20 2a 2f 0a 58 7b 0a | 58 09 5a 56 41 4c 55 45 |s */.X{.|X.ZVALUE|
|00002950| 20 74 6d 70 31 2c 20 74 | 6d 70 32 2c 20 74 6d 70 | tmp1, t|mp2, tmp|
|00002960| 33 3b 0a 58 09 46 55 4c | 4c 20 64 69 67 69 74 3b |3;.X.FUL|L digit;|
|00002970| 0a 58 09 4c 45 4e 20 6c | 65 6e 3b 0a 58 09 69 6e |.X.LEN l|en;.X.in|
|00002980| 74 20 63 76 3b 0a 58 0a | 58 09 69 66 20 28 69 73 |t cv;.X.|X.if (is|
|00002990| 6e 65 67 28 7a 33 29 20 | 7c 7c 20 69 73 7a 65 72 |neg(z3) ||| iszer|
|000029a0| 6f 28 7a 33 29 29 0a 58 | 09 09 65 72 72 6f 72 28 |o(z3)).X|..error(|
|000029b0| 22 4e 6f 6e 2d 70 6f 73 | 69 74 69 76 65 20 6d 6f |"Non-pos|itive mo|
|000029c0| 64 75 6c 75 73 20 69 6e | 20 7a 63 6d 70 6d 6f 64 |dulus in| zcmpmod|
|000029d0| 22 29 3b 0a 58 09 69 66 | 20 28 69 73 74 77 6f 28 |");.X.if| (istwo(|
|000029e0| 7a 33 29 29 0a 58 09 09 | 72 65 74 75 72 6e 20 28 |z3)).X..|return (|
|000029f0| 28 28 7a 31 2e 76 5b 30 | 5d 20 2b 20 7a 32 2e 76 |((z1.v[0|] + z2.v|
|00002a00| 5b 30 5d 29 20 26 20 30 | 78 31 29 20 21 3d 20 30 |[0]) & 0|x1) != 0|
|00002a10| 29 3b 0a 58 0a 58 09 2f | 2a 0a 58 09 20 2a 20 49 |);.X.X./|*.X. * I|
|00002a20| 66 20 74 68 65 20 74 77 | 6f 20 6e 75 6d 62 65 72 |f the tw|o number|
|00002a30| 73 20 61 72 65 20 65 71 | 75 61 6c 2c 20 74 68 65 |s are eq|ual, the|
|00002a40| 6e 20 74 68 65 69 72 20 | 6d 6f 64 73 20 61 72 65 |n their |mods are|
|00002a50| 20 65 71 75 61 6c 2e 0a | 58 09 20 2a 2f 0a 58 09 | equal..|X. */.X.|
|00002a60| 69 66 20 28 28 7a 31 2e | 73 69 67 6e 20 3d 3d 20 |if ((z1.|sign == |
|00002a70| 7a 32 2e 73 69 67 6e 29 | 20 26 26 20 28 7a 31 2e |z2.sign)| && (z1.|
|00002a80| 6c 65 6e 20 3d 3d 20 7a | 32 2e 6c 65 6e 29 20 26 |len == z|2.len) &|
|00002a90| 26 0a 58 09 09 28 7a 31 | 2e 76 5b 30 5d 20 3d 3d |&.X..(z1|.v[0] ==|
|00002aa0| 20 7a 32 2e 76 5b 30 5d | 29 20 26 26 20 28 7a 63 | z2.v[0]|) && (zc|
|00002ab0| 6d 70 28 7a 31 2c 20 7a | 32 29 20 3d 3d 20 30 29 |mp(z1, z|2) == 0)|
|00002ac0| 29 0a 58 09 09 09 72 65 | 74 75 72 6e 20 46 41 4c |).X...re|turn FAL|
|00002ad0| 53 45 3b 0a 58 0a 58 09 | 2f 2a 0a 58 09 20 2a 20 |SE;.X.X.|/*.X. * |
|00002ae0| 49 66 20 62 6f 74 68 20 | 6e 75 6d 62 65 72 73 20 |If both |numbers |
|00002af0| 61 72 65 20 6e 65 67 61 | 74 69 76 65 2c 20 74 68 |are nega|tive, th|
|00002b00| 65 6e 20 77 65 20 63 61 | 6e 20 6d 61 6b 65 20 74 |en we ca|n make t|
|00002b10| 68 65 6d 20 70 6f 73 69 | 74 69 76 65 2e 0a 58 09 |hem posi|tive..X.|
|00002b20| 20 2a 2f 0a 58 09 69 66 | 20 28 69 73 6e 65 67 28 | */.X.if| (isneg(|
|00002b30| 7a 31 29 20 26 26 20 69 | 73 6e 65 67 28 7a 32 29 |z1) && i|sneg(z2)|
|00002b40| 29 20 7b 0a 58 09 09 7a | 31 2e 73 69 67 6e 20 3d |) {.X..z|1.sign =|
|00002b50| 20 30 3b 0a 58 09 09 7a | 32 2e 73 69 67 6e 20 3d | 0;.X..z|2.sign =|
|00002b60| 20 30 3b 0a 58 09 7d 0a | 58 0a 58 09 2f 2a 0a 58 | 0;.X.}.|X.X./*.X|
|00002b70| 09 20 2a 20 46 6f 72 20 | 73 6d 61 6c 6c 20 6e 65 |. * For |small ne|
|00002b80| 67 61 74 69 76 65 20 6e | 75 6d 62 65 72 73 2c 20 |gative n|umbers, |
|00002b90| 6d 61 6b 65 20 74 68 65 | 6d 20 70 6f 73 69 74 69 |make the|m positi|
|00002ba0| 76 65 20 62 65 66 6f 72 | 65 20 63 6f 6d 70 61 72 |ve befor|e compar|
|00002bb0| 69 6e 67 2e 0a 58 09 20 | 2a 20 49 6e 20 61 6e 79 |ing..X. |* In any|
|00002bc0| 20 63 61 73 65 2c 20 74 | 68 65 20 72 65 73 75 6c | case, t|he resul|
|00002bd0| 74 69 6e 67 20 6e 75 6d | 62 65 72 73 20 61 72 65 |ting num|bers are|
|00002be0| 20 69 6e 20 74 6d 70 31 | 20 61 6e 64 20 74 6d 70 | in tmp1| and tmp|
|00002bf0| 32 2e 0a 58 09 20 2a 2f | 0a 58 09 74 6d 70 31 20 |2..X. */|.X.tmp1 |
|00002c00| 3d 20 7a 31 3b 0a 58 09 | 74 6d 70 32 20 3d 20 7a |= z1;.X.|tmp2 = z|
|00002c10| 32 3b 0a 58 09 6c 65 6e | 20 3d 20 7a 33 2e 6c 65 |2;.X.len| = z3.le|
|00002c20| 6e 3b 0a 58 09 64 69 67 | 69 74 20 3d 20 7a 33 2e |n;.X.dig|it = z3.|
|00002c30| 76 5b 6c 65 6e 20 2d 20 | 31 5d 3b 0a 58 0a 58 09 |v[len - |1];.X.X.|
|00002c40| 69 66 20 28 69 73 6e 65 | 67 28 7a 31 29 20 26 26 |if (isne|g(z1) &&|
|00002c50| 20 28 28 7a 31 2e 6c 65 | 6e 20 3c 20 6c 65 6e 29 | ((z1.le|n < len)|
|00002c60| 20 7c 7c 0a 58 09 09 28 | 28 7a 31 2e 6c 65 6e 20 | ||.X..(|(z1.len |
|00002c70| 3d 3d 20 6c 65 6e 29 20 | 26 26 20 28 7a 31 2e 76 |== len) |&& (z1.v|
|00002c80| 5b 7a 31 2e 6c 65 6e 20 | 2d 20 31 5d 20 3c 20 64 |[z1.len |- 1] < d|
|00002c90| 69 67 69 74 29 29 29 29 | 0a 58 09 09 09 7a 61 64 |igit))))|.X...zad|
|00002ca0| 64 28 7a 31 2c 20 7a 33 | 2c 20 26 74 6d 70 31 29 |d(z1, z3|, &tmp1)|
|00002cb0| 3b 0a 58 0a 58 09 69 66 | 20 28 69 73 6e 65 67 28 |;.X.X.if| (isneg(|
|00002cc0| 7a 32 29 20 26 26 20 28 | 28 7a 32 2e 6c 65 6e 20 |z2) && (|(z2.len |
|00002cd0| 3c 20 6c 65 6e 29 20 7c | 7c 0a 58 09 09 28 28 7a |< len) |||.X..((z|
|00002ce0| 32 2e 6c 65 6e 20 3d 3d | 20 6c 65 6e 29 20 26 26 |2.len ==| len) &&|
|00002cf0| 20 28 7a 32 2e 76 5b 7a | 32 2e 6c 65 6e 20 2d 20 | (z2.v[z|2.len - |
|00002d00| 31 5d 20 3c 20 64 69 67 | 69 74 29 29 29 29 0a 58 |1] < dig|it)))).X|
|00002d10| 09 09 09 7a 61 64 64 28 | 7a 32 2c 20 7a 33 2c 20 |...zadd(|z2, z3, |
|00002d20| 26 74 6d 70 32 29 3b 0a | 58 0a 58 09 2f 2a 0a 58 |&tmp2);.|X.X./*.X|
|00002d30| 09 20 2a 20 4e 6f 77 20 | 63 6f 6d 70 61 72 65 20 |. * Now |compare |
|00002d40| 74 68 65 20 74 77 6f 20 | 6e 75 6d 62 65 72 73 20 |the two |numbers |
|00002d50| 66 6f 72 20 65 71 75 61 | 6c 69 74 79 2e 0a 58 09 |for equa|lity..X.|
|00002d60| 20 2a 20 49 66 20 74 68 | 65 79 20 61 72 65 20 65 | * If th|ey are e|
|00002d70| 71 75 61 6c 20 77 65 20 | 61 72 65 20 61 6c 6c 20 |qual we |are all |
|00002d80| 64 6f 6e 65 2e 0a 58 09 | 20 2a 2f 0a 58 09 69 66 |done..X.| */.X.if|
|00002d90| 20 28 7a 63 6d 70 28 74 | 6d 70 31 2c 20 74 6d 70 | (zcmp(t|mp1, tmp|
|00002da0| 32 29 20 3d 3d 20 30 29 | 20 7b 0a 58 09 09 69 66 |2) == 0)| {.X..if|
|00002db0| 20 28 74 6d 70 31 2e 76 | 20 21 3d 20 7a 31 2e 76 | (tmp1.v| != z1.v|
|00002dc0| 29 0a 58 09 09 09 66 72 | 65 65 68 28 74 6d 70 31 |).X...fr|eeh(tmp1|
|00002dd0| 2e 76 29 3b 0a 58 09 09 | 69 66 20 28 74 6d 70 32 |.v);.X..|if (tmp2|
|00002de0| 2e 76 20 21 3d 20 7a 32 | 2e 76 29 0a 58 09 09 09 |.v != z2|.v).X...|
|00002df0| 66 72 65 65 68 28 74 6d | 70 32 2e 76 29 3b 0a 58 |freeh(tm|p2.v);.X|
|00002e00| 09 09 72 65 74 75 72 6e | 20 46 41 4c 53 45 3b 0a |..return| FALSE;.|
|00002e10| 58 09 7d 0a 58 0a 58 09 | 2f 2a 0a 58 09 20 2a 20 |X.}.X.X.|/*.X. * |
|00002e20| 54 68 65 79 20 61 72 65 | 20 6e 6f 74 20 69 64 65 |They are| not ide|
|00002e30| 6e 74 69 63 61 6c 2e 20 | 20 4e 6f 77 20 69 66 20 |ntical. | Now if |
|00002e40| 62 6f 74 68 20 6e 75 6d | 62 65 72 73 20 61 72 65 |both num|bers are|
|00002e50| 20 70 6f 73 69 74 69 76 | 65 0a 58 09 20 2a 20 61 | positiv|e.X. * a|
|00002e60| 6e 64 20 6c 65 73 73 20 | 74 68 61 6e 20 74 68 65 |nd less |than the|
|00002e70| 20 6d 6f 64 75 6c 75 73 | 2c 20 74 68 65 6e 20 74 | modulus|, then t|
|00002e80| 68 65 79 20 61 72 65 20 | 64 65 66 69 6e 69 74 65 |hey are |definite|
|00002e90| 6c 79 20 6e 6f 74 20 65 | 71 75 61 6c 2e 0a 58 09 |ly not e|qual..X.|
|00002ea0| 20 2a 2f 0a 58 09 69 66 | 20 28 28 74 6d 70 31 2e | */.X.if| ((tmp1.|
|00002eb0| 73 69 67 6e 20 3d 3d 20 | 74 6d 70 32 2e 73 69 67 |sign == |tmp2.sig|
|00002ec0| 6e 29 20 26 26 0a 58 09 | 09 28 28 74 6d 70 31 2e |n) &&.X.|.((tmp1.|
|00002ed0| 6c 65 6e 20 3c 20 6c 65 | 6e 29 20 7c 7c 20 28 7a |len < le|n) || (z|
|00002ee0| 72 65 6c 28 74 6d 70 31 | 2c 20 7a 33 29 20 3c 20 |rel(tmp1|, z3) < |
|00002ef0| 30 29 29 20 26 26 0a 58 | 09 09 28 28 74 6d 70 32 |0)) &&.X|..((tmp2|
|00002f00| 2e 6c 65 6e 20 3c 20 6c | 65 6e 29 20 7c 7c 20 28 |.len < l|en) || (|
|00002f10| 7a 72 65 6c 28 74 6d 70 | 32 2c 20 7a 33 29 20 3c |zrel(tmp|2, z3) <|
|00002f20| 20 30 29 29 29 0a 58 09 | 7b 0a 58 09 09 69 66 20 | 0))).X.|{.X..if |
|00002f30| 28 74 6d 70 31 2e 76 20 | 21 3d 20 7a 31 2e 76 29 |(tmp1.v |!= z1.v)|
|00002f40| 0a 58 09 09 09 66 72 65 | 65 68 28 74 6d 70 31 2e |.X...fre|eh(tmp1.|
|00002f50| 76 29 3b 0a 58 09 09 69 | 66 20 28 74 6d 70 32 2e |v);.X..i|f (tmp2.|
|00002f60| 76 20 21 3d 20 7a 32 2e | 76 29 0a 58 09 09 09 66 |v != z2.|v).X...f|
|00002f70| 72 65 65 68 28 74 6d 70 | 32 2e 76 29 3b 0a 58 09 |reeh(tmp|2.v);.X.|
|00002f80| 09 72 65 74 75 72 6e 20 | 54 52 55 45 3b 0a 58 09 |.return |TRUE;.X.|
|00002f90| 7d 0a 58 0a 58 09 2f 2a | 0a 58 09 20 2a 20 45 69 |}.X.X./*|.X. * Ei|
|00002fa0| 74 68 65 72 20 6f 6e 65 | 20 6f 66 20 74 68 65 20 |ther one| of the |
|00002fb0| 6e 75 6d 62 65 72 73 20 | 69 73 20 6e 65 67 61 74 |numbers |is negat|
|00002fc0| 69 76 65 20 6f 72 20 69 | 73 20 6c 61 72 67 65 2e |ive or i|s large.|
|00002fd0| 0a 58 09 20 2a 20 53 6f | 20 64 6f 20 74 68 65 20 |.X. * So| do the |
|00002fe0| 73 74 61 6e 64 61 72 64 | 20 74 68 69 6e 67 20 61 |standard| thing a|
|00002ff0| 6e 64 20 73 75 62 74 72 | 61 63 74 20 74 68 65 20 |nd subtr|act the |
|00003000| 74 77 6f 20 6e 75 6d 62 | 65 72 73 2e 0a 58 09 20 |two numb|ers..X. |
|00003010| 2a 20 54 68 65 6e 20 74 | 68 65 79 20 61 72 65 20 |* Then t|hey are |
|00003020| 65 71 75 61 6c 20 69 66 | 20 74 68 65 20 72 65 73 |equal if| the res|
|00003030| 75 6c 74 20 69 73 20 30 | 20 28 6d 6f 64 20 7a 33 |ult is 0| (mod z3|
|00003040| 29 2e 0a 58 09 20 2a 2f | 0a 58 09 7a 73 75 62 28 |)..X. */|.X.zsub(|
|00003050| 74 6d 70 31 2c 20 74 6d | 70 32 2c 20 26 74 6d 70 |tmp1, tm|p2, &tmp|
|00003060| 33 29 3b 0a 58 09 69 66 | 20 28 74 6d 70 31 2e 76 |3);.X.if| (tmp1.v|
|00003070| 20 21 3d 20 7a 31 2e 76 | 29 0a 58 09 09 66 72 65 | != z1.v|).X..fre|
|00003080| 65 68 28 74 6d 70 31 2e | 76 29 3b 0a 58 09 69 66 |eh(tmp1.|v);.X.if|
|00003090| 20 28 74 6d 70 32 2e 76 | 20 21 3d 20 7a 32 2e 76 | (tmp2.v| != z2.v|
|000030a0| 29 0a 58 09 09 66 72 65 | 65 68 28 74 6d 70 32 2e |).X..fre|eh(tmp2.|
|000030b0| 76 29 3b 0a 58 0a 58 09 | 2f 2a 0a 58 09 20 2a 20 |v);.X.X.|/*.X. * |
|000030c0| 43 6f 6d 70 61 72 65 20 | 74 68 65 20 72 65 73 75 |Compare |the resu|
|000030d0| 6c 74 20 77 69 74 68 20 | 74 68 65 20 6d 6f 64 75 |lt with |the modu|
|000030e0| 6c 75 73 20 74 6f 20 73 | 65 65 20 69 66 20 69 74 |lus to s|ee if it|
|000030f0| 20 69 73 20 65 71 75 61 | 6c 20 74 6f 0a 58 09 20 | is equa|l to.X. |
|00003100| 2a 20 6f 72 20 6c 65 73 | 73 20 74 68 61 6e 20 74 |* or les|s than t|
|00003110| 68 65 20 6d 6f 64 75 6c | 75 73 2e 20 20 49 66 20 |he modul|us.  If |
|00003120| 73 6f 2c 20 77 65 20 6b | 6e 6f 77 20 74 68 65 20 |so, we k|now the |
|00003130| 6d 6f 64 20 72 65 73 75 | 6c 74 2e 0a 58 09 20 2a |mod resu|lt..X. *|
|00003140| 2f 0a 58 09 74 6d 70 33 | 2e 73 69 67 6e 20 3d 20 |/.X.tmp3|.sign = |
|00003150| 30 3b 0a 58 09 63 76 20 | 3d 20 7a 72 65 6c 28 74 |0;.X.cv |= zrel(t|
|00003160| 6d 70 33 2c 20 7a 33 29 | 3b 0a 58 09 69 66 20 28 |mp3, z3)|;.X.if (|
|00003170| 63 76 20 3d 3d 20 30 29 | 20 7b 0a 58 09 09 66 72 |cv == 0)| {.X..fr|
|00003180| 65 65 68 28 74 6d 70 33 | 2e 76 29 3b 0a 58 09 09 |eeh(tmp3|.v);.X..|
|00003190| 72 65 74 75 72 6e 20 46 | 41 4c 53 45 3b 0a 58 09 |return F|ALSE;.X.|
|000031a0| 7d 0a 58 09 69 66 20 28 | 63 76 20 3c 20 30 29 20 |}.X.if (|cv < 0) |
|000031b0| 7b 0a 58 09 09 66 72 65 | 65 68 28 74 6d 70 33 2e |{.X..fre|eh(tmp3.|
|000031c0| 76 29 3b 0a 58 09 09 72 | 65 74 75 72 6e 20 54 52 |v);.X..r|eturn TR|
|000031d0| 55 45 3b 0a 58 09 7d 0a | 58 0a 58 09 2f 2a 0a 58 |UE;.X.}.|X.X./*.X|
|000031e0| 09 20 2a 20 57 65 20 61 | 72 65 20 66 6f 72 63 65 |. * We a|re force|
|000031f0| 64 20 74 6f 20 61 63 74 | 75 61 6c 6c 79 20 64 6f |d to act|ually do|
|00003200| 20 74 68 65 20 64 69 76 | 69 73 69 6f 6e 2e 0a 58 | the div|ision..X|
|00003210| 09 20 2a 20 54 68 65 20 | 6e 75 6d 62 65 72 73 20 |. * The |numbers |
|00003220| 61 72 65 20 63 6f 6e 67 | 72 75 65 6e 74 20 69 66 |are cong|ruent if|
|00003230| 20 74 68 65 20 72 65 73 | 75 6c 74 20 69 73 20 7a | the res|ult is z|
|00003240| 65 72 6f 2e 0a 58 09 20 | 2a 2f 0a 58 09 7a 6d 6f |ero..X. |*/.X.zmo|
|00003250| 64 28 74 6d 70 33 2c 20 | 7a 33 2c 20 26 74 6d 70 |d(tmp3, |z3, &tmp|
|00003260| 31 29 3b 0a 58 09 66 72 | 65 65 68 28 74 6d 70 33 |1);.X.fr|eeh(tmp3|
|00003270| 2e 76 29 3b 0a 58 09 69 | 66 20 28 69 73 7a 65 72 |.v);.X.i|f (iszer|
|00003280| 6f 28 74 6d 70 31 29 29 | 20 7b 0a 58 09 09 66 72 |o(tmp1))| {.X..fr|
|00003290| 65 65 68 28 74 6d 70 31 | 2e 76 29 3b 0a 58 09 09 |eeh(tmp1|.v);.X..|
|000032a0| 72 65 74 75 72 6e 20 46 | 41 4c 53 45 3b 0a 58 09 |return F|ALSE;.X.|
|000032b0| 7d 20 65 6c 73 65 20 7b | 0a 58 09 09 66 72 65 65 |} else {|.X..free|
|000032c0| 68 28 74 6d 70 31 2e 76 | 29 3b 0a 58 09 09 72 65 |h(tmp1.v|);.X..re|
|000032d0| 74 75 72 6e 20 54 52 55 | 45 3b 0a 58 09 7d 0a 58 |turn TRU|E;.X.}.X|
|000032e0| 7d 0a 58 0a 58 0a 58 2f | 2a 0a 58 20 2a 20 43 6f |}.X.X.X/|*.X * Co|
|000032f0| 6d 70 75 74 65 20 74 68 | 65 20 72 65 73 75 6c 74 |mpute th|e result|
|00003300| 20 6f 66 20 72 61 69 73 | 69 6e 67 20 6f 6e 65 20 | of rais|ing one |
|00003310| 6e 75 6d 62 65 72 20 74 | 6f 20 61 20 70 6f 77 65 |number t|o a powe|
|00003320| 72 20 6d 6f 64 75 6c 6f | 20 61 6e 6f 74 68 65 72 |r modulo| another|
|00003330| 20 6e 75 6d 62 65 72 2e | 0a 58 20 2a 20 54 68 61 | number.|.X * Tha|
|00003340| 74 20 69 73 2c 20 74 68 | 69 73 20 63 6f 6d 70 75 |t is, th|is compu|
|00003350| 74 65 73 3a 20 20 61 5e | 62 20 28 6d 6f 64 75 6c |tes:  a^|b (modul|
|00003360| 6f 20 63 29 2e 0a 58 20 | 2a 20 54 68 69 73 20 63 |o c)..X |* This c|
|00003370| 61 6c 63 75 6c 61 74 65 | 73 20 74 68 65 20 72 65 |alculate|s the re|
|00003380| 73 75 6c 74 20 62 79 20 | 65 78 61 6d 69 6e 69 6e |sult by |examinin|
|00003390| 67 20 74 68 65 20 70 6f | 77 65 72 20 50 4f 57 42 |g the po|wer POWB|
|000033a0| 49 54 53 20 62 69 74 73 | 20 61 74 20 61 20 74 69 |ITS bits| at a ti|
|000033b0| 6d 65 2c 0a 58 20 2a 20 | 75 73 69 6e 67 20 61 20 |me,.X * |using a |
|000033c0| 73 6d 61 6c 6c 20 74 61 | 62 6c 65 20 6f 66 20 50 |small ta|ble of P|
|000033d0| 4f 57 4e 55 4d 53 20 6c | 6f 77 20 70 6f 77 65 72 |OWNUMS l|ow power|
|000033e0| 73 20 74 6f 20 63 61 6c | 63 75 6c 61 74 65 20 70 |s to cal|culate p|
|000033f0| 6f 77 65 72 73 20 66 6f | 72 20 74 68 6f 73 65 20 |owers fo|r those |
|00003400| 62 69 74 73 2c 0a 58 20 | 2a 20 61 6e 64 20 72 65 |bits,.X |* and re|
|00003410| 70 65 61 74 65 64 20 73 | 71 75 61 72 69 6e 67 20 |peated s|quaring |
|00003420| 61 6e 64 20 6d 75 6c 74 | 69 70 6c 79 69 6e 67 20 |and mult|iplying |
|00003430| 62 79 20 74 68 65 20 70 | 61 72 74 69 61 6c 20 70 |by the p|artial p|
|00003440| 6f 77 65 72 73 20 74 6f | 20 67 65 6e 65 72 61 74 |owers to| generat|
|00003450| 65 0a 58 20 2a 20 74 68 | 65 20 63 6f 6d 70 6c 65 |e.X * th|e comple|
|00003460| 74 65 20 70 6f 77 65 72 | 2e 20 20 49 66 20 74 68 |te power|.  If th|
|00003470| 65 20 70 6f 77 65 72 20 | 62 65 69 6e 67 20 72 61 |e power |being ra|
|00003480| 69 73 65 64 20 74 6f 20 | 69 73 20 68 69 67 68 20 |ised to |is high |
|00003490| 65 6e 6f 75 67 68 2c 20 | 74 68 65 6e 0a 58 20 2a |enough, |then.X *|
|000034a0| 20 74 68 69 73 20 75 73 | 65 73 20 74 68 65 20 52 | this us|es the R|
|000034b0| 45 44 43 20 61 6c 67 6f | 72 69 74 68 6d 20 74 6f |EDC algo|rithm to|
|000034c0| 20 61 76 6f 69 64 20 64 | 6f 69 6e 67 20 6d 61 6e | avoid d|oing man|
|000034d0| 79 20 64 69 76 69 73 69 | 6f 6e 73 2e 20 20 57 68 |y divisi|ons.  Wh|
|000034e0| 65 6e 20 75 73 69 6e 67 | 0a 58 20 2a 20 52 45 44 |en using|.X * RED|
|000034f0| 43 2c 20 6d 75 6c 74 69 | 70 6c 65 20 63 61 6c 6c |C, multi|ple call|
|00003500| 73 20 74 6f 20 74 68 69 | 73 20 72 6f 75 74 69 6e |s to thi|s routin|
|00003510| 65 20 75 73 69 6e 67 20 | 74 68 65 20 73 61 6d 65 |e using |the same|
|00003520| 20 6d 6f 64 75 6c 75 73 | 20 77 69 6c 6c 20 62 65 | modulus| will be|
|00003530| 0a 58 20 2a 20 73 6c 69 | 67 68 74 6c 79 20 66 61 |.X * sli|ghtly fa|
|00003540| 73 74 65 72 2e 0a 58 20 | 2a 2f 0a 58 76 6f 69 64 |ster..X |*/.Xvoid|
|00003550| 0a 58 7a 70 6f 77 65 72 | 6d 6f 64 28 7a 31 2c 20 |.Xzpower|mod(z1, |
|00003560| 7a 32 2c 20 7a 33 2c 20 | 72 65 73 29 0a 58 09 5a |z2, z3, |res).X.Z|
|00003570| 56 41 4c 55 45 20 7a 31 | 2c 20 7a 32 2c 20 7a 33 |VALUE z1|, z2, z3|
|00003580| 2c 20 2a 72 65 73 3b 0a | 58 7b 0a 58 09 48 41 4c |, *res;.|X{.X.HAL|
|00003590| 46 20 2a 68 70 3b 09 09 | 2f 2a 20 70 6f 69 6e 74 |F *hp;..|/* point|
|000035a0| 65 72 20 74 6f 20 63 75 | 72 72 65 6e 74 20 77 6f |er to cu|rrent wo|
|000035b0| 72 64 20 6f 66 20 74 68 | 65 20 70 6f 77 65 72 20 |rd of th|e power |
|000035c0| 2a 2f 0a 58 09 52 45 44 | 43 20 2a 72 70 3b 09 09 |*/.X.RED|C *rp;..|
|000035d0| 2f 2a 20 52 45 44 43 20 | 69 6e 66 6f 72 6d 61 74 |/* REDC |informat|
|000035e0| 69 6f 6e 20 74 6f 20 62 | 65 20 75 73 65 64 20 2a |ion to b|e used *|
|000035f0| 2f 0a 58 09 5a 56 41 4c | 55 45 20 2a 70 70 3b 09 |/.X.ZVAL|UE *pp;.|
|00003600| 09 2f 2a 20 70 6f 69 6e | 74 65 72 20 74 6f 20 6c |./* poin|ter to l|
|00003610| 6f 77 20 70 6f 77 65 72 | 20 74 61 62 6c 65 20 2a |ow power| table *|
|00003620| 2f 0a 58 09 5a 56 41 4c | 55 45 20 61 6e 73 2c 20 |/.X.ZVAL|UE ans, |
|00003630| 74 65 6d 70 3b 09 2f 2a | 20 63 61 6c 63 75 6c 61 |temp;./*| calcula|
|00003640| 74 69 6f 6e 20 76 61 6c | 75 65 73 20 2a 2f 0a 58 |tion val|ues */.X|
|00003650| 09 5a 56 41 4c 55 45 20 | 6d 6f 64 70 6f 77 3b 09 |.ZVALUE |modpow;.|
|00003660| 09 2f 2a 20 63 75 72 72 | 65 6e 74 20 73 6d 61 6c |./* curr|ent smal|
|00003670| 6c 20 70 6f 77 65 72 20 | 2a 2f 0a 58 09 5a 56 41 |l power |*/.X.ZVA|
|00003680| 4c 55 45 20 6c 6f 77 70 | 6f 77 65 72 73 5b 50 4f |LUE lowp|owers[PO|
|00003690| 57 4e 55 4d 53 5d 3b 09 | 2f 2a 20 6c 6f 77 20 70 |WNUMS];.|/* low p|
|000036a0| 6f 77 65 72 73 20 2a 2f | 0a 58 09 69 6e 74 20 73 |owers */|.X.int s|
|000036b0| 69 67 6e 3b 09 09 2f 2a | 20 6f 72 69 67 69 6e 61 |ign;../*| origina|
|000036c0| 6c 20 73 69 67 6e 20 6f | 66 20 6e 75 6d 62 65 72 |l sign o|f number|
|000036d0| 20 2a 2f 0a 58 09 69 6e | 74 20 63 75 72 73 68 69 | */.X.in|t curshi|
|000036e0| 66 74 3b 09 09 2f 2a 20 | 73 68 69 66 74 20 76 61 |ft;../* |shift va|
|000036f0| 6c 75 65 20 66 6f 72 20 | 77 6f 72 64 20 6f 66 20 |lue for |word of |
|00003700| 70 6f 77 65 72 20 2a 2f | 0a 58 09 48 41 4c 46 20 |power */|.X.HALF |
|00003710| 63 75 72 68 61 6c 66 3b | 09 09 2f 2a 20 63 75 72 |curhalf;|../* cur|
|00003720| 72 65 6e 74 20 77 6f 72 | 64 20 6f 66 20 70 6f 77 |rent wor|d of pow|
|00003730| 65 72 20 2a 2f 0a 58 09 | 75 6e 73 69 67 6e 65 64 |er */.X.|unsigned|
|00003740| 20 69 6e 74 20 63 75 72 | 70 6f 77 3b 09 2f 2a 20 | int cur|pow;./* |
|00003750| 63 75 72 72 65 6e 74 20 | 6c 6f 77 20 70 6f 77 65 |current |low powe|
|00003760| 72 20 2a 2f 0a 58 09 75 | 6e 73 69 67 6e 65 64 20 |r */.X.u|nsigned |
|00003770| 69 6e 74 20 63 75 72 62 | 69 74 3b 09 2f 2a 20 63 |int curb|it;./* c|
|00003780| 75 72 72 65 6e 74 20 62 | 69 74 20 6f 66 20 6c 6f |urrent b|it of lo|
|00003790| 77 20 70 6f 77 65 72 20 | 2a 2f 0a 58 09 69 6e 74 |w power |*/.X.int|
|000037a0| 20 69 3b 0a 58 0a 58 09 | 69 66 20 28 69 73 6e 65 | i;.X.X.|if (isne|
|000037b0| 67 28 7a 33 29 20 7c 7c | 20 69 73 7a 65 72 6f 28 |g(z3) ||| iszero(|
|000037c0| 7a 33 29 29 0a 58 09 09 | 65 72 72 6f 72 28 22 4e |z3)).X..|error("N|
|000037d0| 6f 6e 2d 70 6f 73 69 74 | 69 76 65 20 6d 6f 64 75 |on-posit|ive modu|
|000037e0| 6c 75 73 20 69 6e 20 7a | 70 6f 77 65 72 6d 6f 64 |lus in z|powermod|
|000037f0| 22 29 3b 0a 58 09 69 66 | 20 28 69 73 6e 65 67 28 |");.X.if| (isneg(|
|00003800| 7a 32 29 29 0a 58 09 09 | 65 72 72 6f 72 28 22 4e |z2)).X..|error("N|
|00003810| 65 67 61 74 69 76 65 20 | 70 6f 77 65 72 20 69 6e |egative |power in|
|00003820| 20 7a 70 6f 77 65 72 6d | 6f 64 22 29 3b 0a 58 0a | zpowerm|od");.X.|
|00003830| 58 09 73 69 67 6e 20 3d | 20 7a 31 2e 73 69 67 6e |X.sign =| z1.sign|
|00003840| 3b 0a 58 09 7a 31 2e 73 | 69 67 6e 20 3d 20 30 3b |;.X.z1.s|ign = 0;|
|00003850| 0a 58 0a 58 09 2f 2a 0a | 58 09 20 2a 20 43 68 65 |.X.X./*.|X. * Che|
|00003860| 63 6b 20 65 61 73 79 20 | 63 61 73 65 73 20 66 69 |ck easy |cases fi|
|00003870| 72 73 74 2e 0a 58 09 20 | 2a 2f 0a 58 09 69 66 20 |rst..X. |*/.X.if |
|00003880| 28 69 73 7a 65 72 6f 28 | 7a 31 29 20 7c 7c 20 69 |(iszero(|z1) || i|
|00003890| 73 75 6e 69 74 28 7a 33 | 29 29 20 7b 09 09 2f 2a |sunit(z3|)) {../*|
|000038a0| 20 30 5e 78 20 6f 72 20 | 6d 6f 64 20 31 20 2a 2f | 0^x or |mod 1 */|
|000038b0| 0a 58 09 09 2a 72 65 73 | 20 3d 20 5f 7a 65 72 6f |.X..*res| = _zero|
|000038c0| 5f 3b 0a 58 09 09 72 65 | 74 75 72 6e 3b 0a 58 09 |_;.X..re|turn;.X.|
|000038d0| 7d 0a 58 09 69 66 20 28 | 69 73 74 77 6f 28 7a 33 |}.X.if (|istwo(z3|
|000038e0| 29 29 20 7b 09 09 09 2f | 2a 20 6d 6f 64 20 32 20 |)) {.../|* mod 2 |
|000038f0| 2a 2f 0a 58 09 09 69 66 | 20 28 69 73 6f 64 64 28 |*/.X..if| (isodd(|
|00003900| 7a 31 29 29 0a 58 09 09 | 09 2a 72 65 73 20 3d 20 |z1)).X..|.*res = |
|00003910| 5f 6f 6e 65 5f 3b 0a 58 | 09 09 65 6c 73 65 0a 58 |_one_;.X|..else.X|
|00003920| 09 09 09 2a 72 65 73 20 | 3d 20 5f 7a 65 72 6f 5f |...*res |= _zero_|
|00003930| 3b 0a 58 09 09 72 65 74 | 75 72 6e 3b 0a 58 09 7d |;.X..ret|urn;.X.}|
|00003940| 0a 58 09 69 66 20 28 69 | 73 75 6e 69 74 28 7a 31 |.X.if (i|sunit(z1|
|00003950| 29 20 26 26 20 28 21 73 | 69 67 6e 20 7c 7c 20 69 |) && (!s|ign || i|
|00003960| 73 65 76 65 6e 28 7a 32 | 29 29 29 20 7b 0a 58 09 |seven(z2|))) {.X.|
|00003970| 09 2f 2a 20 31 5e 78 20 | 6f 72 20 28 2d 31 29 5e |./* 1^x |or (-1)^|
|00003980| 28 32 78 29 20 2a 2f 0a | 58 09 09 2a 72 65 73 20 |(2x) */.|X..*res |
|00003990| 3d 20 5f 6f 6e 65 5f 3b | 0a 58 09 09 72 65 74 75 |= _one_;|.X..retu|
|000039a0| 72 6e 3b 0a 58 09 7d 0a | 58 0a 58 09 2f 2a 0a 58 |rn;.X.}.|X.X./*.X|
|000039b0| 09 20 2a 20 4e 6f 72 6d | 61 6c 69 7a 65 20 74 68 |. * Norm|alize th|
|000039c0| 65 20 6e 75 6d 62 65 72 | 20 62 65 69 6e 67 20 72 |e number| being r|
|000039d0| 61 69 73 65 64 20 74 6f | 20 62 65 20 6e 6f 6e 2d |aised to| be non-|
|000039e0| 6e 65 67 61 74 69 76 65 | 20 61 6e 64 20 74 6f 20 |negative| and to |
|000039f0| 6c 69 65 0a 58 09 20 2a | 20 77 69 74 68 69 6e 20 |lie.X. *| within |
|00003a00| 74 68 65 20 6d 6f 64 75 | 6c 6f 20 72 61 6e 67 65 |the modu|lo range|
|00003a10| 2e 20 20 54 68 65 6e 20 | 63 68 65 63 6b 20 66 6f |.  Then |check fo|
|00003a20| 72 20 7a 65 72 6f 20 6f | 72 20 6f 6e 65 20 73 70 |r zero o|r one sp|
|00003a30| 65 63 69 61 6c 6c 79 2e | 0a 58 09 20 2a 2f 0a 58 |ecially.|.X. */.X|
|00003a40| 09 7a 6d 6f 64 28 7a 31 | 2c 20 7a 33 2c 20 26 74 |.zmod(z1|, z3, &t|
|00003a50| 65 6d 70 29 3b 0a 58 09 | 69 66 20 28 69 73 7a 65 |emp);.X.|if (isze|
|00003a60| 72 6f 28 74 65 6d 70 29 | 29 20 7b 0a 58 09 09 66 |ro(temp)|) {.X..f|
|00003a70| 72 65 65 68 28 74 65 6d | 70 2e 76 29 3b 0a 58 09 |reeh(tem|p.v);.X.|
|00003a80| 09 2a 72 65 73 20 3d 20 | 5f 7a 65 72 6f 5f 3b 0a |.*res = |_zero_;.|
|00003a90| 58 09 09 72 65 74 75 72 | 6e 3b 0a 58 09 7d 0a 58 |X..retur|n;.X.}.X|
|00003aa0| 09 7a 31 20 3d 20 74 65 | 6d 70 3b 0a 58 09 69 66 |.z1 = te|mp;.X.if|
|00003ab0| 20 28 73 69 67 6e 29 20 | 7b 0a 58 09 09 7a 73 75 | (sign) |{.X..zsu|
|00003ac0| 62 28 7a 33 2c 20 7a 31 | 2c 20 26 74 65 6d 70 29 |b(z3, z1|, &temp)|
|00003ad0| 3b 0a 58 09 09 66 72 65 | 65 68 28 7a 31 2e 76 29 |;.X..fre|eh(z1.v)|
|00003ae0| 3b 0a 58 09 09 7a 31 20 | 3d 20 74 65 6d 70 3b 0a |;.X..z1 |= temp;.|
|00003af0| 58 09 7d 0a 58 09 69 66 | 20 28 69 73 75 6e 69 74 |X.}.X.if| (isunit|
|00003b00| 28 7a 31 29 29 20 7b 0a | 58 09 09 66 72 65 65 68 |(z1)) {.|X..freeh|
|00003b10| 28 7a 31 2e 76 29 3b 0a | 58 09 09 2a 72 65 73 20 |(z1.v);.|X..*res |
|00003b20| 3d 20 5f 6f 6e 65 5f 3b | 0a 58 09 09 72 65 74 75 |= _one_;|.X..retu|
|00003b30| 72 6e 3b 0a 58 09 7d 0a | 58 0a 58 09 2f 2a 0a 58 |rn;.X.}.|X.X./*.X|
|00003b40| 09 20 2a 20 49 66 20 74 | 68 65 20 6d 6f 64 75 6c |. * If t|he modul|
|00003b50| 75 73 20 69 73 20 6f 64 | 64 2c 20 6c 61 72 67 65 |us is od|d, large|
|00003b60| 20 65 6e 6f 75 67 68 2c | 20 69 73 20 6e 6f 74 20 | enough,| is not |
|00003b70| 6f 6e 65 20 6c 65 73 73 | 20 74 68 61 6e 20 61 6e |one less| than an|
|00003b80| 0a 58 09 20 2a 20 65 78 | 61 63 74 20 70 6f 77 65 |.X. * ex|act powe|
|00003b90| 72 20 6f 66 20 74 77 6f | 2c 20 61 6e 64 20 69 66 |r of two|, and if|
|00003ba0| 20 74 68 65 20 70 6f 77 | 65 72 20 69 73 20 6c 61 | the pow|er is la|
|00003bb0| 72 67 65 20 65 6e 6f 75 | 67 68 2c 20 74 68 65 6e |rge enou|gh, then|
|00003bc0| 20 75 73 65 0a 58 09 20 | 2a 20 74 68 65 20 52 45 | use.X. |* the RE|
|00003bd0| 44 43 20 61 6c 67 6f 72 | 69 74 68 6d 2e 20 20 54 |DC algor|ithm.  T|
|00003be0| 68 65 20 73 69 7a 65 20 | 77 68 65 72 65 20 74 68 |he size |where th|
|00003bf0| 69 73 20 69 73 20 64 6f | 6e 65 20 69 73 20 63 6f |is is do|ne is co|
|00003c00| 6e 66 69 67 75 72 61 62 | 6c 65 2e 0a 58 09 20 2a |nfigurab|le..X. *|
|00003c10| 2f 0a 58 09 69 66 20 28 | 28 7a 32 2e 6c 65 6e 20 |/.X.if (|(z2.len |
|00003c20| 3e 20 31 29 20 26 26 20 | 28 7a 33 2e 6c 65 6e 20 |> 1) && |(z3.len |
|00003c30| 3e 3d 20 5f 70 6f 77 32 | 5f 29 20 26 26 20 69 73 |>= _pow2|_) && is|
|00003c40| 6f 64 64 28 7a 33 29 0a | 58 09 09 26 26 20 21 7a |odd(z3).|X..&& !z|
|00003c50| 69 73 61 6c 6c 62 69 74 | 73 28 7a 33 29 29 0a 58 |isallbit|s(z3)).X|
|00003c60| 09 7b 0a 58 09 09 69 66 | 20 28 70 6f 77 65 72 6d |.{.X..if| (powerm|
|00003c70| 6f 64 72 65 64 63 20 26 | 26 20 7a 63 6d 70 28 70 |odredc &|& zcmp(p|
|00003c80| 6f 77 65 72 6d 6f 64 72 | 65 64 63 2d 3e 6d 6f 64 |owermodr|edc->mod|
|00003c90| 2c 20 7a 33 29 29 20 7b | 0a 58 09 09 09 7a 72 65 |, z3)) {|.X...zre|
|00003ca0| 64 63 66 72 65 65 28 70 | 6f 77 65 72 6d 6f 64 72 |dcfree(p|owermodr|
|00003cb0| 65 64 63 29 3b 0a 58 09 | 09 09 70 6f 77 65 72 6d |edc);.X.|..powerm|
|00003cc0| 6f 64 72 65 64 63 20 3d | 20 4e 55 4c 4c 3b 0a 58 |odredc =| NULL;.X|
|00003cd0| 09 09 7d 0a 58 09 09 69 | 66 20 28 70 6f 77 65 72 |..}.X..i|f (power|
|00003ce0| 6d 6f 64 72 65 64 63 20 | 3d 3d 20 4e 55 4c 4c 29 |modredc |== NULL)|
|00003cf0| 0a 58 09 09 09 70 6f 77 | 65 72 6d 6f 64 72 65 64 |.X...pow|ermodred|
|00003d00| 63 20 3d 20 7a 72 65 64 | 63 61 6c 6c 6f 63 28 7a |c = zred|calloc(z|
|00003d10| 33 29 3b 0a 58 09 09 72 | 70 20 3d 20 70 6f 77 65 |3);.X..r|p = powe|
|00003d20| 72 6d 6f 64 72 65 64 63 | 3b 0a 58 09 09 7a 72 65 |rmodredc|;.X..zre|
|00003d30| 64 63 65 6e 63 6f 64 65 | 28 72 70 2c 20 7a 31 2c |dcencode|(rp, z1,|
|00003d40| 20 26 74 65 6d 70 29 3b | 0a 58 09 09 7a 72 65 64 | &temp);|.X..zred|
|00003d50| 63 70 6f 77 65 72 28 72 | 70 2c 20 74 65 6d 70 2c |cpower(r|p, temp,|
|00003d60| 20 7a 32 2c 20 26 7a 31 | 29 3b 0a 58 09 09 66 72 | z2, &z1|);.X..fr|
|00003d70| 65 65 68 28 74 65 6d 70 | 2e 76 29 3b 0a 58 09 09 |eeh(temp|.v);.X..|
|00003d80| 7a 72 65 64 63 64 65 63 | 6f 64 65 28 72 70 2c 20 |zredcdec|ode(rp, |
|00003d90| 7a 31 2c 20 72 65 73 29 | 3b 0a 58 09 09 66 72 65 |z1, res)|;.X..fre|
|00003da0| 65 68 28 7a 31 2e 76 29 | 3b 0a 58 09 09 72 65 74 |eh(z1.v)|;.X..ret|
|00003db0| 75 72 6e 3b 0a 58 09 7d | 0a 58 0a 58 09 2f 2a 0a |urn;.X.}|.X.X./*.|
|00003dc0| 58 09 20 2a 20 4d 6f 64 | 75 6c 75 73 20 6f 72 20 |X. * Mod|ulus or |
|00003dd0| 70 6f 77 65 72 20 69 73 | 20 73 6d 61 6c 6c 20 65 |power is| small e|
|00003de0| 6e 6f 75 67 68 20 74 6f | 20 70 65 72 66 6f 72 6d |nough to| perform|
|00003df0| 20 74 68 65 20 70 6f 77 | 65 72 20 72 61 69 73 69 | the pow|er raisi|
|00003e00| 6e 67 0a 58 09 20 2a 20 | 64 69 72 65 63 74 6c 79 |ng.X. * |directly|
|00003e10| 2e 20 20 49 6e 69 74 69 | 61 6c 69 7a 65 20 74 68 |.  Initi|alize th|
|00003e20| 65 20 74 61 62 6c 65 20 | 6f 66 20 70 6f 77 65 72 |e table |of power|
|00003e30| 73 2e 0a 58 09 20 2a 2f | 0a 58 09 66 6f 72 20 28 |s..X. */|.X.for (|
|00003e40| 70 70 20 3d 20 26 6c 6f | 77 70 6f 77 65 72 73 5b |pp = &lo|wpowers[|
|00003e50| 32 5d 3b 20 70 70 20 3c | 20 26 6c 6f 77 70 6f 77 |2]; pp <| &lowpow|
|00003e60| 65 72 73 5b 50 4f 57 4e | 55 4d 53 5d 3b 20 70 70 |ers[POWN|UMS]; pp|
|00003e70| 2b 2b 29 0a 58 09 09 70 | 70 2d 3e 6c 65 6e 20 3d |++).X..p|p->len =|
|00003e80| 20 30 3b 0a 58 09 6c 6f | 77 70 6f 77 65 72 73 5b | 0;.X.lo|wpowers[|
|00003e90| 30 5d 20 3d 20 5f 6f 6e | 65 5f 3b 0a 58 09 6c 6f |0] = _on|e_;.X.lo|
|00003ea0| 77 70 6f 77 65 72 73 5b | 31 5d 20 3d 20 7a 31 3b |wpowers[|1] = z1;|
|00003eb0| 0a 58 09 61 6e 73 20 3d | 20 5f 6f 6e 65 5f 3b 0a |.X.ans =| _one_;.|
|00003ec0| 58 0a 58 09 68 70 20 3d | 20 26 7a 32 2e 76 5b 7a |X.X.hp =| &z2.v[z|
|00003ed0| 32 2e 6c 65 6e 20 2d 20 | 31 5d 3b 0a 58 09 63 75 |2.len - |1];.X.cu|
|00003ee0| 72 68 61 6c 66 20 3d 20 | 2a 68 70 3b 0a 58 09 63 |rhalf = |*hp;.X.c|
|00003ef0| 75 72 73 68 69 66 74 20 | 3d 20 42 41 53 45 42 20 |urshift |= BASEB |
|00003f00| 2d 20 50 4f 57 42 49 54 | 53 3b 0a 58 09 77 68 69 |- POWBIT|S;.X.whi|
|00003f10| 6c 65 20 28 63 75 72 73 | 68 69 66 74 20 26 26 20 |le (curs|hift && |
|00003f20| 28 28 63 75 72 68 61 6c | 66 20 3e 3e 20 63 75 72 |((curhal|f >> cur|
|00003f30| 73 68 69 66 74 29 20 3d | 3d 20 30 29 29 0a 58 09 |shift) =|= 0)).X.|
|00003f40| 09 63 75 72 73 68 69 66 | 74 20 2d 3d 20 50 4f 57 |.curshif|t -= POW|
|00003f50| 42 49 54 53 3b 0a 58 0a | 58 09 2f 2a 0a 58 09 20 |BITS;.X.|X./*.X. |
|00003f60| 2a 20 43 61 6c 63 75 6c | 61 74 65 20 74 68 65 20 |* Calcul|ate the |
|00003f70| 72 65 73 75 6c 74 20 62 | 79 20 65 78 61 6d 69 6e |result b|y examin|
|00003f80| 69 6e 67 20 74 68 65 20 | 70 6f 77 65 72 20 50 4f |ing the |power PO|
|00003f90| 57 42 49 54 53 20 62 69 | 74 73 20 61 74 20 61 20 |WBITS bi|ts at a |
|00003fa0| 74 69 6d 65 2c 0a 58 09 | 20 2a 20 61 6e 64 20 75 |time,.X.| * and u|
|00003fb0| 73 65 20 74 68 65 20 74 | 61 62 6c 65 20 6f 66 20 |se the t|able of |
|00003fc0| 6c 6f 77 20 70 6f 77 65 | 72 73 20 61 74 20 65 61 |low powe|rs at ea|
|00003fd0| 63 68 20 69 74 65 72 61 | 74 69 6f 6e 2e 0a 58 09 |ch itera|tion..X.|
|00003fe0| 20 2a 2f 0a 58 09 66 6f | 72 20 28 3b 3b 29 20 7b | */.X.fo|r (;;) {|
|00003ff0| 0a 58 09 09 63 75 72 70 | 6f 77 20 3d 20 28 63 75 |.X..curp|ow = (cu|
|00004000| 72 68 61 6c 66 20 3e 3e | 20 63 75 72 73 68 69 66 |rhalf >>| curshif|
|00004010| 74 29 20 26 20 28 50 4f | 57 4e 55 4d 53 20 2d 20 |t) & (PO|WNUMS - |
|00004020| 31 29 3b 0a 58 09 09 70 | 70 20 3d 20 26 6c 6f 77 |1);.X..p|p = &low|
|00004030| 70 6f 77 65 72 73 5b 63 | 75 72 70 6f 77 5d 3b 0a |powers[c|urpow];.|
|00004040| 58 0a 58 09 09 2f 2a 0a | 58 09 09 20 2a 20 49 66 |X.X../*.|X.. * If|
|00004050| 20 74 68 65 20 73 6d 61 | 6c 6c 20 70 6f 77 65 72 | the sma|ll power|
|00004060| 20 69 73 20 6e 6f 74 20 | 79 65 74 20 73 61 76 65 | is not |yet save|
|00004070| 64 20 69 6e 20 74 68 65 | 20 74 61 62 6c 65 2c 20 |d in the| table, |
|00004080| 74 68 65 6e 0a 58 09 09 | 20 2a 20 63 61 6c 63 75 |then.X..| * calcu|
|00004090| 6c 61 74 65 20 69 74 20 | 61 6e 64 20 72 65 6d 65 |late it |and reme|
|000040a0| 6d 62 65 72 20 69 74 20 | 69 6e 20 74 68 65 20 74 |mber it |in the t|
|000040b0| 61 62 6c 65 20 66 6f 72 | 20 66 75 74 75 72 65 20 |able for| future |
|000040c0| 75 73 65 2e 0a 58 09 09 | 20 2a 2f 0a 58 09 09 69 |use..X..| */.X..i|
|000040d0| 66 20 28 70 70 2d 3e 6c | 65 6e 20 3d 3d 20 30 29 |f (pp->l|en == 0)|
|000040e0| 20 7b 0a 58 09 09 09 69 | 66 20 28 63 75 72 70 6f | {.X...i|f (curpo|
|000040f0| 77 20 26 20 30 78 31 29 | 0a 58 09 09 09 09 7a 63 |w & 0x1)|.X....zc|
|00004100| 6f 70 79 28 7a 31 2c 20 | 26 6d 6f 64 70 6f 77 29 |opy(z1, |&modpow)|
|00004110| 3b 0a 58 09 09 09 65 6c | 73 65 0a 58 09 09 09 09 |;.X...el|se.X....|
|00004120| 6d 6f 64 70 6f 77 20 3d | 20 5f 6f 6e 65 5f 3b 0a |modpow =| _one_;.|
|00004130| 58 0a 58 09 09 09 66 6f | 72 20 28 63 75 72 62 69 |X.X...fo|r (curbi|
|00004140| 74 20 3d 20 30 78 32 3b | 20 63 75 72 62 69 74 20 |t = 0x2;| curbit |
|00004150| 3c 3d 20 63 75 72 70 6f | 77 3b 20 63 75 72 62 69 |<= curpo|w; curbi|
|00004160| 74 20 2a 3d 20 32 29 20 | 7b 0a 58 09 09 09 09 70 |t *= 2) |{.X....p|
|00004170| 70 20 3d 20 26 6c 6f 77 | 70 6f 77 65 72 73 5b 63 |p = &low|powers[c|
|00004180| 75 72 62 69 74 5d 3b 0a | 58 09 09 09 09 69 66 20 |urbit];.|X....if |
|00004190| 28 70 70 2d 3e 6c 65 6e | 20 3d 3d 20 30 29 20 7b |(pp->len| == 0) {|
|000041a0| 0a 58 09 09 09 09 09 7a | 73 71 75 61 72 65 28 6c |.X.....z|square(l|
|000041b0| 6f 77 70 6f 77 65 72 73 | 5b 63 75 72 62 69 74 2f |owpowers|[curbit/|
|000041c0| 32 5d 2c 20 26 74 65 6d | 70 29 3b 0a 58 09 09 09 |2], &tem|p);.X...|
|000041d0| 09 09 7a 6d 6f 64 28 74 | 65 6d 70 2c 20 7a 33 2c |..zmod(t|emp, z3,|
|000041e0| 20 70 70 29 3b 0a 58 09 | 09 09 09 09 66 72 65 65 | pp);.X.|....free|
|000041f0| 68 28 74 65 6d 70 2e 76 | 29 3b 0a 58 09 09 09 09 |h(temp.v|);.X....|
|00004200| 7d 0a 58 09 09 09 09 69 | 66 20 28 63 75 72 62 69 |}.X....i|f (curbi|
|00004210| 74 20 26 20 63 75 72 70 | 6f 77 29 20 7b 0a 58 09 |t & curp|ow) {.X.|
|00004220| 09 09 09 09 7a 6d 75 6c | 28 2a 70 70 2c 20 6d 6f |....zmul|(*pp, mo|
|00004230| 64 70 6f 77 2c 20 26 74 | 65 6d 70 29 3b 0a 58 09 |dpow, &t|emp);.X.|
|00004240| 09 09 09 09 66 72 65 65 | 68 28 6d 6f 64 70 6f 77 |....free|h(modpow|
|00004250| 2e 76 29 3b 0a 58 09 09 | 09 09 09 7a 6d 6f 64 28 |.v);.X..|...zmod(|
|00004260| 74 65 6d 70 2c 20 7a 33 | 2c 20 26 6d 6f 64 70 6f |temp, z3|, &modpo|
|00004270| 77 29 3b 0a 58 09 09 09 | 09 09 66 72 65 65 68 28 |w);.X...|..freeh(|
|00004280| 74 65 6d 70 2e 76 29 3b | 0a 58 09 09 09 09 7d 0a |temp.v);|.X....}.|
|00004290| 58 09 09 09 7d 0a 58 09 | 09 09 70 70 20 3d 20 26 |X...}.X.|..pp = &|
|000042a0| 6c 6f 77 70 6f 77 65 72 | 73 5b 63 75 72 70 6f 77 |lowpower|s[curpow|
|000042b0| 5d 3b 0a 58 09 09 09 2a | 70 70 20 3d 20 6d 6f 64 |];.X...*|pp = mod|
|000042c0| 70 6f 77 3b 0a 58 09 09 | 7d 0a 58 0a 58 09 09 2f |pow;.X..|}.X.X../|
|000042d0| 2a 0a 58 09 09 20 2a 20 | 49 66 20 74 68 65 20 70 |*.X.. * |If the p|
|000042e0| 6f 77 65 72 20 69 73 20 | 6e 6f 6e 7a 65 72 6f 2c |ower is |nonzero,|
|000042f0| 20 74 68 65 6e 20 61 63 | 63 75 6d 75 6c 61 74 65 | then ac|cumulate|
|00004300| 20 74 68 65 20 73 6d 61 | 6c 6c 20 70 6f 77 65 72 | the sma|ll power|
|00004310| 0a 58 09 09 20 2a 20 69 | 6e 74 6f 20 74 68 65 20 |.X.. * i|nto the |
|00004320| 72 65 73 75 6c 74 2e 0a | 58 09 09 20 2a 2f 0a 58 |result..|X.. */.X|
|00004330| 09 09 69 66 20 28 63 75 | 72 70 6f 77 29 20 7b 0a |..if (cu|rpow) {.|
|00004340| 58 09 09 09 7a 6d 75 6c | 28 61 6e 73 2c 20 2a 70 |X...zmul|(ans, *p|
|00004350| 70 2c 20 26 74 65 6d 70 | 29 3b 0a 58 09 09 09 66 |p, &temp|);.X...f|
|00004360| 72 65 65 68 28 61 6e 73 | 2e 76 29 3b 0a 58 09 09 |reeh(ans|.v);.X..|
|00004370| 09 7a 6d 6f 64 28 74 65 | 6d 70 2c 20 7a 33 2c 20 |.zmod(te|mp, z3, |
|00004380| 26 61 6e 73 29 3b 0a 58 | 09 09 09 66 72 65 65 68 |&ans);.X|...freeh|
|00004390| 28 74 65 6d 70 2e 76 29 | 3b 0a 58 09 09 7d 0a 58 |(temp.v)|;.X..}.X|
|000043a0| 0a 58 09 09 2f 2a 0a 58 | 09 09 20 2a 20 53 65 6c |.X../*.X|.. * Sel|
|000043b0| 65 63 74 20 74 68 65 20 | 6e 65 78 74 20 50 4f 57 |ect the |next POW|
|000043c0| 42 49 54 53 20 62 69 74 | 73 20 6f 66 20 74 68 65 |BITS bit|s of the|
|000043d0| 20 70 6f 77 65 72 2c 20 | 69 66 20 74 68 65 72 65 | power, |if there|
|000043e0| 20 69 73 0a 58 09 09 20 | 2a 20 61 6e 79 20 6d 6f | is.X.. |* any mo|
|000043f0| 72 65 20 74 6f 20 67 65 | 6e 65 72 61 74 65 2e 0a |re to ge|nerate..|
|00004400| 58 09 09 20 2a 2f 0a 58 | 09 09 63 75 72 73 68 69 |X.. */.X|..curshi|
|00004410| 66 74 20 2d 3d 20 50 4f | 57 42 49 54 53 3b 0a 58 |ft -= PO|WBITS;.X|
|00004420| 09 09 69 66 20 28 63 75 | 72 73 68 69 66 74 20 3c |..if (cu|rshift <|
|00004430| 20 30 29 20 7b 0a 58 09 | 09 09 69 66 20 28 68 70 | 0) {.X.|..if (hp|
|00004440| 2d 2d 20 3d 3d 20 7a 32 | 2e 76 29 0a 58 09 09 09 |-- == z2|.v).X...|
|00004450| 09 62 72 65 61 6b 3b 0a | 58 09 09 09 63 75 72 68 |.break;.|X...curh|
|00004460| 61 6c 66 20 3d 20 2a 68 | 70 3b 0a 58 09 09 09 63 |alf = *h|p;.X...c|
|00004470| 75 72 73 68 69 66 74 20 | 3d 20 42 41 53 45 42 20 |urshift |= BASEB |
|00004480| 2d 20 50 4f 57 42 49 54 | 53 3b 0a 58 09 09 7d 0a |- POWBIT|S;.X..}.|
|00004490| 58 0a 58 09 09 2f 2a 0a | 58 09 09 20 2a 20 53 71 |X.X../*.|X.. * Sq|
|000044a0| 75 61 72 65 20 74 68 65 | 20 72 65 73 75 6c 74 20 |uare the| result |
|000044b0| 50 4f 57 42 49 54 53 20 | 74 69 6d 65 73 20 74 6f |POWBITS |times to|
|000044c0| 20 6d 61 6b 65 20 72 6f | 6f 6d 20 66 6f 72 20 74 | make ro|om for t|
|000044d0| 68 65 20 6e 65 78 74 0a | 58 09 09 20 2a 20 63 68 |he next.|X.. * ch|
|000044e0| 75 6e 6b 20 6f 66 20 62 | 69 74 73 2e 0a 58 09 09 |unk of b|its..X..|
|000044f0| 20 2a 2f 0a 58 09 09 66 | 6f 72 20 28 69 20 3d 20 | */.X..f|or (i = |
|00004500| 30 3b 20 69 20 3c 20 50 | 4f 57 42 49 54 53 3b 20 |0; i < P|OWBITS; |
|00004510| 69 2b 2b 29 20 7b 0a 58 | 09 09 09 7a 73 71 75 61 |i++) {.X|...zsqua|
|00004520| 72 65 28 61 6e 73 2c 20 | 26 74 65 6d 70 29 3b 0a |re(ans, |&temp);.|
|00004530| 58 09 09 09 66 72 65 65 | 68 28 61 6e 73 2e 76 29 |X...free|h(ans.v)|
|00004540| 3b 0a 58 09 09 09 7a 6d | 6f 64 28 74 65 6d 70 2c |;.X...zm|od(temp,|
|00004550| 20 7a 33 2c 20 26 61 6e | 73 29 3b 0a 58 09 09 09 | z3, &an|s);.X...|
|00004560| 66 72 65 65 68 28 74 65 | 6d 70 2e 76 29 3b 0a 58 |freeh(te|mp.v);.X|
|00004570| 09 09 7d 0a 58 09 7d 0a | 58 0a 58 09 66 6f 72 20 |..}.X.}.|X.X.for |
|00004580| 28 70 70 20 3d 20 26 6c | 6f 77 70 6f 77 65 72 73 |(pp = &l|owpowers|
|00004590| 5b 32 5d 3b 20 70 70 20 | 3c 20 26 6c 6f 77 70 6f |[2]; pp |< &lowpo|
|000045a0| 77 65 72 73 5b 50 4f 57 | 4e 55 4d 53 5d 3b 20 70 |wers[POW|NUMS]; p|
|000045b0| 70 2b 2b 29 20 7b 0a 58 | 09 09 69 66 20 28 70 70 |p++) {.X|..if (pp|
|000045c0| 2d 3e 6c 65 6e 29 0a 58 | 09 09 09 66 72 65 65 68 |->len).X|...freeh|
|000045d0| 28 70 70 2d 3e 76 29 3b | 0a 58 09 7d 0a 58 09 2a |(pp->v);|.X.}.X.*|
|000045e0| 72 65 73 20 3d 20 61 6e | 73 3b 0a 58 7d 0a 58 0a |res = an|s;.X}.X.|
|000045f0| 58 0a 58 2f 2a 0a 58 20 | 2a 20 49 6e 69 74 69 61 |X.X/*.X |* Initia|
|00004600| 6c 69 7a 65 20 74 68 65 | 20 52 45 44 43 20 61 6c |lize the| REDC al|
|00004610| 67 6f 72 69 74 68 6d 20 | 66 6f 72 20 61 20 70 61 |gorithm |for a pa|
|00004620| 72 74 69 63 75 6c 61 72 | 20 6d 6f 64 75 6c 75 73 |rticular| modulus|
|00004630| 2c 0a 58 20 2a 20 72 65 | 74 75 72 6e 69 6e 67 20 |,.X * re|turning |
|00004640| 61 20 70 6f 69 6e 74 65 | 72 20 74 6f 20 61 20 73 |a pointe|r to a s|
|00004650| 74 72 75 63 74 75 72 65 | 20 74 68 61 74 20 69 73 |tructure| that is|
|00004660| 20 75 73 65 64 20 66 6f | 72 20 6f 74 68 65 72 0a | used fo|r other.|
|00004670| 58 20 2a 20 52 45 44 43 | 20 63 61 6c 6c 73 2e 20 |X * REDC| calls. |
|00004680| 20 41 6e 20 65 72 72 6f | 72 20 69 73 20 67 65 6e | An erro|r is gen|
|00004690| 65 72 61 74 65 64 20 69 | 66 20 74 68 65 20 73 74 |erated i|f the st|
|000046a0| 72 75 63 74 75 72 65 20 | 63 61 6e 6e 6f 74 0a 58 |ructure |cannot.X|
|000046b0| 20 2a 20 62 65 20 61 6c | 6c 6f 63 61 74 65 64 2e | * be al|located.|
|000046c0| 20 20 54 68 65 20 6d 6f | 64 75 6c 75 73 20 6d 75 |  The mo|dulus mu|
|000046d0| 73 74 20 62 65 20 6f 64 | 64 20 61 6e 64 20 70 6f |st be od|d and po|
|000046e0| 73 69 74 69 76 65 2e 0a | 58 20 2a 2f 0a 58 52 45 |sitive..|X */.XRE|
|000046f0| 44 43 20 2a 0a 58 7a 72 | 65 64 63 61 6c 6c 6f 63 |DC *.Xzr|edcalloc|
|00004700| 28 7a 31 29 0a 58 09 5a | 56 41 4c 55 45 20 7a 31 |(z1).X.Z|VALUE z1|
|00004710| 3b 09 09 2f 2a 20 6d 6f | 64 75 6c 75 73 20 74 6f |;../* mo|dulus to|
|00004720| 20 69 6e 69 74 69 61 6c | 69 7a 65 20 66 6f 72 20 | initial|ize for |
|00004730| 2a 2f 0a 58 7b 0a 58 09 | 52 45 44 43 20 2a 72 70 |*/.X{.X.|REDC *rp|
|00004740| 3b 09 09 2f 2a 20 52 45 | 44 43 20 69 6e 66 6f 72 |;../* RE|DC infor|
|00004750| 6d 61 74 69 6f 6e 20 2a | 2f 0a 58 09 5a 56 41 4c |mation *|/.X.ZVAL|
|00004760| 55 45 20 74 6d 70 3b 0a | 58 09 6c 6f 6e 67 20 62 |UE tmp;.|X.long b|
|00004770| 69 74 3b 0a 58 0a 58 09 | 69 66 20 28 69 73 65 76 |it;.X.X.|if (isev|
|00004780| 65 6e 28 7a 31 29 20 7c | 7c 20 69 73 6e 65 67 28 |en(z1) ||| isneg(|
|00004790| 7a 31 29 29 0a 58 09 09 | 65 72 72 6f 72 28 22 52 |z1)).X..|error("R|
|000047a0| 45 44 43 20 72 65 71 75 | 69 72 65 73 20 70 6f 73 |EDC requ|ires pos|
|000047b0| 69 74 69 76 65 20 6f 64 | 64 20 6d 6f 64 75 6c 75 |itive od|d modulu|
|000047c0| 73 22 29 3b 0a 58 0a 58 | 09 72 70 20 3d 20 28 52 |s");.X.X|.rp = (R|
|000047d0| 45 44 43 20 2a 29 20 6d | 61 6c 6c 6f 63 28 73 69 |EDC *) m|alloc(si|
|000047e0| 7a 65 6f 66 28 52 45 44 | 43 29 29 3b 0a 58 09 69 |zeof(RED|C));.X.i|
|000047f0| 66 20 28 72 70 20 3d 3d | 20 4e 55 4c 4c 29 0a 58 |f (rp ==| NULL).X|
|00004800| 09 09 65 72 72 6f 72 28 | 22 43 61 6e 6e 6f 74 20 |..error(|"Cannot |
|00004810| 61 6c 6c 6f 63 61 74 65 | 20 52 45 44 43 20 73 74 |allocate| REDC st|
|00004820| 72 75 63 74 75 72 65 22 | 29 3b 0a 58 0a 58 09 2f |ructure"|);.X.X./|
|00004830| 2a 0a 58 09 20 2a 20 52 | 6f 75 6e 64 20 75 70 20 |*.X. * R|ound up |
|00004840| 74 68 65 20 62 69 6e 61 | 72 79 20 6d 6f 64 75 6c |the bina|ry modul|
|00004850| 75 73 20 74 6f 20 74 68 | 65 20 6e 65 78 74 20 70 |us to th|e next p|
|00004860| 6f 77 65 72 20 6f 66 20 | 74 77 6f 0a 58 09 20 2a |ower of |two.X. *|
|00004870| 20 77 68 69 63 68 20 69 | 73 20 61 74 20 61 20 77 | which i|s at a w|
|00004880| 6f 72 64 20 62 6f 75 6e | 64 61 72 79 2e 20 20 54 |ord boun|dary.  T|
|00004890| 68 65 6e 20 74 68 65 20 | 73 68 69 66 74 20 61 6e |hen the |shift an|
|000048a0| 64 20 6d 6f 64 75 6c 6f | 0a 58 09 20 2a 20 6f 70 |d modulo|.X. * op|
|000048b0| 65 72 61 74 69 6f 6e 73 | 20 6d 6f 64 20 74 68 65 |erations| mod the|
|000048c0| 20 62 69 6e 61 72 79 20 | 6d 6f 64 75 6c 75 73 20 | binary |modulus |
|000048d0| 63 61 6e 20 62 65 20 64 | 6f 6e 65 20 76 65 72 79 |can be d|one very|
|000048e0| 20 63 68 65 61 70 6c 79 | 2e 0a 58 09 20 2a 20 43 | cheaply|..X. * C|
|000048f0| 61 6c 63 75 6c 61 74 65 | 20 74 68 65 20 52 45 44 |alculate| the RED|
|00004900| 43 20 66 6f 72 6d 61 74 | 20 66 6f 72 20 74 68 65 |C format| for the|
|00004910| 20 6e 75 6d 62 65 72 20 | 31 20 66 6f 72 20 66 75 | number |1 for fu|
|00004920| 74 75 72 65 20 75 73 65 | 2e 0a 58 09 20 2a 2f 0a |ture use|..X. */.|
|00004930| 58 09 62 69 74 20 3d 20 | 7a 68 69 67 68 62 69 74 |X.bit = |zhighbit|
|00004940| 28 7a 31 29 20 2b 20 31 | 3b 0a 58 09 69 66 20 28 |(z1) + 1|;.X.if (|
|00004950| 62 69 74 20 25 20 42 41 | 53 45 42 29 0a 58 09 09 |bit % BA|SEB).X..|
|00004960| 62 69 74 20 2b 3d 20 28 | 42 41 53 45 42 20 2d 20 |bit += (|BASEB - |
|00004970| 28 62 69 74 20 25 20 42 | 41 53 45 42 29 29 3b 0a |(bit % B|ASEB));.|
|00004980| 58 09 7a 63 6f 70 79 28 | 7a 31 2c 20 26 72 70 2d |X.zcopy(|z1, &rp-|
|00004990| 3e 6d 6f 64 29 3b 0a 58 | 09 7a 62 69 74 76 61 6c |>mod);.X|.zbitval|
|000049a0| 75 65 28 62 69 74 2c 20 | 26 74 6d 70 29 3b 0a 58 |ue(bit, |&tmp);.X|
|000049b0| 09 7a 31 2e 73 69 67 6e | 20 3d 20 31 3b 0a 58 09 |.z1.sign| = 1;.X.|
|000049c0| 28 76 6f 69 64 29 20 7a | 6d 6f 64 69 6e 76 28 7a |(void) z|modinv(z|
|000049d0| 31 2c 20 74 6d 70 2c 20 | 26 72 70 2d 3e 69 6e 76 |1, tmp, |&rp->inv|
|000049e0| 29 3b 0a 58 09 7a 6d 6f | 64 28 74 6d 70 2c 20 72 |);.X.zmo|d(tmp, r|
|000049f0| 70 2d 3e 6d 6f 64 2c 20 | 26 72 70 2d 3e 6f 6e 65 |p->mod, |&rp->one|
|00004a00| 29 3b 0a 58 09 66 72 65 | 65 68 28 74 6d 70 2e 76 |);.X.fre|eh(tmp.v|
|00004a10| 29 3b 0a 58 09 72 70 2d | 3e 6c 65 6e 20 3d 20 62 |);.X.rp-|>len = b|
|00004a20| 69 74 20 2f 20 42 41 53 | 45 42 3b 0a 58 09 72 65 |it / BAS|EB;.X.re|
|00004a30| 74 75 72 6e 20 72 70 3b | 0a 58 7d 0a 58 0a 58 0a |turn rp;|.X}.X.X.|
|00004a40| 58 2f 2a 0a 58 20 2a 20 | 46 72 65 65 20 61 6e 79 |X/*.X * |Free any|
|00004a50| 20 6e 75 6d 62 65 72 73 | 20 61 73 73 6f 63 69 61 | numbers| associa|
|00004a60| 74 65 64 20 77 69 74 68 | 20 74 68 65 20 73 70 65 |ted with| the spe|
|00004a70| 63 69 66 69 65 64 20 52 | 45 44 43 20 73 74 72 75 |cified R|EDC stru|
|00004a80| 63 74 75 72 65 2c 0a 58 | 20 2a 20 61 6e 64 20 74 |cture,.X| * and t|
|00004a90| 68 65 6e 20 74 68 65 20 | 52 45 44 43 20 73 74 72 |hen the |REDC str|
|00004aa0| 75 63 74 75 72 65 20 69 | 74 73 65 6c 66 2e 0a 58 |ucture i|tself..X|
|00004ab0| 20 2a 2f 0a 58 76 6f 69 | 64 0a 58 7a 72 65 64 63 | */.Xvoi|d.Xzredc|
|00004ac0| 66 72 65 65 28 72 70 29 | 0a 58 09 52 45 44 43 20 |free(rp)|.X.REDC |
|00004ad0| 2a 72 70 3b 09 09 2f 2a | 20 52 45 44 43 20 69 6e |*rp;../*| REDC in|
|00004ae0| 66 6f 72 6d 61 74 69 6f | 6e 20 74 6f 20 62 65 20 |formatio|n to be |
|00004af0| 63 6c 65 61 72 65 64 20 | 2a 2f 0a 58 7b 0a 58 09 |cleared |*/.X{.X.|
|00004b00| 66 72 65 65 68 28 72 70 | 2d 3e 6d 6f 64 2e 76 29 |freeh(rp|->mod.v)|
|00004b10| 3b 0a 58 09 66 72 65 65 | 68 28 72 70 2d 3e 69 6e |;.X.free|h(rp->in|
|00004b20| 76 2e 76 29 3b 0a 58 09 | 66 72 65 65 68 28 72 70 |v.v);.X.|freeh(rp|
|00004b30| 2d 3e 6f 6e 65 2e 76 29 | 3b 0a 58 09 66 72 65 65 |->one.v)|;.X.free|
|00004b40| 28 72 70 29 3b 0a 58 7d | 0a 58 0a 58 0a 58 2f 2a |(rp);.X}|.X.X.X/*|
|00004b50| 0a 58 20 2a 20 43 6f 6e | 76 65 72 74 20 61 20 6e |.X * Con|vert a n|
|00004b60| 6f 72 6d 61 6c 20 6e 75 | 6d 62 65 72 20 69 6e 74 |ormal nu|mber int|
|00004b70| 6f 20 74 68 65 20 73 70 | 65 63 69 66 69 65 64 20 |o the sp|ecified |
|00004b80| 52 45 44 43 20 66 6f 72 | 6d 61 74 2e 0a 58 20 2a |REDC for|mat..X *|
|00004b90| 20 54 68 65 20 6e 75 6d | 62 65 72 20 74 6f 20 62 | The num|ber to b|
|00004ba0| 65 20 63 6f 6e 76 65 72 | 74 65 64 20 63 61 6e 20 |e conver|ted can |
|00004bb0| 62 65 20 6e 65 67 61 74 | 69 76 65 20 6f 72 20 6f |be negat|ive or o|
|00004bc0| 75 74 20 6f 66 20 6d 6f | 64 75 6c 6f 20 72 61 6e |ut of mo|dulo ran|
|00004bd0| 67 65 2e 0a 58 20 2a 20 | 54 68 65 20 72 65 73 75 |ge..X * |The resu|
|00004be0| 6c 74 69 6e 67 20 6e 75 | 6d 62 65 72 20 63 61 6e |lting nu|mber can|
|00004bf0| 20 62 65 20 75 73 65 64 | 20 66 6f 72 20 6d 75 6c | be used| for mul|
|00004c00| 74 69 70 6c 79 69 6e 67 | 2c 20 61 64 64 69 6e 67 |tiplying|, adding|
|00004c10| 2c 20 73 75 62 74 72 61 | 63 74 69 6e 67 2c 0a 58 |, subtra|cting,.X|
|00004c20| 20 2a 20 6f 72 20 63 6f | 6d 70 61 72 69 6e 67 20 | * or co|mparing |
|00004c30| 77 69 74 68 20 61 6e 79 | 20 6f 74 68 65 72 20 73 |with any| other s|
|00004c40| 75 63 68 20 63 6f 6e 76 | 65 72 74 65 64 20 6e 75 |uch conv|erted nu|
|00004c50| 6d 62 65 72 73 2c 20 61 | 73 20 69 66 20 74 68 65 |mbers, a|s if the|
|00004c60| 20 6e 75 6d 62 65 72 73 | 0a 58 20 2a 20 77 65 72 | numbers|.X * wer|
|00004c70| 65 20 62 65 69 6e 67 20 | 63 61 6c 63 75 6c 61 74 |e being |calculat|
|00004c80| 65 64 20 6d 6f 64 75 6c | 6f 20 74 68 65 20 6e 75 |ed modul|o the nu|
|00004c90| 6d 62 65 72 20 77 68 69 | 63 68 20 69 6e 69 74 69 |mber whi|ch initi|
|00004ca0| 61 6c 69 7a 65 64 20 74 | 68 65 20 52 45 44 43 0a |alized t|he REDC.|
|00004cb0| 58 20 2a 20 69 6e 66 6f | 72 6d 61 74 69 6f 6e 2e |X * info|rmation.|
|00004cc0| 20 20 57 68 65 6e 20 74 | 68 65 20 66 69 6e 61 6c |  When t|he final|
|00004cd0| 20 76 61 6c 75 65 20 69 | 73 20 75 6e 63 6f 6e 76 | value i|s unconv|
|00004ce0| 65 72 74 65 64 2c 20 74 | 68 65 20 72 65 73 75 6c |erted, t|he resul|
|00004cf0| 74 20 69 73 20 74 68 65 | 0a 58 20 2a 20 73 61 6d |t is the|.X * sam|
|00004d00| 65 20 61 73 20 69 66 20 | 74 68 65 20 75 73 75 61 |e as if |the usua|
|00004d10| 6c 20 6f 70 65 72 61 74 | 69 6f 6e 73 20 77 65 72 |l operat|ions wer|
|00004d20| 65 20 64 6f 6e 65 20 77 | 69 74 68 20 74 68 65 20 |e done w|ith the |
|00004d30| 6f 72 69 67 69 6e 61 6c | 20 6e 75 6d 62 65 72 73 |original| numbers|
|00004d40| 2e 0a 58 20 2a 2f 0a 58 | 76 6f 69 64 0a 58 7a 72 |..X */.X|void.Xzr|
|00004d50| 65 64 63 65 6e 63 6f 64 | 65 28 72 70 2c 20 7a 31 |edcencod|e(rp, z1|
|00004d60| 2c 20 72 65 73 29 0a 58 | 09 52 45 44 43 20 2a 72 |, res).X|.REDC *r|
|00004d70| 70 3b 09 09 2f 2a 20 52 | 45 44 43 20 69 6e 66 6f |p;../* R|EDC info|
|00004d80| 72 6d 61 74 69 6f 6e 20 | 2a 2f 0a 58 09 5a 56 41 |rmation |*/.X.ZVA|
|00004d90| 4c 55 45 20 7a 31 3b 09 | 09 2f 2a 20 6e 75 6d 62 |LUE z1;.|./* numb|
|00004da0| 65 72 20 74 6f 20 62 65 | 20 63 6f 6e 76 65 72 74 |er to be| convert|
|00004db0| 65 64 20 2a 2f 0a 58 09 | 5a 56 41 4c 55 45 20 2a |ed */.X.|ZVALUE *|
|00004dc0| 72 65 73 3b 09 09 2f 2a | 20 72 65 74 75 72 6e 65 |res;../*| returne|
|00004dd0| 64 20 63 6f 6e 76 65 72 | 74 65 64 20 6e 75 6d 62 |d conver|ted numb|
|00004de0| 65 72 20 2a 2f 0a 58 7b | 0a 58 09 5a 56 41 4c 55 |er */.X{|.X.ZVALU|
|00004df0| 45 20 74 6d 70 31 2c 20 | 74 6d 70 32 3b 0a 58 0a |E tmp1, |tmp2;.X.|
|00004e00| 58 09 2f 2a 0a 58 09 20 | 2a 20 48 61 6e 64 6c 65 |X./*.X. |* Handle|
|00004e10| 20 74 68 65 20 63 61 73 | 65 73 20 30 2c 20 31 2c | the cas|es 0, 1,|
|00004e20| 20 2d 31 2c 20 61 6e 64 | 20 32 20 73 70 65 63 69 | -1, and| 2 speci|
|00004e30| 61 6c 6c 79 20 73 69 6e | 63 65 20 74 68 65 73 65 |ally sin|ce these|
|00004e40| 20 61 72 65 0a 58 09 20 | 2a 20 65 61 73 79 20 74 | are.X. |* easy t|
|00004e50| 6f 20 63 61 6c 63 75 6c | 61 74 65 2e 20 20 5a 65 |o calcul|ate.  Ze|
|00004e60| 72 6f 20 74 72 61 6e 73 | 66 6f 72 6d 73 20 74 6f |ro trans|forms to|
|00004e70| 20 7a 65 72 6f 2c 20 61 | 6e 64 20 74 68 65 20 6f | zero, a|nd the o|
|00004e80| 74 68 65 72 73 0a 58 09 | 20 2a 20 63 61 6e 20 62 |thers.X.| * can b|
|00004e90| 65 20 6f 62 74 61 69 6e | 65 64 20 66 72 6f 6d 20 |e obtain|ed from |
|00004ea0| 74 68 65 20 70 72 65 63 | 6f 6d 70 75 74 65 64 20 |the prec|omputed |
|00004eb0| 52 45 44 43 20 66 6f 72 | 6d 61 74 20 66 6f 72 20 |REDC for|mat for |
|00004ec0| 31 20 73 69 6e 63 65 0a | 58 09 20 2a 20 61 64 64 |1 since.|X. * add|
|00004ed0| 69 74 69 6f 6e 20 61 6e | 64 20 73 75 62 74 72 61 |ition an|d subtra|
|00004ee0| 63 74 69 6f 6e 20 61 63 | 74 20 6e 6f 72 6d 61 6c |ction ac|t normal|
|00004ef0| 6c 79 20 66 6f 72 20 52 | 45 44 43 20 66 6f 72 6d |ly for R|EDC form|
|00004f00| 61 74 20 6e 75 6d 62 65 | 72 73 2e 0a 58 09 20 2a |at numbe|rs..X. *|
|00004f10| 2f 0a 58 09 69 66 20 28 | 69 73 7a 65 72 6f 28 7a |/.X.if (|iszero(z|
|00004f20| 31 29 29 20 7b 0a 58 09 | 09 2a 72 65 73 20 3d 20 |1)) {.X.|.*res = |
|00004f30| 5f 7a 65 72 6f 5f 3b 0a | 58 09 09 72 65 74 75 72 |_zero_;.|X..retur|
|00004f40| 6e 3b 0a 58 09 7d 0a 58 | 09 69 66 20 28 69 73 6f |n;.X.}.X|.if (iso|
|00004f50| 6e 65 28 7a 31 29 29 20 | 7b 0a 58 09 09 7a 63 6f |ne(z1)) |{.X..zco|
|00004f60| 70 79 28 72 70 2d 3e 6f | 6e 65 2c 20 72 65 73 29 |py(rp->o|ne, res)|
|00004f70| 3b 0a 58 09 09 72 65 74 | 75 72 6e 3b 0a 58 09 7d |;.X..ret|urn;.X.}|
|00004f80| 0a 58 09 69 66 20 28 69 | 73 75 6e 69 74 28 7a 31 |.X.if (i|sunit(z1|
|00004f90| 29 29 20 7b 0a 58 09 09 | 7a 73 75 62 28 72 70 2d |)) {.X..|zsub(rp-|
|00004fa0| 3e 6d 6f 64 2c 20 72 70 | 2d 3e 6f 6e 65 2c 20 72 |>mod, rp|->one, r|
|00004fb0| 65 73 29 3b 0a 58 09 09 | 72 65 74 75 72 6e 3b 0a |es);.X..|return;.|
|00004fc0| 58 09 7d 0a 58 09 69 66 | 20 28 69 73 74 77 6f 28 |X.}.X.if| (istwo(|
|00004fd0| 7a 31 29 29 20 7b 0a 58 | 09 09 7a 61 64 64 28 72 |z1)) {.X|..zadd(r|
|00004fe0| 70 2d 3e 6f 6e 65 2c 20 | 72 70 2d 3e 6f 6e 65 2c |p->one, |rp->one,|
|00004ff0| 20 26 74 6d 70 31 29 3b | 0a 58 09 09 69 66 20 28 | &tmp1);|.X..if (|
|00005000| 7a 72 65 6c 28 74 6d 70 | 31 2c 20 72 70 2d 3e 6d |zrel(tmp|1, rp->m|
|00005010| 6f 64 29 20 3c 20 30 29 | 20 7b 0a 58 09 09 09 2a |od) < 0)| {.X...*|
|00005020| 72 65 73 20 3d 20 74 6d | 70 31 3b 0a 58 09 09 09 |res = tm|p1;.X...|
|00005030| 72 65 74 75 72 6e 3b 0a | 58 09 09 7d 0a 58 09 09 |return;.|X..}.X..|
|00005040| 7a 73 75 62 28 74 6d 70 | 31 2c 20 72 70 2d 3e 6d |zsub(tmp|1, rp->m|
|00005050| 6f 64 2c 20 72 65 73 29 | 3b 0a 58 09 09 66 72 65 |od, res)|;.X..fre|
|00005060| 65 68 28 74 6d 70 31 2e | 76 29 3b 0a 58 09 09 72 |eh(tmp1.|v);.X..r|
|00005070| 65 74 75 72 6e 3b 0a 58 | 09 7d 0a 58 0a 58 09 2f |eturn;.X|.}.X.X./|
|00005080| 2a 0a 58 09 20 2a 20 4e | 6f 74 20 61 20 74 72 69 |*.X. * N|ot a tri|
|00005090| 76 69 61 6c 20 6e 75 6d | 62 65 72 20 74 6f 20 63 |vial num|ber to c|
|000050a0| 6f 6e 76 65 72 74 2c 20 | 73 6f 20 64 6f 20 74 68 |onvert, |so do th|
|000050b0| 65 20 66 75 6c 6c 20 74 | 72 61 6e 73 66 6f 72 6d |e full t|ransform|
|000050c0| 61 74 69 6f 6e 2e 0a 58 | 09 20 2a 20 43 6f 6e 76 |ation..X|. * Conv|
|000050d0| 65 72 74 20 6e 65 67 61 | 74 69 76 65 20 6e 75 6d |ert nega|tive num|
|000050e0| 62 65 72 73 20 74 6f 20 | 70 6f 73 69 74 69 76 65 |bers to |positive|
|000050f0| 20 6e 75 6d 62 65 72 73 | 20 62 65 66 6f 72 65 20 | numbers| before |
|00005100| 63 6f 6e 76 65 72 74 69 | 6e 67 2e 0a 58 09 20 2a |converti|ng..X. *|
|00005110| 2f 0a 58 09 74 6d 70 31 | 2e 6c 65 6e 20 3d 20 30 |/.X.tmp1|.len = 0|
|00005120| 3b 0a 58 09 69 66 20 28 | 69 73 6e 65 67 28 7a 31 |;.X.if (|isneg(z1|
|00005130| 29 29 20 7b 0a 58 09 09 | 7a 6d 6f 64 28 7a 31 2c |)) {.X..|zmod(z1,|
|00005140| 20 72 70 2d 3e 6d 6f 64 | 2c 20 26 74 6d 70 31 29 | rp->mod|, &tmp1)|
|00005150| 3b 0a 58 09 09 7a 31 20 | 3d 20 74 6d 70 31 3b 0a |;.X..z1 |= tmp1;.|
|00005160| 58 09 7d 0a 58 09 7a 73 | 68 69 66 74 28 7a 31 2c |X.}.X.zs|hift(z1,|
|00005170| 20 72 70 2d 3e 6c 65 6e | 20 2a 20 42 41 53 45 42 | rp->len| * BASEB|
|00005180| 2c 20 26 74 6d 70 32 29 | 3b 0a 58 09 69 66 20 28 |, &tmp2)|;.X.if (|
|00005190| 74 6d 70 31 2e 6c 65 6e | 29 0a 58 09 09 66 72 65 |tmp1.len|).X..fre|
|000051a0| 65 68 28 74 6d 70 31 2e | 76 29 3b 0a 58 09 7a 6d |eh(tmp1.|v);.X.zm|
|000051b0| 6f 64 28 74 6d 70 32 2c | 20 72 70 2d 3e 6d 6f 64 |od(tmp2,| rp->mod|
|000051c0| 2c 20 72 65 73 29 3b 0a | 58 09 66 72 65 65 68 28 |, res);.|X.freeh(|
|000051d0| 74 6d 70 32 2e 76 29 3b | 0a 58 7d 0a 58 0a 58 0a |tmp2.v);|.X}.X.X.|
|000051e0| 58 2f 2a 0a 58 20 2a 20 | 54 68 65 20 52 45 44 43 |X/*.X * |The REDC|
|000051f0| 20 61 6c 67 6f 72 69 74 | 68 6d 20 75 73 65 64 20 | algorit|hm used |
|00005200| 74 6f 20 63 6f 6e 76 65 | 72 74 20 6e 75 6d 62 65 |to conve|rt numbe|
|00005210| 72 73 20 6f 75 74 20 6f | 66 20 52 45 44 43 20 66 |rs out o|f REDC f|
|00005220| 6f 72 6d 61 74 20 61 6e | 64 20 61 6c 73 6f 0a 58 |ormat an|d also.X|
|00005230| 20 2a 20 75 73 65 64 20 | 61 66 74 65 72 20 6d 75 | * used |after mu|
|00005240| 6c 74 69 70 6c 69 63 61 | 74 69 6f 6e 20 6f 66 20 |ltiplica|tion of |
|00005250| 74 77 6f 20 52 45 44 43 | 20 6e 75 6d 62 65 72 73 |two REDC| numbers|
|00005260| 2e 20 20 55 73 69 6e 67 | 20 74 68 69 73 20 72 6f |.  Using| this ro|
|00005270| 75 74 69 6e 65 0a 58 20 | 2a 20 61 76 6f 69 64 73 |utine.X |* avoids|
|00005280| 20 61 6e 79 20 64 69 76 | 69 64 65 73 2c 20 72 65 | any div|ides, re|
|00005290| 70 6c 61 63 69 6e 67 20 | 74 68 65 20 64 69 76 69 |placing |the divi|
|000052a0| 64 65 20 62 79 20 74 77 | 6f 20 6d 75 6c 74 69 70 |de by tw|o multip|
|000052b0| 6c 69 63 61 74 69 6f 6e | 73 2e 0a 58 20 2a 20 49 |lication|s..X * I|
|000052c0| 66 20 74 68 65 20 6e 75 | 6d 62 65 72 73 20 61 72 |f the nu|mbers ar|
|000052d0| 65 20 76 65 72 79 20 6c | 61 72 67 65 2c 20 74 68 |e very l|arge, th|
|000052e0| 65 6e 20 74 68 65 73 65 | 20 74 77 6f 20 6d 75 6c |en these| two mul|
|000052f0| 74 69 70 6c 69 65 73 20 | 77 69 6c 6c 20 62 65 0a |tiplies |will be.|
|00005300| 58 20 2a 20 71 75 69 63 | 6b 65 72 20 74 68 61 6e |X * quic|ker than|
|00005310| 20 74 68 65 20 64 69 76 | 69 64 65 2c 20 73 69 6e | the div|ide, sin|
|00005320| 63 65 20 64 69 76 69 64 | 69 6e 67 20 69 73 20 68 |ce divid|ing is h|
|00005330| 61 72 64 65 72 20 74 68 | 61 6e 20 6d 75 6c 74 69 |arder th|an multi|
|00005340| 70 6c 79 69 6e 67 2e 0a | 58 20 2a 2f 0a 58 76 6f |plying..|X */.Xvo|
|00005350| 69 64 0a 58 7a 72 65 64 | 63 64 65 63 6f 64 65 28 |id.Xzred|cdecode(|
|00005360| 72 70 2c 20 7a 31 2c 20 | 72 65 73 29 0a 58 09 52 |rp, z1, |res).X.R|
|00005370| 45 44 43 20 2a 72 70 3b | 09 09 2f 2a 20 52 45 44 |EDC *rp;|../* RED|
|00005380| 43 20 69 6e 66 6f 72 6d | 61 74 69 6f 6e 20 2a 2f |C inform|ation */|
|00005390| 0a 58 09 5a 56 41 4c 55 | 45 20 7a 31 3b 09 09 2f |.X.ZVALU|E z1;../|
|000053a0| 2a 20 6e 75 6d 62 65 72 | 20 74 6f 20 62 65 20 74 |* number| to be t|
|000053b0| 72 61 6e 73 66 6f 72 6d | 65 64 20 2a 2f 0a 58 09 |ransform|ed */.X.|
|000053c0| 5a 56 41 4c 55 45 20 2a | 72 65 73 3b 09 09 2f 2a |ZVALUE *|res;../*|
|000053d0| 20 72 65 74 75 72 6e 65 | 64 20 74 72 61 6e 73 66 | returne|d transf|
|000053e0| 6f 72 6d 65 64 20 6e 75 | 6d 62 65 72 20 2a 2f 0a |ormed nu|mber */.|
|000053f0| 58 7b 0a 58 09 5a 56 41 | 4c 55 45 20 74 6d 70 31 |X{.X.ZVA|LUE tmp1|
|00005400| 2c 20 74 6d 70 32 3b 09 | 2f 2a 20 74 65 6d 70 6f |, tmp2;.|/* tempo|
|00005410| 72 61 72 69 65 73 20 2a | 2f 0a 58 09 48 41 4c 46 |raries *|/.X.HALF|
|00005420| 20 2a 68 70 3b 09 09 2f | 2a 20 73 61 76 65 64 20 | *hp;../|* saved |
|00005430| 70 6f 69 6e 74 65 72 20 | 74 6f 20 74 6d 70 32 20 |pointer |to tmp2 |
|00005440| 76 61 6c 75 65 20 2a 2f | 0a 58 0a 58 09 69 66 20 |value */|.X.X.if |
|00005450| 28 69 73 6e 65 67 28 7a | 31 29 29 0a 58 09 09 65 |(isneg(z|1)).X..e|
|00005460| 72 72 6f 72 28 22 4e 65 | 67 61 74 69 76 65 20 6e |rror("Ne|gative n|
|00005470| 75 6d 62 65 72 20 66 6f | 72 20 7a 72 65 64 63 22 |umber fo|r zredc"|
|00005480| 29 3b 0a 58 0a 58 09 2f | 2a 0a 58 09 20 2a 20 43 |);.X.X./|*.X. * C|
|00005490| 68 65 63 6b 20 66 69 72 | 73 74 20 66 6f 72 20 74 |heck fir|st for t|
|000054a0| 68 65 20 73 70 65 63 69 | 61 6c 20 76 61 6c 75 65 |he speci|al value|
|000054b0| 73 20 66 6f 72 20 30 20 | 61 6e 64 20 31 20 74 68 |s for 0 |and 1 th|
|000054c0| 61 74 20 61 72 65 20 65 | 61 73 79 2e 0a 58 09 20 |at are e|asy..X. |
|000054d0| 2a 2f 0a 58 09 69 66 20 | 28 69 73 7a 65 72 6f 28 |*/.X.if |(iszero(|
|000054e0| 7a 31 29 29 20 7b 0a 58 | 09 09 2a 72 65 73 20 3d |z1)) {.X|..*res =|
|000054f0| 20 5f 7a 65 72 6f 5f 3b | 0a 58 09 09 72 65 74 75 | _zero_;|.X..retu|
|00005500| 72 6e 3b 0a 58 09 7d 0a | 58 09 69 66 20 28 28 7a |rn;.X.}.|X.if ((z|
|00005510| 31 2e 6c 65 6e 20 3d 3d | 20 72 70 2d 3e 6f 6e 65 |1.len ==| rp->one|
|00005520| 2e 6c 65 6e 29 20 26 26 | 20 28 7a 31 2e 76 5b 30 |.len) &&| (z1.v[0|
|00005530| 5d 20 3d 3d 20 72 70 2d | 3e 6f 6e 65 2e 76 5b 30 |] == rp-|>one.v[0|
|00005540| 5d 29 20 26 26 0a 58 09 | 09 28 7a 63 6d 70 28 7a |]) &&.X.|.(zcmp(z|
|00005550| 31 2c 20 72 70 2d 3e 6f | 6e 65 29 20 3d 3d 20 30 |1, rp->o|ne) == 0|
|00005560| 29 29 20 7b 0a 58 09 09 | 09 2a 72 65 73 20 3d 20 |)) {.X..|.*res = |
|00005570| 5f 6f 6e 65 5f 3b 0a 58 | 09 09 09 72 65 74 75 72 |_one_;.X|...retur|
|00005580| 6e 3b 0a 58 09 7d 0a 58 | 0a 58 09 2f 2a 0a 58 09 |n;.X.}.X|.X./*.X.|
|00005590| 20 2a 20 46 69 72 73 74 | 20 63 61 6c 63 75 6c 61 | * First| calcula|
|000055a0| 74 65 20 74 68 65 20 66 | 6f 6c 6c 6f 77 69 6e 67 |te the f|ollowing|
|000055b0| 3a 0a 58 09 20 2a 20 09 | 74 6d 70 32 20 3d 20 28 |:.X. * .|tmp2 = (|
|000055c0| 28 7a 31 20 25 20 32 5e | 62 69 74 6e 75 6d 29 20 |(z1 % 2^|bitnum) |
|000055d0| 2a 20 69 6e 76 29 20 25 | 20 32 5e 62 69 74 6e 75 |* inv) %| 2^bitnu|
|000055e0| 6d 2e 0a 58 09 20 2a 20 | 54 68 65 20 6d 6f 64 20 |m..X. * |The mod |
|000055f0| 6f 70 65 72 61 74 69 6f | 6e 73 20 63 61 6e 20 62 |operatio|ns can b|
|00005600| 65 20 64 6f 6e 65 20 77 | 69 74 68 20 6e 6f 20 77 |e done w|ith no w|
|00005610| 6f 72 6b 20 73 69 6e 63 | 65 20 74 68 65 20 62 69 |ork sinc|e the bi|
|00005620| 74 0a 58 09 20 2a 20 6e | 75 6d 62 65 72 20 77 61 |t.X. * n|umber wa|
|00005630| 73 20 73 65 6c 65 63 74 | 65 64 20 61 73 20 61 20 |s select|ed as a |
|00005640| 6d 75 6c 74 69 70 6c 65 | 20 6f 66 20 74 68 65 20 |multiple| of the |
|00005650| 77 6f 72 64 20 73 69 7a | 65 2e 20 20 4a 75 73 74 |word siz|e.  Just|
|00005660| 0a 58 09 20 2a 20 72 65 | 64 75 63 65 20 74 68 65 |.X. * re|duce the|
|00005670| 20 73 69 7a 65 73 20 6f | 66 20 74 68 65 20 6e 75 | sizes o|f the nu|
|00005680| 6d 62 65 72 73 20 61 73 | 20 72 65 71 75 69 72 65 |mbers as| require|
|00005690| 64 2e 0a 58 09 20 2a 2f | 0a 58 09 74 6d 70 31 20 |d..X. */|.X.tmp1 |
|000056a0| 3d 20 7a 31 3b 0a 58 09 | 69 66 20 28 74 6d 70 31 |= z1;.X.|if (tmp1|
|000056b0| 2e 6c 65 6e 20 3e 20 72 | 70 2d 3e 6c 65 6e 29 0a |.len > r|p->len).|
|000056c0| 58 09 09 74 6d 70 31 2e | 6c 65 6e 20 3d 20 72 70 |X..tmp1.|len = rp|
|000056d0| 2d 3e 6c 65 6e 3b 0a 58 | 09 7a 6d 75 6c 28 74 6d |->len;.X|.zmul(tm|
|000056e0| 70 31 2c 20 72 70 2d 3e | 69 6e 76 2c 20 26 74 6d |p1, rp->|inv, &tm|
|000056f0| 70 32 29 3b 0a 58 09 69 | 66 20 28 74 6d 70 32 2e |p2);.X.i|f (tmp2.|
|00005700| 6c 65 6e 20 3e 20 72 70 | 2d 3e 6c 65 6e 29 0a 58 |len > rp|->len).X|
|00005710| 09 09 74 6d 70 32 2e 6c | 65 6e 20 3d 20 72 70 2d |..tmp2.l|en = rp-|
|00005720| 3e 6c 65 6e 3b 0a 58 0a | 58 09 2f 2a 0a 58 09 20 |>len;.X.|X./*.X. |
|00005730| 2a 20 4e 65 78 74 20 63 | 61 6c 63 75 6c 61 74 65 |* Next c|alculate|
|00005740| 20 74 68 65 20 66 6f 6c | 6c 6f 77 69 6e 67 3a 0a | the fol|lowing:.|
|00005750| 58 09 20 2a 09 72 65 73 | 20 3d 20 28 7a 31 20 2b |X. *.res| = (z1 +|
|00005760| 20 74 6d 70 32 20 2a 20 | 6d 6f 64 75 6c 75 73 29 | tmp2 * |modulus)|
|00005770| 20 2f 20 32 5e 62 69 74 | 6e 75 6d 0a 58 09 20 2a | / 2^bit|num.X. *|
|00005780| 20 54 68 65 20 64 69 76 | 69 73 69 6f 6e 20 62 79 | The div|ision by|
|00005790| 20 61 20 70 6f 77 65 72 | 20 6f 66 20 32 20 69 73 | a power| of 2 is|
|000057a0| 20 61 6c 77 61 79 73 20 | 65 78 61 63 74 2c 20 61 | always |exact, a|
|000057b0| 6e 64 20 72 65 71 75 69 | 72 65 73 20 6e 6f 0a 58 |nd requi|res no.X|
|000057c0| 09 20 2a 20 77 6f 72 6b | 2e 20 20 4a 75 73 74 20 |. * work|.  Just |
|000057d0| 61 64 6a 75 73 74 20 74 | 68 65 20 61 64 64 72 65 |adjust t|he addre|
|000057e0| 73 73 20 61 6e 64 20 6c | 65 6e 67 74 68 20 6f 66 |ss and l|ength of|
|000057f0| 20 74 68 65 20 6e 75 6d | 62 65 72 20 74 6f 20 64 | the num|ber to d|
|00005800| 6f 0a 58 09 20 2a 20 74 | 68 65 20 64 69 76 69 64 |o.X. * t|he divid|
|00005810| 65 2c 20 62 75 74 20 73 | 61 76 65 20 74 68 65 20 |e, but s|ave the |
|00005820| 6f 72 69 67 69 6e 61 6c | 20 70 6f 69 6e 74 65 72 |original| pointer|
|00005830| 20 66 6f 72 20 66 72 65 | 65 69 6e 67 20 6c 61 74 | for fre|eing lat|
|00005840| 65 72 2e 0a 58 09 20 2a | 2f 0a 58 09 7a 6d 75 6c |er..X. *|/.X.zmul|
|00005850| 28 74 6d 70 32 2c 20 72 | 70 2d 3e 6d 6f 64 2c 20 |(tmp2, r|p->mod, |
|00005860| 26 74 6d 70 31 29 3b 0a | 58 09 66 72 65 65 68 28 |&tmp1);.|X.freeh(|
|00005870| 74 6d 70 32 2e 76 29 3b | 0a 58 09 7a 61 64 64 28 |tmp2.v);|.X.zadd(|
|00005880| 7a 31 2c 20 74 6d 70 31 | 2c 20 26 74 6d 70 32 29 |z1, tmp1|, &tmp2)|
|00005890| 3b 0a 58 09 66 72 65 65 | 68 28 74 6d 70 31 2e 76 |;.X.free|h(tmp1.v|
|000058a0| 29 3b 0a 58 09 68 70 20 | 3d 20 74 6d 70 32 2e 76 |);.X.hp |= tmp2.v|
|000058b0| 3b 0a 58 09 69 66 20 28 | 74 6d 70 32 2e 6c 65 6e |;.X.if (|tmp2.len|
|000058c0| 20 3c 3d 20 72 70 2d 3e | 6c 65 6e 29 20 7b 0a 58 | <= rp->|len) {.X|
|000058d0| 09 09 66 72 65 65 68 28 | 68 70 29 3b 0a 58 09 09 |..freeh(|hp);.X..|
|000058e0| 2a 72 65 73 20 3d 20 5f | 7a 65 72 6f 5f 3b 0a 58 |*res = _|zero_;.X|
|000058f0| 09 09 72 65 74 75 72 6e | 3b 0a 58 09 7d 0a 58 09 |..return|;.X.}.X.|
|00005900| 74 6d 70 32 2e 76 20 2b | 3d 20 72 70 2d 3e 6c 65 |tmp2.v +|= rp->le|
|00005910| 6e 3b 0a 58 09 74 6d 70 | 32 2e 6c 65 6e 20 2d 3d |n;.X.tmp|2.len -=|
|00005920| 20 72 70 2d 3e 6c 65 6e | 3b 0a 58 0a 58 09 2f 2a | rp->len|;.X.X./*|
|00005930| 0a 58 09 20 2a 20 46 69 | 6e 61 6c 6c 79 20 64 6f |.X. * Fi|nally do|
|00005940| 20 61 20 66 69 6e 61 6c | 20 6d 6f 64 75 6c 6f 20 | a final| modulo |
|00005950| 62 79 20 61 20 73 69 6d | 70 6c 65 20 73 75 62 74 |by a sim|ple subt|
|00005960| 72 61 63 74 69 6f 6e 20 | 69 66 20 6e 65 63 65 73 |raction |if neces|
|00005970| 73 61 72 79 2e 0a 58 09 | 20 2a 20 54 68 69 73 20 |sary..X.| * This |
|00005980| 69 73 20 61 6c 6c 20 74 | 68 61 74 20 69 73 20 6e |is all t|hat is n|
|00005990| 65 65 64 65 64 20 62 65 | 63 61 75 73 65 20 74 68 |eeded be|cause th|
|000059a0| 65 20 70 72 65 76 69 6f | 75 73 20 63 61 6c 63 75 |e previo|us calcu|
|000059b0| 6c 61 74 69 6f 6e 20 69 | 73 0a 58 09 20 2a 20 67 |lation i|s.X. * g|
|000059c0| 75 61 72 61 6e 74 65 65 | 64 20 74 6f 20 61 6c 77 |uarantee|d to alw|
|000059d0| 61 79 73 20 62 65 20 6c | 65 73 73 20 74 68 61 6e |ays be l|ess than|
|000059e0| 20 74 77 69 63 65 20 74 | 68 65 20 6d 6f 64 75 6c | twice t|he modul|
|000059f0| 75 73 2e 0a 58 09 20 2a | 2f 0a 58 09 69 66 20 28 |us..X. *|/.X.if (|
|00005a00| 7a 72 65 6c 28 74 6d 70 | 32 2c 20 72 70 2d 3e 6d |zrel(tmp|2, rp->m|
|00005a10| 6f 64 29 20 3c 20 30 29 | 0a 58 09 09 7a 63 6f 70 |od) < 0)|.X..zcop|
|00005a20| 79 28 74 6d 70 32 2c 20 | 72 65 73 29 3b 0a 58 09 |y(tmp2, |res);.X.|
|00005a30| 65 6c 73 65 0a 58 09 09 | 7a 73 75 62 28 74 6d 70 |else.X..|zsub(tmp|
|00005a40| 32 2c 20 72 70 2d 3e 6d | 6f 64 2c 20 72 65 73 29 |2, rp->m|od, res)|
|00005a50| 3b 0a 58 09 66 72 65 65 | 68 28 68 70 29 3b 0a 58 |;.X.free|h(hp);.X|
|00005a60| 7d 0a 58 0a 58 0a 58 2f | 2a 0a 58 20 2a 20 4d 75 |}.X.X.X/|*.X * Mu|
|00005a70| 6c 74 69 70 6c 79 20 74 | 77 6f 20 6e 75 6d 62 65 |ltiply t|wo numbe|
|00005a80| 72 73 20 69 6e 20 52 45 | 44 43 20 66 6f 72 6d 61 |rs in RE|DC forma|
|00005a90| 74 20 74 6f 67 65 74 68 | 65 72 20 70 72 6f 64 75 |t togeth|er produ|
|00005aa0| 63 69 6e 67 20 61 20 72 | 65 73 75 6c 74 20 61 6c |cing a r|esult al|
|00005ab0| 73 6f 0a 58 20 2a 20 69 | 6e 20 52 45 44 43 20 66 |so.X * i|n REDC f|
|00005ac0| 6f 72 6d 61 74 2e 20 20 | 49 66 20 74 68 65 20 72 |ormat.  |If the r|
|00005ad0| 65 73 75 6c 74 20 69 73 | 20 63 6f 6e 76 65 72 74 |esult is| convert|
|00005ae0| 65 64 20 62 61 63 6b 20 | 74 6f 20 61 20 6e 6f 72 |ed back |to a nor|
|00005af0| 6d 61 6c 20 6e 75 6d 62 | 65 72 2c 0a 58 20 2a 20 |mal numb|er,.X * |
|00005b00| 74 68 65 6e 20 74 68 65 | 20 72 65 73 75 6c 74 20 |then the| result |
|00005b10| 69 73 20 74 68 65 20 73 | 61 6d 65 20 61 73 20 74 |is the s|ame as t|
|00005b20| 68 65 20 6d 6f 64 75 6c | 6f 27 64 20 6d 75 6c 74 |he modul|o'd mult|
|00005b30| 69 70 6c 69 63 61 74 69 | 6f 6e 20 6f 66 20 74 68 |iplicati|on of th|
|00005b40| 65 0a 58 20 2a 20 6f 72 | 69 67 69 6e 61 6c 20 6e |e.X * or|iginal n|
|00005b50| 75 6d 62 65 72 73 20 62 | 65 66 6f 72 65 20 74 68 |umbers b|efore th|
|00005b60| 65 79 20 77 65 72 65 20 | 63 6f 6e 76 65 72 74 65 |ey were |converte|
|00005b70| 64 20 74 6f 20 52 45 44 | 43 20 66 6f 72 6d 61 74 |d to RED|C format|
|00005b80| 2e 20 20 54 68 69 73 0a | 58 20 2a 20 63 61 6c 63 |.  This.|X * calc|
|00005b90| 75 6c 61 74 69 6f 6e 20 | 69 73 20 64 6f 6e 65 20 |ulation |is done |
|00005ba0| 69 6e 20 6f 6e 65 20 6f | 66 20 74 77 6f 20 77 61 |in one o|f two wa|
|00005bb0| 79 73 2c 20 64 65 70 65 | 6e 64 69 6e 67 20 6f 6e |ys, depe|nding on|
|00005bc0| 20 74 68 65 20 73 69 7a | 65 20 6f 66 20 74 68 65 | the siz|e of the|
|00005bd0| 0a 58 20 2a 20 6d 6f 64 | 75 6c 75 73 2e 20 20 46 |.X * mod|ulus.  F|
|00005be0| 6f 72 20 6c 61 72 67 65 | 20 6e 75 6d 62 65 72 73 |or large| numbers|
|00005bf0| 2c 20 74 68 65 20 52 45 | 44 43 20 64 65 66 69 6e |, the RE|DC defin|
|00005c00| 69 74 69 6f 6e 20 69 73 | 20 75 73 65 64 20 64 69 |ition is| used di|
|00005c10| 72 65 63 74 6c 79 0a 58 | 20 2a 20 77 68 69 63 68 |rectly.X| * which|
|00005c20| 20 69 6e 76 6f 6c 76 65 | 73 20 74 68 72 65 65 20 | involve|s three |
|00005c30| 6d 75 6c 74 69 70 6c 69 | 65 73 20 6f 76 65 72 61 |multipli|es overa|
|00005c40| 6c 6c 2e 20 20 46 6f 72 | 20 73 6d 61 6c 6c 20 6e |ll.  For| small n|
|00005c50| 75 6d 62 65 72 73 2c 20 | 61 0a 58 20 2a 20 63 6f |umbers, |a.X * co|
|00005c60| 6d 70 6c 69 63 61 74 65 | 64 20 72 6f 75 74 69 6e |mplicate|d routin|
|00005c70| 65 20 69 73 20 75 73 65 | 64 20 77 68 69 63 68 20 |e is use|d which |
|00005c80| 64 6f 65 73 20 74 68 65 | 20 69 6e 64 69 63 61 74 |does the| indicat|
|00005c90| 65 64 20 6d 75 6c 74 69 | 70 6c 69 63 61 74 69 6f |ed multi|plicatio|
|00005ca0| 6e 0a 58 20 2a 20 61 6e | 64 20 74 68 65 20 52 45 |n.X * an|d the RE|
|00005cb0| 44 43 20 61 6c 67 6f 72 | 69 74 68 6d 20 61 74 20 |DC algor|ithm at |
|00005cc0| 74 68 65 20 73 61 6d 65 | 20 74 69 6d 65 20 74 6f |the same| time to|
|00005cd0| 20 70 72 6f 64 75 63 65 | 20 74 68 65 20 72 65 73 | produce| the res|
|00005ce0| 75 6c 74 2e 0a 58 20 2a | 2f 0a 58 76 6f 69 64 0a |ult..X *|/.Xvoid.|
|00005cf0| 58 7a 72 65 64 63 6d 75 | 6c 28 72 70 2c 20 7a 31 |Xzredcmu|l(rp, z1|
|00005d00| 2c 20 7a 32 2c 20 72 65 | 73 29 0a 58 09 52 45 44 |, z2, re|s).X.RED|
|00005d10| 43 20 2a 72 70 3b 09 09 | 2f 2a 20 52 45 44 43 20 |C *rp;..|/* REDC |
|00005d20| 69 6e 66 6f 72 6d 61 74 | 69 6f 6e 20 2a 2f 0a 58 |informat|ion */.X|
|00005d30| 09 5a 56 41 4c 55 45 20 | 7a 31 3b 09 09 2f 2a 20 |.ZVALUE |z1;../* |
|00005d40| 66 69 72 73 74 20 52 45 | 44 43 20 6e 75 6d 62 65 |first RE|DC numbe|
|00005d50| 72 20 74 6f 20 62 65 20 | 6d 75 6c 74 69 70 6c 69 |r to be |multipli|
|00005d60| 65 64 20 2a 2f 0a 58 09 | 5a 56 41 4c 55 45 20 7a |ed */.X.|ZVALUE z|
|00005d70| 32 3b 09 09 2f 2a 20 73 | 65 63 6f 6e 64 20 52 45 |2;../* s|econd RE|
|00005d80| 44 43 20 6e 75 6d 62 65 | 72 20 74 6f 20 62 65 20 |DC numbe|r to be |
|00005d90| 6d 75 6c 74 69 70 6c 69 | 65 64 20 2a 2f 0a 58 09 |multipli|ed */.X.|
|00005da0| 5a 56 41 4c 55 45 20 2a | 72 65 73 3b 09 09 2f 2a |ZVALUE *|res;../*|
|00005db0| 20 72 65 73 75 6c 74 69 | 6e 67 20 52 45 44 43 20 | resulti|ng REDC |
|00005dc0| 6e 75 6d 62 65 72 20 2a | 2f 0a 58 7b 0a 58 09 46 |number *|/.X{.X.F|
|00005dd0| 55 4c 4c 20 6d 75 6c 62 | 3b 0a 58 09 46 55 4c 4c |ULL mulb|;.X.FULL|
|00005de0| 20 6d 75 6c 6e 3b 0a 58 | 09 48 41 4c 46 20 2a 68 | muln;.X|.HALF *h|
|00005df0| 31 3b 0a 58 09 48 41 4c | 46 20 2a 68 32 3b 0a 58 |1;.X.HAL|F *h2;.X|
|00005e00| 09 48 41 4c 46 20 2a 68 | 33 3b 0a 58 09 48 41 4c |.HALF *h|3;.X.HAL|
|00005e10| 46 20 2a 68 64 3b 0a 58 | 09 48 41 4c 46 20 4e 69 |F *hd;.X|.HALF Ni|
|00005e20| 6e 76 3b 0a 58 09 48 41 | 4c 46 20 74 6f 70 64 69 |nv;.X.HA|LF topdi|
|00005e30| 67 69 74 3b 0a 58 09 4c | 45 4e 20 6d 6f 64 6c 65 |git;.X.L|EN modle|
|00005e40| 6e 3b 0a 58 09 4c 45 4e | 20 6c 65 6e 3b 0a 58 09 |n;.X.LEN| len;.X.|
|00005e50| 4c 45 4e 20 6c 65 6e 32 | 3b 0a 58 09 53 49 55 4e |LEN len2|;.X.SIUN|
|00005e60| 49 4f 4e 20 73 69 76 61 | 6c 31 3b 0a 58 09 53 49 |ION siva|l1;.X.SI|
|00005e70| 55 4e 49 4f 4e 20 73 69 | 76 61 6c 32 3b 0a 58 09 |UNION si|val2;.X.|
|00005e80| 53 49 55 4e 49 4f 4e 20 | 73 69 76 61 6c 33 3b 0a |SIUNION |sival3;.|
|00005e90| 58 09 53 49 55 4e 49 4f | 4e 20 63 61 72 72 79 3b |X.SIUNIO|N carry;|
|00005ea0| 0a 58 09 5a 56 41 4c 55 | 45 20 74 6d 70 3b 0a 58 |.X.ZVALU|E tmp;.X|
|00005eb0| 0a 58 09 69 66 20 28 69 | 73 6e 65 67 28 7a 31 29 |.X.if (i|sneg(z1)|
|00005ec0| 20 7c 7c 20 28 7a 31 2e | 6c 65 6e 20 3e 20 72 70 | || (z1.|len > rp|
|00005ed0| 2d 3e 6d 6f 64 2e 6c 65 | 6e 29 20 7c 7c 0a 58 09 |->mod.le|n) ||.X.|
|00005ee0| 09 69 73 6e 65 67 28 7a | 32 29 20 7c 7c 20 28 7a |.isneg(z|2) || (z|
|00005ef0| 32 2e 6c 65 6e 20 3e 20 | 72 70 2d 3e 6d 6f 64 2e |2.len > |rp->mod.|
|00005f00| 6c 65 6e 29 29 0a 58 09 | 09 09 65 72 72 6f 72 28 |len)).X.|..error(|
|00005f10| 22 4e 65 67 61 74 69 76 | 65 20 6f 72 20 74 6f 6f |"Negativ|e or too|
|00005f20| 20 6c 61 72 67 65 20 6e | 75 6d 62 65 72 20 69 6e | large n|umber in|
|00005f30| 20 7a 72 65 64 63 6d 75 | 6c 22 29 3b 0a 58 0a 58 | zredcmu|l");.X.X|
|00005f40| 09 2f 2a 0a 58 09 20 2a | 20 43 68 65 63 6b 20 66 |./*.X. *| Check f|
|00005f50| 6f 72 20 73 70 65 63 69 | 61 6c 20 76 61 6c 75 65 |or speci|al value|
|00005f60| 73 20 77 68 69 63 68 20 | 77 65 20 65 61 73 69 6c |s which |we easil|
|00005f70| 79 20 6b 6e 6f 77 20 74 | 68 65 20 61 6e 73 77 65 |y know t|he answe|
|00005f80| 72 2e 0a 58 09 20 2a 2f | 0a 58 09 69 66 20 28 69 |r..X. */|.X.if (i|
|00005f90| 73 7a 65 72 6f 28 7a 31 | 29 20 7c 7c 20 69 73 7a |szero(z1|) || isz|
|00005fa0| 65 72 6f 28 7a 32 29 29 | 20 7b 0a 58 09 09 2a 72 |ero(z2))| {.X..*r|
|00005fb0| 65 73 20 3d 20 5f 7a 65 | 72 6f 5f 3b 0a 58 09 09 |es = _ze|ro_;.X..|
|00005fc0| 72 65 74 75 72 6e 3b 0a | 58 09 7d 0a 58 0a 58 09 |return;.|X.}.X.X.|
|00005fd0| 69 66 20 28 28 7a 31 2e | 6c 65 6e 20 3d 3d 20 72 |if ((z1.|len == r|
|00005fe0| 70 2d 3e 6f 6e 65 2e 6c | 65 6e 29 20 26 26 20 28 |p->one.l|en) && (|
|00005ff0| 7a 31 2e 76 5b 30 5d 20 | 3d 3d 20 72 70 2d 3e 6f |z1.v[0] |== rp->o|
|00006000| 6e 65 2e 76 5b 30 5d 29 | 20 26 26 0a 58 09 09 28 |ne.v[0])| &&.X..(|
|00006010| 7a 63 6d 70 28 7a 31 2c | 20 72 70 2d 3e 6f 6e 65 |zcmp(z1,| rp->one|
|00006020| 29 20 3d 3d 20 30 29 29 | 20 7b 0a 58 09 09 09 7a |) == 0))| {.X...z|
|00006030| 63 6f 70 79 28 7a 32 2c | 20 72 65 73 29 3b 0a 58 |copy(z2,| res);.X|
|00006040| 09 09 09 72 65 74 75 72 | 6e 3b 0a 58 09 7d 0a 58 |...retur|n;.X.}.X|
|00006050| 0a 58 09 69 66 20 28 28 | 7a 32 2e 6c 65 6e 20 3d |.X.if ((|z2.len =|
|00006060| 3d 20 72 70 2d 3e 6f 6e | 65 2e 6c 65 6e 29 20 26 |= rp->on|e.len) &|
|00006070| 26 20 28 7a 32 2e 76 5b | 30 5d 20 3d 3d 20 72 70 |& (z2.v[|0] == rp|
|00006080| 2d 3e 6f 6e 65 2e 76 5b | 30 5d 29 20 26 26 0a 58 |->one.v[|0]) &&.X|
|00006090| 09 09 28 7a 63 6d 70 28 | 7a 32 2c 20 72 70 2d 3e |..(zcmp(|z2, rp->|
|000060a0| 6f 6e 65 29 20 3d 3d 20 | 30 29 29 20 7b 0a 58 09 |one) == |0)) {.X.|
|000060b0| 09 09 7a 63 6f 70 79 28 | 7a 31 2c 20 72 65 73 29 |..zcopy(|z1, res)|
|000060c0| 3b 0a 58 09 09 09 72 65 | 74 75 72 6e 3b 0a 58 09 |;.X...re|turn;.X.|
|000060d0| 7d 0a 58 0a 58 09 2f 2a | 0a 58 09 20 2a 20 49 66 |}.X.X./*|.X. * If|
|000060e0| 20 74 68 65 20 73 69 7a | 65 20 6f 66 20 74 68 65 | the siz|e of the|
|000060f0| 20 6d 6f 64 75 6c 75 73 | 20 69 73 20 6c 61 72 67 | modulus| is larg|
|00006100| 65 2c 20 74 68 65 6e 20 | 6a 75 73 74 20 64 6f 20 |e, then |just do |
|00006110| 74 68 65 20 6d 75 6c 74 | 69 70 6c 79 2c 0a 58 09 |the mult|iply,.X.|
|00006120| 20 2a 20 66 6f 6c 6c 6f | 77 65 64 20 62 79 20 74 | * follo|wed by t|
|00006130| 68 65 20 74 77 6f 20 6d | 75 6c 74 69 70 6c 69 65 |he two m|ultiplie|
|00006140| 73 20 63 6f 6e 74 61 69 | 6e 65 64 20 69 6e 20 74 |s contai|ned in t|
|00006150| 68 65 20 52 45 44 43 20 | 72 6f 75 74 69 6e 65 2e |he REDC |routine.|
|00006160| 0a 58 09 20 2a 20 54 68 | 69 73 20 77 69 6c 6c 20 |.X. * Th|is will |
|00006170| 62 65 20 71 75 69 63 6b | 65 72 20 74 68 61 6e 20 |be quick|er than |
|00006180| 64 69 72 65 63 74 6c 79 | 20 64 6f 69 6e 67 20 74 |directly| doing t|
|00006190| 68 65 20 52 45 44 43 20 | 63 61 6c 63 75 6c 61 74 |he REDC |calculat|
|000061a0| 69 6f 6e 0a 58 09 20 2a | 20 62 65 63 61 75 73 65 |ion.X. *| because|
|000061b0| 20 6f 66 20 74 68 65 20 | 4f 28 4e 5e 31 2e 35 38 | of the |O(N^1.58|
|000061c0| 35 29 20 73 70 65 65 64 | 20 6f 66 20 74 68 65 20 |5) speed| of the |
|000061d0| 6d 75 6c 74 69 70 6c 69 | 65 73 2e 20 20 54 68 65 |multipli|es.  The|
|000061e0| 20 73 69 7a 65 0a 58 09 | 20 2a 20 6f 66 20 74 68 | size.X.| * of th|
|000061f0| 65 20 6e 75 6d 62 65 72 | 20 77 68 69 63 68 20 74 |e number| which t|
|00006200| 68 69 73 20 69 73 20 64 | 6f 6e 65 20 69 73 20 63 |his is d|one is c|
|00006210| 6f 6e 66 69 67 75 72 61 | 62 6c 65 2e 0a 58 09 20 |onfigura|ble..X. |
|00006220| 2a 2f 0a 58 09 69 66 20 | 28 72 70 2d 3e 6d 6f 64 |*/.X.if |(rp->mod|
|00006230| 2e 6c 65 6e 20 3e 3d 20 | 5f 72 65 64 63 32 5f 29 |.len >= |_redc2_)|
|00006240| 20 7b 0a 58 09 09 7a 6d | 75 6c 28 7a 31 2c 20 7a | {.X..zm|ul(z1, z|
|00006250| 32 2c 20 26 74 6d 70 29 | 3b 0a 58 09 09 7a 72 65 |2, &tmp)|;.X..zre|
|00006260| 64 63 64 65 63 6f 64 65 | 28 72 70 2c 20 74 6d 70 |dcdecode|(rp, tmp|
|00006270| 2c 20 72 65 73 29 3b 0a | 58 09 09 66 72 65 65 68 |, res);.|X..freeh|
|00006280| 28 74 6d 70 2e 76 29 3b | 0a 58 09 09 72 65 74 75 |(tmp.v);|.X..retu|
|00006290| 72 6e 3b 0a 58 09 7d 0a | 58 0a 58 09 2f 2a 0a 58 |rn;.X.}.|X.X./*.X|
|000062a0| 09 20 2a 20 54 68 65 20 | 6e 75 6d 62 65 72 20 69 |. * The |number i|
|000062b0| 73 20 73 6d 61 6c 6c 20 | 65 6e 6f 75 67 68 20 74 |s small |enough t|
|000062c0| 6f 20 63 61 6c 63 75 6c | 61 74 65 20 62 79 20 64 |o calcul|ate by d|
|000062d0| 6f 69 6e 67 20 74 68 65 | 20 4f 28 4e 5e 32 29 20 |oing the| O(N^2) |
|000062e0| 52 45 44 43 0a 58 09 20 | 2a 20 61 6c 67 6f 72 69 |REDC.X. |* algori|
|000062f0| 74 68 6d 20 64 69 72 65 | 63 74 6c 79 2e 20 20 54 |thm dire|ctly.  T|
|00006300| 68 69 73 20 61 6c 67 6f | 72 69 74 68 6d 20 70 65 |his algo|rithm pe|
|00006310| 72 66 6f 72 6d 73 20 74 | 68 65 20 6d 75 6c 74 69 |rforms t|he multi|
|00006320| 70 6c 69 63 61 74 69 6f | 6e 20 61 6e 64 0a 58 09 |plicatio|n and.X.|
|00006330| 20 2a 20 74 68 65 20 72 | 65 64 75 63 74 69 6f 6e | * the r|eduction|
|00006340| 20 61 74 20 74 68 65 20 | 73 61 6d 65 20 74 69 6d | at the |same tim|
|00006350| 65 2e 20 20 4e 6f 74 69 | 63 65 20 74 68 65 20 6f |e.  Noti|ce the o|
|00006360| 62 73 63 75 72 65 20 66 | 61 63 74 73 20 74 68 61 |bscure f|acts tha|
|00006370| 74 0a 58 09 20 2a 20 6f | 6e 6c 79 20 74 68 65 20 |t.X. * o|nly the |
|00006380| 6c 6f 77 65 73 74 20 77 | 6f 72 64 20 6f 66 20 74 |lowest w|ord of t|
|00006390| 68 65 20 69 6e 76 65 72 | 73 65 20 76 61 6c 75 65 |he inver|se value|
|000063a0| 20 69 73 20 75 73 65 64 | 2c 20 61 6e 64 20 74 68 | is used|, and th|
|000063b0| 61 74 0a 58 09 20 2a 20 | 74 68 65 72 65 20 69 73 |at.X. * |there is|
|000063c0| 20 6e 6f 20 73 68 69 66 | 74 69 6e 67 20 6f 66 20 | no shif|ting of |
|000063d0| 74 68 65 20 70 61 72 74 | 69 61 6c 20 70 72 6f 64 |the part|ial prod|
|000063e0| 75 63 74 73 20 61 73 20 | 74 68 65 72 65 20 69 73 |ucts as |there is|
|000063f0| 20 69 6e 20 61 0a 58 09 | 20 2a 20 6e 6f 72 6d 61 | in a.X.| * norma|
+--------+-------------------------+-------------------------+--------+--------+