home *** CD-ROM | disk | FTP | other *** search
open in: 
MacOS 8.1
|
Win98
|
DOS
browse contents |
view JSON data
|
view as text
This file was processed as: SHell self-extracting ARchive
(archive/shar).
| Confidence | Program | Detection | Match Type | Support
|
|---|
| 100%
| dexvert
| SHell self-extracting ARchive (archive/shar)
| magic
| Supported |
| 1%
| dexvert
| Text File (text/txt)
| fallback
| Supported |
| 100%
| file
| ASCII text
| default
| |
| 100%
| checkBytes
| Printable ASCII
| default
| |
| 100%
| perlTextCheck
| Likely Text (Perl)
| default
| |
| 100%
| siegfried
| fmt/329 Shell Archive Format
| default
| |
| 100%
| detectItEasy
| Format: plain text[LF]
| default (weak)
|
|
hex view+--------+-------------------------+-------------------------+--------+--------+
|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 34 |Subject:| v26i044|
|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 | 38 2f 32 31 0a 53 65 6e |r, Part1|8/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 34 |ume 26, |Issue 44|
|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 38 0a 0a 23 21 20 2f |alc/part|18..#! /|
|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 38 20 28 6f 66 20 | 32 31 29 2e 22 0a 23 20 | 18 (of |21).".# |
|000002e0| 43 6f 6e 74 65 6e 74 73 | 3a 20 20 7a 66 75 6e 63 |Contents|: zfunc|
|000002f0| 2e 63 0a 23 20 57 72 61 | 70 70 65 64 20 62 79 20 |.c.# Wra|pped by |
|00000300| 64 62 65 6c 6c 40 65 6c | 6d 20 6f 6e 20 54 75 65 |dbell@el|m on Tue|
|00000310| 20 46 65 62 20 32 35 20 | 31 35 3a 32 31 3a 31 36 | Feb 25 |15:21:16|
|00000320| 20 31 39 39 32 0a 50 41 | 54 48 3d 2f 62 69 6e 3a | 1992.PA|TH=/bin:|
|00000330| 2f 75 73 72 2f 62 69 6e | 3a 2f 75 73 72 2f 75 63 |/usr/bin|:/usr/uc|
|00000340| 62 20 3b 20 65 78 70 6f | 72 74 20 50 41 54 48 0a |b ; expo|rt PATH.|
|00000350| 69 66 20 74 65 73 74 20 | 2d 66 20 27 7a 66 75 6e |if test |-f 'zfun|
|00000360| 63 2e 63 27 20 2d 61 20 | 22 24 7b 31 7d 22 20 21 |c.c' -a |"${1}" !|
|00000370| 3d 20 22 2d 63 22 20 3b | 20 74 68 65 6e 20 0a 20 |= "-c" ;| then . |
|00000380| 20 65 63 68 6f 20 73 68 | 61 72 3a 20 57 69 6c 6c | echo sh|ar: Will|
|00000390| 20 6e 6f 74 20 63 6c 6f | 62 62 65 72 20 65 78 69 | not clo|bber exi|
|000003a0| 73 74 69 6e 67 20 66 69 | 6c 65 20 5c 22 27 7a 66 |sting fi|le \"'zf|
|000003b0| 75 6e 63 2e 63 27 5c 22 | 0a 65 6c 73 65 0a 65 63 |unc.c'\"|.else.ec|
|000003c0| 68 6f 20 73 68 61 72 3a | 20 45 78 74 72 61 63 74 |ho shar:| Extract|
|000003d0| 69 6e 67 20 5c 22 27 7a | 66 75 6e 63 2e 63 27 5c |ing \"'z|func.c'\|
|000003e0| 22 20 5c 28 33 34 33 34 | 35 20 63 68 61 72 61 63 |" \(3434|5 charac|
|000003f0| 74 65 72 73 5c 29 0a 73 | 65 64 20 22 73 2f 5e 58 |ters\).s|ed "s/^X|
|00000400| 2f 2f 22 20 3e 27 7a 66 | 75 6e 63 2e 63 27 20 3c |//" >'zf|unc.c' <|
|00000410| 3c 27 45 4e 44 5f 4f 46 | 5f 46 49 4c 45 27 0a 58 |<'END_OF|_FILE'.X|
|00000420| 2f 2a 0a 58 20 2a 20 43 | 6f 70 79 72 69 67 68 74 |/*.X * C|opyright|
|00000430| 20 28 63 29 20 31 39 39 | 32 20 44 61 76 69 64 20 | (c) 199|2 David |
|00000440| 49 2e 20 42 65 6c 6c 0a | 58 20 2a 20 50 65 72 6d |I. Bell.|X * Perm|
|00000450| 69 73 73 69 6f 6e 20 69 | 73 20 67 72 61 6e 74 65 |ission i|s grante|
|00000460| 64 20 74 6f 20 75 73 65 | 2c 20 64 69 73 74 72 69 |d to use|, distri|
|00000470| 62 75 74 65 2c 20 6f 72 | 20 6d 6f 64 69 66 79 20 |bute, or| modify |
|00000480| 74 68 69 73 20 73 6f 75 | 72 63 65 2c 0a 58 20 2a |this sou|rce,.X *|
|00000490| 20 70 72 6f 76 69 64 65 | 64 20 74 68 61 74 20 74 | provide|d that t|
|000004a0| 68 69 73 20 63 6f 70 79 | 72 69 67 68 74 20 6e 6f |his copy|right no|
|000004b0| 74 69 63 65 20 72 65 6d | 61 69 6e 73 20 69 6e 74 |tice rem|ains int|
|000004c0| 61 63 74 2e 0a 58 20 2a | 0a 58 20 2a 20 45 78 74 |act..X *|.X * Ext|
|000004d0| 65 6e 64 65 64 20 70 72 | 65 63 69 73 69 6f 6e 20 |ended pr|ecision |
|000004e0| 69 6e 74 65 67 72 61 6c | 20 61 72 69 74 68 6d 65 |integral| arithme|
|000004f0| 74 69 63 20 6e 6f 6e 2d | 70 72 69 6d 69 74 69 76 |tic non-|primitiv|
|00000500| 65 20 72 6f 75 74 69 6e | 65 73 0a 58 20 2a 2f 0a |e routin|es.X */.|
|00000510| 58 0a 58 23 69 6e 63 6c | 75 64 65 20 22 6d 61 74 |X.X#incl|ude "mat|
|00000520| 68 2e 68 22 0a 58 0a 58 | 73 74 61 74 69 63 20 5a |h.h".X.X|static Z|
|00000530| 56 41 4c 55 45 20 70 72 | 69 6d 65 70 72 6f 64 3b |VALUE pr|imeprod;|
|00000540| 09 09 2f 2a 20 70 72 6f | 64 75 63 74 20 6f 66 20 |../* pro|duct of |
|00000550| 70 72 69 6d 65 73 20 75 | 6e 64 65 72 20 31 30 30 |primes u|nder 100|
|00000560| 20 2a 2f 0a 58 5a 56 41 | 4c 55 45 20 5f 74 65 6e | */.XZVA|LUE _ten|
|00000570| 70 6f 77 65 72 73 5f 5b | 33 32 5d 3b 09 09 09 2f |powers_[|32];.../|
|00000580| 2a 20 74 61 62 6c 65 20 | 6f 66 20 31 30 5e 32 5e |* table |of 10^2^|
|00000590| 6e 20 2a 2f 0a 58 0a 58 | 23 69 66 20 30 0a 58 73 |n */.X.X|#if 0.Xs|
|000005a0| 74 61 74 69 63 20 63 68 | 61 72 20 2a 61 62 6f 72 |tatic ch|ar *abor|
|000005b0| 74 6d 73 67 20 3d 20 22 | 43 61 6c 63 75 6c 61 74 |tmsg = "|Calculat|
|000005c0| 69 6f 6e 20 61 62 6f 72 | 74 65 64 22 3b 0a 58 73 |ion abor|ted";.Xs|
|000005d0| 74 61 74 69 63 20 63 68 | 61 72 20 2a 6d 65 6d 6d |tatic ch|ar *memm|
|000005e0| 73 67 20 3d 20 22 4e 6f | 74 20 65 6e 6f 75 67 68 |sg = "No|t enough|
|000005f0| 20 6d 65 6d 6f 72 79 22 | 3b 0a 58 23 65 6e 64 69 | memory"|;.X#endi|
|00000600| 66 0a 58 0a 58 0a 58 2f | 2a 0a 58 20 2a 20 43 6f |f.X.X.X/|*.X * Co|
|00000610| 6d 70 75 74 65 20 74 68 | 65 20 66 61 63 74 6f 72 |mpute th|e factor|
|00000620| 69 61 6c 20 6f 66 20 61 | 20 6e 75 6d 62 65 72 2e |ial of a| number.|
|00000630| 0a 58 20 2a 2f 0a 58 76 | 6f 69 64 0a 58 7a 66 61 |.X */.Xv|oid.Xzfa|
|00000640| 63 74 28 7a 2c 20 64 65 | 73 74 29 0a 58 09 5a 56 |ct(z, de|st).X.ZV|
|00000650| 41 4c 55 45 20 7a 2c 20 | 2a 64 65 73 74 3b 0a 58 |ALUE z, |*dest;.X|
|00000660| 7b 0a 58 09 6c 6f 6e 67 | 20 70 74 77 6f 3b 09 09 |{.X.long| ptwo;..|
|00000670| 2f 2a 20 63 6f 75 6e 74 | 20 6f 66 20 70 6f 77 65 |/* count| of powe|
|00000680| 72 73 20 6f 66 20 74 77 | 6f 20 2a 2f 0a 58 09 6c |rs of tw|o */.X.l|
|00000690| 6f 6e 67 20 6e 3b 09 09 | 09 2f 2a 20 63 75 72 72 |ong n;..|./* curr|
|000006a0| 65 6e 74 20 6d 75 6c 74 | 69 70 6c 69 63 61 74 69 |ent mult|iplicati|
|000006b0| 6f 6e 20 76 61 6c 75 65 | 20 2a 2f 0a 58 09 6c 6f |on value| */.X.lo|
|000006c0| 6e 67 20 6d 3b 09 09 09 | 2f 2a 20 72 65 64 75 63 |ng m;...|/* reduc|
|000006d0| 65 64 20 6d 75 6c 74 69 | 70 6c 69 63 61 74 69 6f |ed multi|plicatio|
|000006e0| 6e 20 76 61 6c 75 65 20 | 2a 2f 0a 58 09 6c 6f 6e |n value |*/.X.lon|
|000006f0| 67 20 6d 75 6c 3b 09 09 | 2f 2a 20 63 6f 6c 6c 65 |g mul;..|/* colle|
|00000700| 63 74 65 64 20 76 61 6c | 75 65 20 74 6f 20 6d 75 |cted val|ue to mu|
|00000710| 6c 74 69 70 6c 79 20 62 | 79 20 2a 2f 0a 58 09 5a |ltiply b|y */.X.Z|
|00000720| 56 41 4c 55 45 20 72 65 | 73 2c 20 74 65 6d 70 3b |VALUE re|s, temp;|
|00000730| 0a 58 0a 58 09 69 66 20 | 28 69 73 6e 65 67 28 7a |.X.X.if |(isneg(z|
|00000740| 29 29 0a 58 09 09 65 72 | 72 6f 72 28 22 4e 65 67 |)).X..er|ror("Neg|
|00000750| 61 74 69 76 65 20 61 72 | 67 75 6d 65 6e 74 20 66 |ative ar|gument f|
|00000760| 6f 72 20 66 61 63 74 6f | 72 69 61 6c 22 29 3b 0a |or facto|rial");.|
|00000770| 58 09 69 66 20 28 69 73 | 62 69 67 28 7a 29 29 0a |X.if (is|big(z)).|
|00000780| 58 09 09 65 72 72 6f 72 | 28 22 56 65 72 79 20 6c |X..error|("Very l|
|00000790| 61 72 67 65 20 66 61 63 | 74 6f 72 69 61 6c 22 29 |arge fac|torial")|
|000007a0| 3b 0a 58 09 6e 20 3d 20 | 28 69 73 74 69 6e 79 28 |;.X.n = |(istiny(|
|000007b0| 7a 29 20 3f 20 7a 31 74 | 6f 6c 28 7a 29 20 3a 20 |z) ? z1t|ol(z) : |
|000007c0| 7a 32 74 6f 6c 28 7a 29 | 29 3b 0a 58 09 70 74 77 |z2tol(z)|);.X.ptw|
|000007d0| 6f 20 3d 20 30 3b 0a 58 | 09 6d 75 6c 20 3d 20 31 |o = 0;.X|.mul = 1|
|000007e0| 3b 0a 58 09 72 65 73 20 | 3d 20 5f 6f 6e 65 5f 3b |;.X.res |= _one_;|
|000007f0| 0a 58 09 2f 2a 0a 58 09 | 20 2a 20 4d 75 6c 74 69 |.X./*.X.| * Multi|
|00000800| 70 6c 79 20 6e 75 6d 62 | 65 72 73 20 74 6f 67 65 |ply numb|ers toge|
|00000810| 74 68 65 72 2c 20 62 75 | 74 20 73 71 75 65 65 7a |ther, bu|t squeez|
|00000820| 65 20 6f 75 74 20 61 6c | 6c 20 70 6f 77 65 72 73 |e out al|l powers|
|00000830| 20 6f 66 20 74 77 6f 2e | 0a 58 09 20 2a 20 57 65 | of two.|.X. * We|
|00000840| 20 77 69 6c 6c 20 70 75 | 74 20 74 68 65 6d 20 62 | will pu|t them b|
|00000850| 61 63 6b 20 69 6e 20 61 | 74 20 74 68 65 20 65 6e |ack in a|t the en|
|00000860| 64 2e 20 20 41 6c 73 6f | 20 63 6f 6c 6c 65 63 74 |d. Also| collect|
|00000870| 20 6d 75 6c 74 69 70 6c | 65 0a 58 09 20 2a 20 6e | multipl|e.X. * n|
|00000880| 75 6d 62 65 72 73 20 74 | 6f 67 65 74 68 65 72 20 |umbers t|ogether |
|00000890| 75 6e 74 69 6c 20 74 68 | 65 72 65 20 69 73 20 61 |until th|ere is a|
|000008a0| 20 72 69 73 6b 20 6f 66 | 20 6f 76 65 72 66 6c 6f | risk of| overflo|
|000008b0| 77 2e 0a 58 09 20 2a 2f | 0a 58 09 66 6f 72 20 28 |w..X. */|.X.for (|
|000008c0| 3b 20 6e 20 3e 20 31 3b | 20 6e 2d 2d 29 20 7b 0a |; n > 1;| n--) {.|
|000008d0| 58 09 09 66 6f 72 20 28 | 6d 20 3d 20 6e 3b 20 28 |X..for (|m = n; (|
|000008e0| 28 6d 20 26 20 30 78 31 | 29 20 3d 3d 20 30 29 3b |(m & 0x1|) == 0);|
|000008f0| 20 6d 20 3e 3e 3d 20 31 | 29 0a 58 09 09 09 70 74 | m >>= 1|).X...pt|
|00000900| 77 6f 2b 2b 3b 0a 58 09 | 09 6d 75 6c 20 2a 3d 20 |wo++;.X.|.mul *= |
|00000910| 6d 3b 0a 58 09 09 69 66 | 20 28 6d 75 6c 20 3c 20 |m;.X..if| (mul < |
|00000920| 42 41 53 45 31 2f 32 29 | 0a 58 09 09 09 63 6f 6e |BASE1/2)|.X...con|
|00000930| 74 69 6e 75 65 3b 0a 58 | 09 09 7a 6d 75 6c 69 28 |tinue;.X|..zmuli(|
|00000940| 72 65 73 2c 20 6d 75 6c | 2c 20 26 74 65 6d 70 29 |res, mul|, &temp)|
|00000950| 3b 0a 58 09 09 66 72 65 | 65 68 28 72 65 73 2e 76 |;.X..fre|eh(res.v|
|00000960| 29 3b 0a 58 09 09 72 65 | 73 20 3d 20 74 65 6d 70 |);.X..re|s = temp|
|00000970| 3b 0a 58 09 09 6d 75 6c | 20 3d 20 31 3b 0a 58 09 |;.X..mul| = 1;.X.|
|00000980| 7d 0a 58 09 2f 2a 0a 58 | 09 20 2a 20 4d 75 6c 74 |}.X./*.X|. * Mult|
|00000990| 69 70 6c 79 20 62 79 20 | 74 68 65 20 72 65 6d 61 |iply by |the rema|
|000009a0| 69 6e 69 6e 67 20 76 61 | 6c 75 65 2c 20 74 68 65 |ining va|lue, the|
|000009b0| 6e 20 73 63 61 6c 65 20 | 72 65 73 75 6c 74 20 62 |n scale |result b|
|000009c0| 79 0a 58 09 20 2a 20 74 | 68 65 20 70 72 6f 70 65 |y.X. * t|he prope|
|000009d0| 72 20 70 6f 77 65 72 20 | 6f 66 20 74 77 6f 2e 0a |r power |of two..|
|000009e0| 58 09 20 2a 2f 0a 58 09 | 69 66 20 28 6d 75 6c 20 |X. */.X.|if (mul |
|000009f0| 3e 20 31 29 20 7b 0a 58 | 09 09 7a 6d 75 6c 69 28 |> 1) {.X|..zmuli(|
|00000a00| 72 65 73 2c 20 6d 75 6c | 2c 20 26 74 65 6d 70 29 |res, mul|, &temp)|
|00000a10| 3b 0a 58 09 09 66 72 65 | 65 68 28 72 65 73 2e 76 |;.X..fre|eh(res.v|
|00000a20| 29 3b 0a 58 09 09 72 65 | 73 20 3d 20 74 65 6d 70 |);.X..re|s = temp|
|00000a30| 3b 0a 58 09 7d 0a 58 09 | 7a 73 68 69 66 74 28 72 |;.X.}.X.|zshift(r|
|00000a40| 65 73 2c 20 70 74 77 6f | 2c 20 26 74 65 6d 70 29 |es, ptwo|, &temp)|
|00000a50| 3b 0a 58 09 66 72 65 65 | 68 28 72 65 73 2e 76 29 |;.X.free|h(res.v)|
|00000a60| 3b 0a 58 09 2a 64 65 73 | 74 20 3d 20 74 65 6d 70 |;.X.*des|t = temp|
|00000a70| 3b 0a 58 7d 0a 58 0a 58 | 0a 58 2f 2a 0a 58 20 2a |;.X}.X.X|.X/*.X *|
|00000a80| 20 43 6f 6d 70 75 74 65 | 20 74 68 65 20 70 72 6f | Compute| the pro|
|00000a90| 64 75 63 74 20 6f 66 20 | 74 68 65 20 70 72 69 6d |duct of |the prim|
|00000aa0| 65 73 20 75 70 20 74 6f | 20 74 68 65 20 73 70 65 |es up to| the spe|
|00000ab0| 63 69 66 69 65 64 20 6e | 75 6d 62 65 72 2e 0a 58 |cified n|umber..X|
|00000ac0| 20 2a 2f 0a 58 76 6f 69 | 64 0a 58 7a 70 66 61 63 | */.Xvoi|d.Xzpfac|
|00000ad0| 74 28 7a 2c 20 64 65 73 | 74 29 0a 58 09 5a 56 41 |t(z, des|t).X.ZVA|
|00000ae0| 4c 55 45 20 7a 2c 20 2a | 64 65 73 74 3b 0a 58 7b |LUE z, *|dest;.X{|
|00000af0| 0a 58 09 6c 6f 6e 67 20 | 6e 3b 09 09 09 2f 2a 20 |.X.long |n;.../* |
|00000b00| 6c 69 6d 69 74 69 6e 67 | 20 6e 75 6d 62 65 72 20 |limiting| number |
|00000b10| 74 6f 20 6d 75 6c 74 69 | 70 6c 79 20 62 79 20 2a |to multi|ply by *|
|00000b20| 2f 0a 58 09 6c 6f 6e 67 | 20 70 3b 09 09 09 2f 2a |/.X.long| p;.../*|
|00000b30| 20 63 75 72 72 65 6e 74 | 20 70 72 69 6d 65 20 2a | current| prime *|
|00000b40| 2f 0a 58 09 6c 6f 6e 67 | 20 69 3b 09 09 09 2f 2a |/.X.long| i;.../*|
|00000b50| 20 74 65 73 74 20 76 61 | 6c 75 65 20 2a 2f 0a 58 | test va|lue */.X|
|00000b60| 09 6c 6f 6e 67 20 6d 75 | 6c 3b 09 09 2f 2a 20 63 |.long mu|l;../* c|
|00000b70| 6f 6c 6c 65 63 74 65 64 | 20 76 61 6c 75 65 20 74 |ollected| value t|
|00000b80| 6f 20 6d 75 6c 74 69 70 | 6c 79 20 62 79 20 2a 2f |o multip|ly by */|
|00000b90| 0a 58 09 5a 56 41 4c 55 | 45 20 72 65 73 2c 20 74 |.X.ZVALU|E res, t|
|00000ba0| 65 6d 70 3b 0a 58 0a 58 | 09 69 66 20 28 69 73 6e |emp;.X.X|.if (isn|
|00000bb0| 65 67 28 7a 29 29 0a 58 | 09 09 65 72 72 6f 72 28 |eg(z)).X|..error(|
|00000bc0| 22 4e 65 67 61 74 69 76 | 65 20 61 72 67 75 6d 65 |"Negativ|e argume|
|00000bd0| 6e 74 20 66 6f 72 20 66 | 61 63 74 6f 72 69 61 6c |nt for f|actorial|
|00000be0| 22 29 3b 0a 58 09 69 66 | 20 28 69 73 62 69 67 28 |");.X.if| (isbig(|
|00000bf0| 7a 29 29 0a 58 09 09 65 | 72 72 6f 72 28 22 56 65 |z)).X..e|rror("Ve|
|00000c00| 72 79 20 6c 61 72 67 65 | 20 66 61 63 74 6f 72 69 |ry large| factori|
|00000c10| 61 6c 22 29 3b 0a 58 09 | 6e 20 3d 20 28 69 73 74 |al");.X.|n = (ist|
|00000c20| 69 6e 79 28 7a 29 20 3f | 20 7a 31 74 6f 6c 28 7a |iny(z) ?| z1tol(z|
|00000c30| 29 20 3a 20 7a 32 74 6f | 6c 28 7a 29 29 3b 0a 58 |) : z2to|l(z));.X|
|00000c40| 09 2f 2a 0a 58 09 20 2a | 20 4d 75 6c 74 69 70 6c |./*.X. *| Multipl|
|00000c50| 79 20 62 79 20 74 68 65 | 20 70 72 69 6d 65 73 20 |y by the| primes |
|00000c60| 69 6e 20 6f 72 64 65 72 | 2c 20 63 6f 6c 6c 65 63 |in order|, collec|
|00000c70| 74 69 6e 67 20 6d 75 6c | 74 69 70 6c 65 20 6e 75 |ting mul|tiple nu|
|00000c80| 6d 62 65 72 73 0a 58 09 | 20 2a 20 74 6f 67 65 74 |mbers.X.| * toget|
|00000c90| 68 65 72 20 75 6e 74 69 | 6c 20 74 68 65 72 65 20 |her unti|l there |
|00000ca0| 69 73 20 61 20 63 68 61 | 6e 67 65 20 6f 66 20 6f |is a cha|nge of o|
|00000cb0| 76 65 72 66 6c 6f 77 2e | 0a 58 09 20 2a 2f 0a 58 |verflow.|.X. */.X|
|00000cc0| 09 6d 75 6c 20 3d 20 31 | 20 2b 20 28 6e 20 3e 20 |.mul = 1| + (n > |
|00000cd0| 31 29 3b 0a 58 09 72 65 | 73 20 3d 20 5f 6f 6e 65 |1);.X.re|s = _one|
|00000ce0| 5f 3b 0a 58 09 66 6f 72 | 20 28 70 20 3d 20 33 3b |_;.X.for| (p = 3;|
|00000cf0| 20 70 20 3c 3d 20 6e 3b | 20 70 20 2b 3d 20 32 29 | p <= n;| p += 2)|
|00000d00| 20 7b 0a 58 09 09 66 6f | 72 20 28 69 20 3d 20 33 | {.X..fo|r (i = 3|
|00000d10| 3b 20 28 69 20 2a 20 69 | 29 20 3c 3d 20 70 3b 20 |; (i * i|) <= p; |
|00000d20| 69 20 2b 3d 20 32 29 20 | 7b 0a 58 09 09 09 69 66 |i += 2) |{.X...if|
|00000d30| 20 28 28 70 20 25 20 69 | 29 20 3d 3d 20 30 29 0a | ((p % i|) == 0).|
|00000d40| 58 09 09 09 09 67 6f 74 | 6f 20 6e 65 78 74 3b 0a |X....got|o next;.|
|00000d50| 58 09 09 7d 0a 58 09 09 | 6d 75 6c 20 2a 3d 20 70 |X..}.X..|mul *= p|
|00000d60| 3b 0a 58 09 09 69 66 20 | 28 6d 75 6c 20 3c 20 42 |;.X..if |(mul < B|
|00000d70| 41 53 45 31 2f 32 29 0a | 58 09 09 09 63 6f 6e 74 |ASE1/2).|X...cont|
|00000d80| 69 6e 75 65 3b 0a 58 09 | 09 7a 6d 75 6c 69 28 72 |inue;.X.|.zmuli(r|
|00000d90| 65 73 2c 20 6d 75 6c 2c | 20 26 74 65 6d 70 29 3b |es, mul,| &temp);|
|00000da0| 0a 58 09 09 66 72 65 65 | 68 28 72 65 73 2e 76 29 |.X..free|h(res.v)|
|00000db0| 3b 0a 58 09 09 72 65 73 | 20 3d 20 74 65 6d 70 3b |;.X..res| = temp;|
|00000dc0| 0a 58 09 09 6d 75 6c 20 | 3d 20 31 3b 0a 58 6e 65 |.X..mul |= 1;.Xne|
|00000dd0| 78 74 3a 20 3b 0a 58 09 | 7d 0a 58 09 2f 2a 0a 58 |xt: ;.X.|}.X./*.X|
|00000de0| 09 20 2a 20 4d 75 6c 74 | 69 70 6c 79 20 62 79 20 |. * Mult|iply by |
|00000df0| 74 68 65 20 66 69 6e 61 | 6c 20 76 61 6c 75 65 20 |the fina|l value |
|00000e00| 69 66 20 61 6e 79 2e 0a | 58 09 20 2a 2f 0a 58 09 |if any..|X. */.X.|
|00000e10| 69 66 20 28 6d 75 6c 20 | 3e 20 31 29 20 7b 0a 58 |if (mul |> 1) {.X|
|00000e20| 09 09 7a 6d 75 6c 69 28 | 72 65 73 2c 20 6d 75 6c |..zmuli(|res, mul|
|00000e30| 2c 20 26 74 65 6d 70 29 | 3b 0a 58 09 09 66 72 65 |, &temp)|;.X..fre|
|00000e40| 65 68 28 72 65 73 2e 76 | 29 3b 0a 58 09 09 72 65 |eh(res.v|);.X..re|
|00000e50| 73 20 3d 20 74 65 6d 70 | 3b 0a 58 09 7d 0a 58 09 |s = temp|;.X.}.X.|
|00000e60| 2a 64 65 73 74 20 3d 20 | 72 65 73 3b 0a 58 7d 0a |*dest = |res;.X}.|
|00000e70| 58 0a 58 0a 58 2f 2a 0a | 58 20 2a 20 43 6f 6d 70 |X.X.X/*.|X * Comp|
|00000e80| 75 74 65 20 74 68 65 20 | 6c 65 61 73 74 20 63 6f |ute the |least co|
|00000e90| 6d 6d 6f 6e 20 6d 75 6c | 74 69 70 6c 65 20 6f 66 |mmon mul|tiple of|
|00000ea0| 20 61 6c 6c 20 74 68 65 | 20 6e 75 6d 62 65 72 73 | all the| numbers|
|00000eb0| 20 75 70 20 74 6f 20 74 | 68 65 0a 58 20 2a 20 73 | up to t|he.X * s|
|00000ec0| 70 65 63 69 66 69 65 64 | 20 6e 75 6d 62 65 72 2e |pecified| number.|
|00000ed0| 0a 58 20 2a 2f 0a 58 76 | 6f 69 64 0a 58 7a 6c 63 |.X */.Xv|oid.Xzlc|
|00000ee0| 6d 66 61 63 74 28 7a 2c | 20 64 65 73 74 29 0a 58 |mfact(z,| dest).X|
|00000ef0| 09 5a 56 41 4c 55 45 20 | 7a 2c 20 2a 64 65 73 74 |.ZVALUE |z, *dest|
|00000f00| 3b 0a 58 7b 0a 58 09 6c | 6f 6e 67 20 6e 3b 09 09 |;.X{.X.l|ong n;..|
|00000f10| 09 2f 2a 20 6c 69 6d 69 | 74 69 6e 67 20 6e 75 6d |./* limi|ting num|
|00000f20| 62 65 72 20 74 6f 20 6d | 75 6c 74 69 70 6c 79 20 |ber to m|ultiply |
|00000f30| 62 79 20 2a 2f 0a 58 09 | 6c 6f 6e 67 20 70 3b 09 |by */.X.|long p;.|
|00000f40| 09 09 2f 2a 20 63 75 72 | 72 65 6e 74 20 70 72 69 |../* cur|rent pri|
|00000f50| 6d 65 20 2a 2f 0a 58 09 | 6c 6f 6e 67 20 70 70 3b |me */.X.|long pp;|
|00000f60| 09 09 2f 2a 20 70 6f 77 | 65 72 20 6f 66 20 70 72 |../* pow|er of pr|
|00000f70| 69 6d 65 20 2a 2f 0a 58 | 09 6c 6f 6e 67 20 69 3b |ime */.X|.long i;|
|00000f80| 09 09 09 2f 2a 20 74 65 | 73 74 20 76 61 6c 75 65 |.../* te|st value|
|00000f90| 20 2a 2f 0a 58 09 5a 56 | 41 4c 55 45 20 72 65 73 | */.X.ZV|ALUE res|
|00000fa0| 2c 20 74 65 6d 70 3b 0a | 58 0a 58 09 69 66 20 28 |, temp;.|X.X.if (|
|00000fb0| 69 73 6e 65 67 28 7a 29 | 20 7c 7c 20 69 73 7a 65 |isneg(z)| || isze|
|00000fc0| 72 6f 28 7a 29 29 0a 58 | 09 09 65 72 72 6f 72 28 |ro(z)).X|..error(|
|00000fd0| 22 4e 6f 6e 2d 70 6f 73 | 69 74 69 76 65 20 61 72 |"Non-pos|itive ar|
|00000fe0| 67 75 6d 65 6e 74 20 66 | 6f 72 20 6c 63 6d 66 61 |gument f|or lcmfa|
|00000ff0| 63 74 22 29 3b 0a 58 09 | 69 66 20 28 69 73 62 69 |ct");.X.|if (isbi|
|00001000| 67 28 7a 29 29 0a 58 09 | 09 65 72 72 6f 72 28 22 |g(z)).X.|.error("|
|00001010| 56 65 72 79 20 6c 61 72 | 67 65 20 6c 63 6d 66 61 |Very lar|ge lcmfa|
|00001020| 63 74 22 29 3b 0a 58 09 | 6e 20 3d 20 28 69 73 74 |ct");.X.|n = (ist|
|00001030| 69 6e 79 28 7a 29 20 3f | 20 7a 31 74 6f 6c 28 7a |iny(z) ?| z1tol(z|
|00001040| 29 20 3a 20 7a 32 74 6f | 6c 28 7a 29 29 3b 0a 58 |) : z2to|l(z));.X|
|00001050| 09 2f 2a 0a 58 09 20 2a | 20 4d 75 6c 74 69 70 6c |./*.X. *| Multipl|
|00001060| 79 20 62 79 20 70 6f 77 | 65 72 73 20 6f 66 20 74 |y by pow|ers of t|
|00001070| 68 65 20 6e 65 63 65 73 | 73 61 72 79 20 6f 64 64 |he neces|sary odd|
|00001080| 20 70 72 69 6d 65 73 20 | 69 6e 20 6f 72 64 65 72 | primes |in order|
|00001090| 2e 0a 58 09 20 2a 20 54 | 68 65 20 70 6f 77 65 72 |..X. * T|he power|
|000010a0| 20 66 6f 72 20 65 61 63 | 68 20 70 72 69 6d 65 20 | for eac|h prime |
|000010b0| 69 73 20 74 68 65 20 68 | 69 67 68 65 73 74 20 6f |is the h|ighest o|
|000010c0| 6e 65 20 77 68 69 63 68 | 20 69 73 20 6e 6f 74 0a |ne which| is not.|
|000010d0| 58 09 20 2a 20 6d 6f 72 | 65 20 74 68 61 6e 20 74 |X. * mor|e than t|
|000010e0| 68 65 20 73 70 65 63 69 | 66 69 65 64 20 6e 75 6d |he speci|fied num|
|000010f0| 62 65 72 2e 0a 58 09 20 | 2a 2f 0a 58 09 72 65 73 |ber..X. |*/.X.res|
|00001100| 20 3d 20 5f 6f 6e 65 5f | 3b 0a 58 09 66 6f 72 20 | = _one_|;.X.for |
|00001110| 28 70 20 3d 20 33 3b 20 | 70 20 3c 3d 20 6e 3b 20 |(p = 3; |p <= n; |
|00001120| 70 20 2b 3d 20 32 29 20 | 7b 0a 58 09 09 66 6f 72 |p += 2) |{.X..for|
|00001130| 20 28 69 20 3d 20 33 3b | 20 28 69 20 2a 20 69 29 | (i = 3;| (i * i)|
|00001140| 20 3c 3d 20 70 3b 20 69 | 20 2b 3d 20 32 29 20 7b | <= p; i| += 2) {|
|00001150| 0a 58 09 09 09 69 66 20 | 28 28 70 20 25 20 69 29 |.X...if |((p % i)|
|00001160| 20 3d 3d 20 30 29 0a 58 | 09 09 09 09 67 6f 74 6f | == 0).X|....goto|
|00001170| 20 6e 65 78 74 3b 0a 58 | 09 09 7d 0a 58 09 09 69 | next;.X|..}.X..i|
|00001180| 20 3d 20 70 3b 0a 58 09 | 09 77 68 69 6c 65 20 28 | = p;.X.|.while (|
|00001190| 69 20 3c 3d 20 6e 29 20 | 7b 0a 58 09 09 09 70 70 |i <= n) |{.X...pp|
|000011a0| 20 3d 20 69 3b 0a 58 09 | 09 09 69 20 2a 3d 20 70 | = i;.X.|..i *= p|
|000011b0| 3b 0a 58 09 09 7d 0a 58 | 09 09 7a 6d 75 6c 69 28 |;.X..}.X|..zmuli(|
|000011c0| 72 65 73 2c 20 70 70 2c | 20 26 74 65 6d 70 29 3b |res, pp,| &temp);|
|000011d0| 0a 58 09 09 66 72 65 65 | 68 28 72 65 73 2e 76 29 |.X..free|h(res.v)|
|000011e0| 3b 0a 58 09 09 72 65 73 | 20 3d 20 74 65 6d 70 3b |;.X..res| = temp;|
|000011f0| 0a 58 6e 65 78 74 3a 20 | 3b 0a 58 09 7d 0a 58 09 |.Xnext: |;.X.}.X.|
|00001200| 2f 2a 0a 58 09 20 2a 20 | 46 69 6e 69 73 68 20 62 |/*.X. * |Finish b|
|00001210| 79 20 73 63 61 6c 69 6e | 67 20 62 79 20 74 68 65 |y scalin|g by the|
|00001220| 20 6e 65 63 65 73 73 61 | 72 79 20 70 6f 77 65 72 | necessa|ry power|
|00001230| 20 6f 66 20 74 77 6f 2e | 0a 58 09 20 2a 2f 0a 58 | of two.|.X. */.X|
|00001240| 09 7a 73 68 69 66 74 28 | 72 65 73 2c 20 7a 68 69 |.zshift(|res, zhi|
|00001250| 67 68 62 69 74 28 7a 29 | 2c 20 64 65 73 74 29 3b |ghbit(z)|, dest);|
|00001260| 0a 58 09 66 72 65 65 68 | 28 72 65 73 2e 76 29 3b |.X.freeh|(res.v);|
|00001270| 0a 58 7d 0a 58 0a 58 0a | 58 2f 2a 0a 58 20 2a 20 |.X}.X.X.|X/*.X * |
|00001280| 43 6f 6d 70 75 74 65 20 | 74 68 65 20 70 65 72 6d |Compute |the perm|
|00001290| 75 61 74 69 6f 6e 20 66 | 75 6e 63 74 69 6f 6e 20 |uation f|unction |
|000012a0| 20 4d 21 20 2f 20 28 4d | 20 2d 20 4e 29 21 2e 0a | M! / (M| - N)!..|
|000012b0| 58 20 2a 2f 0a 58 76 6f | 69 64 0a 58 7a 70 65 72 |X */.Xvo|id.Xzper|
|000012c0| 6d 28 7a 31 2c 20 7a 32 | 2c 20 72 65 73 29 0a 58 |m(z1, z2|, res).X|
|000012d0| 09 5a 56 41 4c 55 45 20 | 7a 31 2c 20 7a 32 2c 20 |.ZVALUE |z1, z2, |
|000012e0| 2a 72 65 73 3b 0a 58 7b | 0a 58 09 6c 6f 6e 67 20 |*res;.X{|.X.long |
|000012f0| 63 6f 75 6e 74 3b 0a 58 | 09 5a 56 41 4c 55 45 20 |count;.X|.ZVALUE |
|00001300| 63 75 72 2c 20 74 6d 70 | 2c 20 61 6e 73 3b 0a 58 |cur, tmp|, ans;.X|
|00001310| 0a 58 09 69 66 20 28 69 | 73 6e 65 67 28 7a 31 29 |.X.if (i|sneg(z1)|
|00001320| 20 7c 7c 20 69 73 6e 65 | 67 28 7a 32 29 29 0a 58 | || isne|g(z2)).X|
|00001330| 09 09 65 72 72 6f 72 28 | 22 4e 65 67 61 74 69 76 |..error(|"Negativ|
|00001340| 65 20 61 72 67 75 6d 65 | 6e 74 20 66 6f 72 20 70 |e argume|nt for p|
|00001350| 65 72 6d 75 74 61 74 69 | 6f 6e 22 29 3b 0a 58 09 |ermutati|on");.X.|
|00001360| 69 66 20 28 7a 72 65 6c | 28 7a 31 2c 20 7a 32 29 |if (zrel|(z1, z2)|
|00001370| 20 3c 20 30 29 0a 58 09 | 09 65 72 72 6f 72 28 22 | < 0).X.|.error("|
|00001380| 53 65 63 6f 6e 64 20 61 | 72 67 20 6c 61 72 67 65 |Second a|rg large|
|00001390| 72 20 74 68 61 6e 20 66 | 69 72 73 74 20 69 6e 20 |r than f|irst in |
|000013a0| 70 65 72 6d 75 74 61 74 | 69 6f 6e 22 29 3b 0a 58 |permutat|ion");.X|
|000013b0| 09 69 66 20 28 69 73 62 | 69 67 28 7a 32 29 29 0a |.if (isb|ig(z2)).|
|000013c0| 58 09 09 65 72 72 6f 72 | 28 22 56 65 72 79 20 6c |X..error|("Very l|
|000013d0| 61 72 67 65 20 70 65 72 | 6d 75 74 61 74 69 6f 6e |arge per|mutation|
|000013e0| 22 29 3b 0a 58 09 63 6f | 75 6e 74 20 3d 20 28 69 |");.X.co|unt = (i|
|000013f0| 73 74 69 6e 79 28 7a 32 | 29 20 3f 20 7a 31 74 6f |stiny(z2|) ? z1to|
|00001400| 6c 28 7a 32 29 20 3a 20 | 7a 32 74 6f 6c 28 7a 32 |l(z2) : |z2tol(z2|
|00001410| 29 29 3b 0a 58 09 7a 63 | 6f 70 79 28 7a 31 2c 20 |));.X.zc|opy(z1, |
|00001420| 26 61 6e 73 29 3b 0a 58 | 09 7a 73 75 62 28 7a 31 |&ans);.X|.zsub(z1|
|00001430| 2c 20 5f 6f 6e 65 5f 2c | 20 26 63 75 72 29 3b 0a |, _one_,| &cur);.|
|00001440| 58 09 77 68 69 6c 65 20 | 28 2d 2d 63 6f 75 6e 74 |X.while |(--count|
|00001450| 20 3e 20 30 29 20 7b 0a | 58 09 09 7a 6d 75 6c 28 | > 0) {.|X..zmul(|
|00001460| 61 6e 73 2c 20 63 75 72 | 2c 20 26 74 6d 70 29 3b |ans, cur|, &tmp);|
|00001470| 0a 58 09 09 66 72 65 65 | 68 28 61 6e 73 2e 76 29 |.X..free|h(ans.v)|
|00001480| 3b 0a 58 09 09 61 6e 73 | 20 3d 20 74 6d 70 3b 0a |;.X..ans| = tmp;.|
|00001490| 58 09 09 7a 73 75 62 28 | 63 75 72 2c 20 5f 6f 6e |X..zsub(|cur, _on|
|000014a0| 65 5f 2c 20 26 74 6d 70 | 29 3b 0a 58 09 09 66 72 |e_, &tmp|);.X..fr|
|000014b0| 65 65 68 28 63 75 72 2e | 76 29 3b 0a 58 09 09 63 |eeh(cur.|v);.X..c|
|000014c0| 75 72 20 3d 20 74 6d 70 | 3b 0a 58 09 7d 0a 58 09 |ur = tmp|;.X.}.X.|
|000014d0| 66 72 65 65 68 28 63 75 | 72 2e 76 29 3b 0a 58 09 |freeh(cu|r.v);.X.|
|000014e0| 2a 72 65 73 20 3d 20 61 | 6e 73 3b 0a 58 7d 0a 58 |*res = a|ns;.X}.X|
|000014f0| 0a 58 0a 58 2f 2a 0a 58 | 20 2a 20 43 6f 6d 70 75 |.X.X/*.X| * Compu|
|00001500| 74 65 20 74 68 65 20 63 | 6f 6d 62 69 6e 61 74 6f |te the c|ombinato|
|00001510| 72 69 61 6c 20 66 75 6e | 63 74 69 6f 6e 20 20 4d |rial fun|ction M|
|00001520| 21 20 2f 20 28 20 4e 21 | 20 2a 20 28 4d 20 2d 20 |! / ( N!| * (M - |
|00001530| 4e 29 21 20 29 2e 0a 58 | 20 2a 2f 0a 58 76 6f 69 |N)! )..X| */.Xvoi|
|00001540| 64 0a 58 7a 63 6f 6d 62 | 28 7a 31 2c 20 7a 32 2c |d.Xzcomb|(z1, z2,|
|00001550| 20 72 65 73 29 0a 58 09 | 5a 56 41 4c 55 45 20 7a | res).X.|ZVALUE z|
|00001560| 31 2c 20 7a 32 2c 20 2a | 72 65 73 3b 0a 58 7b 0a |1, z2, *|res;.X{.|
|00001570| 58 09 5a 56 41 4c 55 45 | 20 61 6e 73 3b 0a 58 09 |X.ZVALUE| ans;.X.|
|00001580| 5a 56 41 4c 55 45 20 6d | 75 6c 2c 20 64 69 76 2c |ZVALUE m|ul, div,|
|00001590| 20 74 65 6d 70 3b 0a 58 | 09 46 55 4c 4c 20 63 6f | temp;.X|.FULL co|
|000015a0| 75 6e 74 2c 20 69 3b 0a | 58 09 48 41 4c 46 20 64 |unt, i;.|X.HALF d|
|000015b0| 68 5b 32 5d 3b 0a 58 0a | 58 09 69 66 20 28 69 73 |h[2];.X.|X.if (is|
|000015c0| 6e 65 67 28 7a 31 29 20 | 7c 7c 20 69 73 6e 65 67 |neg(z1) ||| isneg|
|000015d0| 28 7a 32 29 29 0a 58 09 | 09 65 72 72 6f 72 28 22 |(z2)).X.|.error("|
|000015e0| 4e 65 67 61 74 69 76 65 | 20 61 72 67 75 6d 65 6e |Negative| argumen|
|000015f0| 74 20 66 6f 72 20 63 6f | 6d 62 69 6e 61 74 6f 72 |t for co|mbinator|
|00001600| 69 61 6c 22 29 3b 0a 58 | 09 7a 73 75 62 28 7a 31 |ial");.X|.zsub(z1|
|00001610| 2c 20 7a 32 2c 20 26 74 | 65 6d 70 29 3b 0a 58 09 |, z2, &t|emp);.X.|
|00001620| 69 66 20 28 69 73 6e 65 | 67 28 74 65 6d 70 29 29 |if (isne|g(temp))|
|00001630| 20 7b 0a 58 09 09 66 72 | 65 65 68 28 74 65 6d 70 | {.X..fr|eeh(temp|
|00001640| 2e 76 29 3b 0a 58 09 09 | 65 72 72 6f 72 28 22 53 |.v);.X..|error("S|
|00001650| 65 63 6f 6e 64 20 61 72 | 67 20 6c 61 72 67 65 72 |econd ar|g larger|
|00001660| 20 74 68 61 6e 20 66 69 | 72 73 74 20 66 6f 72 20 | than fi|rst for |
|00001670| 63 6f 6d 62 69 6e 61 74 | 6f 72 69 61 6c 22 29 3b |combinat|orial");|
|00001680| 0a 58 09 7d 0a 58 09 69 | 66 20 28 69 73 62 69 67 |.X.}.X.i|f (isbig|
|00001690| 28 7a 32 29 20 26 26 20 | 69 73 62 69 67 28 74 65 |(z2) && |isbig(te|
|000016a0| 6d 70 29 29 20 7b 0a 58 | 09 09 66 72 65 65 68 28 |mp)) {.X|..freeh(|
|000016b0| 74 65 6d 70 2e 76 29 3b | 0a 58 09 09 65 72 72 6f |temp.v);|.X..erro|
|000016c0| 72 28 22 56 65 72 79 20 | 6c 61 72 67 65 20 63 6f |r("Very |large co|
|000016d0| 6d 62 69 6e 61 74 6f 72 | 69 61 6c 22 29 3b 0a 58 |mbinator|ial");.X|
|000016e0| 09 7d 0a 58 09 63 6f 75 | 6e 74 20 3d 20 28 69 73 |.}.X.cou|nt = (is|
|000016f0| 74 69 6e 79 28 7a 32 29 | 20 3f 20 7a 31 74 6f 6c |tiny(z2)| ? z1tol|
|00001700| 28 7a 32 29 20 3a 20 7a | 32 74 6f 6c 28 7a 32 29 |(z2) : z|2tol(z2)|
|00001710| 29 3b 0a 58 09 69 20 3d | 20 28 69 73 74 69 6e 79 |);.X.i =| (istiny|
|00001720| 28 74 65 6d 70 29 20 3f | 20 7a 31 74 6f 6c 28 74 |(temp) ?| z1tol(t|
|00001730| 65 6d 70 29 20 3a 20 7a | 32 74 6f 6c 28 74 65 6d |emp) : z|2tol(tem|
|00001740| 70 29 29 3b 0a 58 09 69 | 66 20 28 69 73 62 69 67 |p));.X.i|f (isbig|
|00001750| 28 7a 32 29 20 7c 7c 20 | 28 21 69 73 62 69 67 28 |(z2) || |(!isbig(|
|00001760| 74 65 6d 70 29 20 26 26 | 20 28 69 20 3c 20 63 6f |temp) &&| (i < co|
|00001770| 75 6e 74 29 29 29 0a 58 | 09 09 63 6f 75 6e 74 20 |unt))).X|..count |
|00001780| 3d 20 69 3b 0a 58 09 66 | 72 65 65 68 28 74 65 6d |= i;.X.f|reeh(tem|
|00001790| 70 2e 76 29 3b 0a 58 09 | 6d 75 6c 20 3d 20 7a 31 |p.v);.X.|mul = z1|
|000017a0| 3b 0a 58 09 64 69 76 2e | 73 69 67 6e 20 3d 20 30 |;.X.div.|sign = 0|
|000017b0| 3b 0a 58 09 64 69 76 2e | 76 20 3d 20 64 68 3b 0a |;.X.div.|v = dh;.|
|000017c0| 58 09 61 6e 73 20 3d 20 | 5f 6f 6e 65 5f 3b 0a 58 |X.ans = |_one_;.X|
|000017d0| 09 66 6f 72 20 28 69 20 | 3d 20 31 3b 20 69 20 3c |.for (i |= 1; i <|
|000017e0| 3d 20 63 6f 75 6e 74 3b | 20 69 2b 2b 29 20 7b 0a |= count;| i++) {.|
|000017f0| 58 09 09 64 68 5b 30 5d | 20 3d 20 69 20 26 20 42 |X..dh[0]| = i & B|
|00001800| 41 53 45 31 3b 0a 58 09 | 09 64 68 5b 31 5d 20 3d |ASE1;.X.|.dh[1] =|
|00001810| 20 69 20 2f 20 42 41 53 | 45 3b 0a 58 09 09 64 69 | i / BAS|E;.X..di|
|00001820| 76 2e 6c 65 6e 20 3d 20 | 31 20 2b 20 28 64 68 5b |v.len = |1 + (dh[|
|00001830| 31 5d 20 21 3d 20 30 29 | 3b 0a 58 09 09 7a 6d 75 |1] != 0)|;.X..zmu|
|00001840| 6c 28 61 6e 73 2c 20 6d | 75 6c 2c 20 26 74 65 6d |l(ans, m|ul, &tem|
|00001850| 70 29 3b 0a 58 09 09 66 | 72 65 65 68 28 61 6e 73 |p);.X..f|reeh(ans|
|00001860| 2e 76 29 3b 0a 58 09 09 | 7a 71 75 6f 28 74 65 6d |.v);.X..|zquo(tem|
|00001870| 70 2c 20 64 69 76 2c 20 | 26 61 6e 73 29 3b 0a 58 |p, div, |&ans);.X|
|00001880| 09 09 66 72 65 65 68 28 | 74 65 6d 70 2e 76 29 3b |..freeh(|temp.v);|
|00001890| 0a 58 09 09 7a 73 75 62 | 28 6d 75 6c 2c 20 5f 6f |.X..zsub|(mul, _o|
|000018a0| 6e 65 5f 2c 20 26 74 65 | 6d 70 29 3b 0a 58 09 09 |ne_, &te|mp);.X..|
|000018b0| 69 66 20 28 6d 75 6c 2e | 76 20 21 3d 20 7a 31 2e |if (mul.|v != z1.|
|000018c0| 76 29 0a 58 09 09 09 66 | 72 65 65 68 28 6d 75 6c |v).X...f|reeh(mul|
|000018d0| 2e 76 29 3b 0a 58 09 09 | 6d 75 6c 20 3d 20 74 65 |.v);.X..|mul = te|
|000018e0| 6d 70 3b 0a 58 09 7d 0a | 58 09 69 66 20 28 6d 75 |mp;.X.}.|X.if (mu|
|000018f0| 6c 2e 76 20 21 3d 20 7a | 31 2e 76 29 0a 58 09 09 |l.v != z|1.v).X..|
|00001900| 66 72 65 65 68 28 6d 75 | 6c 2e 76 29 3b 0a 58 09 |freeh(mu|l.v);.X.|
|00001910| 2a 72 65 73 20 3d 20 61 | 6e 73 3b 0a 58 7d 0a 58 |*res = a|ns;.X}.X|
|00001920| 0a 58 0a 58 2f 2a 0a 58 | 20 2a 20 50 65 72 66 6f |.X.X/*.X| * Perfo|
|00001930| 72 6d 20 61 20 70 72 6f | 62 61 62 69 6c 69 73 74 |rm a pro|babilist|
|00001940| 69 63 20 70 72 69 6d 61 | 6c 69 74 79 20 74 65 73 |ic prima|lity tes|
|00001950| 74 20 28 61 6c 67 6f 72 | 69 74 68 6d 20 50 20 69 |t (algor|ithm P i|
|00001960| 6e 20 4b 6e 75 74 68 29 | 2e 0a 58 20 2a 20 52 65 |n Knuth)|..X * Re|
|00001970| 74 75 72 6e 73 20 46 41 | 4c 53 45 20 69 66 20 64 |turns FA|LSE if d|
|00001980| 65 66 69 6e 69 74 65 6c | 79 20 6e 6f 74 20 70 72 |efinitel|y not pr|
|00001990| 69 6d 65 2c 20 6f 72 20 | 54 52 55 45 20 69 66 20 |ime, or |TRUE if |
|000019a0| 70 72 6f 62 61 62 6c 79 | 20 70 72 69 6d 65 2e 0a |probably| prime..|
|000019b0| 58 20 2a 20 43 6f 75 6e | 74 20 64 65 74 65 72 6d |X * Coun|t determ|
|000019c0| 69 6e 65 73 20 68 6f 77 | 20 6d 61 6e 79 20 74 69 |ines how| many ti|
|000019d0| 6d 65 73 20 74 6f 20 63 | 68 65 63 6b 20 66 6f 72 |mes to c|heck for|
|000019e0| 20 70 72 69 6d 61 6c 69 | 74 79 2e 0a 58 20 2a 20 | primali|ty..X * |
|000019f0| 54 68 65 20 63 68 61 6e | 63 65 20 6f 66 20 61 20 |The chan|ce of a |
|00001a00| 6e 6f 6e 2d 70 72 69 6d | 65 20 70 61 73 73 69 6e |non-prim|e passin|
|00001a10| 67 20 74 68 69 73 20 74 | 65 73 74 20 69 73 20 6c |g this t|est is l|
|00001a20| 65 73 73 20 74 68 61 6e | 20 28 31 2f 34 29 5e 63 |ess than| (1/4)^c|
|00001a30| 6f 75 6e 74 2e 0a 58 20 | 2a 20 46 6f 72 20 65 78 |ount..X |* For ex|
|00001a40| 61 6d 70 6c 65 2c 20 61 | 20 63 6f 75 6e 74 20 6f |ample, a| count o|
|00001a50| 66 20 31 30 30 20 66 61 | 69 6c 73 20 66 6f 72 20 |f 100 fa|ils for |
|00001a60| 6f 6e 6c 79 20 31 20 69 | 6e 20 31 30 5e 36 30 20 |only 1 i|n 10^60 |
|00001a70| 6e 75 6d 62 65 72 73 2e | 0a 58 20 2a 2f 0a 58 42 |numbers.|.X */.XB|
|00001a80| 4f 4f 4c 0a 58 7a 70 72 | 69 6d 65 74 65 73 74 28 |OOL.Xzpr|imetest(|
|00001a90| 7a 2c 20 63 6f 75 6e 74 | 29 0a 58 09 5a 56 41 4c |z, count|).X.ZVAL|
|00001aa0| 55 45 20 7a 3b 09 09 2f | 2a 20 6e 75 6d 62 65 72 |UE z;../|* number|
|00001ab0| 20 74 6f 20 74 65 73 74 | 20 66 6f 72 20 70 72 69 | to test| for pri|
|00001ac0| 6d 65 6e 65 73 73 20 2a | 2f 0a 58 09 6c 6f 6e 67 |meness *|/.X.long|
|00001ad0| 20 63 6f 75 6e 74 3b 0a | 58 7b 0a 58 09 6c 6f 6e | count;.|X{.X.lon|
|00001ae0| 67 20 69 6a 2c 20 69 6b | 2c 20 69 78 3b 0a 58 09 |g ij, ik|, ix;.X.|
|00001af0| 5a 56 41 4c 55 45 20 7a | 6d 31 2c 20 7a 31 2c 20 |ZVALUE z|m1, z1, |
|00001b00| 7a 32 2c 20 7a 33 2c 20 | 7a 74 6d 70 3b 0a 58 09 |z2, z3, |ztmp;.X.|
|00001b10| 48 41 4c 46 20 76 61 6c | 5b 32 5d 3b 0a 58 0a 58 |HALF val|[2];.X.X|
|00001b20| 09 7a 2e 73 69 67 6e 20 | 3d 20 30 3b 0a 58 09 69 |.z.sign |= 0;.X.i|
|00001b30| 66 20 28 69 73 65 76 65 | 6e 28 7a 29 29 09 09 2f |f (iseve|n(z))../|
|00001b40| 2a 20 69 66 20 65 76 65 | 6e 2c 20 6e 6f 74 20 70 |* if eve|n, not p|
|00001b50| 72 69 6d 65 20 69 66 20 | 6e 6f 74 20 32 20 2a 2f |rime if |not 2 */|
|00001b60| 0a 58 09 09 72 65 74 75 | 72 6e 20 28 69 73 74 77 |.X..retu|rn (istw|
|00001b70| 6f 28 7a 29 20 21 3d 20 | 30 29 3b 0a 58 09 2f 2a |o(z) != |0);.X./*|
|00001b80| 0a 58 09 20 2a 20 53 65 | 65 20 69 66 20 74 68 65 |.X. * Se|e if the|
|00001b90| 20 6e 75 6d 62 65 72 20 | 69 73 20 73 6d 61 6c 6c | number |is small|
|00001ba0| 2c 20 61 6e 64 20 69 73 | 20 65 69 74 68 65 72 20 |, and is| either |
|00001bb0| 61 20 73 6d 61 6c 6c 20 | 70 72 69 6d 65 2c 0a 58 |a small |prime,.X|
|00001bc0| 09 20 2a 20 6f 72 20 69 | 73 20 64 69 76 69 73 69 |. * or i|s divisi|
|00001bd0| 62 6c 65 20 62 79 20 61 | 20 73 6d 61 6c 6c 20 70 |ble by a| small p|
|00001be0| 72 69 6d 65 2e 0a 58 09 | 20 2a 2f 0a 58 09 69 66 |rime..X.| */.X.if|
|00001bf0| 20 28 69 73 74 69 6e 79 | 28 7a 29 20 26 26 20 28 | (istiny|(z) && (|
|00001c00| 2a 7a 2e 76 20 3c 3d 20 | 28 48 41 4c 46 29 28 31 |*z.v <= |(HALF)(1|
|00001c10| 30 31 2a 31 30 31 2d 31 | 29 29 29 20 7b 0a 58 09 |01*101-1|))) {.X.|
|00001c20| 09 69 78 20 3d 20 2a 7a | 2e 76 3b 0a 58 09 09 66 |.ix = *z|.v;.X..f|
|00001c30| 6f 72 20 28 69 6b 20 3d | 20 33 3b 20 28 69 6b 20 |or (ik =| 3; (ik |
|00001c40| 3c 3d 20 39 37 29 20 26 | 26 20 28 28 69 6b 20 2a |<= 97) &|& ((ik *|
|00001c50| 20 69 6b 29 20 3c 3d 20 | 69 78 29 3b 20 69 6b 20 | ik) <= |ix); ik |
|00001c60| 2b 3d 20 32 29 0a 58 09 | 09 09 69 66 20 28 28 69 |+= 2).X.|..if ((i|
|00001c70| 78 20 25 20 69 6b 29 20 | 3d 3d 20 30 29 0a 58 09 |x % ik) |== 0).X.|
|00001c80| 09 09 09 72 65 74 75 72 | 6e 20 46 41 4c 53 45 3b |...retur|n FALSE;|
|00001c90| 0a 58 09 09 72 65 74 75 | 72 6e 20 54 52 55 45 3b |.X..retu|rn TRUE;|
|00001ca0| 0a 58 09 7d 0a 58 09 2f | 2a 0a 58 09 20 2a 20 53 |.X.}.X./|*.X. * S|
|00001cb0| 65 65 20 69 66 20 74 68 | 65 20 6e 75 6d 62 65 72 |ee if th|e number|
|00001cc0| 20 69 73 20 64 69 76 69 | 73 69 62 6c 65 20 62 79 | is divi|sible by|
|00001cd0| 20 6f 6e 65 20 6f 66 20 | 74 68 65 20 70 72 69 6d | one of |the prim|
|00001ce0| 65 73 20 33 2c 20 35 2c | 0a 58 09 20 2a 20 37 2c |es 3, 5,|.X. * 7,|
|00001cf0| 20 31 31 2c 20 6f 72 20 | 31 33 2e 20 20 54 68 69 | 11, or |13. Thi|
|00001d00| 73 20 69 73 20 61 20 76 | 65 72 79 20 65 61 73 79 |s is a v|ery easy|
|00001d10| 20 63 68 65 63 6b 2e 0a | 58 09 20 2a 2f 0a 58 09 | check..|X. */.X.|
|00001d20| 69 6a 20 3d 20 7a 6d 6f | 64 69 28 7a 2c 20 31 35 |ij = zmo|di(z, 15|
|00001d30| 30 31 35 4c 29 3b 0a 58 | 09 69 66 20 28 21 28 69 |015L);.X|.if (!(i|
|00001d40| 6a 20 25 20 33 29 20 7c | 7c 20 21 28 69 6a 20 25 |j % 3) ||| !(ij %|
|00001d50| 20 35 29 20 7c 7c 20 21 | 28 69 6a 20 25 20 37 29 | 5) || !|(ij % 7)|
|00001d60| 20 7c 7c 20 21 28 69 6a | 20 25 20 31 31 29 20 7c | || !(ij| % 11) ||
|00001d70| 7c 20 21 28 69 6a 20 25 | 20 31 33 29 29 0a 58 09 || !(ij %| 13)).X.|
|00001d80| 09 72 65 74 75 72 6e 20 | 46 41 4c 53 45 3b 0a 58 |.return |FALSE;.X|
|00001d90| 09 2f 2a 0a 58 09 20 2a | 20 43 68 65 63 6b 20 74 |./*.X. *| Check t|
|00001da0| 68 65 20 67 63 64 20 6f | 66 20 74 68 65 20 6e 75 |he gcd o|f the nu|
|00001db0| 6d 62 65 72 20 61 6e 64 | 20 74 68 65 20 70 72 6f |mber and| the pro|
|00001dc0| 64 75 63 74 20 6f 66 20 | 6d 6f 72 65 20 6f 66 20 |duct of |more of |
|00001dd0| 74 68 65 20 66 69 72 73 | 74 0a 58 09 20 2a 20 66 |the firs|t.X. * f|
|00001de0| 65 77 20 6f 64 64 20 70 | 72 69 6d 65 73 2e 20 20 |ew odd p|rimes. |
|00001df0| 57 65 20 6d 75 73 74 20 | 62 75 69 6c 64 20 74 68 |We must |build th|
|00001e00| 65 20 70 72 69 6d 65 20 | 70 72 6f 64 75 63 74 20 |e prime |product |
|00001e10| 6f 6e 20 74 68 65 20 66 | 69 72 73 74 20 63 61 6c |on the f|irst cal|
|00001e20| 6c 2e 0a 58 09 20 2a 2f | 0a 58 09 7a 74 6d 70 2e |l..X. */|.X.ztmp.|
|00001e30| 73 69 67 6e 20 3d 20 30 | 3b 0a 58 09 7a 74 6d 70 |sign = 0|;.X.ztmp|
|00001e40| 2e 6c 65 6e 20 3d 20 31 | 3b 0a 58 09 7a 74 6d 70 |.len = 1|;.X.ztmp|
|00001e50| 2e 76 20 3d 20 76 61 6c | 3b 0a 58 09 69 66 20 28 |.v = val|;.X.if (|
|00001e60| 70 72 69 6d 65 70 72 6f | 64 2e 6c 65 6e 20 3d 3d |primepro|d.len ==|
|00001e70| 20 30 29 20 7b 0a 58 09 | 09 76 61 6c 5b 30 5d 20 | 0) {.X.|.val[0] |
|00001e80| 3d 20 31 30 31 3b 0a 58 | 09 09 7a 70 66 61 63 74 |= 101;.X|..zpfact|
|00001e90| 28 7a 74 6d 70 2c 20 26 | 70 72 69 6d 65 70 72 6f |(ztmp, &|primepro|
|00001ea0| 64 29 3b 0a 58 09 7d 0a | 58 09 7a 67 63 64 28 7a |d);.X.}.|X.zgcd(z|
|00001eb0| 2c 20 70 72 69 6d 65 70 | 72 6f 64 2c 20 26 7a 31 |, primep|rod, &z1|
|00001ec0| 29 3b 0a 58 09 69 66 20 | 28 21 69 73 75 6e 69 74 |);.X.if |(!isunit|
|00001ed0| 28 7a 31 29 29 20 7b 0a | 58 09 09 66 72 65 65 68 |(z1)) {.|X..freeh|
|00001ee0| 28 7a 31 2e 76 29 3b 0a | 58 09 09 72 65 74 75 72 |(z1.v);.|X..retur|
|00001ef0| 6e 20 46 41 4c 53 45 3b | 0a 58 09 7d 0a 58 09 66 |n FALSE;|.X.}.X.f|
|00001f00| 72 65 65 68 28 7a 31 2e | 76 29 3b 0a 58 09 2f 2a |reeh(z1.|v);.X./*|
|00001f10| 0a 58 09 20 2a 20 4e 6f | 74 20 64 69 76 69 73 69 |.X. * No|t divisi|
|00001f20| 62 6c 65 20 62 79 20 61 | 20 73 6d 61 6c 6c 20 70 |ble by a| small p|
|00001f30| 72 69 6d 65 2c 20 73 6f | 20 6f 6e 77 61 72 64 20 |rime, so| onward |
|00001f40| 77 69 74 68 20 74 68 65 | 20 72 65 61 6c 20 74 65 |with the| real te|
|00001f50| 73 74 2e 0a 58 09 20 2a | 20 4d 61 6b 65 20 73 75 |st..X. *| Make su|
|00001f60| 72 65 20 74 68 65 20 63 | 6f 75 6e 74 20 69 73 20 |re the c|ount is |
|00001f70| 6c 69 6d 69 74 65 64 20 | 62 79 20 74 68 65 20 6e |limited |by the n|
|00001f80| 75 6d 62 65 72 20 6f 66 | 20 6f 64 64 20 6e 75 6d |umber of| odd num|
|00001f90| 62 65 72 73 20 62 65 74 | 77 65 65 6e 0a 58 09 20 |bers bet|ween.X. |
|00001fa0| 2a 20 74 68 72 65 65 20 | 61 6e 64 20 74 68 65 20 |* three |and the |
|00001fb0| 6e 75 6d 62 65 72 20 62 | 65 69 6e 67 20 74 65 73 |number b|eing tes|
|00001fc0| 74 65 64 2e 0a 58 09 20 | 2a 2f 0a 58 09 69 78 20 |ted..X. |*/.X.ix |
|00001fd0| 3d 20 28 28 69 73 74 69 | 6e 79 28 7a 29 20 3f 20 |= ((isti|ny(z) ? |
|00001fe0| 7a 31 74 6f 6c 28 7a 29 | 20 3a 20 7a 32 74 6f 6c |z1tol(z)| : z2tol|
|00001ff0| 28 7a 29 20 2d 20 33 29 | 20 2f 20 32 29 3b 0a 58 |(z) - 3)| / 2);.X|
|00002000| 09 69 66 20 28 63 6f 75 | 6e 74 20 3e 20 69 78 29 |.if (cou|nt > ix)|
|00002010| 20 63 6f 75 6e 74 20 3d | 20 69 78 3b 0a 58 09 7a | count =| ix;.X.z|
|00002020| 73 75 62 28 7a 2c 20 5f | 6f 6e 65 5f 2c 20 26 7a |sub(z, _|one_, &z|
|00002030| 6d 31 29 3b 0a 58 09 69 | 6b 20 3d 20 7a 6c 6f 77 |m1);.X.i|k = zlow|
|00002040| 62 69 74 28 7a 6d 31 29 | 3b 0a 58 09 7a 73 68 69 |bit(zm1)|;.X.zshi|
|00002050| 66 74 28 7a 6d 31 2c 20 | 2d 69 6b 2c 20 26 7a 31 |ft(zm1, |-ik, &z1|
|00002060| 29 3b 0a 58 09 2f 2a 0a | 58 09 20 2a 20 4c 6f 6f |);.X./*.|X. * Loo|
|00002070| 70 20 6f 76 65 72 20 76 | 61 72 69 6f 75 73 20 22 |p over v|arious "|
|00002080| 72 61 6e 64 6f 6d 22 20 | 6e 75 6d 62 65 72 73 2c |random" |numbers,|
|00002090| 20 74 65 73 74 69 6e 67 | 20 65 61 63 68 20 6f 6e | testing| each on|
|000020a0| 65 2e 0a 58 09 20 2a 20 | 54 68 65 73 65 20 6e 75 |e..X. * |These nu|
|000020b0| 6d 62 65 72 73 20 61 72 | 65 20 74 68 65 20 6f 64 |mbers ar|e the od|
|000020c0| 64 20 6e 75 6d 62 65 72 | 73 20 73 74 61 72 74 69 |d number|s starti|
|000020d0| 6e 67 20 66 72 6f 6d 20 | 74 68 72 65 65 2e 0a 58 |ng from |three..X|
|000020e0| 09 20 2a 2f 0a 58 09 66 | 6f 72 20 28 69 78 20 3d |. */.X.f|or (ix =|
|000020f0| 20 30 3b 20 69 78 20 3c | 20 63 6f 75 6e 74 3b 20 | 0; ix <| count; |
|00002100| 69 78 2b 2b 29 20 7b 0a | 58 09 09 76 61 6c 5b 30 |ix++) {.|X..val[0|
|00002110| 5d 20 3d 20 28 69 78 20 | 2a 20 32 29 20 2b 20 33 |] = (ix |* 2) + 3|
|00002120| 3b 0a 58 09 09 69 6a 20 | 3d 20 30 3b 0a 58 09 09 |;.X..ij |= 0;.X..|
|00002130| 7a 70 6f 77 65 72 6d 6f | 64 28 7a 74 6d 70 2c 20 |zpowermo|d(ztmp, |
|00002140| 7a 31 2c 20 7a 2c 20 26 | 7a 33 29 3b 0a 58 09 09 |z1, z, &|z3);.X..|
|00002150| 66 6f 72 20 28 3b 3b 29 | 20 7b 0a 58 09 09 09 69 |for (;;)| {.X...i|
|00002160| 66 20 28 69 73 6f 6e 65 | 28 7a 33 29 29 20 7b 0a |f (isone|(z3)) {.|
|00002170| 58 09 09 09 09 69 66 20 | 28 69 6a 29 09 2f 2a 20 |X....if |(ij)./* |
|00002180| 6e 75 6d 62 65 72 20 69 | 73 20 64 65 66 69 6e 69 |number i|s defini|
|00002190| 74 65 6c 79 20 6e 6f 74 | 20 70 72 69 6d 65 20 2a |tely not| prime *|
|000021a0| 2f 0a 58 09 09 09 09 09 | 67 6f 74 6f 20 6e 6f 74 |/.X.....|goto not|
|000021b0| 70 72 69 6d 65 3b 0a 58 | 09 09 09 09 62 72 65 61 |prime;.X|....brea|
|000021c0| 6b 3b 0a 58 09 09 09 7d | 0a 58 09 09 09 69 66 20 |k;.X...}|.X...if |
|000021d0| 28 7a 63 6d 70 28 7a 33 | 2c 20 7a 6d 31 29 20 3d |(zcmp(z3|, zm1) =|
|000021e0| 3d 20 30 29 0a 58 09 09 | 09 09 62 72 65 61 6b 3b |= 0).X..|..break;|
|000021f0| 0a 58 09 09 09 69 66 20 | 28 2b 2b 69 6a 20 3e 3d |.X...if |(++ij >=|
|00002200| 20 69 6b 29 0a 58 09 09 | 09 09 67 6f 74 6f 20 6e | ik).X..|..goto n|
|00002210| 6f 74 70 72 69 6d 65 3b | 09 2f 2a 20 6e 75 6d 62 |otprime;|./* numb|
|00002220| 65 72 20 69 73 20 64 65 | 66 69 6e 69 74 65 6c 79 |er is de|finitely|
|00002230| 20 6e 6f 74 20 70 72 69 | 6d 65 20 2a 2f 0a 58 09 | not pri|me */.X.|
|00002240| 09 09 7a 73 71 75 61 72 | 65 28 7a 33 2c 20 26 7a |..zsquar|e(z3, &z|
|00002250| 32 29 3b 0a 58 09 09 09 | 66 72 65 65 68 28 7a 33 |2);.X...|freeh(z3|
|00002260| 2e 76 29 3b 0a 58 09 09 | 09 7a 6d 6f 64 28 7a 32 |.v);.X..|.zmod(z2|
|00002270| 2c 20 7a 2c 20 26 7a 33 | 29 3b 0a 58 09 09 09 66 |, z, &z3|);.X...f|
|00002280| 72 65 65 68 28 7a 32 2e | 76 29 3b 0a 58 09 09 7d |reeh(z2.|v);.X..}|
|00002290| 0a 58 09 09 66 72 65 65 | 68 28 7a 33 2e 76 29 3b |.X..free|h(z3.v);|
|000022a0| 0a 58 09 7d 0a 58 09 66 | 72 65 65 68 28 7a 6d 31 |.X.}.X.f|reeh(zm1|
|000022b0| 2e 76 29 3b 0a 58 09 66 | 72 65 65 68 28 7a 31 2e |.v);.X.f|reeh(z1.|
|000022c0| 76 29 3b 0a 58 09 72 65 | 74 75 72 6e 20 54 52 55 |v);.X.re|turn TRU|
|000022d0| 45 3b 09 2f 2a 20 6e 75 | 6d 62 65 72 20 6d 69 67 |E;./* nu|mber mig|
|000022e0| 68 74 20 62 65 20 70 72 | 69 6d 65 20 2a 2f 0a 58 |ht be pr|ime */.X|
|000022f0| 0a 58 6e 6f 74 70 72 69 | 6d 65 3a 0a 58 09 66 72 |.Xnotpri|me:.X.fr|
|00002300| 65 65 68 28 7a 33 2e 76 | 29 3b 0a 58 09 66 72 65 |eeh(z3.v|);.X.fre|
|00002310| 65 68 28 7a 6d 31 2e 76 | 29 3b 0a 58 09 66 72 65 |eh(zm1.v|);.X.fre|
|00002320| 65 68 28 7a 31 2e 76 29 | 3b 0a 58 09 72 65 74 75 |eh(z1.v)|;.X.retu|
|00002330| 72 6e 20 46 41 4c 53 45 | 3b 0a 58 7d 0a 58 0a 58 |rn FALSE|;.X}.X.X|
|00002340| 0a 58 2f 2a 0a 58 20 2a | 20 43 6f 6d 70 75 74 65 |.X/*.X *| Compute|
|00002350| 20 74 68 65 20 4a 61 63 | 6f 62 69 20 66 75 6e 63 | the Jac|obi func|
|00002360| 74 69 6f 6e 20 28 70 20 | 2f 20 71 29 20 66 6f 72 |tion (p |/ q) for|
|00002370| 20 6f 64 64 20 71 2e 0a | 58 20 2a 20 49 66 20 71 | odd q..|X * If q|
|00002380| 20 69 73 20 70 72 69 6d | 65 20 74 68 65 6e 20 74 | is prim|e then t|
|00002390| 68 65 20 72 65 73 75 6c | 74 20 69 73 3a 0a 58 20 |he resul|t is:.X |
|000023a0| 2a 09 31 20 69 66 20 70 | 20 3d 3d 20 78 5e 32 20 |*.1 if p| == x^2 |
|000023b0| 28 6d 6f 64 20 71 29 20 | 66 6f 72 20 73 6f 6d 65 |(mod q) |for some|
|000023c0| 20 78 2e 0a 58 20 2a 09 | 2d 31 20 6f 74 68 65 72 | x..X *.|-1 other|
|000023d0| 77 69 73 65 2e 0a 58 20 | 2a 20 49 66 20 71 20 69 |wise..X |* If q i|
|000023e0| 73 20 6e 6f 74 20 70 72 | 69 6d 65 2c 20 74 68 65 |s not pr|ime, the|
|000023f0| 6e 20 74 68 65 20 72 65 | 73 75 6c 74 20 69 73 20 |n the re|sult is |
|00002400| 6e 6f 74 20 6d 65 61 6e | 69 6e 67 66 75 6c 20 69 |not mean|ingful i|
|00002410| 66 20 69 74 20 69 73 20 | 31 2e 0a 58 20 2a 20 54 |f it is |1..X * T|
|00002420| 68 69 73 20 66 75 6e 63 | 74 69 6f 6e 20 72 65 74 |his func|tion ret|
|00002430| 75 72 6e 73 20 30 20 69 | 66 20 71 20 69 73 20 65 |urns 0 i|f q is e|
|00002440| 76 65 6e 20 6f 72 20 71 | 20 3c 20 30 2e 0a 58 20 |ven or q| < 0..X |
|00002450| 2a 2f 0a 58 46 4c 41 47 | 0a 58 7a 6a 61 63 6f 62 |*/.XFLAG|.Xzjacob|
|00002460| 69 28 7a 31 2c 20 7a 32 | 29 0a 58 09 5a 56 41 4c |i(z1, z2|).X.ZVAL|
|00002470| 55 45 20 7a 31 2c 20 7a | 32 3b 0a 58 7b 0a 58 09 |UE z1, z|2;.X{.X.|
|00002480| 5a 56 41 4c 55 45 20 70 | 2c 20 71 2c 20 74 6d 70 |ZVALUE p|, q, tmp|
|00002490| 3b 0a 58 09 6c 6f 6e 67 | 20 6c 6f 77 62 69 74 3b |;.X.long| lowbit;|
|000024a0| 0a 58 09 69 6e 74 20 76 | 61 6c 3b 0a 58 0a 58 09 |.X.int v|al;.X.X.|
|000024b0| 69 66 20 28 69 73 65 76 | 65 6e 28 7a 32 29 20 7c |if (isev|en(z2) ||
|000024c0| 7c 20 69 73 6e 65 67 28 | 7a 32 29 29 0a 58 09 09 || isneg(|z2)).X..|
|000024d0| 72 65 74 75 72 6e 20 30 | 3b 0a 58 09 76 61 6c 20 |return 0|;.X.val |
|000024e0| 3d 20 31 3b 0a 58 09 69 | 66 20 28 69 73 7a 65 72 |= 1;.X.i|f (iszer|
|000024f0| 6f 28 7a 31 29 20 7c 7c | 20 69 73 6f 6e 65 28 7a |o(z1) ||| isone(z|
|00002500| 31 29 29 0a 58 09 09 72 | 65 74 75 72 6e 20 76 61 |1)).X..r|eturn va|
|00002510| 6c 3b 0a 58 09 69 66 20 | 28 69 73 75 6e 69 74 28 |l;.X.if |(isunit(|
|00002520| 7a 31 29 29 20 7b 0a 58 | 09 09 69 66 20 28 28 2a |z1)) {.X|..if ((*|
|00002530| 7a 32 2e 76 20 2d 20 31 | 29 20 26 20 30 78 32 29 |z2.v - 1|) & 0x2)|
|00002540| 0a 58 09 09 09 76 61 6c | 20 3d 20 2d 76 61 6c 3b |.X...val| = -val;|
|00002550| 0a 58 09 09 72 65 74 75 | 72 6e 20 76 61 6c 3b 0a |.X..retu|rn val;.|
|00002560| 58 09 7d 0a 58 09 7a 63 | 6f 70 79 28 7a 31 2c 20 |X.}.X.zc|opy(z1, |
|00002570| 26 70 29 3b 0a 58 09 7a | 63 6f 70 79 28 7a 32 2c |&p);.X.z|copy(z2,|
|00002580| 20 26 71 29 3b 0a 58 09 | 66 6f 72 20 28 3b 3b 29 | &q);.X.|for (;;)|
|00002590| 20 7b 0a 58 09 09 7a 6d | 6f 64 28 70 2c 20 71 2c | {.X..zm|od(p, q,|
|000025a0| 20 26 74 6d 70 29 3b 0a | 58 09 09 66 72 65 65 68 | &tmp);.|X..freeh|
|000025b0| 28 70 2e 76 29 3b 0a 58 | 09 09 70 20 3d 20 74 6d |(p.v);.X|..p = tm|
|000025c0| 70 3b 0a 58 09 09 69 66 | 20 28 69 73 7a 65 72 6f |p;.X..if| (iszero|
|000025d0| 28 70 29 29 20 7b 0a 58 | 09 09 09 66 72 65 65 68 |(p)) {.X|...freeh|
|000025e0| 28 70 2e 76 29 3b 0a 58 | 09 09 09 70 20 3d 20 5f |(p.v);.X|...p = _|
|000025f0| 6f 6e 65 5f 3b 0a 58 09 | 09 7d 0a 58 09 09 69 66 |one_;.X.|.}.X..if|
|00002600| 20 28 69 73 65 76 65 6e | 28 70 29 29 20 7b 0a 58 | (iseven|(p)) {.X|
|00002610| 09 09 09 6c 6f 77 62 69 | 74 20 3d 20 7a 6c 6f 77 |...lowbi|t = zlow|
|00002620| 62 69 74 28 70 29 3b 0a | 58 09 09 09 7a 73 68 69 |bit(p);.|X...zshi|
|00002630| 66 74 28 70 2c 20 2d 6c | 6f 77 62 69 74 2c 20 26 |ft(p, -l|owbit, &|
|00002640| 74 6d 70 29 3b 0a 58 09 | 09 09 66 72 65 65 68 28 |tmp);.X.|..freeh(|
|00002650| 70 2e 76 29 3b 0a 58 09 | 09 09 70 20 3d 20 74 6d |p.v);.X.|..p = tm|
|00002660| 70 3b 0a 58 09 09 09 69 | 66 20 28 28 6c 6f 77 62 |p;.X...i|f ((lowb|
|00002670| 69 74 20 26 20 31 29 20 | 26 26 20 28 28 28 2a 71 |it & 1) |&& (((*q|
|00002680| 2e 76 20 26 20 30 78 37 | 29 20 3d 3d 20 33 29 20 |.v & 0x7|) == 3) |
|00002690| 7c 7c 20 28 28 2a 71 2e | 76 20 26 20 30 78 37 29 ||| ((*q.|v & 0x7)|
|000026a0| 20 3d 3d 20 35 29 29 29 | 0a 58 09 09 09 09 76 61 | == 5)))|.X....va|
|000026b0| 6c 20 3d 20 2d 76 61 6c | 3b 0a 58 09 09 7d 0a 58 |l = -val|;.X..}.X|
|000026c0| 09 09 69 66 20 28 69 73 | 75 6e 69 74 28 70 29 29 |..if (is|unit(p))|
|000026d0| 20 7b 0a 58 09 09 09 66 | 72 65 65 68 28 70 2e 76 | {.X...f|reeh(p.v|
|000026e0| 29 3b 0a 58 09 09 09 66 | 72 65 65 68 28 71 2e 76 |);.X...f|reeh(q.v|
|000026f0| 29 3b 0a 58 09 09 09 72 | 65 74 75 72 6e 20 76 61 |);.X...r|eturn va|
|00002700| 6c 3b 0a 58 09 09 7d 0a | 58 09 09 69 66 20 28 28 |l;.X..}.|X..if ((|
|00002710| 2a 70 2e 76 20 26 20 2a | 71 2e 76 20 26 20 30 78 |*p.v & *|q.v & 0x|
|00002720| 33 29 20 3d 3d 20 33 29 | 0a 58 09 09 09 76 61 6c |3) == 3)|.X...val|
|00002730| 20 3d 20 2d 76 61 6c 3b | 0a 58 09 09 74 6d 70 20 | = -val;|.X..tmp |
|00002740| 3d 20 71 3b 0a 58 09 09 | 71 20 3d 20 70 3b 0a 58 |= q;.X..|q = p;.X|
|00002750| 09 09 70 20 3d 20 74 6d | 70 3b 0a 58 09 7d 0a 58 |..p = tm|p;.X.}.X|
|00002760| 7d 0a 58 0a 58 0a 58 2f | 2a 0a 58 20 2a 20 52 65 |}.X.X.X/|*.X * Re|
|00002770| 74 75 72 6e 20 74 68 65 | 20 46 69 62 6f 6e 61 63 |turn the| Fibonac|
|00002780| 63 69 20 6e 75 6d 62 65 | 72 20 46 28 6e 29 2e 0a |ci numbe|r F(n)..|
|00002790| 58 20 2a 20 54 68 69 73 | 20 69 73 20 65 76 61 6c |X * This| is eval|
|000027a0| 75 61 74 65 64 20 62 79 | 20 72 65 63 75 72 73 69 |uated by| recursi|
|000027b0| 76 65 6c 79 20 75 73 69 | 6e 67 20 74 68 65 20 66 |vely usi|ng the f|
|000027c0| 6f 72 6d 75 6c 61 73 3a | 0a 58 20 2a 09 46 28 32 |ormulas:|.X *.F(2|
|000027d0| 4e 2b 31 29 20 3d 20 46 | 28 4e 2b 31 29 5e 32 20 |N+1) = F|(N+1)^2 |
|000027e0| 2b 20 46 28 4e 29 5e 32 | 0a 58 20 2a 20 61 6e 64 |+ F(N)^2|.X * and|
|000027f0| 0a 58 20 2a 09 46 28 32 | 4e 29 20 3d 20 46 28 4e |.X *.F(2|N) = F(N|
|00002800| 2b 31 29 5e 32 20 2d 20 | 46 28 4e 2d 31 29 5e 32 |+1)^2 - |F(N-1)^2|
|00002810| 0a 58 20 2a 2f 0a 58 76 | 6f 69 64 0a 58 7a 66 69 |.X */.Xv|oid.Xzfi|
|00002820| 62 28 7a 2c 20 72 65 73 | 29 0a 58 09 5a 56 41 4c |b(z, res|).X.ZVAL|
|00002830| 55 45 20 7a 2c 20 2a 72 | 65 73 3b 0a 58 7b 0a 58 |UE z, *r|es;.X{.X|
|00002840| 09 75 6e 73 69 67 6e 65 | 64 20 6c 6f 6e 67 20 69 |.unsigne|d long i|
|00002850| 3b 0a 58 09 6c 6f 6e 67 | 20 6e 3b 0a 58 09 69 6e |;.X.long| n;.X.in|
|00002860| 74 20 73 69 67 6e 3b 0a | 58 09 5a 56 41 4c 55 45 |t sign;.|X.ZVALUE|
|00002870| 20 66 6e 6d 31 2c 20 66 | 6e 2c 20 66 6e 70 31 3b | fnm1, f|n, fnp1;|
|00002880| 09 09 2f 2a 20 63 6f 6e | 73 65 63 75 74 69 76 65 |../* con|secutive|
|00002890| 20 66 69 62 6f 6e 61 63 | 63 69 20 76 61 6c 75 65 | fibonac|ci value|
|000028a0| 73 20 2a 2f 0a 58 09 5a | 56 41 4c 55 45 20 74 31 |s */.X.Z|VALUE t1|
|000028b0| 2c 20 74 32 2c 20 74 33 | 3b 0a 58 0a 58 09 69 66 |, t2, t3|;.X.X.if|
|000028c0| 20 28 69 73 62 69 67 28 | 7a 29 29 0a 58 09 09 65 | (isbig(|z)).X..e|
|000028d0| 72 72 6f 72 28 22 56 65 | 72 79 20 6c 61 72 67 65 |rror("Ve|ry large|
|000028e0| 20 46 69 62 6f 6e 61 63 | 63 69 20 6e 75 6d 62 65 | Fibonac|ci numbe|
|000028f0| 72 22 29 3b 0a 58 09 6e | 20 3d 20 28 69 73 74 69 |r");.X.n| = (isti|
|00002900| 6e 79 28 7a 29 20 3f 20 | 7a 31 74 6f 6c 28 7a 29 |ny(z) ? |z1tol(z)|
|00002910| 20 3a 20 7a 32 74 6f 6c | 28 7a 29 29 3b 0a 58 09 | : z2tol|(z));.X.|
|00002920| 69 66 20 28 6e 20 3d 3d | 20 30 29 20 7b 0a 58 09 |if (n ==| 0) {.X.|
|00002930| 09 2a 72 65 73 20 3d 20 | 5f 7a 65 72 6f 5f 3b 0a |.*res = |_zero_;.|
|00002940| 58 09 09 72 65 74 75 72 | 6e 3b 0a 58 09 7d 0a 58 |X..retur|n;.X.}.X|
|00002950| 09 73 69 67 6e 20 3d 20 | 7a 2e 73 69 67 6e 20 26 |.sign = |z.sign &|
|00002960| 26 20 28 28 6e 20 26 20 | 30 78 31 29 20 3d 3d 20 |& ((n & |0x1) == |
|00002970| 30 29 3b 0a 58 09 69 66 | 20 28 6e 20 3c 3d 20 32 |0);.X.if| (n <= 2|
|00002980| 29 20 7b 0a 58 09 09 2a | 72 65 73 20 3d 20 5f 6f |) {.X..*|res = _o|
|00002990| 6e 65 5f 3b 0a 58 09 09 | 72 65 73 2d 3e 73 69 67 |ne_;.X..|res->sig|
|000029a0| 6e 20 3d 20 28 42 4f 4f | 4c 29 73 69 67 6e 3b 0a |n = (BOO|L)sign;.|
|000029b0| 58 09 09 72 65 74 75 72 | 6e 3b 0a 58 09 7d 0a 58 |X..retur|n;.X.}.X|
|000029c0| 09 69 20 3d 20 54 4f 50 | 46 55 4c 4c 3b 0a 58 09 |.i = TOP|FULL;.X.|
|000029d0| 77 68 69 6c 65 20 28 28 | 69 20 26 20 6e 29 20 3d |while ((|i & n) =|
|000029e0| 3d 20 30 29 0a 58 09 09 | 69 20 3e 3e 3d 20 31 4c |= 0).X..|i >>= 1L|
|000029f0| 3b 0a 58 09 69 20 3e 3e | 3d 20 31 4c 3b 0a 58 09 |;.X.i >>|= 1L;.X.|
|00002a00| 66 6e 6d 31 20 3d 20 5f | 7a 65 72 6f 5f 3b 0a 58 |fnm1 = _|zero_;.X|
|00002a10| 09 66 6e 20 3d 20 5f 6f | 6e 65 5f 3b 0a 58 09 66 |.fn = _o|ne_;.X.f|
|00002a20| 6e 70 31 20 3d 20 5f 6f | 6e 65 5f 3b 0a 58 09 77 |np1 = _o|ne_;.X.w|
|00002a30| 68 69 6c 65 20 28 69 29 | 20 7b 0a 58 09 09 7a 73 |hile (i)| {.X..zs|
|00002a40| 71 75 61 72 65 28 66 6e | 6d 31 2c 20 26 74 31 29 |quare(fn|m1, &t1)|
|00002a50| 3b 0a 58 09 09 7a 73 71 | 75 61 72 65 28 66 6e 2c |;.X..zsq|uare(fn,|
|00002a60| 20 26 74 32 29 3b 0a 58 | 09 09 7a 73 71 75 61 72 | &t2);.X|..zsquar|
|00002a70| 65 28 66 6e 70 31 2c 20 | 26 74 33 29 3b 0a 58 09 |e(fnp1, |&t3);.X.|
|00002a80| 09 66 72 65 65 68 28 66 | 6e 6d 31 2e 76 29 3b 0a |.freeh(f|nm1.v);.|
|00002a90| 58 09 09 66 72 65 65 68 | 28 66 6e 2e 76 29 3b 0a |X..freeh|(fn.v);.|
|00002aa0| 58 09 09 66 72 65 65 68 | 28 66 6e 70 31 2e 76 29 |X..freeh|(fnp1.v)|
|00002ab0| 3b 0a 58 09 09 7a 61 64 | 64 28 74 32 2c 20 74 33 |;.X..zad|d(t2, t3|
|00002ac0| 2c 20 26 66 6e 70 31 29 | 3b 0a 58 09 09 7a 73 75 |, &fnp1)|;.X..zsu|
|00002ad0| 62 28 74 33 2c 20 74 31 | 2c 20 26 66 6e 29 3b 0a |b(t3, t1|, &fn);.|
|00002ae0| 58 09 09 66 72 65 65 68 | 28 74 31 2e 76 29 3b 0a |X..freeh|(t1.v);.|
|00002af0| 58 09 09 66 72 65 65 68 | 28 74 32 2e 76 29 3b 0a |X..freeh|(t2.v);.|
|00002b00| 58 09 09 66 72 65 65 68 | 28 74 33 2e 76 29 3b 0a |X..freeh|(t3.v);.|
|00002b10| 58 09 09 69 66 20 28 69 | 20 26 20 6e 29 20 7b 0a |X..if (i| & n) {.|
|00002b20| 58 09 09 09 66 6e 6d 31 | 20 3d 20 66 6e 3b 0a 58 |X...fnm1| = fn;.X|
|00002b30| 09 09 09 66 6e 20 3d 20 | 66 6e 70 31 3b 0a 58 09 |...fn = |fnp1;.X.|
|00002b40| 09 09 7a 61 64 64 28 66 | 6e 6d 31 2c 20 66 6e 2c |..zadd(f|nm1, fn,|
|00002b50| 20 26 66 6e 70 31 29 3b | 0a 58 09 09 7d 20 65 6c | &fnp1);|.X..} el|
|00002b60| 73 65 0a 58 09 09 09 7a | 73 75 62 28 66 6e 70 31 |se.X...z|sub(fnp1|
|00002b70| 2c 20 66 6e 2c 20 26 66 | 6e 6d 31 29 3b 0a 58 09 |, fn, &f|nm1);.X.|
|00002b80| 09 69 20 3e 3e 3d 20 31 | 4c 3b 0a 58 09 7d 0a 58 |.i >>= 1|L;.X.}.X|
|00002b90| 09 66 72 65 65 68 28 66 | 6e 6d 31 2e 76 29 3b 0a |.freeh(f|nm1.v);.|
|00002ba0| 58 09 66 72 65 65 68 28 | 66 6e 70 31 2e 76 29 3b |X.freeh(|fnp1.v);|
|00002bb0| 0a 58 09 2a 72 65 73 20 | 3d 20 66 6e 3b 0a 58 09 |.X.*res |= fn;.X.|
|00002bc0| 72 65 73 2d 3e 73 69 67 | 6e 20 3d 20 28 42 4f 4f |res->sig|n = (BOO|
|00002bd0| 4c 29 73 69 67 6e 3b 0a | 58 7d 0a 58 0a 58 0a 58 |L)sign;.|X}.X.X.X|
|00002be0| 2f 2a 0a 58 20 2a 20 43 | 6f 6d 70 75 74 65 20 74 |/*.X * C|ompute t|
|00002bf0| 68 65 20 72 65 73 75 6c | 74 20 6f 66 20 72 61 69 |he resul|t of rai|
|00002c00| 73 69 6e 67 20 6f 6e 65 | 20 6e 75 6d 62 65 72 20 |sing one| number |
|00002c10| 74 6f 20 74 68 65 20 70 | 6f 77 65 72 20 6f 66 20 |to the p|ower of |
|00002c20| 61 6e 6f 74 68 65 72 0a | 58 20 2a 20 54 68 65 20 |another.|X * The |
|00002c30| 73 65 63 6f 6e 64 20 6e | 75 6d 62 65 72 20 69 73 |second n|umber is|
|00002c40| 20 61 73 73 75 6d 65 64 | 20 74 6f 20 62 65 20 6e | assumed| to be n|
|00002c50| 6f 6e 2d 6e 65 67 61 74 | 69 76 65 2e 0a 58 20 2a |on-negat|ive..X *|
|00002c60| 20 49 74 20 63 61 6e 6e | 6f 74 20 62 65 20 74 6f | It cann|ot be to|
|00002c70| 6f 20 6c 61 72 67 65 20 | 65 78 63 65 70 74 20 66 |o large |except f|
|00002c80| 6f 72 20 74 72 69 76 69 | 61 6c 20 63 61 73 65 73 |or trivi|al cases|
|00002c90| 2e 0a 58 20 2a 2f 0a 58 | 76 6f 69 64 0a 58 7a 70 |..X */.X|void.Xzp|
|00002ca0| 6f 77 69 28 7a 31 2c 20 | 7a 32 2c 20 72 65 73 29 |owi(z1, |z2, res)|
|00002cb0| 0a 58 09 5a 56 41 4c 55 | 45 20 7a 31 2c 20 7a 32 |.X.ZVALU|E z1, z2|
|00002cc0| 2c 20 2a 72 65 73 3b 0a | 58 7b 0a 58 09 69 6e 74 |, *res;.|X{.X.int|
|00002cd0| 20 73 69 67 6e 3b 09 09 | 2f 2a 20 66 69 6e 61 6c | sign;..|/* final|
|00002ce0| 20 73 69 67 6e 20 6f 66 | 20 6e 75 6d 62 65 72 20 | sign of| number |
|00002cf0| 2a 2f 0a 58 09 75 6e 73 | 69 67 6e 65 64 20 6c 6f |*/.X.uns|igned lo|
|00002d00| 6e 67 20 70 6f 77 65 72 | 3b 09 2f 2a 20 70 6f 77 |ng power|;./* pow|
|00002d10| 65 72 20 74 6f 20 72 61 | 69 73 65 20 74 6f 20 2a |er to ra|ise to *|
|00002d20| 2f 0a 58 09 75 6e 73 69 | 67 6e 65 64 20 6c 6f 6e |/.X.unsi|gned lon|
|00002d30| 67 20 62 69 74 3b 09 2f | 2a 20 63 75 72 72 65 6e |g bit;./|* curren|
|00002d40| 74 20 62 69 74 20 76 61 | 6c 75 65 20 2a 2f 0a 58 |t bit va|lue */.X|
|00002d50| 09 6c 6f 6e 67 20 74 77 | 6f 73 3b 09 09 2f 2a 20 |.long tw|os;../* |
|00002d60| 63 6f 75 6e 74 20 6f 66 | 20 74 69 6d 65 73 20 32 |count of| times 2|
|00002d70| 20 69 73 20 69 6e 20 72 | 65 73 75 6c 74 20 2a 2f | is in r|esult */|
|00002d80| 0a 58 09 5a 56 41 4c 55 | 45 20 61 6e 73 2c 20 74 |.X.ZVALU|E ans, t|
|00002d90| 65 6d 70 3b 0a 58 0a 58 | 09 73 69 67 6e 20 3d 20 |emp;.X.X|.sign = |
|00002da0| 28 7a 31 2e 73 69 67 6e | 20 26 26 20 69 73 6f 64 |(z1.sign| && isod|
|00002db0| 64 28 7a 32 29 29 3b 0a | 58 09 7a 31 2e 73 69 67 |d(z2));.|X.z1.sig|
|00002dc0| 6e 20 3d 20 30 3b 0a 58 | 09 7a 32 2e 73 69 67 6e |n = 0;.X|.z2.sign|
|00002dd0| 20 3d 20 30 3b 0a 58 09 | 69 66 20 28 69 73 7a 65 | = 0;.X.|if (isze|
|00002de0| 72 6f 28 7a 32 29 29 20 | 7b 09 2f 2a 20 6e 75 6d |ro(z2)) |{./* num|
|00002df0| 62 65 72 20 72 61 69 73 | 65 64 20 74 6f 20 70 6f |ber rais|ed to po|
|00002e00| 77 65 72 20 30 20 2a 2f | 0a 58 09 09 69 66 20 28 |wer 0 */|.X..if (|
|00002e10| 69 73 7a 65 72 6f 28 7a | 31 29 29 0a 58 09 09 09 |iszero(z|1)).X...|
|00002e20| 65 72 72 6f 72 28 22 5a | 65 72 6f 20 72 61 69 73 |error("Z|ero rais|
|00002e30| 65 64 20 74 6f 20 7a 65 | 72 6f 20 70 6f 77 65 72 |ed to ze|ro power|
|00002e40| 22 29 3b 0a 58 09 09 2a | 72 65 73 20 3d 20 5f 6f |");.X..*|res = _o|
|00002e50| 6e 65 5f 3b 0a 58 09 09 | 72 65 74 75 72 6e 3b 0a |ne_;.X..|return;.|
|00002e60| 58 09 7d 0a 58 09 69 66 | 20 28 69 73 6c 65 6f 6e |X.}.X.if| (isleon|
|00002e70| 65 28 7a 31 29 29 20 7b | 09 2f 2a 20 30 2c 20 31 |e(z1)) {|./* 0, 1|
|00002e80| 2c 20 6f 72 20 2d 31 20 | 72 61 69 73 65 64 20 74 |, or -1 |raised t|
|00002e90| 6f 20 61 20 70 6f 77 65 | 72 20 2a 2f 0a 58 09 09 |o a powe|r */.X..|
|00002ea0| 61 6e 73 20 3d 20 5f 6f | 6e 65 5f 3b 0a 58 09 09 |ans = _o|ne_;.X..|
|00002eb0| 61 6e 73 2e 73 69 67 6e | 20 3d 20 28 42 4f 4f 4c |ans.sign| = (BOOL|
|00002ec0| 29 73 69 67 6e 3b 0a 58 | 09 09 69 66 20 28 2a 7a |)sign;.X|..if (*z|
|00002ed0| 31 2e 76 20 3d 3d 20 30 | 29 0a 58 09 09 09 61 6e |1.v == 0|).X...an|
|00002ee0| 73 20 3d 20 5f 7a 65 72 | 6f 5f 3b 0a 58 09 09 2a |s = _zer|o_;.X..*|
|00002ef0| 72 65 73 20 3d 20 61 6e | 73 3b 0a 58 09 09 72 65 |res = an|s;.X..re|
|00002f00| 74 75 72 6e 3b 0a 58 09 | 7d 0a 58 09 69 66 20 28 |turn;.X.|}.X.if (|
|00002f10| 69 73 62 69 67 28 7a 32 | 29 29 0a 58 09 09 65 72 |isbig(z2|)).X..er|
|00002f20| 72 6f 72 28 22 52 61 69 | 73 69 6e 67 20 74 6f 20 |ror("Rai|sing to |
|00002f30| 76 65 72 79 20 6c 61 72 | 67 65 20 70 6f 77 65 72 |very lar|ge power|
|00002f40| 22 29 3b 0a 58 09 70 6f | 77 65 72 20 3d 20 28 69 |");.X.po|wer = (i|
|00002f50| 73 74 69 6e 79 28 7a 32 | 29 20 3f 20 7a 31 74 6f |stiny(z2|) ? z1to|
|00002f60| 6c 28 7a 32 29 20 3a 20 | 7a 32 74 6f 6c 28 7a 32 |l(z2) : |z2tol(z2|
|00002f70| 29 29 3b 0a 58 09 69 66 | 20 28 69 73 74 77 6f 28 |));.X.if| (istwo(|
|00002f80| 7a 31 29 29 20 7b 09 2f | 2a 20 74 77 6f 20 72 61 |z1)) {./|* two ra|
|00002f90| 69 73 65 64 20 74 6f 20 | 61 20 70 6f 77 65 72 20 |ised to |a power |
|00002fa0| 2a 2f 0a 58 09 09 7a 62 | 69 74 76 61 6c 75 65 28 |*/.X..zb|itvalue(|
|00002fb0| 28 6c 6f 6e 67 29 20 70 | 6f 77 65 72 2c 20 72 65 |(long) p|ower, re|
|00002fc0| 73 29 3b 0a 58 09 09 72 | 65 74 75 72 6e 3b 0a 58 |s);.X..r|eturn;.X|
|00002fd0| 09 7d 0a 58 09 2f 2a 0a | 58 09 20 2a 20 53 65 65 |.}.X./*.|X. * See|
|00002fe0| 20 69 66 20 74 68 69 73 | 20 69 73 20 61 20 70 6f | if this| is a po|
|00002ff0| 77 65 72 20 6f 66 20 74 | 65 6e 0a 58 09 20 2a 2f |wer of t|en.X. */|
|00003000| 0a 58 09 69 66 20 28 69 | 73 74 69 6e 79 28 7a 31 |.X.if (i|stiny(z1|
|00003010| 29 20 26 26 20 28 2a 7a | 31 2e 76 20 3d 3d 20 31 |) && (*z|1.v == 1|
|00003020| 30 29 29 20 7b 0a 58 09 | 09 7a 74 65 6e 70 6f 77 |0)) {.X.|.ztenpow|
|00003030| 28 28 6c 6f 6e 67 29 20 | 70 6f 77 65 72 2c 20 72 |((long) |power, r|
|00003040| 65 73 29 3b 0a 58 09 09 | 72 65 73 2d 3e 73 69 67 |es);.X..|res->sig|
|00003050| 6e 20 3d 20 28 42 4f 4f | 4c 29 73 69 67 6e 3b 0a |n = (BOO|L)sign;.|
|00003060| 58 09 09 72 65 74 75 72 | 6e 3b 0a 58 09 7d 0a 58 |X..retur|n;.X.}.X|
|00003070| 09 2f 2a 0a 58 09 20 2a | 20 48 61 6e 64 6c 65 20 |./*.X. *| Handle |
|00003080| 6c 6f 77 20 70 6f 77 65 | 72 73 20 73 70 65 63 69 |low powe|rs speci|
|00003090| 61 6c 6c 79 0a 58 09 20 | 2a 2f 0a 58 09 69 66 20 |ally.X. |*/.X.if |
|000030a0| 28 70 6f 77 65 72 20 3c | 3d 20 34 29 20 7b 0a 58 |(power <|= 4) {.X|
|000030b0| 09 09 73 77 69 74 63 68 | 20 28 28 69 6e 74 29 20 |..switch| ((int) |
|000030c0| 70 6f 77 65 72 29 20 7b | 0a 58 09 09 09 63 61 73 |power) {|.X...cas|
|000030d0| 65 20 31 3a 0a 58 09 09 | 09 09 61 6e 73 2e 6c 65 |e 1:.X..|..ans.le|
|000030e0| 6e 20 3d 20 7a 31 2e 6c | 65 6e 3b 0a 58 09 09 09 |n = z1.l|en;.X...|
|000030f0| 09 61 6e 73 2e 76 20 3d | 20 61 6c 6c 6f 63 28 61 |.ans.v =| alloc(a|
|00003100| 6e 73 2e 6c 65 6e 29 3b | 0a 58 09 09 09 09 63 6f |ns.len);|.X....co|
|00003110| 70 79 76 61 6c 28 7a 31 | 2c 20 61 6e 73 29 3b 0a |pyval(z1|, ans);.|
|00003120| 58 09 09 09 09 61 6e 73 | 2e 73 69 67 6e 20 3d 20 |X....ans|.sign = |
|00003130| 28 42 4f 4f 4c 29 73 69 | 67 6e 3b 0a 58 09 09 09 |(BOOL)si|gn;.X...|
|00003140| 09 2a 72 65 73 20 3d 20 | 61 6e 73 3b 0a 58 09 09 |.*res = |ans;.X..|
|00003150| 09 09 72 65 74 75 72 6e | 3b 0a 58 09 09 09 63 61 |..return|;.X...ca|
|00003160| 73 65 20 32 3a 0a 58 09 | 09 09 09 7a 73 71 75 61 |se 2:.X.|...zsqua|
|00003170| 72 65 28 7a 31 2c 20 72 | 65 73 29 3b 0a 58 09 09 |re(z1, r|es);.X..|
|00003180| 09 09 72 65 74 75 72 6e | 3b 0a 58 09 09 09 63 61 |..return|;.X...ca|
|00003190| 73 65 20 33 3a 0a 58 09 | 09 09 09 7a 73 71 75 61 |se 3:.X.|...zsqua|
|000031a0| 72 65 28 7a 31 2c 20 26 | 74 65 6d 70 29 3b 0a 58 |re(z1, &|temp);.X|
|000031b0| 09 09 09 09 7a 6d 75 6c | 28 7a 31 2c 20 74 65 6d |....zmul|(z1, tem|
|000031c0| 70 2c 20 72 65 73 29 3b | 0a 58 09 09 09 09 66 72 |p, res);|.X....fr|
|000031d0| 65 65 68 28 74 65 6d 70 | 2e 76 29 3b 0a 58 09 09 |eeh(temp|.v);.X..|
|000031e0| 09 09 72 65 73 2d 3e 73 | 69 67 6e 20 3d 20 28 42 |..res->s|ign = (B|
|000031f0| 4f 4f 4c 29 73 69 67 6e | 3b 0a 58 09 09 09 09 72 |OOL)sign|;.X....r|
|00003200| 65 74 75 72 6e 3b 0a 58 | 09 09 09 63 61 73 65 20 |eturn;.X|...case |
|00003210| 34 3a 0a 58 09 09 09 09 | 7a 73 71 75 61 72 65 28 |4:.X....|zsquare(|
|00003220| 7a 31 2c 20 26 74 65 6d | 70 29 3b 0a 58 09 09 09 |z1, &tem|p);.X...|
|00003230| 09 7a 73 71 75 61 72 65 | 28 74 65 6d 70 2c 20 72 |.zsquare|(temp, r|
|00003240| 65 73 29 3b 0a 58 09 09 | 09 09 66 72 65 65 68 28 |es);.X..|..freeh(|
|00003250| 74 65 6d 70 2e 76 29 3b | 0a 58 09 09 09 09 72 65 |temp.v);|.X....re|
|00003260| 74 75 72 6e 3b 0a 58 09 | 09 7d 0a 58 09 7d 0a 58 |turn;.X.|.}.X.}.X|
|00003270| 09 2f 2a 0a 58 09 20 2a | 20 53 68 69 66 74 20 6f |./*.X. *| Shift o|
|00003280| 75 74 20 61 6c 6c 20 70 | 6f 77 65 72 73 20 6f 66 |ut all p|owers of|
|00003290| 20 74 77 6f 73 20 73 6f | 20 74 68 65 20 6d 75 6c | twos so| the mul|
|000032a0| 74 69 70 6c 69 65 73 20 | 61 72 65 20 73 6d 61 6c |tiplies |are smal|
|000032b0| 6c 65 72 2e 0a 58 09 20 | 2a 20 57 65 20 77 69 6c |ler..X. |* We wil|
|000032c0| 6c 20 73 68 69 66 74 20 | 62 61 63 6b 20 74 68 65 |l shift |back the|
|000032d0| 20 72 69 67 68 74 20 61 | 6d 6f 75 6e 74 20 77 68 | right a|mount wh|
|000032e0| 65 6e 20 64 6f 6e 65 2e | 0a 58 09 20 2a 2f 0a 58 |en done.|.X. */.X|
|000032f0| 09 74 77 6f 73 20 3d 20 | 30 3b 0a 58 09 69 66 20 |.twos = |0;.X.if |
|00003300| 28 69 73 65 76 65 6e 28 | 7a 31 29 29 20 7b 0a 58 |(iseven(|z1)) {.X|
|00003310| 09 09 74 77 6f 73 20 3d | 20 7a 6c 6f 77 62 69 74 |..twos =| zlowbit|
|00003320| 28 7a 31 29 3b 0a 58 09 | 09 61 6e 73 2e 76 20 3d |(z1);.X.|.ans.v =|
|00003330| 20 61 6c 6c 6f 63 28 7a | 31 2e 6c 65 6e 29 3b 0a | alloc(z|1.len);.|
|00003340| 58 09 09 61 6e 73 2e 6c | 65 6e 20 3d 20 7a 31 2e |X..ans.l|en = z1.|
|00003350| 6c 65 6e 3b 0a 58 09 09 | 63 6f 70 79 76 61 6c 28 |len;.X..|copyval(|
|00003360| 7a 31 2c 20 61 6e 73 29 | 3b 0a 58 09 09 73 68 69 |z1, ans)|;.X..shi|
|00003370| 66 74 72 28 61 6e 73 2c | 20 74 77 6f 73 29 3b 0a |ftr(ans,| twos);.|
|00003380| 58 09 09 74 72 69 6d 28 | 26 61 6e 73 29 3b 0a 58 |X..trim(|&ans);.X|
|00003390| 09 09 7a 31 20 3d 20 61 | 6e 73 3b 0a 58 09 09 74 |..z1 = a|ns;.X..t|
|000033a0| 77 6f 73 20 2a 3d 20 70 | 6f 77 65 72 3b 0a 58 09 |wos *= p|ower;.X.|
|000033b0| 7d 0a 58 09 2f 2a 0a 58 | 09 20 2a 20 43 6f 6d 70 |}.X./*.X|. * Comp|
|000033c0| 75 74 65 20 74 68 65 20 | 70 6f 77 65 72 20 62 79 |ute the |power by|
|000033d0| 20 73 71 75 61 72 69 6e | 67 20 61 6e 64 20 6d 75 | squarin|g and mu|
|000033e0| 6c 74 69 70 6c 79 69 6e | 67 2e 0a 58 09 20 2a 20 |ltiplyin|g..X. * |
|000033f0| 54 68 69 73 20 75 73 65 | 73 20 74 68 65 20 6c 65 |This use|s the le|
|00003400| 66 74 20 74 6f 20 72 69 | 67 68 74 20 6d 65 74 68 |ft to ri|ght meth|
|00003410| 6f 64 20 6f 66 20 70 6f | 77 65 72 20 72 61 69 73 |od of po|wer rais|
|00003420| 69 6e 67 2e 0a 58 09 20 | 2a 2f 0a 58 09 62 69 74 |ing..X. |*/.X.bit|
|00003430| 20 3d 20 54 4f 50 46 55 | 4c 4c 3b 0a 58 09 77 68 | = TOPFU|LL;.X.wh|
|00003440| 69 6c 65 20 28 28 62 69 | 74 20 26 20 70 6f 77 65 |ile ((bi|t & powe|
|00003450| 72 29 20 3d 3d 20 30 29 | 0a 58 09 09 62 69 74 20 |r) == 0)|.X..bit |
|00003460| 3e 3e 3d 20 31 4c 3b 0a | 58 09 62 69 74 20 3e 3e |>>= 1L;.|X.bit >>|
|00003470| 3d 20 31 4c 3b 0a 58 09 | 7a 73 71 75 61 72 65 28 |= 1L;.X.|zsquare(|
|00003480| 7a 31 2c 20 26 61 6e 73 | 29 3b 0a 58 09 69 66 20 |z1, &ans|);.X.if |
|00003490| 28 62 69 74 20 26 20 70 | 6f 77 65 72 29 20 7b 0a |(bit & p|ower) {.|
|000034a0| 58 09 09 7a 6d 75 6c 28 | 61 6e 73 2c 20 7a 31 2c |X..zmul(|ans, z1,|
|000034b0| 20 26 74 65 6d 70 29 3b | 0a 58 09 09 66 72 65 65 | &temp);|.X..free|
|000034c0| 68 28 61 6e 73 2e 76 29 | 3b 0a 58 09 09 61 6e 73 |h(ans.v)|;.X..ans|
|000034d0| 20 3d 20 74 65 6d 70 3b | 0a 58 09 7d 0a 58 09 62 | = temp;|.X.}.X.b|
|000034e0| 69 74 20 3e 3e 3d 20 31 | 4c 3b 0a 58 09 77 68 69 |it >>= 1|L;.X.whi|
|000034f0| 6c 65 20 28 62 69 74 29 | 20 7b 0a 58 09 09 7a 73 |le (bit)| {.X..zs|
|00003500| 71 75 61 72 65 28 61 6e | 73 2c 20 26 74 65 6d 70 |quare(an|s, &temp|
|00003510| 29 3b 0a 58 09 09 66 72 | 65 65 68 28 61 6e 73 2e |);.X..fr|eeh(ans.|
|00003520| 76 29 3b 0a 58 09 09 61 | 6e 73 20 3d 20 74 65 6d |v);.X..a|ns = tem|
|00003530| 70 3b 0a 58 09 09 69 66 | 20 28 62 69 74 20 26 20 |p;.X..if| (bit & |
|00003540| 70 6f 77 65 72 29 20 7b | 0a 58 09 09 09 7a 6d 75 |power) {|.X...zmu|
|00003550| 6c 28 61 6e 73 2c 20 7a | 31 2c 20 26 74 65 6d 70 |l(ans, z|1, &temp|
|00003560| 29 3b 0a 58 09 09 09 66 | 72 65 65 68 28 61 6e 73 |);.X...f|reeh(ans|
|00003570| 2e 76 29 3b 0a 58 09 09 | 09 61 6e 73 20 3d 20 74 |.v);.X..|.ans = t|
|00003580| 65 6d 70 3b 0a 58 09 09 | 7d 0a 58 09 09 62 69 74 |emp;.X..|}.X..bit|
|00003590| 20 3e 3e 3d 20 31 4c 3b | 0a 58 09 7d 0a 58 09 2f | >>= 1L;|.X.}.X./|
|000035a0| 2a 0a 58 09 20 2a 20 53 | 63 61 6c 65 20 62 61 63 |*.X. * S|cale bac|
|000035b0| 6b 20 75 70 20 62 79 20 | 70 72 6f 70 65 72 20 70 |k up by |proper p|
|000035c0| 6f 77 65 72 20 6f 66 20 | 74 77 6f 0a 58 09 20 2a |ower of |two.X. *|
|000035d0| 2f 0a 58 09 69 66 20 28 | 74 77 6f 73 29 20 7b 0a |/.X.if (|twos) {.|
|000035e0| 58 09 09 7a 73 68 69 66 | 74 28 61 6e 73 2c 20 74 |X..zshif|t(ans, t|
|000035f0| 77 6f 73 2c 20 26 74 65 | 6d 70 29 3b 0a 58 09 09 |wos, &te|mp);.X..|
|00003600| 66 72 65 65 68 28 61 6e | 73 2e 76 29 3b 0a 58 09 |freeh(an|s.v);.X.|
|00003610| 09 61 6e 73 20 3d 20 74 | 65 6d 70 3b 0a 58 09 09 |.ans = t|emp;.X..|
|00003620| 66 72 65 65 68 28 7a 31 | 2e 76 29 3b 0a 58 09 7d |freeh(z1|.v);.X.}|
|00003630| 0a 58 09 61 6e 73 2e 73 | 69 67 6e 20 3d 20 28 42 |.X.ans.s|ign = (B|
|00003640| 4f 4f 4c 29 73 69 67 6e | 3b 0a 58 09 2a 72 65 73 |OOL)sign|;.X.*res|
|00003650| 20 3d 20 61 6e 73 3b 0a | 58 7d 0a 58 0a 58 0a 58 | = ans;.|X}.X.X.X|
|00003660| 2f 2a 0a 58 20 2a 20 43 | 6f 6d 70 75 74 65 20 74 |/*.X * C|ompute t|
|00003670| 65 6e 20 74 6f 20 74 68 | 65 20 73 70 65 63 69 66 |en to th|e specif|
|00003680| 69 65 64 20 70 6f 77 65 | 72 0a 58 20 2a 20 54 68 |ied powe|r.X * Th|
|00003690| 69 73 20 73 61 76 65 73 | 20 73 6f 6d 65 20 77 6f |is saves| some wo|
|000036a0| 72 6b 20 73 69 6e 63 65 | 20 74 68 65 20 73 71 75 |rk since| the squ|
|000036b0| 61 72 65 73 20 6f 66 20 | 74 65 6e 20 61 72 65 20 |ares of |ten are |
|000036c0| 73 61 76 65 64 2e 0a 58 | 20 2a 2f 0a 58 76 6f 69 |saved..X| */.Xvoi|
|000036d0| 64 0a 58 7a 74 65 6e 70 | 6f 77 28 70 6f 77 65 72 |d.Xztenp|ow(power|
|000036e0| 2c 20 72 65 73 29 0a 58 | 09 6c 6f 6e 67 20 70 6f |, res).X|.long po|
|000036f0| 77 65 72 3b 0a 58 09 5a | 56 41 4c 55 45 20 2a 72 |wer;.X.Z|VALUE *r|
|00003700| 65 73 3b 0a 58 7b 0a 58 | 09 6c 6f 6e 67 20 69 3b |es;.X{.X|.long i;|
|00003710| 0a 58 09 5a 56 41 4c 55 | 45 20 61 6e 73 3b 0a 58 |.X.ZVALU|E ans;.X|
|00003720| 09 5a 56 41 4c 55 45 20 | 74 65 6d 70 3b 0a 58 0a |.ZVALUE |temp;.X.|
|00003730| 58 09 69 66 20 28 70 6f | 77 65 72 20 3c 3d 20 30 |X.if (po|wer <= 0|
|00003740| 29 20 7b 0a 58 09 09 2a | 72 65 73 20 3d 20 5f 6f |) {.X..*|res = _o|
|00003750| 6e 65 5f 3b 0a 58 09 09 | 72 65 74 75 72 6e 3b 0a |ne_;.X..|return;.|
|00003760| 58 09 7d 0a 58 09 61 6e | 73 20 3d 20 5f 6f 6e 65 |X.}.X.an|s = _one|
|00003770| 5f 3b 0a 58 09 5f 74 65 | 6e 70 6f 77 65 72 73 5f |_;.X._te|npowers_|
|00003780| 5b 30 5d 20 3d 20 5f 74 | 65 6e 5f 3b 0a 58 09 66 |[0] = _t|en_;.X.f|
|00003790| 6f 72 20 28 69 20 3d 20 | 30 3b 20 70 6f 77 65 72 |or (i = |0; power|
|000037a0| 3b 20 69 2b 2b 29 20 7b | 0a 58 09 09 69 66 20 28 |; i++) {|.X..if (|
|000037b0| 5f 74 65 6e 70 6f 77 65 | 72 73 5f 5b 69 5d 2e 6c |_tenpowe|rs_[i].l|
|000037c0| 65 6e 20 3d 3d 20 30 29 | 0a 58 09 09 09 7a 73 71 |en == 0)|.X...zsq|
|000037d0| 75 61 72 65 28 5f 74 65 | 6e 70 6f 77 65 72 73 5f |uare(_te|npowers_|
|000037e0| 5b 69 2d 31 5d 2c 20 26 | 5f 74 65 6e 70 6f 77 65 |[i-1], &|_tenpowe|
|000037f0| 72 73 5f 5b 69 5d 29 3b | 0a 58 09 09 69 66 20 28 |rs_[i]);|.X..if (|
|00003800| 70 6f 77 65 72 20 26 20 | 30 78 31 29 20 7b 0a 58 |power & |0x1) {.X|
|00003810| 09 09 09 7a 6d 75 6c 28 | 61 6e 73 2c 20 5f 74 65 |...zmul(|ans, _te|
|00003820| 6e 70 6f 77 65 72 73 5f | 5b 69 5d 2c 20 26 74 65 |npowers_|[i], &te|
|00003830| 6d 70 29 3b 0a 58 09 09 | 09 66 72 65 65 68 28 61 |mp);.X..|.freeh(a|
|00003840| 6e 73 2e 76 29 3b 0a 58 | 09 09 09 61 6e 73 20 3d |ns.v);.X|...ans =|
|00003850| 20 74 65 6d 70 3b 0a 58 | 09 09 7d 0a 58 09 09 70 | temp;.X|..}.X..p|
|00003860| 6f 77 65 72 20 2f 3d 20 | 32 3b 0a 58 09 7d 0a 58 |ower /= |2;.X.}.X|
|00003870| 09 2a 72 65 73 20 3d 20 | 61 6e 73 3b 0a 58 7d 0a |.*res = |ans;.X}.|
|00003880| 58 0a 58 0a 58 2f 2a 0a | 58 20 2a 20 43 61 6c 63 |X.X.X/*.|X * Calc|
|00003890| 75 6c 61 74 65 20 6d 6f | 64 75 6c 61 72 20 69 6e |ulate mo|dular in|
|000038a0| 76 65 72 73 65 20 73 75 | 70 70 72 65 73 73 69 6e |verse su|ppressin|
|000038b0| 67 20 75 6e 6e 65 63 65 | 73 73 61 72 79 20 64 69 |g unnece|ssary di|
|000038c0| 76 69 73 69 6f 6e 73 2e | 0a 58 20 2a 20 54 68 69 |visions.|.X * Thi|
|000038d0| 73 20 69 73 20 62 61 73 | 65 64 20 6f 6e 20 74 68 |s is bas|ed on th|
|000038e0| 65 20 45 75 63 6c 69 64 | 69 61 6e 20 61 6c 67 6f |e Euclid|ian algo|
|000038f0| 72 69 74 68 6d 20 66 6f | 72 20 6c 61 72 67 65 20 |rithm fo|r large |
|00003900| 6e 75 6d 62 65 72 73 2e | 0a 58 20 2a 20 28 41 6c |numbers.|.X * (Al|
|00003910| 67 6f 72 69 74 68 6d 20 | 58 20 66 72 6f 6d 20 4b |gorithm |X from K|
|00003920| 6e 75 74 68 20 56 6f 6c | 20 32 2c 20 73 65 63 74 |nuth Vol| 2, sect|
|00003930| 69 6f 6e 20 34 2e 35 2e | 32 2e 20 61 6e 64 20 65 |ion 4.5.|2. and e|
|00003940| 78 65 72 63 69 73 65 20 | 31 37 29 0a 58 20 2a 20 |xercise |17).X * |
|00003950| 52 65 74 75 72 6e 73 20 | 54 52 55 45 20 69 66 20 |Returns |TRUE if |
|00003960| 74 68 65 72 65 20 69 73 | 20 6e 6f 20 73 6f 6c 75 |there is| no solu|
|00003970| 74 69 6f 6e 20 62 65 63 | 61 75 73 65 20 74 68 65 |tion bec|ause the|
|00003980| 20 6e 75 6d 62 65 72 73 | 0a 58 20 2a 20 61 72 65 | numbers|.X * are|
|00003990| 20 6e 6f 74 20 72 65 6c | 61 74 69 76 65 6c 79 20 | not rel|atively |
|000039a0| 70 72 69 6d 65 2e 0a 58 | 20 2a 2f 0a 58 42 4f 4f |prime..X| */.XBOO|
|000039b0| 4c 0a 58 7a 6d 6f 64 69 | 6e 76 28 75 2c 20 76 2c |L.Xzmodi|nv(u, v,|
|000039c0| 20 72 65 73 29 0a 58 09 | 5a 56 41 4c 55 45 20 75 | res).X.|ZVALUE u|
|000039d0| 2c 20 76 3b 0a 58 09 5a | 56 41 4c 55 45 20 2a 72 |, v;.X.Z|VALUE *r|
|000039e0| 65 73 3b 0a 58 7b 0a 58 | 09 46 55 4c 4c 09 71 31 |es;.X{.X|.FULL.q1|
|000039f0| 2c 20 71 32 2c 20 75 69 | 33 2c 20 76 69 33 2c 20 |, q2, ui|3, vi3, |
|00003a00| 75 68 2c 20 76 68 2c 20 | 41 2c 20 42 2c 20 43 2c |uh, vh, |A, B, C,|
|00003a10| 20 44 2c 20 54 3b 0a 58 | 09 5a 56 41 4c 55 45 09 | D, T;.X|.ZVALUE.|
|00003a20| 75 32 2c 20 75 33 2c 20 | 76 32 2c 20 76 33 2c 20 |u2, u3, |v2, v3, |
|00003a30| 71 7a 2c 20 74 6d 70 31 | 2c 20 74 6d 70 32 2c 20 |qz, tmp1|, tmp2, |
|00003a40| 74 6d 70 33 3b 0a 58 0a | 58 09 69 66 20 28 69 73 |tmp3;.X.|X.if (is|
|00003a50| 6e 65 67 28 75 29 20 7c | 7c 20 69 73 6e 65 67 28 |neg(u) ||| isneg(|
|00003a60| 76 29 20 7c 7c 20 28 7a | 72 65 6c 28 75 2c 20 76 |v) || (z|rel(u, v|
|00003a70| 29 20 3e 3d 20 30 29 29 | 0a 58 09 09 7a 6d 6f 64 |) >= 0))|.X..zmod|
|00003a80| 28 75 2c 20 76 2c 20 26 | 76 33 29 3b 0a 58 09 65 |(u, v, &|v3);.X.e|
|00003a90| 6c 73 65 0a 58 09 09 7a | 63 6f 70 79 28 75 2c 20 |lse.X..z|copy(u, |
|00003aa0| 26 76 33 29 3b 0a 58 09 | 7a 63 6f 70 79 28 76 2c |&v3);.X.|zcopy(v,|
|00003ab0| 20 26 75 33 29 3b 0a 58 | 09 75 32 20 3d 20 5f 7a | &u3);.X|.u2 = _z|
|00003ac0| 65 72 6f 5f 3b 0a 58 09 | 76 32 20 3d 20 5f 6f 6e |ero_;.X.|v2 = _on|
|00003ad0| 65 5f 3b 0a 58 0a 58 09 | 2f 2a 0a 58 09 20 2a 20 |e_;.X.X.|/*.X. * |
|00003ae0| 4c 6f 6f 70 20 68 65 72 | 65 20 77 68 69 6c 65 20 |Loop her|e while |
|00003af0| 74 68 65 20 73 69 7a 65 | 20 6f 66 20 74 68 65 20 |the size| of the |
|00003b00| 6e 75 6d 62 65 72 73 20 | 72 65 6d 61 69 6e 20 61 |numbers |remain a|
|00003b10| 62 6f 76 65 0a 58 09 20 | 2a 20 74 68 65 20 73 69 |bove.X. |* the si|
|00003b20| 7a 65 20 6f 66 20 61 20 | 46 55 4c 4c 2e 20 20 54 |ze of a |FULL. T|
|00003b30| 68 72 6f 75 67 68 6f 75 | 74 20 74 68 69 73 20 6c |hroughou|t this l|
|00003b40| 6f 6f 70 20 75 33 20 3e | 3d 20 76 33 2e 0a 58 09 |oop u3 >|= v3..X.|
|00003b50| 20 2a 2f 0a 58 09 77 68 | 69 6c 65 20 28 28 75 33 | */.X.wh|ile ((u3|
|00003b60| 2e 6c 65 6e 20 3e 20 31 | 29 20 26 26 20 21 69 73 |.len > 1|) && !is|
|00003b70| 7a 65 72 6f 28 76 33 29 | 29 20 7b 0a 58 09 09 75 |zero(v3)|) {.X..u|
|00003b80| 68 20 3d 20 28 28 28 46 | 55 4c 4c 29 20 75 33 2e |h = (((F|ULL) u3.|
|00003b90| 76 5b 75 33 2e 6c 65 6e | 20 2d 20 31 5d 29 20 3c |v[u3.len| - 1]) <|
|00003ba0| 3c 20 42 41 53 45 42 29 | 20 2b 20 75 33 2e 76 5b |< BASEB)| + u3.v[|
|00003bb0| 75 33 2e 6c 65 6e 20 2d | 20 32 5d 3b 0a 58 09 09 |u3.len -| 2];.X..|
|00003bc0| 76 68 20 3d 20 30 3b 0a | 58 09 09 69 66 20 28 28 |vh = 0;.|X..if ((|
|00003bd0| 76 33 2e 6c 65 6e 20 2b | 20 31 29 20 3e 3d 20 75 |v3.len +| 1) >= u|
|00003be0| 33 2e 6c 65 6e 29 0a 58 | 09 09 09 76 68 20 3d 20 |3.len).X|...vh = |
|00003bf0| 76 33 2e 76 5b 76 33 2e | 6c 65 6e 20 2d 20 31 5d |v3.v[v3.|len - 1]|
|00003c00| 3b 0a 58 09 09 69 66 20 | 28 76 33 2e 6c 65 6e 20 |;.X..if |(v3.len |
|00003c10| 3d 3d 20 75 33 2e 6c 65 | 6e 29 0a 58 09 09 09 76 |== u3.le|n).X...v|
|00003c20| 68 20 3d 20 28 76 68 20 | 3c 3c 20 42 41 53 45 42 |h = (vh |<< BASEB|
|00003c30| 29 20 2b 20 76 33 2e 76 | 5b 76 33 2e 6c 65 6e 20 |) + v3.v|[v3.len |
|00003c40| 2d 20 32 5d 3b 0a 58 09 | 09 41 20 3d 20 31 3b 0a |- 2];.X.|.A = 1;.|
|00003c50| 58 09 09 42 20 3d 20 30 | 3b 0a 58 09 09 43 20 3d |X..B = 0|;.X..C =|
|00003c60| 20 30 3b 0a 58 09 09 44 | 20 3d 20 31 3b 0a 58 0a | 0;.X..D| = 1;.X.|
|00003c70| 58 09 09 2f 2a 0a 58 09 | 09 20 2a 20 43 61 6c 63 |X../*.X.|. * Calc|
|00003c80| 75 6c 61 74 65 20 73 75 | 63 63 65 73 73 69 76 65 |ulate su|ccessive|
|00003c90| 20 71 75 6f 74 69 65 6e | 74 73 20 6f 66 20 74 68 | quotien|ts of th|
|00003ca0| 65 20 63 6f 6e 74 69 6e | 75 65 64 20 66 72 61 63 |e contin|ued frac|
|00003cb0| 74 69 6f 6e 0a 58 09 09 | 20 2a 20 65 78 70 61 6e |tion.X..| * expan|
|00003cc0| 73 69 6f 6e 20 75 73 69 | 6e 67 20 6f 6e 6c 79 20 |sion usi|ng only |
|00003cd0| 73 69 6e 67 6c 65 20 70 | 72 65 63 69 73 69 6f 6e |single p|recision|
|00003ce0| 20 61 72 69 74 68 6d 65 | 74 69 63 20 75 6e 74 69 | arithme|tic unti|
|00003cf0| 6c 0a 58 09 09 20 2a 20 | 67 72 65 61 74 65 72 20 |l.X.. * |greater |
|00003d00| 70 72 65 63 69 73 69 6f | 6e 20 69 73 20 72 65 71 |precisio|n is req|
|00003d10| 75 69 72 65 64 2e 0a 58 | 09 09 20 2a 2f 0a 58 09 |uired..X|.. */.X.|
|00003d20| 09 77 68 69 6c 65 20 28 | 28 76 68 20 2b 20 43 29 |.while (|(vh + C)|
|00003d30| 20 26 26 20 28 76 68 20 | 2b 20 44 29 29 20 7b 0a | && (vh |+ D)) {.|
|00003d40| 58 09 09 09 71 31 20 3d | 20 28 75 68 20 2b 20 41 |X...q1 =| (uh + A|
|00003d50| 29 20 2f 20 28 76 68 20 | 2b 20 43 29 3b 0a 58 09 |) / (vh |+ C);.X.|
|00003d60| 09 09 71 32 20 3d 20 28 | 75 68 20 2b 20 42 29 20 |..q2 = (|uh + B) |
|00003d70| 2f 20 28 76 68 20 2b 20 | 44 29 3b 0a 58 09 09 09 |/ (vh + |D);.X...|
|00003d80| 69 66 20 28 71 31 20 21 | 3d 20 71 32 29 0a 58 09 |if (q1 !|= q2).X.|
|00003d90| 09 09 09 62 72 65 61 6b | 3b 0a 58 09 09 09 54 20 |...break|;.X...T |
|00003da0| 3d 20 41 20 2d 20 71 31 | 20 2a 20 43 3b 0a 58 09 |= A - q1| * C;.X.|
|00003db0| 09 09 41 20 3d 20 43 3b | 0a 58 09 09 09 43 20 3d |..A = C;|.X...C =|
|00003dc0| 20 54 3b 0a 58 09 09 09 | 54 20 3d 20 42 20 2d 20 | T;.X...|T = B - |
|00003dd0| 71 31 20 2a 20 44 3b 0a | 58 09 09 09 42 20 3d 20 |q1 * D;.|X...B = |
|00003de0| 44 3b 0a 58 09 09 09 44 | 20 3d 20 54 3b 0a 58 09 |D;.X...D| = T;.X.|
|00003df0| 09 09 54 20 3d 20 75 68 | 20 2d 20 71 31 20 2a 20 |..T = uh| - q1 * |
|00003e00| 76 68 3b 0a 58 09 09 09 | 75 68 20 3d 20 76 68 3b |vh;.X...|uh = vh;|
|00003e10| 0a 58 09 09 09 76 68 20 | 3d 20 54 3b 0a 58 09 09 |.X...vh |= T;.X..|
|00003e20| 7d 0a 58 09 0a 58 09 09 | 2f 2a 0a 58 09 09 20 2a |}.X..X..|/*.X.. *|
|00003e30| 20 49 66 20 42 20 69 73 | 20 7a 65 72 6f 2c 20 74 | If B is| zero, t|
|00003e40| 68 65 6e 20 77 65 20 6d | 61 64 65 20 6e 6f 20 70 |hen we m|ade no p|
|00003e50| 72 6f 67 72 65 73 73 20 | 62 65 63 61 75 73 65 0a |rogress |because.|
|00003e60| 58 09 09 20 2a 20 74 68 | 65 20 63 61 6c 63 75 6c |X.. * th|e calcul|
|00003e70| 61 74 69 6f 6e 20 72 65 | 71 75 69 72 65 73 20 61 |ation re|quires a|
|00003e80| 20 76 65 72 79 20 6c 61 | 72 67 65 20 71 75 6f 74 | very la|rge quot|
|00003e90| 69 65 6e 74 2e 0a 58 09 | 09 20 2a 20 53 6f 20 77 |ient..X.|. * So w|
|00003ea0| 65 20 6d 75 73 74 20 64 | 6f 20 74 68 69 73 20 73 |e must d|o this s|
|00003eb0| 74 65 70 20 6f 66 20 74 | 68 65 20 63 61 6c 63 75 |tep of t|he calcu|
|00003ec0| 6c 61 74 69 6f 6e 20 69 | 6e 0a 58 09 09 20 2a 20 |lation i|n.X.. * |
|00003ed0| 66 75 6c 6c 20 70 72 65 | 63 69 73 69 6f 6e 0a 58 |full pre|cision.X|
|00003ee0| 09 09 20 2a 2f 0a 58 09 | 09 69 66 20 28 42 20 3d |.. */.X.|.if (B =|
|00003ef0| 3d 20 30 29 20 7b 0a 58 | 09 09 09 7a 71 75 6f 28 |= 0) {.X|...zquo(|
|00003f00| 75 33 2c 20 76 33 2c 20 | 26 71 7a 29 3b 0a 58 09 |u3, v3, |&qz);.X.|
|00003f10| 09 09 7a 6d 75 6c 28 71 | 7a 2c 20 76 32 2c 20 26 |..zmul(q|z, v2, &|
|00003f20| 74 6d 70 31 29 3b 0a 58 | 09 09 09 7a 73 75 62 28 |tmp1);.X|...zsub(|
|00003f30| 75 32 2c 20 74 6d 70 31 | 2c 20 26 74 6d 70 32 29 |u2, tmp1|, &tmp2)|
|00003f40| 3b 0a 58 09 09 09 66 72 | 65 65 68 28 74 6d 70 31 |;.X...fr|eeh(tmp1|
|00003f50| 2e 76 29 3b 0a 58 09 09 | 09 66 72 65 65 68 28 75 |.v);.X..|.freeh(u|
|00003f60| 32 2e 76 29 3b 0a 58 09 | 09 09 75 32 20 3d 20 76 |2.v);.X.|..u2 = v|
|00003f70| 32 3b 0a 58 09 09 09 76 | 32 20 3d 20 74 6d 70 32 |2;.X...v|2 = tmp2|
|00003f80| 3b 0a 58 09 09 09 7a 6d | 75 6c 28 71 7a 2c 20 76 |;.X...zm|ul(qz, v|
|00003f90| 33 2c 20 26 74 6d 70 31 | 29 3b 0a 58 09 09 09 7a |3, &tmp1|);.X...z|
|00003fa0| 73 75 62 28 75 33 2c 20 | 74 6d 70 31 2c 20 26 74 |sub(u3, |tmp1, &t|
|00003fb0| 6d 70 32 29 3b 0a 58 09 | 09 09 66 72 65 65 68 28 |mp2);.X.|..freeh(|
|00003fc0| 74 6d 70 31 2e 76 29 3b | 0a 58 09 09 09 66 72 65 |tmp1.v);|.X...fre|
|00003fd0| 65 68 28 75 33 2e 76 29 | 3b 0a 58 09 09 09 75 33 |eh(u3.v)|;.X...u3|
|00003fe0| 20 3d 20 76 33 3b 0a 58 | 09 09 09 76 33 20 3d 20 | = v3;.X|...v3 = |
|00003ff0| 74 6d 70 32 3b 0a 58 09 | 09 09 66 72 65 65 68 28 |tmp2;.X.|..freeh(|
|00004000| 71 7a 2e 76 29 3b 0a 58 | 09 09 09 63 6f 6e 74 69 |qz.v);.X|...conti|
|00004010| 6e 75 65 3b 0a 58 09 09 | 7d 0a 58 09 09 2f 2a 0a |nue;.X..|}.X../*.|
|00004020| 58 09 09 20 2a 20 41 70 | 70 6c 79 20 74 68 65 20 |X.. * Ap|ply the |
|00004030| 63 61 6c 63 75 6c 61 74 | 65 64 20 41 2c 42 2c 43 |calculat|ed A,B,C|
|00004040| 2c 44 20 6e 75 6d 62 65 | 72 73 20 74 6f 20 74 68 |,D numbe|rs to th|
|00004050| 65 20 63 75 72 72 65 6e | 74 0a 58 09 09 20 2a 20 |e curren|t.X.. * |
|00004060| 76 61 6c 75 65 73 20 74 | 6f 20 75 70 64 61 74 65 |values t|o update|
|00004070| 20 74 68 65 6d 20 61 73 | 20 69 66 20 74 68 65 20 | them as| if the |
|00004080| 66 75 6c 6c 20 70 72 65 | 63 69 73 69 6f 6e 0a 58 |full pre|cision.X|
|00004090| 09 09 20 2a 20 63 61 6c | 63 75 6c 61 74 69 6f 6e |.. * cal|culation|
|000040a0| 73 20 68 61 64 20 62 65 | 65 6e 20 63 61 72 72 69 |s had be|en carri|
|000040b0| 65 64 20 6f 75 74 2e 0a | 58 09 09 20 2a 2f 0a 58 |ed out..|X.. */.X|
|000040c0| 09 09 7a 6d 75 6c 69 28 | 75 32 2c 20 28 6c 6f 6e |..zmuli(|u2, (lon|
|000040d0| 67 29 20 41 2c 20 26 74 | 6d 70 31 29 3b 0a 58 09 |g) A, &t|mp1);.X.|
|000040e0| 09 7a 6d 75 6c 69 28 76 | 32 2c 20 28 6c 6f 6e 67 |.zmuli(v|2, (long|
|000040f0| 29 20 42 2c 20 26 74 6d | 70 32 29 3b 0a 58 09 09 |) B, &tm|p2);.X..|
|00004100| 7a 61 64 64 28 74 6d 70 | 31 2c 20 74 6d 70 32 2c |zadd(tmp|1, tmp2,|
|00004110| 20 26 74 6d 70 33 29 3b | 0a 58 09 09 66 72 65 65 | &tmp3);|.X..free|
|00004120| 68 28 74 6d 70 31 2e 76 | 29 3b 0a 58 09 09 66 72 |h(tmp1.v|);.X..fr|
|00004130| 65 65 68 28 74 6d 70 32 | 2e 76 29 3b 0a 58 09 09 |eeh(tmp2|.v);.X..|
|00004140| 7a 6d 75 6c 69 28 75 32 | 2c 20 28 6c 6f 6e 67 29 |zmuli(u2|, (long)|
|00004150| 20 43 2c 20 26 74 6d 70 | 31 29 3b 0a 58 09 09 7a | C, &tmp|1);.X..z|
|00004160| 6d 75 6c 69 28 76 32 2c | 20 28 6c 6f 6e 67 29 20 |muli(v2,| (long) |
|00004170| 44 2c 20 26 74 6d 70 32 | 29 3b 0a 58 09 09 66 72 |D, &tmp2|);.X..fr|
|00004180| 65 65 68 28 75 32 2e 76 | 29 3b 0a 58 09 09 66 72 |eeh(u2.v|);.X..fr|
|00004190| 65 65 68 28 76 32 2e 76 | 29 3b 0a 58 09 09 75 32 |eeh(v2.v|);.X..u2|
|000041a0| 20 3d 20 74 6d 70 33 3b | 0a 58 09 09 7a 61 64 64 | = tmp3;|.X..zadd|
|000041b0| 28 74 6d 70 31 2c 20 74 | 6d 70 32 2c 20 26 76 32 |(tmp1, t|mp2, &v2|
|000041c0| 29 3b 0a 58 09 09 66 72 | 65 65 68 28 74 6d 70 31 |);.X..fr|eeh(tmp1|
|000041d0| 2e 76 29 3b 0a 58 09 09 | 66 72 65 65 68 28 74 6d |.v);.X..|freeh(tm|
|000041e0| 70 32 2e 76 29 3b 0a 58 | 09 09 7a 6d 75 6c 69 28 |p2.v);.X|..zmuli(|
|000041f0| 75 33 2c 20 28 6c 6f 6e | 67 29 20 41 2c 20 26 74 |u3, (lon|g) A, &t|
|00004200| 6d 70 31 29 3b 0a 58 09 | 09 7a 6d 75 6c 69 28 76 |mp1);.X.|.zmuli(v|
|00004210| 33 2c 20 28 6c 6f 6e 67 | 29 20 42 2c 20 26 74 6d |3, (long|) B, &tm|
|00004220| 70 32 29 3b 0a 58 09 09 | 7a 61 64 64 28 74 6d 70 |p2);.X..|zadd(tmp|
|00004230| 31 2c 20 74 6d 70 32 2c | 20 26 74 6d 70 33 29 3b |1, tmp2,| &tmp3);|
|00004240| 0a 58 09 09 66 72 65 65 | 68 28 74 6d 70 31 2e 76 |.X..free|h(tmp1.v|
|00004250| 29 3b 0a 58 09 09 66 72 | 65 65 68 28 74 6d 70 32 |);.X..fr|eeh(tmp2|
|00004260| 2e 76 29 3b 0a 58 09 09 | 7a 6d 75 6c 69 28 75 33 |.v);.X..|zmuli(u3|
|00004270| 2c 20 28 6c 6f 6e 67 29 | 20 43 2c 20 26 74 6d 70 |, (long)| C, &tmp|
|00004280| 31 29 3b 0a 58 09 09 7a | 6d 75 6c 69 28 76 33 2c |1);.X..z|muli(v3,|
|00004290| 20 28 6c 6f 6e 67 29 20 | 44 2c 20 26 74 6d 70 32 | (long) |D, &tmp2|
|000042a0| 29 3b 0a 58 09 09 66 72 | 65 65 68 28 75 33 2e 76 |);.X..fr|eeh(u3.v|
|000042b0| 29 3b 0a 58 09 09 66 72 | 65 65 68 28 76 33 2e 76 |);.X..fr|eeh(v3.v|
|000042c0| 29 3b 0a 58 09 09 75 33 | 20 3d 20 74 6d 70 33 3b |);.X..u3| = tmp3;|
|000042d0| 0a 58 09 09 7a 61 64 64 | 28 74 6d 70 31 2c 20 74 |.X..zadd|(tmp1, t|
|000042e0| 6d 70 32 2c 20 26 76 33 | 29 3b 0a 58 09 09 66 72 |mp2, &v3|);.X..fr|
|000042f0| 65 65 68 28 74 6d 70 31 | 2e 76 29 3b 0a 58 09 09 |eeh(tmp1|.v);.X..|
|00004300| 66 72 65 65 68 28 74 6d | 70 32 2e 76 29 3b 0a 58 |freeh(tm|p2.v);.X|
|00004310| 09 7d 0a 58 0a 58 09 2f | 2a 0a 58 09 20 2a 20 48 |.}.X.X./|*.X. * H|
|00004320| 65 72 65 20 77 68 65 6e | 20 74 68 65 20 72 65 6d |ere when| the rem|
|00004330| 61 69 6e 69 6e 67 20 6e | 75 6d 62 65 72 73 20 62 |aining n|umbers b|
|00004340| 65 63 6f 6d 65 20 73 69 | 6e 67 6c 65 20 70 72 65 |ecome si|ngle pre|
|00004350| 63 69 73 69 6f 6e 20 69 | 6e 20 73 69 7a 65 2e 0a |cision i|n size..|
|00004360| 58 09 20 2a 20 46 69 6e | 69 73 68 20 74 68 65 20 |X. * Fin|ish the |
|00004370| 70 72 6f 63 65 64 75 72 | 65 20 75 73 69 6e 67 20 |procedur|e using |
|00004380| 73 69 6e 67 6c 65 20 70 | 72 65 63 69 73 69 6f 6e |single p|recision|
|00004390| 20 63 61 6c 63 75 6c 61 | 74 69 6f 6e 73 2e 0a 58 | calcula|tions..X|
|000043a0| 09 20 2a 2f 0a 58 09 69 | 66 20 28 69 73 7a 65 72 |. */.X.i|f (iszer|
|000043b0| 6f 28 76 33 29 20 26 26 | 20 21 69 73 6f 6e 65 28 |o(v3) &&| !isone(|
|000043c0| 75 33 29 29 20 7b 0a 58 | 09 09 66 72 65 65 68 28 |u3)) {.X|..freeh(|
|000043d0| 75 33 2e 76 29 3b 0a 58 | 09 09 66 72 65 65 68 28 |u3.v);.X|..freeh(|
|000043e0| 76 33 2e 76 29 3b 0a 58 | 09 09 66 72 65 65 68 28 |v3.v);.X|..freeh(|
|000043f0| 75 32 2e 76 29 3b 0a 58 | 09 09 66 72 65 65 68 28 |u2.v);.X|..freeh(|
|00004400| 76 32 2e 76 29 3b 0a 58 | 09 09 72 65 74 75 72 6e |v2.v);.X|..return|
|00004410| 20 54 52 55 45 3b 0a 58 | 09 7d 0a 58 09 75 69 33 | TRUE;.X|.}.X.ui3|
|00004420| 20 3d 20 28 69 73 74 69 | 6e 79 28 75 33 29 20 3f | = (isti|ny(u3) ?|
|00004430| 20 7a 31 74 6f 6c 28 75 | 33 29 20 3a 20 7a 32 74 | z1tol(u|3) : z2t|
|00004440| 6f 6c 28 75 33 29 29 3b | 0a 58 09 76 69 33 20 3d |ol(u3));|.X.vi3 =|
|00004450| 20 28 69 73 74 69 6e 79 | 28 76 33 29 20 3f 20 7a | (istiny|(v3) ? z|
|00004460| 31 74 6f 6c 28 76 33 29 | 20 3a 20 7a 32 74 6f 6c |1tol(v3)| : z2tol|
|00004470| 28 76 33 29 29 3b 0a 58 | 09 66 72 65 65 68 28 75 |(v3));.X|.freeh(u|
|00004480| 33 2e 76 29 3b 0a 58 09 | 66 72 65 65 68 28 76 33 |3.v);.X.|freeh(v3|
|00004490| 2e 76 29 3b 0a 58 09 77 | 68 69 6c 65 20 28 76 69 |.v);.X.w|hile (vi|
|000044a0| 33 29 20 7b 0a 58 09 09 | 71 31 20 3d 20 75 69 33 |3) {.X..|q1 = ui3|
|000044b0| 20 2f 20 76 69 33 3b 0a | 58 09 09 7a 6d 75 6c 69 | / vi3;.|X..zmuli|
|000044c0| 28 76 32 2c 20 28 6c 6f | 6e 67 29 20 71 31 2c 20 |(v2, (lo|ng) q1, |
|000044d0| 26 74 6d 70 31 29 3b 0a | 58 09 09 7a 73 75 62 28 |&tmp1);.|X..zsub(|
|000044e0| 75 32 2c 20 74 6d 70 31 | 2c 20 26 74 6d 70 32 29 |u2, tmp1|, &tmp2)|
|000044f0| 3b 0a 58 09 09 66 72 65 | 65 68 28 74 6d 70 31 2e |;.X..fre|eh(tmp1.|
|00004500| 76 29 3b 0a 58 09 09 66 | 72 65 65 68 28 75 32 2e |v);.X..f|reeh(u2.|
|00004510| 76 29 3b 0a 58 09 09 75 | 32 20 3d 20 76 32 3b 0a |v);.X..u|2 = v2;.|
|00004520| 58 09 09 76 32 20 3d 20 | 74 6d 70 32 3b 0a 58 09 |X..v2 = |tmp2;.X.|
|00004530| 09 71 32 20 3d 20 75 69 | 33 20 2d 20 71 31 20 2a |.q2 = ui|3 - q1 *|
|00004540| 20 76 69 33 3b 0a 58 09 | 09 75 69 33 20 3d 20 76 | vi3;.X.|.ui3 = v|
|00004550| 69 33 3b 0a 58 09 09 76 | 69 33 20 3d 20 71 32 3b |i3;.X..v|i3 = q2;|
|00004560| 0a 58 09 7d 0a 58 09 66 | 72 65 65 68 28 76 32 2e |.X.}.X.f|reeh(v2.|
|00004570| 76 29 3b 0a 58 09 69 66 | 20 28 75 69 33 20 21 3d |v);.X.if| (ui3 !=|
|00004580| 20 31 29 20 7b 0a 58 09 | 09 66 72 65 65 68 28 75 | 1) {.X.|.freeh(u|
|00004590| 32 2e 76 29 3b 0a 58 09 | 09 72 65 74 75 72 6e 20 |2.v);.X.|.return |
|000045a0| 54 52 55 45 3b 0a 58 09 | 7d 0a 58 09 69 66 20 28 |TRUE;.X.|}.X.if (|
|000045b0| 69 73 6e 65 67 28 75 32 | 29 29 20 7b 0a 58 09 09 |isneg(u2|)) {.X..|
|000045c0| 7a 61 64 64 28 76 2c 20 | 75 32 2c 20 72 65 73 29 |zadd(v, |u2, res)|
|000045d0| 3b 0a 58 09 09 66 72 65 | 65 68 28 75 32 2e 76 29 |;.X..fre|eh(u2.v)|
|000045e0| 3b 0a 58 09 09 72 65 74 | 75 72 6e 20 46 41 4c 53 |;.X..ret|urn FALS|
|000045f0| 45 3b 0a 58 09 7d 0a 58 | 09 2a 72 65 73 20 3d 20 |E;.X.}.X|.*res = |
|00004600| 75 32 3b 0a 58 09 72 65 | 74 75 72 6e 20 46 41 4c |u2;.X.re|turn FAL|
|00004610| 53 45 3b 0a 58 7d 0a 58 | 0a 58 0a 58 23 69 66 20 |SE;.X}.X|.X.X#if |
|00004620| 30 0a 58 2f 2a 0a 58 20 | 2a 20 41 70 70 72 6f 78 |0.X/*.X |* Approx|
|00004630| 69 6d 61 74 65 20 74 68 | 65 20 71 75 6f 74 69 65 |imate th|e quotie|
|00004640| 6e 74 20 6f 66 20 74 77 | 6f 20 69 6e 74 65 67 65 |nt of tw|o intege|
|00004650| 72 73 20 62 79 20 61 6e | 6f 74 68 65 72 20 73 65 |rs by an|other se|
|00004660| 74 20 6f 66 20 73 6d 61 | 6c 6c 65 72 0a 58 20 2a |t of sma|ller.X *|
|00004670| 20 69 6e 74 65 67 65 72 | 73 2e 20 20 54 68 69 73 | integer|s. This|
|00004680| 20 75 73 65 73 20 63 6f | 6e 74 69 6e 75 65 64 20 | uses co|ntinued |
|00004690| 66 72 61 63 74 69 6f 6e | 73 20 74 6f 20 64 65 74 |fraction|s to det|
|000046a0| 65 72 6d 69 6e 65 20 74 | 68 65 20 73 6d 61 6c 6c |ermine t|he small|
|000046b0| 65 72 20 73 65 74 2e 0a | 58 20 2a 2f 0a 58 76 6f |er set..|X */.Xvo|
|000046c0| 69 64 0a 58 7a 61 70 70 | 72 6f 78 28 7a 31 2c 20 |id.Xzapp|rox(z1, |
|000046d0| 7a 32 2c 20 72 65 73 31 | 2c 20 72 65 73 32 29 0a |z2, res1|, res2).|
|000046e0| 58 09 5a 56 41 4c 55 45 | 20 7a 31 2c 20 7a 32 2c |X.ZVALUE| z1, z2,|
|000046f0| 20 2a 72 65 73 31 2c 20 | 2a 72 65 73 32 3b 0a 58 | *res1, |*res2;.X|
|00004700| 7b 0a 58 09 69 6e 74 20 | 73 69 67 6e 3b 0a 58 09 |{.X.int |sign;.X.|
|00004710| 5a 56 41 4c 55 45 20 75 | 31 2c 20 76 31 2c 20 75 |ZVALUE u|1, v1, u|
|00004720| 33 2c 20 76 33 2c 20 71 | 2c 20 74 31 2c 20 74 32 |3, v3, q|, t1, t2|
|00004730| 2c 20 74 33 3b 0a 58 0a | 58 09 73 69 67 6e 20 3d |, t3;.X.|X.sign =|
|00004740| 20 28 28 7a 31 2e 73 69 | 67 6e 20 21 3d 20 30 29 | ((z1.si|gn != 0)|
|00004750| 20 5e 20 28 7a 32 2e 73 | 69 67 6e 20 21 3d 20 30 | ^ (z2.s|ign != 0|
|00004760| 29 29 3b 0a 58 09 7a 31 | 2e 73 69 67 6e 20 3d 20 |));.X.z1|.sign = |
|00004770| 30 3b 0a 58 09 7a 32 2e | 73 69 67 6e 20 3d 20 30 |0;.X.z2.|sign = 0|
|00004780| 3b 0a 58 09 76 33 20 3d | 20 7a 32 3b 0a 58 09 75 |;.X.v3 =| z2;.X.u|
|00004790| 33 20 3d 20 7a 31 3b 0a | 58 09 75 31 20 3d 20 5f |3 = z1;.|X.u1 = _|
|000047a0| 6f 6e 65 5f 3b 0a 58 09 | 76 31 20 3d 20 5f 7a 65 |one_;.X.|v1 = _ze|
|000047b0| 72 6f 5f 3b 0a 58 09 77 | 68 69 6c 65 20 28 21 69 |ro_;.X.w|hile (!i|
|000047c0| 73 7a 65 72 6f 28 76 33 | 29 29 20 7b 0a 58 09 09 |szero(v3|)) {.X..|
|000047d0| 7a 64 69 76 28 75 33 2c | 20 76 33 2c 20 26 71 2c |zdiv(u3,| v3, &q,|
|000047e0| 20 26 74 31 29 3b 0a 58 | 09 09 7a 6d 75 6c 28 76 | &t1);.X|..zmul(v|
|000047f0| 31 2c 20 71 2c 20 26 74 | 32 29 3b 0a 58 09 09 7a |1, q, &t|2);.X..z|
|00004800| 73 75 62 28 75 31 2c 20 | 74 32 2c 20 26 74 33 29 |sub(u1, |t2, &t3)|
|00004810| 3b 0a 58 09 09 66 72 65 | 65 68 28 71 2e 76 29 3b |;.X..fre|eh(q.v);|
|00004820| 0a 58 09 09 66 72 65 65 | 68 28 74 32 2e 76 29 3b |.X..free|h(t2.v);|
|00004830| 0a 58 09 09 66 72 65 65 | 68 28 75 31 2e 76 29 3b |.X..free|h(u1.v);|
|00004840| 0a 58 09 09 69 66 20 28 | 28 75 33 2e 76 20 21 3d |.X..if (|(u3.v !=|
|00004850| 20 7a 31 2e 76 29 20 26 | 26 20 28 75 33 2e 76 20 | z1.v) &|& (u3.v |
|00004860| 21 3d 20 7a 32 2e 76 29 | 29 0a 58 09 09 09 66 72 |!= z2.v)|).X...fr|
|00004870| 65 65 68 28 75 33 2e 76 | 29 3b 0a 58 09 09 75 31 |eeh(u3.v|);.X..u1|
|00004880| 20 3d 20 76 31 3b 0a 58 | 09 09 75 33 20 3d 20 76 | = v1;.X|..u3 = v|
|00004890| 33 3b 0a 58 09 09 76 31 | 20 3d 20 74 33 3b 0a 58 |3;.X..v1| = t3;.X|
|000048a0| 09 09 76 33 20 3d 20 74 | 31 3b 0a 58 09 7d 0a 58 |..v3 = t|1;.X.}.X|
|000048b0| 09 69 66 20 28 21 69 73 | 75 6e 69 74 28 75 33 29 |.if (!is|unit(u3)|
|000048c0| 29 0a 58 09 09 65 72 72 | 6f 72 28 22 4e 6f 6e 2d |).X..err|or("Non-|
|000048d0| 72 65 6c 61 74 69 76 6c | 79 20 70 72 69 6d 65 20 |relativl|y prime |
|000048e0| 6e 75 6d 62 65 72 73 20 | 66 6f 72 20 61 70 70 72 |numbers |for appr|
|000048f0| 6f 78 22 29 3b 0a 58 09 | 69 66 20 28 28 75 33 2e |ox");.X.|if ((u3.|
|00004900| 76 20 21 3d 20 7a 31 2e | 76 29 20 26 26 20 28 75 |v != z1.|v) && (u|
|00004910| 33 2e 76 20 21 3d 20 7a | 32 2e 76 29 29 0a 58 09 |3.v != z|2.v)).X.|
|00004920| 09 66 72 65 65 68 28 75 | 33 2e 76 29 3b 0a 58 09 |.freeh(u|3.v);.X.|
|00004930| 69 66 20 28 28 76 33 2e | 76 20 21 3d 20 7a 31 2e |if ((v3.|v != z1.|
|00004940| 76 29 20 26 26 20 28 76 | 33 2e 76 20 21 3d 20 7a |v) && (v|3.v != z|
|00004950| 32 2e 76 29 29 0a 58 09 | 09 66 72 65 65 68 28 76 |2.v)).X.|.freeh(v|
|00004960| 33 2e 76 29 3b 0a 58 09 | 66 72 65 65 68 28 76 31 |3.v);.X.|freeh(v1|
|00004970| 2e 76 29 3b 0a 58 09 7a | 6d 75 6c 28 75 31 2c 20 |.v);.X.z|mul(u1, |
|00004980| 7a 31 2c 20 26 74 31 29 | 3b 0a 58 09 7a 73 75 62 |z1, &t1)|;.X.zsub|
|00004990| 28 74 31 2c 20 5f 6f 6e | 65 5f 2c 20 26 74 32 29 |(t1, _on|e_, &t2)|
|000049a0| 3b 0a 58 09 66 72 65 65 | 68 28 74 31 2e 76 29 3b |;.X.free|h(t1.v);|
|000049b0| 0a 58 09 7a 71 75 6f 28 | 74 32 2c 20 7a 32 2c 20 |.X.zquo(|t2, z2, |
|000049c0| 26 74 31 29 3b 0a 58 09 | 66 72 65 65 68 28 74 32 |&t1);.X.|freeh(t2|
|000049d0| 2e 76 29 3b 0a 58 09 75 | 31 2e 73 69 67 6e 20 3d |.v);.X.u|1.sign =|
|000049e0| 20 28 42 4f 4f 4c 29 73 | 69 67 6e 3b 0a 58 09 74 | (BOOL)s|ign;.X.t|
|000049f0| 31 2e 73 69 67 6e 20 3d | 20 30 3b 0a 58 09 2a 72 |1.sign =| 0;.X.*r|
|00004a00| 65 73 31 20 3d 20 74 31 | 3b 0a 58 09 2a 72 65 73 |es1 = t1|;.X.*res|
|00004a10| 32 20 3d 20 75 31 3b 0a | 58 7d 0a 58 23 65 6e 64 |2 = u1;.|X}.X#end|
|00004a20| 69 66 0a 58 0a 58 0a 58 | 2f 2a 0a 58 20 2a 20 42 |if.X.X.X|/*.X * B|
|00004a30| 69 6e 61 72 79 20 67 63 | 64 20 61 6c 67 6f 72 69 |inary gc|d algori|
|00004a40| 74 68 6d 0a 58 20 2a 20 | 54 68 69 73 20 61 6c 67 |thm.X * |This alg|
|00004a50| 6f 72 69 74 68 6d 20 74 | 61 6b 65 6e 20 66 72 6f |orithm t|aken fro|
|00004a60| 6d 20 4b 6e 75 74 68 0a | 58 20 2a 2f 0a 58 76 6f |m Knuth.|X */.Xvo|
|00004a70| 69 64 0a 58 7a 67 63 64 | 28 7a 31 2c 20 7a 32 2c |id.Xzgcd|(z1, z2,|
|00004a80| 20 72 65 73 29 0a 58 09 | 5a 56 41 4c 55 45 20 7a | res).X.|ZVALUE z|
|00004a90| 31 2c 20 7a 32 2c 20 2a | 72 65 73 3b 0a 58 7b 0a |1, z2, *|res;.X{.|
|00004aa0| 58 09 5a 56 41 4c 55 45 | 20 75 2c 20 76 2c 20 74 |X.ZVALUE| u, v, t|
|00004ab0| 3b 0a 58 09 72 65 67 69 | 73 74 65 72 20 6c 6f 6e |;.X.regi|ster lon|
|00004ac0| 67 20 6a 2c 20 6b 2c 20 | 6f 6c 65 6e 2c 20 6d 61 |g j, k, |olen, ma|
|00004ad0| 73 6b 3b 0a 58 09 72 65 | 67 69 73 74 65 72 20 48 |sk;.X.re|gister H|
|00004ae0| 41 4c 46 20 68 3b 0a 58 | 09 48 41 4c 46 20 2a 6f |ALF h;.X|.HALF *o|
|00004af0| 6c 64 76 31 2c 20 2a 6f | 6c 64 76 32 3b 0a 58 0a |ldv1, *o|ldv2;.X.|
|00004b00| 58 09 2f 2a 0a 58 09 20 | 2a 20 46 69 72 73 74 20 |X./*.X. |* First |
|00004b10| 73 65 65 20 69 66 20 6f | 6e 65 20 6e 75 6d 62 65 |see if o|ne numbe|
|00004b20| 72 20 69 73 20 76 65 72 | 79 20 6d 75 63 68 20 6c |r is ver|y much l|
|00004b30| 61 72 67 65 72 20 74 68 | 61 6e 20 74 68 65 20 6f |arger th|an the o|
|00004b40| 74 68 65 72 2e 0a 58 09 | 20 2a 20 49 66 20 73 6f |ther..X.| * If so|
|00004b50| 2c 20 74 68 65 6e 20 64 | 69 76 69 64 65 20 61 73 |, then d|ivide as|
|00004b60| 20 6e 65 63 65 73 73 61 | 72 79 20 74 6f 20 67 65 | necessa|ry to ge|
|00004b70| 74 20 74 68 65 20 6e 75 | 6d 62 65 72 73 20 6e 65 |t the nu|mbers ne|
|00004b80| 61 72 20 65 61 63 68 20 | 6f 74 68 65 72 2e 0a 58 |ar each |other..X|
|00004b90| 09 20 2a 2f 0a 58 09 7a | 31 2e 73 69 67 6e 20 3d |. */.X.z|1.sign =|
|00004ba0| 20 30 3b 0a 58 09 7a 32 | 2e 73 69 67 6e 20 3d 20 | 0;.X.z2|.sign = |
|00004bb0| 30 3b 0a 58 09 6f 6c 64 | 76 31 20 3d 20 7a 31 2e |0;.X.old|v1 = z1.|
|00004bc0| 76 3b 0a 58 09 6f 6c 64 | 76 32 20 3d 20 7a 32 2e |v;.X.old|v2 = z2.|
|00004bd0| 76 3b 0a 58 09 69 66 20 | 28 7a 31 2e 6c 65 6e 20 |v;.X.if |(z1.len |
|00004be0| 3c 20 7a 32 2e 6c 65 6e | 29 20 7b 0a 58 09 09 74 |< z2.len|) {.X..t|
|00004bf0| 20 3d 20 7a 31 3b 0a 58 | 09 09 7a 31 20 3d 20 7a | = z1;.X|..z1 = z|
|00004c00| 32 3b 0a 58 09 09 7a 32 | 20 3d 20 74 3b 0a 58 09 |2;.X..z2| = t;.X.|
|00004c10| 7d 0a 58 09 77 68 69 6c | 65 20 28 28 7a 31 2e 6c |}.X.whil|e ((z1.l|
|00004c20| 65 6e 20 3e 20 28 7a 32 | 2e 6c 65 6e 20 2b 20 35 |en > (z2|.len + 5|
|00004c30| 29 29 20 26 26 20 21 69 | 73 7a 65 72 6f 28 7a 32 |)) && !i|szero(z2|
|00004c40| 29 29 20 7b 0a 58 09 09 | 7a 6d 6f 64 28 7a 31 2c |)) {.X..|zmod(z1,|
|00004c50| 20 7a 32 2c 20 26 74 29 | 3b 0a 58 09 09 69 66 20 | z2, &t)|;.X..if |
|00004c60| 28 28 7a 31 2e 76 20 21 | 3d 20 6f 6c 64 76 31 29 |((z1.v !|= oldv1)|
|00004c70| 20 26 26 20 28 7a 31 2e | 76 20 21 3d 20 6f 6c 64 | && (z1.|v != old|
|00004c80| 76 32 29 29 0a 58 09 09 | 09 66 72 65 65 68 28 7a |v2)).X..|.freeh(z|
|00004c90| 31 2e 76 29 3b 0a 58 09 | 09 7a 31 20 3d 20 7a 32 |1.v);.X.|.z1 = z2|
|00004ca0| 3b 0a 58 09 09 7a 32 20 | 3d 20 74 3b 0a 58 09 7d |;.X..z2 |= t;.X.}|
|00004cb0| 0a 58 09 2f 2a 0a 58 09 | 20 2a 20 4f 6b 2c 20 6e |.X./*.X.| * Ok, n|
|00004cc0| 6f 77 20 64 6f 20 74 68 | 65 20 62 69 6e 61 72 79 |ow do th|e binary|
|00004cd0| 20 6d 65 74 68 6f 64 20 | 70 72 6f 70 65 72 0a 58 | method |proper.X|
|00004ce0| 09 20 2a 2f 0a 58 09 75 | 2e 6c 65 6e 20 3d 20 7a |. */.X.u|.len = z|
|00004cf0| 31 2e 6c 65 6e 3b 0a 58 | 09 76 2e 6c 65 6e 20 3d |1.len;.X|.v.len =|
|00004d00| 20 7a 32 2e 6c 65 6e 3b | 0a 58 09 75 2e 73 69 67 | z2.len;|.X.u.sig|
|00004d10| 6e 20 3d 20 30 3b 0a 58 | 09 76 2e 73 69 67 6e 20 |n = 0;.X|.v.sign |
|00004d20| 3d 20 30 3b 0a 58 09 69 | 66 20 28 21 7a 74 65 73 |= 0;.X.i|f (!ztes|
|00004d30| 74 28 7a 31 29 29 20 7b | 0a 58 09 09 76 2e 76 20 |t(z1)) {|.X..v.v |
|00004d40| 3d 20 61 6c 6c 6f 63 28 | 76 2e 6c 65 6e 29 3b 0a |= alloc(|v.len);.|
|00004d50| 58 09 09 63 6f 70 79 76 | 61 6c 28 7a 32 2c 20 76 |X..copyv|al(z2, v|
|00004d60| 29 3b 0a 58 09 09 2a 72 | 65 73 20 3d 20 76 3b 0a |);.X..*r|es = v;.|
|00004d70| 58 09 09 67 6f 74 6f 20 | 64 6f 6e 65 3b 0a 58 09 |X..goto |done;.X.|
|00004d80| 7d 0a 58 09 69 66 20 28 | 21 7a 74 65 73 74 28 7a |}.X.if (|!ztest(z|
|00004d90| 32 29 29 20 7b 0a 58 09 | 09 75 2e 76 20 3d 20 61 |2)) {.X.|.u.v = a|
|00004da0| 6c 6c 6f 63 28 75 2e 6c | 65 6e 29 3b 0a 58 09 09 |lloc(u.l|en);.X..|
|00004db0| 63 6f 70 79 76 61 6c 28 | 7a 31 2c 20 75 29 3b 0a |copyval(|z1, u);.|
|00004dc0| 58 09 09 2a 72 65 73 20 | 3d 20 75 3b 0a 58 09 09 |X..*res |= u;.X..|
|00004dd0| 67 6f 74 6f 20 64 6f 6e | 65 3b 0a 58 09 7d 0a 58 |goto don|e;.X.}.X|
|00004de0| 09 75 2e 76 20 3d 20 61 | 6c 6c 6f 63 28 75 2e 6c |.u.v = a|lloc(u.l|
|00004df0| 65 6e 29 3b 0a 58 09 76 | 2e 76 20 3d 20 61 6c 6c |en);.X.v|.v = all|
|00004e00| 6f 63 28 76 2e 6c 65 6e | 29 3b 0a 58 09 63 6f 70 |oc(v.len|);.X.cop|
|00004e10| 79 76 61 6c 28 7a 31 2c | 20 75 29 3b 0a 58 09 63 |yval(z1,| u);.X.c|
|00004e20| 6f 70 79 76 61 6c 28 7a | 32 2c 20 76 29 3b 0a 58 |opyval(z|2, v);.X|
|00004e30| 09 6b 20 3d 20 30 3b 0a | 58 09 77 68 69 6c 65 20 |.k = 0;.|X.while |
|00004e40| 28 75 2e 76 5b 6b 5d 20 | 3d 3d 20 30 20 26 26 20 |(u.v[k] |== 0 && |
|00004e50| 76 2e 76 5b 6b 5d 20 3d | 3d 20 30 29 0a 58 09 09 |v.v[k] =|= 0).X..|
|00004e60| 2b 2b 6b 3b 0a 58 09 6d | 61 73 6b 20 3d 20 30 31 |++k;.X.m|ask = 01|
|00004e70| 3b 0a 58 09 68 20 3d 20 | 75 2e 76 5b 6b 5d 20 7c |;.X.h = |u.v[k] ||
|00004e80| 20 76 2e 76 5b 6b 5d 3b | 0a 58 09 6b 20 2a 3d 20 | v.v[k];|.X.k *= |
|00004e90| 42 41 53 45 42 3b 0a 58 | 09 77 68 69 6c 65 20 28 |BASEB;.X|.while (|
|00004ea0| 21 28 68 20 26 20 6d 61 | 73 6b 29 29 20 7b 0a 58 |!(h & ma|sk)) {.X|
|00004eb0| 09 09 6d 61 73 6b 20 3c | 3c 3d 20 31 3b 0a 58 09 |..mask <|<= 1;.X.|
|00004ec0| 09 6b 2b 2b 3b 0a 58 09 | 7d 0a 58 09 73 68 69 66 |.k++;.X.|}.X.shif|
|00004ed0| 74 72 28 75 2c 20 6b 29 | 3b 0a 58 09 73 68 69 66 |tr(u, k)|;.X.shif|
|00004ee0| 74 72 28 76 2c 20 6b 29 | 3b 0a 58 09 74 72 69 6d |tr(v, k)|;.X.trim|
|00004ef0| 28 26 75 29 3b 0a 58 09 | 74 72 69 6d 28 26 76 29 |(&u);.X.|trim(&v)|
|00004f00| 3b 0a 58 09 69 66 20 28 | 69 73 6f 64 64 28 75 29 |;.X.if (|isodd(u)|
|00004f10| 29 20 7b 0a 58 09 09 74 | 2e 76 20 3d 20 61 6c 6c |) {.X..t|.v = all|
|00004f20| 6f 63 28 76 2e 6c 65 6e | 29 3b 0a 58 09 09 74 2e |oc(v.len|);.X..t.|
|00004f30| 6c 65 6e 20 3d 20 76 2e | 6c 65 6e 3b 0a 58 09 09 |len = v.|len;.X..|
|00004f40| 63 6f 70 79 76 61 6c 28 | 76 2c 20 74 29 3b 0a 58 |copyval(|v, t);.X|
|00004f50| 09 09 74 2e 73 69 67 6e | 20 3d 20 21 76 2e 73 69 |..t.sign| = !v.si|
|00004f60| 67 6e 3b 0a 58 09 7d 20 | 65 6c 73 65 20 7b 0a 58 |gn;.X.} |else {.X|
|00004f70| 09 09 74 2e 76 20 3d 20 | 61 6c 6c 6f 63 28 75 2e |..t.v = |alloc(u.|
|00004f80| 6c 65 6e 29 3b 0a 58 09 | 09 74 2e 6c 65 6e 20 3d |len);.X.|.t.len =|
|00004f90| 20 75 2e 6c 65 6e 3b 0a | 58 09 09 63 6f 70 79 76 | u.len;.|X..copyv|
|00004fa0| 61 6c 28 75 2c 20 74 29 | 3b 0a 58 09 09 74 2e 73 |al(u, t)|;.X..t.s|
|00004fb0| 69 67 6e 20 3d 20 75 2e | 73 69 67 6e 3b 0a 58 09 |ign = u.|sign;.X.|
|00004fc0| 7d 0a 58 09 77 68 69 6c | 65 20 28 7a 74 65 73 74 |}.X.whil|e (ztest|
|00004fd0| 28 74 29 29 20 7b 0a 58 | 09 09 6a 20 3d 20 30 3b |(t)) {.X|..j = 0;|
|00004fe0| 0a 58 09 09 77 68 69 6c | 65 20 28 21 74 2e 76 5b |.X..whil|e (!t.v[|
|00004ff0| 6a 5d 29 0a 58 09 09 09 | 2b 2b 6a 3b 0a 58 09 09 |j]).X...|++j;.X..|
|00005000| 6d 61 73 6b 20 3d 20 30 | 31 3b 0a 58 09 09 68 20 |mask = 0|1;.X..h |
|00005010| 3d 20 74 2e 76 5b 6a 5d | 3b 0a 58 09 09 6a 20 2a |= t.v[j]|;.X..j *|
|00005020| 3d 20 42 41 53 45 42 3b | 0a 58 09 09 77 68 69 6c |= BASEB;|.X..whil|
|00005030| 65 20 28 21 28 68 20 26 | 20 6d 61 73 6b 29 29 20 |e (!(h &| mask)) |
|00005040| 7b 0a 58 09 09 09 6d 61 | 73 6b 20 3c 3c 3d 20 31 |{.X...ma|sk <<= 1|
|00005050| 3b 0a 58 09 09 09 6a 2b | 2b 3b 0a 58 09 09 7d 0a |;.X...j+|+;.X..}.|
|00005060| 58 09 09 73 68 69 66 74 | 72 28 74 2c 20 6a 29 3b |X..shift|r(t, j);|
|00005070| 0a 58 09 09 74 72 69 6d | 28 26 74 29 3b 0a 58 09 |.X..trim|(&t);.X.|
|00005080| 09 69 66 20 28 7a 74 65 | 73 74 28 74 29 20 3e 20 |.if (zte|st(t) > |
|00005090| 30 29 20 7b 0a 58 09 09 | 09 66 72 65 65 68 28 75 |0) {.X..|.freeh(u|
|000050a0| 2e 76 29 3b 0a 58 09 09 | 09 75 20 3d 20 74 3b 0a |.v);.X..|.u = t;.|
|000050b0| 58 09 09 7d 20 65 6c 73 | 65 20 7b 0a 58 09 09 09 |X..} els|e {.X...|
|000050c0| 66 72 65 65 68 28 76 2e | 76 29 3b 0a 58 09 09 09 |freeh(v.|v);.X...|
|000050d0| 76 20 3d 20 74 3b 0a 58 | 09 09 09 76 2e 73 69 67 |v = t;.X|...v.sig|
|000050e0| 6e 20 3d 20 21 74 2e 73 | 69 67 6e 3b 0a 58 09 09 |n = !t.s|ign;.X..|
|000050f0| 7d 0a 58 09 09 7a 73 75 | 62 28 75 2c 20 76 2c 20 |}.X..zsu|b(u, v, |
|00005100| 26 74 29 3b 0a 58 09 7d | 0a 58 09 66 72 65 65 68 |&t);.X.}|.X.freeh|
|00005110| 28 74 2e 76 29 3b 0a 58 | 09 66 72 65 65 68 28 76 |(t.v);.X|.freeh(v|
|00005120| 2e 76 29 3b 0a 58 09 69 | 66 20 28 6b 29 20 7b 0a |.v);.X.i|f (k) {.|
|00005130| 58 09 09 6f 6c 65 6e 20 | 3d 20 75 2e 6c 65 6e 3b |X..olen |= u.len;|
|00005140| 0a 58 09 09 75 2e 6c 65 | 6e 20 2b 3d 20 6b 20 2f |.X..u.le|n += k /|
|00005150| 20 42 41 53 45 42 20 2b | 20 31 3b 0a 58 09 09 75 | BASEB +| 1;.X..u|
|00005160| 2e 76 20 3d 20 28 48 41 | 4c 46 20 2a 29 20 72 65 |.v = (HA|LF *) re|
|00005170| 61 6c 6c 6f 63 28 75 2e | 76 2c 20 75 2e 6c 65 6e |alloc(u.|v, u.len|
|00005180| 20 2a 20 73 69 7a 65 6f | 66 28 48 41 4c 46 29 29 | * sizeo|f(HALF))|
|00005190| 3b 0a 58 09 09 77 68 69 | 6c 65 20 28 6f 6c 65 6e |;.X..whi|le (olen|
|000051a0| 20 21 3d 20 75 2e 6c 65 | 6e 29 0a 58 09 09 09 75 | != u.le|n).X...u|
|000051b0| 2e 76 5b 6f 6c 65 6e 2b | 2b 5d 20 3d 20 30 3b 0a |.v[olen+|+] = 0;.|
|000051c0| 58 09 09 73 68 69 66 74 | 6c 28 75 2c 20 6b 29 3b |X..shift|l(u, k);|
|000051d0| 0a 58 09 7d 0a 58 09 74 | 72 69 6d 28 26 75 29 3b |.X.}.X.t|rim(&u);|
|000051e0| 0a 58 09 2a 72 65 73 20 | 3d 20 75 3b 0a 58 0a 58 |.X.*res |= u;.X.X|
|000051f0| 64 6f 6e 65 3a 0a 58 09 | 69 66 20 28 28 7a 31 2e |done:.X.|if ((z1.|
|00005200| 76 20 21 3d 20 6f 6c 64 | 76 31 29 20 26 26 20 28 |v != old|v1) && (|
|00005210| 7a 31 2e 76 20 21 3d 20 | 6f 6c 64 76 32 29 29 0a |z1.v != |oldv2)).|
|00005220| 58 09 09 66 72 65 65 68 | 28 7a 31 2e 76 29 3b 0a |X..freeh|(z1.v);.|
|00005230| 58 09 69 66 20 28 28 7a | 32 2e 76 20 21 3d 20 6f |X.if ((z|2.v != o|
|00005240| 6c 64 76 31 29 20 26 26 | 20 28 7a 32 2e 76 20 21 |ldv1) &&| (z2.v !|
|00005250| 3d 20 6f 6c 64 76 32 29 | 29 0a 58 09 09 66 72 65 |= oldv2)|).X..fre|
|00005260| 65 68 28 7a 32 2e 76 29 | 3b 0a 58 7d 0a 58 0a 58 |eh(z2.v)|;.X}.X.X|
|00005270| 0a 58 2f 2a 0a 58 20 2a | 20 43 6f 6d 70 75 74 65 |.X/*.X *| Compute|
|00005280| 20 74 68 65 20 6c 63 6d | 20 6f 66 20 74 77 6f 20 | the lcm| of two |
|00005290| 69 6e 74 65 67 65 72 73 | 20 28 6c 65 61 73 74 20 |integers| (least |
|000052a0| 63 6f 6d 6d 6f 6e 20 6d | 75 6c 74 69 70 6c 65 29 |common m|ultiple)|
|000052b0| 2e 0a 58 20 2a 20 54 68 | 69 73 20 69 73 20 64 6f |..X * Th|is is do|
|000052c0| 6e 65 20 75 73 69 6e 67 | 20 74 68 65 20 66 6f 72 |ne using| the for|
|000052d0| 6d 75 6c 61 3a 20 20 67 | 63 64 28 61 2c 62 29 20 |mula: g|cd(a,b) |
|000052e0| 2a 20 6c 63 6d 28 61 2c | 62 29 20 3d 20 61 20 2a |* lcm(a,|b) = a *|
|000052f0| 20 62 2e 0a 58 20 2a 2f | 0a 58 76 6f 69 64 0a 58 | b..X */|.Xvoid.X|
|00005300| 7a 6c 63 6d 28 7a 31 2c | 20 7a 32 2c 20 72 65 73 |zlcm(z1,| z2, res|
|00005310| 29 0a 58 09 5a 56 41 4c | 55 45 20 7a 31 2c 20 7a |).X.ZVAL|UE z1, z|
|00005320| 32 2c 20 2a 72 65 73 3b | 0a 58 7b 0a 58 09 5a 56 |2, *res;|.X{.X.ZV|
|00005330| 41 4c 55 45 20 74 65 6d | 70 31 2c 20 74 65 6d 70 |ALUE tem|p1, temp|
|00005340| 32 3b 0a 58 0a 58 09 7a | 67 63 64 28 7a 31 2c 20 |2;.X.X.z|gcd(z1, |
|00005350| 7a 32 2c 20 26 74 65 6d | 70 31 29 3b 0a 58 09 7a |z2, &tem|p1);.X.z|
|00005360| 71 75 6f 28 7a 31 2c 20 | 74 65 6d 70 31 2c 20 26 |quo(z1, |temp1, &|
|00005370| 74 65 6d 70 32 29 3b 0a | 58 09 66 72 65 65 68 28 |temp2);.|X.freeh(|
|00005380| 74 65 6d 70 31 2e 76 29 | 3b 0a 58 09 7a 6d 75 6c |temp1.v)|;.X.zmul|
|00005390| 28 74 65 6d 70 32 2c 20 | 7a 32 2c 20 72 65 73 29 |(temp2, |z2, res)|
|000053a0| 3b 0a 58 09 66 72 65 65 | 68 28 74 65 6d 70 32 2e |;.X.free|h(temp2.|
|000053b0| 76 29 3b 0a 58 7d 0a 58 | 0a 58 0a 58 2f 2a 0a 58 |v);.X}.X|.X.X/*.X|
|000053c0| 20 2a 20 52 65 74 75 72 | 6e 20 77 68 65 74 68 65 | * Retur|n whethe|
|000053d0| 72 20 6f 72 20 6e 6f 74 | 20 74 77 6f 20 6e 75 6d |r or not| two num|
|000053e0| 62 65 72 73 20 61 72 65 | 20 72 65 6c 61 74 69 76 |bers are| relativ|
|000053f0| 65 6c 79 20 70 72 69 6d | 65 20 74 6f 20 65 61 63 |ely prim|e to eac|
|00005400| 68 20 6f 74 68 65 72 2e | 0a 58 20 2a 2f 0a 58 42 |h other.|.X */.XB|
|00005410| 4f 4f 4c 0a 58 7a 72 65 | 6c 70 72 69 6d 65 28 7a |OOL.Xzre|lprime(z|
|00005420| 31 2c 20 7a 32 29 0a 58 | 09 5a 56 41 4c 55 45 20 |1, z2).X|.ZVALUE |
|00005430| 7a 31 2c 20 7a 32 3b 09 | 09 09 2f 2a 20 6e 75 6d |z1, z2;.|../* num|
|00005440| 62 65 72 73 20 74 6f 20 | 62 65 20 74 65 73 74 65 |bers to |be teste|
|00005450| 64 20 2a 2f 0a 58 7b 0a | 58 09 46 55 4c 4c 20 72 |d */.X{.|X.FULL r|
|00005460| 65 6d 31 2c 20 72 65 6d | 32 3b 09 09 2f 2a 20 72 |em1, rem|2;../* r|
|00005470| 65 6d 61 69 6e 64 65 72 | 73 20 2a 2f 0a 58 09 5a |emainder|s */.X.Z|
|00005480| 56 41 4c 55 45 20 72 65 | 6d 3b 0a 58 09 42 4f 4f |VALUE re|m;.X.BOO|
|00005490| 4c 20 72 65 73 75 6c 74 | 3b 0a 58 0a 58 09 7a 31 |L result|;.X.X.z1|
|000054a0| 2e 73 69 67 6e 20 3d 20 | 30 3b 0a 58 09 7a 32 2e |.sign = |0;.X.z2.|
|000054b0| 73 69 67 6e 20 3d 20 30 | 3b 0a 58 09 69 66 20 28 |sign = 0|;.X.if (|
|000054c0| 69 73 65 76 65 6e 28 7a | 31 29 20 26 26 20 69 73 |iseven(z|1) && is|
|000054d0| 65 76 65 6e 28 7a 32 29 | 29 09 2f 2a 20 66 61 6c |even(z2)|)./* fal|
|000054e0| 73 65 20 69 66 20 62 6f | 74 68 20 65 76 65 6e 20 |se if bo|th even |
|000054f0| 2a 2f 0a 58 09 09 72 65 | 74 75 72 6e 20 46 41 4c |*/.X..re|turn FAL|
|00005500| 53 45 3b 0a 58 09 69 66 | 20 28 69 73 75 6e 69 74 |SE;.X.if| (isunit|
|00005510| 28 7a 31 29 20 7c 7c 20 | 69 73 75 6e 69 74 28 7a |(z1) || |isunit(z|
|00005520| 32 29 29 09 2f 2a 20 74 | 72 75 65 20 69 66 20 65 |2))./* t|rue if e|
|00005530| 69 74 68 65 72 20 69 73 | 20 61 20 75 6e 69 74 20 |ither is| a unit |
|00005540| 2a 2f 0a 58 09 09 72 65 | 74 75 72 6e 20 54 52 55 |*/.X..re|turn TRU|
|00005550| 45 3b 0a 58 09 69 66 20 | 28 69 73 7a 65 72 6f 28 |E;.X.if |(iszero(|
|00005560| 7a 31 29 20 7c 7c 20 69 | 73 7a 65 72 6f 28 7a 32 |z1) || i|szero(z2|
|00005570| 29 29 09 2f 2a 20 66 61 | 6c 73 65 20 69 66 20 65 |))./* fa|lse if e|
|00005580| 69 74 68 65 72 20 69 73 | 20 7a 65 72 6f 20 2a 2f |ither is| zero */|
|00005590| 0a 58 09 09 72 65 74 75 | 72 6e 20 46 41 4c 53 45 |.X..retu|rn FALSE|
|000055a0| 3b 0a 58 09 69 66 20 28 | 69 73 74 77 6f 28 7a 31 |;.X.if (|istwo(z1|
|000055b0| 29 20 7c 7c 20 69 73 74 | 77 6f 28 7a 32 29 29 09 |) || ist|wo(z2)).|
|000055c0| 2f 2a 20 74 72 75 65 20 | 69 66 20 65 69 74 68 65 |/* true |if eithe|
|000055d0| 72 20 69 73 20 74 77 6f | 20 2a 2f 0a 58 09 09 72 |r is two| */.X..r|
|000055e0| 65 74 75 72 6e 20 54 52 | 55 45 3b 0a 58 09 2f 2a |eturn TR|UE;.X./*|
|000055f0| 0a 58 09 20 2a 20 54 72 | 79 20 72 65 64 75 63 69 |.X. * Tr|y reduci|
|00005600| 6e 67 20 65 61 63 68 20 | 6e 75 6d 62 65 72 20 62 |ng each |number b|
|00005610| 79 20 74 68 65 20 70 72 | 6f 64 75 63 74 20 6f 66 |y the pr|oduct of|
|00005620| 20 74 68 65 20 66 69 72 | 73 74 20 66 65 77 20 6f | the fir|st few o|
|00005630| 64 64 20 70 72 69 6d 65 | 73 0a 58 09 20 2a 20 74 |dd prime|s.X. * t|
|00005640| 6f 20 73 65 65 20 69 66 | 20 61 6e 79 20 6f 66 20 |o see if| any of |
|00005650| 74 68 65 6d 20 61 72 65 | 20 61 20 63 6f 6d 6d 6f |them are| a commo|
|00005660| 6e 20 66 61 63 74 6f 72 | 2e 0a 58 09 20 2a 2f 0a |n factor|..X. */.|
|00005670| 58 09 72 65 6d 31 20 3d | 20 7a 6d 6f 64 69 28 7a |X.rem1 =| zmodi(z|
|00005680| 31 2c 20 33 4c 20 2a 20 | 35 20 2a 20 37 20 2a 20 |1, 3L * |5 * 7 * |
|00005690| 31 31 20 2a 20 31 33 29 | 3b 0a 58 09 72 65 6d 32 |11 * 13)|;.X.rem2|
|000056a0| 20 3d 20 7a 6d 6f 64 69 | 28 7a 32 2c 20 33 4c 20 | = zmodi|(z2, 3L |
|000056b0| 2a 20 35 20 2a 20 37 20 | 2a 20 31 31 20 2a 20 31 |* 5 * 7 |* 11 * 1|
|000056c0| 33 29 3b 0a 58 09 69 66 | 20 28 28 28 72 65 6d 31 |3);.X.if| (((rem1|
|000056d0| 20 25 20 33 29 20 3d 3d | 20 30 29 20 26 26 20 28 | % 3) ==| 0) && (|
|000056e0| 28 72 65 6d 32 20 25 20 | 33 29 20 3d 3d 20 30 29 |(rem2 % |3) == 0)|
|000056f0| 29 0a 58 09 09 72 65 74 | 75 72 6e 20 46 41 4c 53 |).X..ret|urn FALS|
|00005700| 45 3b 0a 58 09 69 66 20 | 28 28 28 72 65 6d 31 20 |E;.X.if |(((rem1 |
|00005710| 25 20 35 29 20 3d 3d 20 | 30 29 20 26 26 20 28 28 |% 5) == |0) && ((|
|00005720| 72 65 6d 32 20 25 20 35 | 29 20 3d 3d 20 30 29 29 |rem2 % 5|) == 0))|
|00005730| 0a 58 09 09 72 65 74 75 | 72 6e 20 46 41 4c 53 45 |.X..retu|rn FALSE|
|00005740| 3b 0a 58 09 69 66 20 28 | 28 28 72 65 6d 31 20 25 |;.X.if (|((rem1 %|
|00005750| 20 37 29 20 3d 3d 20 30 | 29 20 26 26 20 28 28 72 | 7) == 0|) && ((r|
|00005760| 65 6d 32 20 25 20 37 29 | 20 3d 3d 20 30 29 29 0a |em2 % 7)| == 0)).|
|00005770| 58 09 09 72 65 74 75 72 | 6e 20 46 41 4c 53 45 3b |X..retur|n FALSE;|
|00005780| 0a 58 09 69 66 20 28 28 | 28 72 65 6d 31 20 25 20 |.X.if ((|(rem1 % |
|00005790| 31 31 29 20 3d 3d 20 30 | 29 20 26 26 20 28 28 72 |11) == 0|) && ((r|
|000057a0| 65 6d 32 20 25 20 31 31 | 29 20 3d 3d 20 30 29 29 |em2 % 11|) == 0))|
|000057b0| 0a 58 09 09 72 65 74 75 | 72 6e 20 46 41 4c 53 45 |.X..retu|rn FALSE|
|000057c0| 3b 0a 58 09 69 66 20 28 | 28 28 72 65 6d 31 20 25 |;.X.if (|((rem1 %|
|000057d0| 20 31 33 29 20 3d 3d 20 | 30 29 20 26 26 20 28 28 | 13) == |0) && ((|
|000057e0| 72 65 6d 32 20 25 20 31 | 33 29 20 3d 3d 20 30 29 |rem2 % 1|3) == 0)|
|000057f0| 29 0a 58 09 09 72 65 74 | 75 72 6e 20 46 41 4c 53 |).X..ret|urn FALS|
|00005800| 45 3b 0a 58 09 2f 2a 0a | 58 09 20 2a 20 54 72 79 |E;.X./*.|X. * Try|
|00005810| 20 61 20 6e 65 77 20 62 | 61 74 63 68 20 6f 66 20 | a new b|atch of |
|00005820| 70 72 69 6d 65 73 20 6e | 6f 77 0a 58 09 20 2a 2f |primes n|ow.X. */|
|00005830| 0a 58 09 72 65 6d 31 20 | 3d 20 7a 6d 6f 64 69 28 |.X.rem1 |= zmodi(|
|00005840| 7a 31 2c 20 31 37 4c 20 | 2a 20 31 39 20 2a 20 32 |z1, 17L |* 19 * 2|
|00005850| 33 29 3b 0a 58 09 72 65 | 6d 32 20 3d 20 7a 6d 6f |3);.X.re|m2 = zmo|
|00005860| 64 69 28 7a 32 2c 20 31 | 37 4c 20 2a 20 31 39 20 |di(z2, 1|7L * 19 |
|00005870| 2a 20 32 33 29 3b 0a 58 | 09 69 66 20 28 28 28 72 |* 23);.X|.if (((r|
|00005880| 65 6d 31 20 25 20 31 37 | 29 20 3d 3d 20 30 29 20 |em1 % 17|) == 0) |
|00005890| 26 26 20 28 28 72 65 6d | 32 20 25 20 31 37 29 20 |&& ((rem|2 % 17) |
|000058a0| 3d 3d 20 30 29 29 0a 58 | 09 09 72 65 74 75 72 6e |== 0)).X|..return|
|000058b0| 20 46 41 4c 53 45 3b 0a | 58 09 69 66 20 28 28 28 | FALSE;.|X.if (((|
|000058c0| 72 65 6d 31 20 25 20 31 | 39 29 20 3d 3d 20 30 29 |rem1 % 1|9) == 0)|
|000058d0| 20 26 26 20 28 28 72 65 | 6d 32 20 25 20 31 39 29 | && ((re|m2 % 19)|
|000058e0| 20 3d 3d 20 30 29 29 0a | 58 09 09 72 65 74 75 72 | == 0)).|X..retur|
|000058f0| 6e 20 46 41 4c 53 45 3b | 0a 58 09 69 66 20 28 28 |n FALSE;|.X.if ((|
|00005900| 28 72 65 6d 31 20 25 20 | 32 33 29 20 3d 3d 20 30 |(rem1 % |23) == 0|
|00005910| 29 20 26 26 20 28 28 72 | 65 6d 32 20 25 20 32 33 |) && ((r|em2 % 23|
|00005920| 29 20 3d 3d 20 30 29 29 | 0a 58 09 09 72 65 74 75 |) == 0))|.X..retu|
|00005930| 72 6e 20 46 41 4c 53 45 | 3b 0a 58 09 2f 2a 0a 58 |rn FALSE|;.X./*.X|
|00005940| 09 20 2a 20 59 75 6b 2c | 20 77 65 20 6d 75 73 74 |. * Yuk,| we must|
|00005950| 20 61 63 74 75 61 6c 6c | 79 20 63 6f 6d 70 75 74 | actuall|y comput|
|00005960| 65 20 74 68 65 20 67 63 | 64 20 74 6f 20 6b 6e 6f |e the gc|d to kno|
|00005970| 77 20 74 68 65 20 61 6e | 73 77 65 72 0a 58 09 20 |w the an|swer.X. |
|00005980| 2a 2f 0a 58 09 7a 67 63 | 64 28 7a 31 2c 20 7a 32 |*/.X.zgc|d(z1, z2|
|00005990| 2c 20 26 72 65 6d 29 3b | 0a 58 09 72 65 73 75 6c |, &rem);|.X.resul|
|000059a0| 74 20 3d 20 69 73 75 6e | 69 74 28 72 65 6d 29 3b |t = isun|it(rem);|
|000059b0| 0a 58 09 66 72 65 65 68 | 28 72 65 6d 2e 76 29 3b |.X.freeh|(rem.v);|
|000059c0| 0a 58 09 72 65 74 75 72 | 6e 20 72 65 73 75 6c 74 |.X.retur|n result|
|000059d0| 3b 0a 58 7d 0a 58 0a 58 | 0a 58 2f 2a 0a 58 20 2a |;.X}.X.X|.X/*.X *|
|000059e0| 20 43 6f 6d 70 75 74 65 | 20 74 68 65 20 6c 6f 67 | Compute| the log|
|000059f0| 20 6f 66 20 6f 6e 65 20 | 6e 75 6d 62 65 72 20 62 | of one |number b|
|00005a00| 61 73 65 20 61 6e 6f 74 | 68 65 72 2c 20 74 6f 20 |ase anot|her, to |
|00005a10| 74 68 65 20 63 6c 6f 73 | 65 73 74 20 69 6e 74 65 |the clos|est inte|
|00005a20| 67 65 72 2e 0a 58 20 2a | 20 54 68 69 73 20 69 73 |ger..X *| This is|
|00005a30| 20 74 68 65 20 6c 61 72 | 67 65 73 74 20 69 6e 74 | the lar|gest int|
|00005a40| 65 67 65 72 20 77 68 69 | 63 68 20 77 68 65 6e 20 |eger whi|ch when |
|00005a50| 74 68 65 20 73 65 63 6f | 6e 64 20 6e 75 6d 62 65 |the seco|nd numbe|
|00005a60| 72 20 69 73 20 72 61 69 | 73 65 64 20 74 6f 20 69 |r is rai|sed to i|
|00005a70| 74 2c 0a 58 20 2a 20 74 | 68 65 20 72 65 73 75 6c |t,.X * t|he resul|
|00005a80| 74 69 6e 67 20 76 61 6c | 75 65 20 69 73 20 6c 65 |ting val|ue is le|
|00005a90| 73 73 20 74 68 61 6e 20 | 6f 72 20 65 71 75 61 6c |ss than |or equal|
|00005aa0| 20 74 6f 20 74 68 65 20 | 66 69 72 73 74 20 6e 75 | to the |first nu|
|00005ab0| 6d 62 65 72 2e 0a 58 20 | 2a 20 45 78 61 6d 70 6c |mber..X |* Exampl|
|00005ac0| 65 3a 20 20 7a 6c 6f 67 | 28 31 32 33 34 35 36 2c |e: zlog|(123456,|
|00005ad0| 20 31 30 29 20 3d 20 35 | 2e 0a 58 20 2a 2f 0a 58 | 10) = 5|..X */.X|
|00005ae0| 6c 6f 6e 67 0a 58 7a 6c | 6f 67 28 7a 31 2c 20 7a |long.Xzl|og(z1, z|
|00005af0| 32 29 0a 58 09 5a 56 41 | 4c 55 45 20 7a 31 2c 20 |2).X.ZVA|LUE z1, |
|00005b00| 7a 32 3b 0a 58 7b 0a 58 | 09 72 65 67 69 73 74 65 |z2;.X{.X|.registe|
|00005b10| 72 20 5a 56 41 4c 55 45 | 20 2a 7a 70 3b 09 09 2f |r ZVALUE| *zp;../|
|00005b20| 2a 20 63 75 72 72 65 6e | 74 20 73 71 75 61 72 65 |* curren|t square|
|00005b30| 20 2a 2f 0a 58 09 6c 6f | 6e 67 20 70 6f 77 65 72 | */.X.lo|ng power|
|00005b40| 3b 09 09 09 2f 2a 20 63 | 75 72 72 65 6e 74 20 70 |;.../* c|urrent p|
|00005b50| 6f 77 65 72 20 2a 2f 0a | 58 09 6c 6f 6e 67 20 77 |ower */.|X.long w|
|00005b60| 6f 72 74 68 3b 09 09 09 | 2f 2a 20 77 6f 72 74 68 |orth;...|/* worth|
|00005b70| 20 6f 66 20 63 75 72 72 | 65 6e 74 20 73 71 75 61 | of curr|ent squa|
|00005b80| 72 65 20 2a 2f 0a 58 09 | 5a 56 41 4c 55 45 20 76 |re */.X.|ZVALUE v|
|00005b90| 61 6c 3b 09 09 09 2f 2a | 20 63 75 72 72 65 6e 74 |al;.../*| current|
|00005ba0| 20 76 61 6c 75 65 20 6f | 66 20 70 6f 77 65 72 20 | value o|f power |
|00005bb0| 2a 2f 0a 58 09 5a 56 41 | 4c 55 45 20 74 65 6d 70 |*/.X.ZVA|LUE temp|
|00005bc0| 3b 09 09 09 2f 2a 20 74 | 65 6d 70 6f 72 61 72 79 |;.../* t|emporary|
|00005bd0| 20 2a 2f 0a 58 09 5a 56 | 41 4c 55 45 20 73 71 75 | */.X.ZV|ALUE squ|
|00005be0| 61 72 65 73 5b 33 32 5d | 3b 09 09 2f 2a 20 74 61 |ares[32]|;../* ta|
|00005bf0| 62 6c 65 20 6f 66 20 73 | 71 75 61 72 65 73 20 6f |ble of s|quares o|
|00005c00| 66 20 62 61 73 65 20 2a | 2f 0a 58 0a 58 09 2f 2a |f base *|/.X.X./*|
|00005c10| 0a 58 09 20 2a 20 4d 61 | 6b 65 20 73 75 72 65 20 |.X. * Ma|ke sure |
|00005c20| 74 68 61 74 20 74 68 65 | 20 6e 75 6d 62 65 72 73 |that the| numbers|
|00005c30| 20 61 72 65 20 6e 6f 6e | 7a 65 72 6f 20 61 6e 64 | are non|zero and|
|00005c40| 20 74 68 65 20 62 61 73 | 65 20 69 73 20 67 72 65 | the bas|e is gre|
|00005c50| 61 74 65 72 20 74 68 61 | 6e 20 6f 6e 65 2e 0a 58 |ater tha|n one..X|
|00005c60| 09 20 2a 2f 0a 58 09 69 | 66 20 28 69 73 6e 65 67 |. */.X.i|f (isneg|
|00005c70| 28 7a 31 29 20 7c 7c 20 | 69 73 7a 65 72 6f 28 7a |(z1) || |iszero(z|
|00005c80| 31 29 20 7c 7c 20 69 73 | 6e 65 67 28 7a 32 29 20 |1) || is|neg(z2) |
|00005c90| 7c 7c 20 69 73 6c 65 6f | 6e 65 28 7a 32 29 29 0a ||| isleo|ne(z2)).|
|00005ca0| 58 09 09 65 72 72 6f 72 | 28 22 42 61 64 20 61 72 |X..error|("Bad ar|
|00005cb0| 67 75 6d 65 6e 74 73 20 | 66 6f 72 20 6c 6f 67 22 |guments |for log"|
|00005cc0| 29 3b 0a 58 09 2f 2a 0a | 58 09 20 2a 20 52 65 6a |);.X./*.|X. * Rej|
|00005cd0| 65 63 74 20 74 72 69 76 | 69 61 6c 20 63 61 73 65 |ect triv|ial case|
|00005ce0| 73 2e 0a 58 09 20 2a 2f | 0a 58 09 69 66 20 28 7a |s..X. */|.X.if (z|
|00005cf0| 31 2e 6c 65 6e 20 3c 20 | 7a 32 2e 6c 65 6e 29 0a |1.len < |z2.len).|
|00005d00| 58 09 09 72 65 74 75 72 | 6e 20 30 3b 0a 58 09 69 |X..retur|n 0;.X.i|
|00005d10| 66 20 28 28 7a 31 2e 6c | 65 6e 20 3d 3d 20 7a 32 |f ((z1.l|en == z2|
|00005d20| 2e 6c 65 6e 29 20 26 26 | 20 28 7a 31 2e 76 5b 7a |.len) &&| (z1.v[z|
|00005d30| 31 2e 6c 65 6e 2d 31 5d | 20 3c 20 7a 32 2e 76 5b |1.len-1]| < z2.v[|
|00005d40| 7a 32 2e 6c 65 6e 2d 31 | 5d 29 29 0a 58 09 09 72 |z2.len-1|])).X..r|
|00005d50| 65 74 75 72 6e 20 30 3b | 0a 58 09 70 6f 77 65 72 |eturn 0;|.X.power|
|00005d60| 20 3d 20 7a 72 65 6c 28 | 7a 31 2c 20 7a 32 29 3b | = zrel(|z1, z2);|
|00005d70| 0a 58 09 69 66 20 28 70 | 6f 77 65 72 20 3c 3d 20 |.X.if (p|ower <= |
|00005d80| 30 29 0a 58 09 09 72 65 | 74 75 72 6e 20 28 70 6f |0).X..re|turn (po|
|00005d90| 77 65 72 20 2b 20 31 29 | 3b 0a 58 09 2f 2a 0a 58 |wer + 1)|;.X./*.X|
|00005da0| 09 20 2a 20 48 61 6e 64 | 6c 65 20 61 6e 79 20 70 |. * Hand|le any p|
|00005db0| 6f 77 65 72 20 6f 66 20 | 74 77 6f 20 73 70 65 63 |ower of |two spec|
|00005dc0| 69 61 6c 2e 0a 58 09 20 | 2a 2f 0a 58 09 69 66 20 |ial..X. |*/.X.if |
|00005dd0| 28 7a 69 73 6f 6e 65 62 | 69 74 28 7a 32 29 29 0a |(zisoneb|it(z2)).|
|00005de0| 58 09 09 72 65 74 75 72 | 6e 20 28 7a 68 69 67 68 |X..retur|n (zhigh|
|00005df0| 62 69 74 28 7a 31 29 20 | 2f 20 7a 6c 6f 77 62 69 |bit(z1) |/ zlowbi|
|00005e00| 74 28 7a 32 29 29 3b 0a | 58 09 2f 2a 0a 58 09 20 |t(z2));.|X./*.X. |
|00005e10| 2a 20 48 61 6e 64 6c 65 | 20 62 61 73 65 20 31 30 |* Handle| base 10|
|00005e20| 20 73 70 65 63 69 61 6c | 0a 58 09 20 2a 2f 0a 58 | special|.X. */.X|
|00005e30| 09 69 66 20 28 28 7a 32 | 2e 6c 65 6e 20 3d 3d 20 |.if ((z2|.len == |
|00005e40| 31 29 20 26 26 20 28 2a | 7a 32 2e 76 20 3d 3d 20 |1) && (*|z2.v == |
|00005e50| 31 30 29 29 0a 58 09 09 | 72 65 74 75 72 6e 20 7a |10)).X..|return z|
|00005e60| 6c 6f 67 31 30 28 7a 31 | 29 3b 0a 58 09 2f 2a 0a |log10(z1|);.X./*.|
|00005e70| 58 09 20 2a 20 4e 6f 77 | 20 6c 6f 6f 70 20 62 79 |X. * Now| loop by|
|00005e80| 20 73 71 75 61 72 69 6e | 67 20 74 68 65 20 62 61 | squarin|g the ba|
|00005e90| 73 65 20 65 61 63 68 20 | 74 69 6d 65 2c 20 61 6e |se each |time, an|
|00005ea0| 64 20 73 65 65 20 77 68 | 65 74 68 65 72 20 6f 72 |d see wh|ether or|
|00005eb0| 0a 58 09 20 2a 20 6e 6f | 74 20 65 61 63 68 20 73 |.X. * no|t each s|
|00005ec0| 75 63 63 65 73 73 69 76 | 65 20 73 71 75 61 72 65 |uccessiv|e square|
|00005ed0| 20 69 73 20 73 74 69 6c | 6c 20 73 6d 61 6c 6c 65 | is stil|l smalle|
|00005ee0| 72 20 74 68 61 6e 20 74 | 68 65 20 6e 75 6d 62 65 |r than t|he numbe|
|00005ef0| 72 2e 0a 58 09 20 2a 2f | 0a 58 09 77 6f 72 74 68 |r..X. */|.X.worth|
|00005f00| 20 3d 20 31 3b 0a 58 09 | 7a 70 20 3d 20 26 73 71 | = 1;.X.|zp = &sq|
|00005f10| 75 61 72 65 73 5b 30 5d | 3b 0a 58 09 2a 7a 70 20 |uares[0]|;.X.*zp |
|00005f20| 3d 20 7a 32 3b 0a 58 09 | 77 68 69 6c 65 20 28 28 |= z2;.X.|while ((|
|00005f30| 28 7a 70 2d 3e 6c 65 6e | 20 2a 20 32 29 20 2d 20 |(zp->len| * 2) - |
|00005f40| 31 29 20 3c 3d 20 7a 31 | 2e 6c 65 6e 29 20 7b 09 |1) <= z1|.len) {.|
|00005f50| 2f 2a 20 77 68 69 6c 65 | 20 73 71 75 61 72 65 20 |/* while| square |
|00005f60| 6e 6f 74 20 74 6f 6f 20 | 6c 61 72 67 65 20 2a 2f |not too |large */|
|00005f70| 0a 58 09 09 7a 73 71 75 | 61 72 65 28 2a 7a 70 2c |.X..zsqu|are(*zp,|
|00005f80| 20 7a 70 20 2b 20 31 29 | 3b 0a 58 09 09 7a 70 2b | zp + 1)|;.X..zp+|
|00005f90| 2b 3b 0a 58 09 09 77 6f | 72 74 68 20 2a 3d 20 32 |+;.X..wo|rth *= 2|
|00005fa0| 3b 0a 58 09 7d 0a 58 09 | 2f 2a 0a 58 09 20 2a 20 |;.X.}.X.|/*.X. * |
|00005fb0| 4e 6f 77 20 62 61 63 6b | 20 64 6f 77 6e 20 74 68 |Now back| down th|
|00005fc0| 65 20 73 71 75 61 72 65 | 73 2c 20 61 6e 64 20 6d |e square|s, and m|
|00005fd0| 75 6c 74 69 70 6c 79 20 | 74 68 65 6d 20 74 6f 67 |ultiply |them tog|
|00005fe0| 65 74 68 65 72 20 74 6f | 20 73 65 65 0a 58 09 20 |ether to| see.X. |
|00005ff0| 2a 20 65 78 61 63 74 6c | 79 20 68 6f 77 20 6d 61 |* exactl|y how ma|
|00006000| 6e 79 20 74 69 6d 65 73 | 20 74 68 65 20 62 61 73 |ny times| the bas|
|00006010| 65 20 63 61 6e 20 62 65 | 20 72 61 69 73 65 64 20 |e can be| raised |
|00006020| 62 79 2e 0a 58 09 20 2a | 2f 0a 58 09 76 61 6c 20 |by..X. *|/.X.val |
|00006030| 3d 20 5f 6f 6e 65 5f 3b | 0a 58 09 70 6f 77 65 72 |= _one_;|.X.power|
|00006040| 20 3d 20 30 3b 0a 58 09 | 66 6f 72 20 28 3b 20 7a | = 0;.X.|for (; z|
|00006050| 70 20 3e 3d 20 73 71 75 | 61 72 65 73 3b 20 7a 70 |p >= squ|ares; zp|
|00006060| 2d 2d 2c 20 77 6f 72 74 | 68 20 2f 3d 20 32 29 20 |--, wort|h /= 2) |
|00006070| 7b 0a 58 09 09 69 66 20 | 28 28 76 61 6c 2e 6c 65 |{.X..if |((val.le|
|00006080| 6e 20 2b 20 7a 70 2d 3e | 6c 65 6e 20 2d 20 31 29 |n + zp->|len - 1)|
|00006090| 20 3c 3d 20 7a 31 2e 6c | 65 6e 29 20 7b 0a 58 09 | <= z1.l|en) {.X.|
|000060a0| 09 09 7a 6d 75 6c 28 76 | 61 6c 2c 20 2a 7a 70 2c |..zmul(v|al, *zp,|
|000060b0| 20 26 74 65 6d 70 29 3b | 0a 58 09 09 09 69 66 20 | &temp);|.X...if |
|000060c0| 28 7a 72 65 6c 28 7a 31 | 2c 20 74 65 6d 70 29 20 |(zrel(z1|, temp) |
|000060d0| 3e 3d 20 30 29 20 7b 0a | 58 09 09 09 09 66 72 65 |>= 0) {.|X....fre|
|000060e0| 65 68 28 76 61 6c 2e 76 | 29 3b 0a 58 09 09 09 09 |eh(val.v|);.X....|
|000060f0| 76 61 6c 20 3d 20 74 65 | 6d 70 3b 0a 58 09 09 09 |val = te|mp;.X...|
|00006100| 09 70 6f 77 65 72 20 2b | 3d 20 77 6f 72 74 68 3b |.power +|= worth;|
|00006110| 0a 58 09 09 09 7d 20 65 | 6c 73 65 0a 58 09 09 09 |.X...} e|lse.X...|
|00006120| 09 66 72 65 65 68 28 74 | 65 6d 70 2e 76 29 3b 0a |.freeh(t|emp.v);.|
|00006130| 58 09 09 7d 0a 58 09 09 | 69 66 20 28 7a 70 20 21 |X..}.X..|if (zp !|
|00006140| 3d 20 73 71 75 61 72 65 | 73 29 0a 58 09 09 09 66 |= square|s).X...f|
|00006150| 72 65 65 68 28 7a 70 2d | 3e 76 29 3b 0a 58 09 7d |reeh(zp-|>v);.X.}|
|00006160| 0a 58 09 72 65 74 75 72 | 6e 20 70 6f 77 65 72 3b |.X.retur|n power;|
|00006170| 0a 58 7d 0a 58 0a 58 0a | 58 2f 2a 0a 58 20 2a 20 |.X}.X.X.|X/*.X * |
|00006180| 52 65 74 75 72 6e 20 74 | 68 65 20 69 6e 74 65 67 |Return t|he integ|
|00006190| 72 61 6c 20 6c 6f 67 20 | 62 61 73 65 20 31 30 20 |ral log |base 10 |
|000061a0| 6f 66 20 61 20 6e 75 6d | 62 65 72 2e 0a 58 20 2a |of a num|ber..X *|
|000061b0| 2f 0a 58 6c 6f 6e 67 0a | 58 7a 6c 6f 67 31 30 28 |/.Xlong.|Xzlog10(|
|000061c0| 7a 29 0a 58 09 5a 56 41 | 4c 55 45 20 7a 3b 0a 58 |z).X.ZVA|LUE z;.X|
|000061d0| 7b 0a 58 09 72 65 67 69 | 73 74 65 72 20 5a 56 41 |{.X.regi|ster ZVA|
|000061e0| 4c 55 45 20 2a 7a 70 3b | 09 09 2f 2a 20 63 75 72 |LUE *zp;|../* cur|
|000061f0| 72 65 6e 74 20 73 71 75 | 61 72 65 20 2a 2f 0a 58 |rent squ|are */.X|
|00006200| 09 6c 6f 6e 67 20 70 6f | 77 65 72 3b 09 09 09 2f |.long po|wer;.../|
|00006210| 2a 20 63 75 72 72 65 6e | 74 20 70 6f 77 65 72 20 |* curren|t power |
|00006220| 2a 2f 0a 58 09 6c 6f 6e | 67 20 77 6f 72 74 68 3b |*/.X.lon|g worth;|
|00006230| 09 09 09 2f 2a 20 77 6f | 72 74 68 20 6f 66 20 63 |.../* wo|rth of c|
|00006240| 75 72 72 65 6e 74 20 73 | 71 75 61 72 65 20 2a 2f |urrent s|quare */|
|00006250| 0a 58 09 5a 56 41 4c 55 | 45 20 76 61 6c 3b 09 09 |.X.ZVALU|E val;..|
|00006260| 09 2f 2a 20 63 75 72 72 | 65 6e 74 20 76 61 6c 75 |./* curr|ent valu|
|00006270| 65 20 6f 66 20 70 6f 77 | 65 72 20 2a 2f 0a 58 09 |e of pow|er */.X.|
|00006280| 5a 56 41 4c 55 45 20 74 | 65 6d 70 3b 09 09 09 2f |ZVALUE t|emp;.../|
|00006290| 2a 20 74 65 6d 70 6f 72 | 61 72 79 20 2a 2f 0a 58 |* tempor|ary */.X|
|000062a0| 0a 58 09 69 66 20 28 21 | 69 73 70 6f 73 28 7a 29 |.X.if (!|ispos(z)|
|000062b0| 29 0a 58 09 09 65 72 72 | 6f 72 28 22 4e 6f 6e 2d |).X..err|or("Non-|
|000062c0| 70 6f 73 69 74 69 76 65 | 20 6e 75 6d 62 65 72 20 |positive| number |
|000062d0| 66 6f 72 20 6c 6f 67 31 | 30 22 29 3b 0a 58 09 2f |for log1|0");.X./|
|000062e0| 2a 0a 58 09 20 2a 20 4c | 6f 6f 70 20 62 79 20 73 |*.X. * L|oop by s|
|000062f0| 71 75 61 72 69 6e 67 20 | 74 68 65 20 62 61 73 65 |quaring |the base|
|00006300| 20 65 61 63 68 20 74 69 | 6d 65 2c 20 61 6e 64 20 | each ti|me, and |
|00006310| 73 65 65 20 77 68 65 74 | 68 65 72 20 6f 72 0a 58 |see whet|her or.X|
|00006320| 09 20 2a 20 6e 6f 74 20 | 65 61 63 68 20 73 75 63 |. * not |each suc|
|00006330| 63 65 73 73 69 76 65 20 | 73 71 75 61 72 65 20 69 |cessive |square i|
|00006340| 73 20 73 74 69 6c 6c 20 | 73 6d 61 6c 6c 65 72 20 |s still |smaller |
|00006350| 74 68 61 6e 20 74 68 65 | 20 6e 75 6d 62 65 72 2e |than the| number.|
|00006360| 0a 58 09 20 2a 2f 0a 58 | 09 77 6f 72 74 68 20 3d |.X. */.X|.worth =|
|00006370| 20 31 3b 0a 58 09 7a 70 | 20 3d 20 26 5f 74 65 6e | 1;.X.zp| = &_ten|
|00006380| 70 6f 77 65 72 73 5f 5b | 30 5d 3b 0a 58 09 2a 7a |powers_[|0];.X.*z|
|00006390| 70 20 3d 20 5f 74 65 6e | 5f 3b 0a 58 09 77 68 69 |p = _ten|_;.X.whi|
|000063a0| 6c 65 20 28 28 28 7a 70 | 2d 3e 6c 65 6e 20 2a 20 |le (((zp|->len * |
|000063b0| 32 29 20 2d 20 31 29 20 | 3c 3d 20 7a 2e 6c 65 6e |2) - 1) |<= z.len|
|000063c0| 29 20 7b 09 2f 2a 20 77 | 68 69 6c 65 20 73 71 75 |) {./* w|hile squ|
|000063d0| 61 72 65 20 6e 6f 74 20 | 74 6f 6f 20 6c 61 72 67 |are not |too larg|
|000063e0| 65 20 2a 2f 0a 58 09 09 | 69 66 20 28 7a 70 5b 31 |e */.X..|if (zp[1|
|000063f0| 5d 2e 6c 65 6e 20 3d 3d | 20 30 29 0a 58 09 09 09 |].len ==| 0).X...|
+--------+-------------------------+-------------------------+--------+--------+
Only 25.0 KB of data is shown above.
