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

< prev

next >

MacOS 8.1

Win98

DOS

browse contents

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 4e 65 77 73 67 72 6f 75 | 70 73 3a 20 63 6f 6d 70 |Newsgrou|ps: comp|
|00000010| 2e 73 6f 75 72 63 65 73 | 2e 75 6e 69 78 0a 46 72 |.sources|.unix.Fr|
|00000020| 6f 6d 3a 20 64 62 65 6c | 6c 40 70 64 61 63 74 2e |om: dbel|l@pdact.|
|00000030| 70 64 2e 6e 65 63 69 73 | 61 2e 6f 7a 2e 61 75 20 |pd.necis|a.oz.au |
|00000040| 28 44 61 76 69 64 20 49 | 2e 20 42 65 6c 6c 29 0a |(David I|. Bell).|
|00000050| 53 75 62 6a 65 63 74 3a | 20 76 32 36 69 30 33 37 |Subject:| v26i037|
|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 | 31 2f 32 31 0a 53 65 6e |r, Part1|1/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 33 37 |ume 26, |Issue 37|
|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 31 0a 0a 23 21 20 2f |alc/part|11..#! /|
|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 31 20 28 6f 66 20 | 32 31 29 2e 22 0a 23 20 | 11 (of |21).".# |
|000002e0| 43 6f 6e 74 65 6e 74 73 | 3a 20 20 7a 6d 75 6c 2e |Contents|:  zmul.|
|000002f0| 63 0a 23 20 57 72 61 70 | 70 65 64 20 62 79 20 64 |c.# Wrap|ped by d|
|00000300| 62 65 6c 6c 40 65 6c 6d | 20 6f 6e 20 54 75 65 20 |bell@elm| on Tue |
|00000310| 46 65 62 20 32 35 20 31 | 35 3a 32 31 3a 30 39 20 |Feb 25 1|5:21:09 |
|00000320| 31 39 39 32 0a 50 41 54 | 48 3d 2f 62 69 6e 3a 2f |1992.PAT|H=/bin:/|
|00000330| 75 73 72 2f 62 69 6e 3a | 2f 75 73 72 2f 75 63 62 |usr/bin:|/usr/ucb|
|00000340| 20 3b 20 65 78 70 6f 72 | 74 20 50 41 54 48 0a 69 | ; expor|t PATH.i|
|00000350| 66 20 74 65 73 74 20 2d | 66 20 27 7a 6d 75 6c 2e |f test -|f 'zmul.|
|00000360| 63 27 20 2d 61 20 22 24 | 7b 31 7d 22 20 21 3d 20 |c' -a "$|{1}" != |
|00000370| 22 2d 63 22 20 3b 20 74 | 68 65 6e 20 0a 20 20 65 |"-c" ; t|hen .  e|
|00000380| 63 68 6f 20 73 68 61 72 | 3a 20 57 69 6c 6c 20 6e |cho shar|: Will n|
|00000390| 6f 74 20 63 6c 6f 62 62 | 65 72 20 65 78 69 73 74 |ot clobb|er exist|
|000003a0| 69 6e 67 20 66 69 6c 65 | 20 5c 22 27 7a 6d 75 6c |ing file| \"'zmul|
|000003b0| 2e 63 27 5c 22 0a 65 6c | 73 65 0a 65 63 68 6f 20 |.c'\".el|se.echo |
|000003c0| 73 68 61 72 3a 20 45 78 | 74 72 61 63 74 69 6e 67 |shar: Ex|tracting|
|000003d0| 20 5c 22 27 7a 6d 75 6c | 2e 63 27 5c 22 20 5c 28 | \"'zmul|.c'\" \(|
|000003e0| 32 37 33 38 39 20 63 68 | 61 72 61 63 74 65 72 73 |27389 ch|aracters|
|000003f0| 5c 29 0a 73 65 64 20 22 | 73 2f 5e 58 2f 2f 22 20 |\).sed "|s/^X//" |
|00000400| 3e 27 7a 6d 75 6c 2e 63 | 27 20 3c 3c 27 45 4e 44 |>'zmul.c|' <<'END|
|00000410| 5f 4f 46 5f 46 49 4c 45 | 27 0a 58 2f 2a 0a 58 20 |_OF_FILE|'.X/*.X |
|00000420| 2a 20 43 6f 70 79 72 69 | 67 68 74 20 28 63 29 20 |* Copyri|ght (c) |
|00000430| 31 39 39 32 20 44 61 76 | 69 64 20 49 2e 20 42 65 |1992 Dav|id I. Be|
|00000440| 6c 6c 0a 58 20 2a 20 50 | 65 72 6d 69 73 73 69 6f |ll.X * P|ermissio|
|00000450| 6e 20 69 73 20 67 72 61 | 6e 74 65 64 20 74 6f 20 |n is gra|nted to |
|00000460| 75 73 65 2c 20 64 69 73 | 74 72 69 62 75 74 65 2c |use, dis|tribute,|
|00000470| 20 6f 72 20 6d 6f 64 69 | 66 79 20 74 68 69 73 20 | or modi|fy this |
|00000480| 73 6f 75 72 63 65 2c 0a | 58 20 2a 20 70 72 6f 76 |source,.|X * prov|
|00000490| 69 64 65 64 20 74 68 61 | 74 20 74 68 69 73 20 63 |ided tha|t this c|
|000004a0| 6f 70 79 72 69 67 68 74 | 20 6e 6f 74 69 63 65 20 |opyright| notice |
|000004b0| 72 65 6d 61 69 6e 73 20 | 69 6e 74 61 63 74 2e 0a |remains |intact..|
|000004c0| 58 20 2a 0a 58 20 2a 20 | 46 61 73 74 65 72 20 74 |X *.X * |Faster t|
|000004d0| 68 61 6e 20 75 73 75 61 | 6c 20 6d 75 6c 74 69 70 |han usua|l multip|
|000004e0| 6c 79 69 6e 67 20 61 6e | 64 20 73 71 75 61 72 69 |lying an|d squari|
|000004f0| 6e 67 20 72 6f 75 74 69 | 6e 65 73 2e 0a 58 20 2a |ng routi|nes..X *|
|00000500| 20 54 68 65 20 61 6c 67 | 6f 72 69 74 68 6d 20 75 | The alg|orithm u|
|00000510| 73 65 64 20 69 73 20 74 | 68 65 20 72 65 61 73 6f |sed is t|he reaso|
|00000520| 6e 61 62 6c 79 20 73 69 | 6d 70 6c 65 20 6f 6e 65 |nably si|mple one|
|00000530| 20 66 72 6f 6d 20 4b 6e | 75 74 68 2c 20 76 6f 6c | from Kn|uth, vol|
|00000540| 75 6d 65 20 32 2c 0a 58 | 20 2a 20 73 65 63 74 69 |ume 2,.X| * secti|
|00000550| 6f 6e 20 34 2e 33 2e 33 | 2e 20 20 54 68 65 73 65 |on 4.3.3|.  These|
|00000560| 20 72 65 63 75 72 73 69 | 76 65 20 72 6f 75 74 69 | recursi|ve routi|
|00000570| 6e 65 73 20 61 72 65 20 | 6f 66 20 73 70 65 65 64 |nes are |of speed|
|00000580| 20 4f 28 4e 5e 31 2e 35 | 38 35 29 0a 58 20 2a 20 | O(N^1.5|85).X * |
|00000590| 69 6e 73 74 65 61 64 20 | 6f 66 20 4f 28 4e 5e 32 |instead |of O(N^2|
|000005a0| 29 2e 20 20 54 68 65 20 | 75 73 75 61 6c 20 6d 75 |).  The |usual mu|
|000005b0| 6c 74 69 70 6c 69 63 61 | 74 69 6f 6e 20 61 6e 64 |ltiplica|tion and|
|000005c0| 20 28 61 6c 6d 6f 73 74 | 20 75 73 75 61 6c 29 20 | (almost| usual) |
|000005d0| 73 71 75 61 72 69 6e 67 | 0a 58 20 2a 20 61 6c 67 |squaring|.X * alg|
|000005e0| 6f 72 69 74 68 6d 73 20 | 61 72 65 20 75 73 65 64 |orithms |are used|
|000005f0| 20 66 6f 72 20 73 6d 61 | 6c 6c 20 6e 75 6d 62 65 | for sma|ll numbe|
|00000600| 72 73 2e 20 20 4f 6e 20 | 61 20 33 38 36 20 77 69 |rs.  On |a 386 wi|
|00000610| 74 68 20 69 74 73 20 63 | 6f 6d 70 69 6c 65 72 2c |th its c|ompiler,|
|00000620| 20 74 68 65 0a 58 20 2a | 20 74 77 6f 20 61 6c 67 | the.X *| two alg|
|00000630| 6f 72 69 74 68 6d 73 20 | 61 72 65 20 65 71 75 61 |orithms |are equa|
|00000640| 6c 20 69 6e 20 73 70 65 | 65 64 20 61 74 20 61 62 |l in spe|ed at ab|
|00000650| 6f 75 74 20 31 30 30 20 | 64 65 63 69 6d 61 6c 20 |out 100 |decimal |
|00000660| 64 69 67 69 74 73 2e 0a | 58 20 2a 2f 0a 58 0a 58 |digits..|X */.X.X|
|00000670| 23 69 6e 63 6c 75 64 65 | 20 3c 73 74 64 69 6f 2e |#include| <stdio.|
|00000680| 68 3e 0a 58 23 69 6e 63 | 6c 75 64 65 20 22 6d 61 |h>.X#inc|lude "ma|
|00000690| 74 68 2e 68 22 0a 58 0a | 58 0a 58 4c 45 4e 20 5f |th.h".X.|X.XLEN _|
|000006a0| 6d 75 6c 32 5f 20 3d 20 | 4d 55 4c 5f 41 4c 47 32 |mul2_ = |MUL_ALG2|
|000006b0| 3b 09 09 2f 2a 20 73 69 | 7a 65 20 6f 66 20 6e 75 |;../* si|ze of nu|
|000006c0| 6d 62 65 72 20 74 6f 20 | 75 73 65 20 6d 75 6c 74 |mber to |use mult|
|000006d0| 69 70 6c 79 20 61 6c 67 | 6f 72 69 74 68 6d 20 32 |iply alg|orithm 2|
|000006e0| 20 2a 2f 0a 58 4c 45 4e | 20 5f 73 71 32 5f 20 3d | */.XLEN| _sq2_ =|
|000006f0| 20 53 51 5f 41 4c 47 32 | 3b 09 09 2f 2a 20 73 69 | SQ_ALG2|;../* si|
|00000700| 7a 65 20 6f 66 20 6e 75 | 6d 62 65 72 20 74 6f 20 |ze of nu|mber to |
|00000710| 75 73 65 20 73 71 75 61 | 72 65 20 61 6c 67 6f 72 |use squa|re algor|
|00000720| 69 74 68 6d 20 32 20 2a | 2f 0a 58 0a 58 0a 58 23 |ithm 2 *|/.X.X.X#|
|00000730| 69 66 20 30 0a 58 73 74 | 61 74 69 63 20 63 68 61 |if 0.Xst|atic cha|
|00000740| 72 20 2a 61 62 6f 72 74 | 6d 73 67 20 3d 20 22 43 |r *abort|msg = "C|
|00000750| 61 6c 63 75 6c 61 74 69 | 6f 6e 20 61 62 6f 72 74 |alculati|on abort|
|00000760| 65 64 22 3b 0a 58 73 74 | 61 74 69 63 20 63 68 61 |ed";.Xst|atic cha|
|00000770| 72 20 2a 6d 65 6d 6d 73 | 67 20 3d 20 22 4e 6f 74 |r *memms|g = "Not|
|00000780| 20 65 6e 6f 75 67 68 20 | 6d 65 6d 6f 72 79 22 3b | enough |memory";|
|00000790| 0a 58 23 65 6e 64 69 66 | 0a 58 0a 58 73 74 61 74 |.X#endif|.X.Xstat|
|000007a0| 69 63 20 48 41 4c 46 20 | 2a 74 65 6d 70 62 75 66 |ic HALF |*tempbuf|
|000007b0| 3b 09 09 2f 2a 20 74 65 | 6d 70 6f 72 61 72 79 20 |;../* te|mporary |
|000007c0| 62 75 66 66 65 72 20 66 | 6f 72 20 6d 75 6c 74 69 |buffer f|or multi|
|000007d0| 70 6c 79 20 61 6e 64 20 | 73 71 75 61 72 65 20 2a |ply and |square *|
|000007e0| 2f 0a 58 0a 58 73 74 61 | 74 69 63 20 4c 45 4e 20 |/.X.Xsta|tic LEN |
|000007f0| 64 6f 6d 75 6c 20 70 72 | 6f 74 6f 28 28 48 41 4c |domul pr|oto((HAL|
|00000800| 46 20 2a 76 31 2c 20 4c | 45 4e 20 73 69 7a 65 31 |F *v1, L|EN size1|
|00000810| 2c 20 48 41 4c 46 20 2a | 76 32 2c 20 4c 45 4e 20 |, HALF *|v2, LEN |
|00000820| 73 69 7a 65 32 2c 20 48 | 41 4c 46 20 2a 61 6e 73 |size2, H|ALF *ans|
|00000830| 29 29 3b 0a 58 73 74 61 | 74 69 63 20 4c 45 4e 20 |));.Xsta|tic LEN |
|00000840| 64 6f 73 71 75 61 72 65 | 20 70 72 6f 74 6f 28 28 |dosquare| proto((|
|00000850| 48 41 4c 46 20 2a 76 70 | 2c 20 4c 45 4e 20 73 69 |HALF *vp|, LEN si|
|00000860| 7a 65 2c 20 48 41 4c 46 | 20 2a 61 6e 73 29 29 3b |ze, HALF| *ans));|
|00000870| 0a 58 0a 58 0a 58 2f 2a | 0a 58 20 2a 20 4d 75 6c |.X.X.X/*|.X * Mul|
|00000880| 74 69 70 6c 79 20 74 77 | 6f 20 6e 75 6d 62 65 72 |tiply tw|o number|
|00000890| 73 20 75 73 69 6e 67 20 | 74 68 65 20 66 6f 6c 6c |s using |the foll|
|000008a0| 6f 77 69 6e 67 20 66 6f | 72 6d 75 6c 61 20 72 65 |owing fo|rmula re|
|000008b0| 63 75 72 73 69 76 65 6c | 79 3a 0a 58 20 2a 09 28 |cursivel|y:.X *.(|
|000008c0| 41 2a 53 2b 42 29 2a 28 | 43 2a 53 2b 44 29 20 3d |A*S+B)*(|C*S+D) =|
|000008d0| 20 28 53 5e 32 2b 53 29 | 2a 41 2a 43 20 2b 20 53 | (S^2+S)|*A*C + S|
|000008e0| 2a 28 41 2d 42 29 2a 28 | 44 2d 43 29 20 2b 20 28 |*(A-B)*(|D-C) + (|
|000008f0| 53 2b 31 29 2a 42 2a 44 | 0a 58 20 2a 20 77 68 65 |S+1)*B*D|.X * whe|
|00000900| 72 65 20 53 20 69 73 20 | 61 20 70 6f 77 65 72 20 |re S is |a power |
|00000910| 6f 66 20 32 5e 31 36 2c | 20 61 6e 64 20 73 6f 20 |of 2^16,| and so |
|00000920| 6d 75 6c 74 69 70 6c 69 | 65 73 20 62 79 20 69 74 |multipli|es by it|
|00000930| 20 61 72 65 20 73 68 69 | 66 74 73 2c 20 61 6e 64 | are shi|fts, and|
|00000940| 0a 58 20 2a 20 41 2c 42 | 2c 43 2c 44 20 61 72 65 |.X * A,B|,C,D are|
|00000950| 20 74 68 65 20 6c 65 66 | 74 20 61 6e 64 20 72 69 | the lef|t and ri|
|00000960| 67 68 74 20 68 61 6c 66 | 73 20 6f 66 20 74 68 65 |ght half|s of the|
|00000970| 20 6e 75 6d 62 65 72 73 | 20 74 6f 20 62 65 20 6d | numbers| to be m|
|00000980| 75 6c 74 69 70 6c 69 65 | 64 2e 0a 58 20 2a 2f 0a |ultiplie|d..X */.|
|00000990| 58 76 6f 69 64 0a 58 7a | 6d 75 6c 28 7a 31 2c 20 |Xvoid.Xz|mul(z1, |
|000009a0| 7a 32 2c 20 72 65 73 29 | 0a 58 09 5a 56 41 4c 55 |z2, res)|.X.ZVALU|
|000009b0| 45 20 7a 31 2c 20 7a 32 | 3b 09 09 2f 2a 20 6e 75 |E z1, z2|;../* nu|
|000009c0| 6d 62 65 72 73 20 74 6f | 20 6d 75 6c 74 69 70 6c |mbers to| multipl|
|000009d0| 79 20 2a 2f 0a 58 09 5a | 56 41 4c 55 45 20 2a 72 |y */.X.Z|VALUE *r|
|000009e0| 65 73 3b 09 09 2f 2a 20 | 72 65 73 75 6c 74 20 6f |es;../* |result o|
|000009f0| 66 20 6d 75 6c 74 69 70 | 6c 69 63 61 74 69 6f 6e |f multip|lication|
|00000a00| 20 2a 2f 0a 58 7b 0a 58 | 09 4c 45 4e 20 6c 65 6e | */.X{.X|.LEN len|
|00000a10| 3b 09 09 2f 2a 20 73 69 | 7a 65 20 6f 66 20 61 72 |;../* si|ze of ar|
|00000a20| 72 61 79 20 2a 2f 0a 58 | 0a 58 09 69 66 20 28 69 |ray */.X|.X.if (i|
|00000a30| 73 7a 65 72 6f 28 7a 31 | 29 20 7c 7c 20 69 73 7a |szero(z1|) || isz|
|00000a40| 65 72 6f 28 7a 32 29 29 | 20 7b 0a 58 09 09 2a 72 |ero(z2))| {.X..*r|
|00000a50| 65 73 20 3d 20 5f 7a 65 | 72 6f 5f 3b 0a 58 09 09 |es = _ze|ro_;.X..|
|00000a60| 72 65 74 75 72 6e 3b 0a | 58 09 7d 0a 58 09 69 66 |return;.|X.}.X.if|
|00000a70| 20 28 69 73 75 6e 69 74 | 28 7a 31 29 29 20 7b 0a | (isunit|(z1)) {.|
|00000a80| 58 09 09 7a 63 6f 70 79 | 28 7a 32 2c 20 72 65 73 |X..zcopy|(z2, res|
|00000a90| 29 3b 0a 58 09 09 72 65 | 73 2d 3e 73 69 67 6e 20 |);.X..re|s->sign |
|00000aa0| 3d 20 28 7a 31 2e 73 69 | 67 6e 20 21 3d 20 7a 32 |= (z1.si|gn != z2|
|00000ab0| 2e 73 69 67 6e 29 3b 0a | 58 09 09 72 65 74 75 72 |.sign);.|X..retur|
|00000ac0| 6e 3b 0a 58 09 7d 0a 58 | 09 69 66 20 28 69 73 75 |n;.X.}.X|.if (isu|
|00000ad0| 6e 69 74 28 7a 32 29 29 | 20 7b 0a 58 09 09 7a 63 |nit(z2))| {.X..zc|
|00000ae0| 6f 70 79 28 7a 31 2c 20 | 72 65 73 29 3b 0a 58 09 |opy(z1, |res);.X.|
|00000af0| 09 72 65 73 2d 3e 73 69 | 67 6e 20 3d 20 28 7a 31 |.res->si|gn = (z1|
|00000b00| 2e 73 69 67 6e 20 21 3d | 20 7a 32 2e 73 69 67 6e |.sign !=| z2.sign|
|00000b10| 29 3b 0a 58 09 09 72 65 | 74 75 72 6e 3b 0a 58 09 |);.X..re|turn;.X.|
|00000b20| 7d 0a 58 0a 58 09 2f 2a | 0a 58 09 20 2a 20 41 6c |}.X.X./*|.X. * Al|
|00000b30| 6c 6f 63 61 74 65 20 61 | 20 74 65 6d 70 6f 72 61 |locate a| tempora|
|00000b40| 72 79 20 62 75 66 66 65 | 72 20 66 6f 72 20 74 68 |ry buffe|r for th|
|00000b50| 65 20 72 65 63 75 72 73 | 69 6f 6e 20 6c 65 76 65 |e recurs|ion leve|
|00000b60| 6c 73 20 74 6f 20 75 73 | 65 2e 0a 58 09 20 2a 20 |ls to us|e..X. * |
|00000b70| 41 6e 20 61 72 72 61 79 | 20 6e 65 65 64 73 20 74 |An array| needs t|
|00000b80| 6f 20 62 65 20 61 6c 6c | 6f 63 61 74 65 64 20 6c |o be all|ocated l|
|00000b90| 61 72 67 65 20 65 6e 6f | 75 67 68 20 66 6f 72 20 |arge eno|ugh for |
|00000ba0| 61 6c 6c 20 6f 66 20 74 | 68 65 0a 58 09 20 2a 20 |all of t|he.X. * |
|00000bb0| 74 65 6d 70 6f 72 61 72 | 79 20 72 65 73 75 6c 74 |temporar|y result|
|00000bc0| 73 20 74 6f 20 66 69 74 | 20 69 6e 2e 20 20 54 68 |s to fit| in.  Th|
|00000bd0| 69 73 20 73 69 7a 65 20 | 69 73 20 61 62 6f 75 74 |is size |is about|
|00000be0| 20 74 77 69 63 65 20 74 | 68 65 20 73 69 7a 65 0a | twice t|he size.|
|00000bf0| 58 09 20 2a 20 6f 66 20 | 74 68 65 20 6c 61 72 67 |X. * of |the larg|
|00000c00| 65 73 74 20 6f 72 69 67 | 69 6e 61 6c 20 6e 75 6d |est orig|inal num|
|00000c10| 62 65 72 2c 20 73 69 6e | 63 65 20 65 61 63 68 20 |ber, sin|ce each |
|00000c20| 72 65 63 75 72 73 69 6f | 6e 20 6c 65 76 65 6c 20 |recursio|n level |
|00000c30| 75 73 65 73 0a 58 09 20 | 2a 20 74 68 65 20 73 69 |uses.X. |* the si|
|00000c40| 7a 65 20 6f 66 20 69 74 | 73 20 67 69 76 65 6e 20 |ze of it|s given |
|00000c50| 6e 75 6d 62 65 72 2c 20 | 61 6e 64 20 77 68 6f 73 |number, |and whos|
|00000c60| 65 20 73 69 7a 65 20 69 | 73 20 31 2f 32 20 74 68 |e size i|s 1/2 th|
|00000c70| 65 20 73 69 7a 65 20 6f | 66 0a 58 09 20 2a 20 74 |e size o|f.X. * t|
|00000c80| 68 65 20 70 72 65 76 69 | 6f 75 73 20 6c 65 76 65 |he previ|ous leve|
|00000c90| 6c 2e 20 20 54 68 65 20 | 73 75 6d 20 6f 66 20 74 |l.  The |sum of t|
|00000ca0| 68 65 20 69 6e 66 69 6e | 69 74 65 20 73 65 72 69 |he infin|ite seri|
|00000cb0| 65 73 20 69 73 20 32 2e | 0a 58 09 20 2a 20 41 64 |es is 2.|.X. * Ad|
|00000cc0| 64 20 73 6f 6d 65 20 65 | 78 74 72 61 20 77 6f 72 |d some e|xtra wor|
|00000cd0| 64 73 20 62 65 63 61 75 | 73 65 20 6f 66 20 72 6f |ds becau|se of ro|
|00000ce0| 75 6e 64 69 6e 67 20 77 | 68 65 6e 20 64 69 76 69 |unding w|hen divi|
|00000cf0| 64 69 6e 67 20 62 79 20 | 32 0a 58 09 20 2a 20 61 |ding by |2.X. * a|
|00000d00| 6e 64 20 61 6c 73 6f 20 | 62 65 63 61 75 73 65 20 |nd also |because |
|00000d10| 6f 66 20 74 68 65 20 65 | 78 74 72 61 20 77 6f 72 |of the e|xtra wor|
|00000d20| 64 20 74 68 61 74 20 65 | 61 63 68 20 6d 75 6c 74 |d that e|ach mult|
|00000d30| 69 70 6c 79 20 6e 65 65 | 64 73 2e 0a 58 09 20 2a |iply nee|ds..X. *|
|00000d40| 2f 0a 58 09 6c 65 6e 20 | 3d 20 7a 31 2e 6c 65 6e |/.X.len |= z1.len|
|00000d50| 3b 0a 58 09 69 66 20 28 | 6c 65 6e 20 3c 20 7a 32 |;.X.if (|len < z2|
|00000d60| 2e 6c 65 6e 29 0a 58 09 | 09 6c 65 6e 20 3d 20 7a |.len).X.|.len = z|
|00000d70| 32 2e 6c 65 6e 3b 0a 58 | 09 6c 65 6e 20 3d 20 32 |2.len;.X|.len = 2|
|00000d80| 20 2a 20 6c 65 6e 20 2b | 20 36 34 3b 0a 58 09 74 | * len +| 64;.X.t|
|00000d90| 65 6d 70 62 75 66 20 3d | 20 7a 61 6c 6c 6f 63 74 |empbuf =| zalloct|
|00000da0| 65 6d 70 28 6c 65 6e 29 | 3b 0a 58 0a 58 09 72 65 |emp(len)|;.X.X.re|
|00000db0| 73 2d 3e 73 69 67 6e 20 | 3d 20 28 7a 31 2e 73 69 |s->sign |= (z1.si|
|00000dc0| 67 6e 20 21 3d 20 7a 32 | 2e 73 69 67 6e 29 3b 0a |gn != z2|.sign);.|
|00000dd0| 58 09 72 65 73 2d 3e 76 | 20 3d 20 61 6c 6c 6f 63 |X.res->v| = alloc|
|00000de0| 28 7a 31 2e 6c 65 6e 20 | 2b 20 7a 32 2e 6c 65 6e |(z1.len |+ z2.len|
|00000df0| 20 2b 20 31 29 3b 0a 58 | 09 72 65 73 2d 3e 6c 65 | + 1);.X|.res->le|
|00000e00| 6e 20 3d 20 64 6f 6d 75 | 6c 28 7a 31 2e 76 2c 20 |n = domu|l(z1.v, |
|00000e10| 7a 31 2e 6c 65 6e 2c 20 | 7a 32 2e 76 2c 20 7a 32 |z1.len, |z2.v, z2|
|00000e20| 2e 6c 65 6e 2c 20 72 65 | 73 2d 3e 76 29 3b 0a 58 |.len, re|s->v);.X|
|00000e30| 7d 0a 58 0a 58 0a 58 2f | 2a 0a 58 20 2a 20 52 65 |}.X.X.X/|*.X * Re|
|00000e40| 63 75 72 73 69 76 65 20 | 72 6f 75 74 69 6e 65 20 |cursive |routine |
|00000e50| 74 6f 20 6d 75 6c 74 69 | 70 6c 79 20 74 77 6f 20 |to multi|ply two |
|00000e60| 6e 75 6d 62 65 72 73 20 | 62 79 20 73 70 6c 69 74 |numbers |by split|
|00000e70| 74 69 6e 67 20 74 68 65 | 6d 20 75 70 20 69 6e 74 |ting the|m up int|
|00000e80| 6f 0a 58 20 2a 20 74 77 | 6f 20 6e 75 6d 62 65 72 |o.X * tw|o number|
|00000e90| 73 20 6f 66 20 68 61 6c | 66 20 74 68 65 20 73 69 |s of hal|f the si|
|00000ea0| 7a 65 2c 20 61 6e 64 20 | 75 73 69 6e 67 20 74 68 |ze, and |using th|
|00000eb0| 65 20 72 65 73 75 6c 74 | 73 20 6f 66 20 6d 75 6c |e result|s of mul|
|00000ec0| 74 69 70 6c 79 69 6e 67 | 20 74 68 65 0a 58 20 2a |tiplying| the.X *|
|00000ed0| 20 73 75 62 70 69 65 63 | 65 73 2e 20 20 54 68 65 | subpiec|es.  The|
|00000ee0| 20 72 65 73 75 6c 74 20 | 69 73 20 70 6c 61 63 65 | result |is place|
|00000ef0| 64 20 69 6e 20 74 68 65 | 20 69 6e 64 69 63 61 74 |d in the| indicat|
|00000f00| 65 64 20 61 72 72 61 79 | 2c 20 77 68 69 63 68 20 |ed array|, which |
|00000f10| 6d 75 73 74 20 62 65 0a | 58 20 2a 20 6c 61 72 67 |must be.|X * larg|
|00000f20| 65 20 65 6e 6f 75 67 68 | 20 66 6f 72 20 74 68 65 |e enough| for the|
|00000f30| 20 72 65 73 75 6c 74 20 | 70 6c 75 73 20 6f 6e 65 | result |plus one|
|00000f40| 20 65 78 74 72 61 20 77 | 6f 72 64 20 28 73 69 7a | extra w|ord (siz|
|00000f50| 65 31 20 2b 20 73 69 7a | 65 32 20 2b 20 31 29 2e |e1 + siz|e2 + 1).|
|00000f60| 0a 58 20 2a 20 52 65 74 | 75 72 6e 73 20 74 68 65 |.X * Ret|urns the|
|00000f70| 20 61 63 74 75 61 6c 20 | 73 69 7a 65 20 6f 66 20 | actual |size of |
|00000f80| 74 68 65 20 72 65 73 75 | 6c 74 20 77 69 74 68 20 |the resu|lt with |
|00000f90| 6c 65 61 64 69 6e 67 20 | 7a 65 72 6f 65 73 20 73 |leading |zeroes s|
|00000fa0| 74 72 69 70 70 65 64 2e | 0a 58 20 2a 20 54 68 69 |tripped.|.X * Thi|
|00000fb0| 73 20 61 6c 73 6f 20 75 | 73 65 73 20 61 20 74 65 |s also u|ses a te|
|00000fc0| 6d 70 6f 72 61 72 79 20 | 61 72 72 61 79 20 77 68 |mporary |array wh|
|00000fd0| 69 63 68 20 6d 75 73 74 | 20 62 65 20 74 77 69 63 |ich must| be twic|
|00000fe0| 65 20 61 73 20 6c 61 72 | 67 65 20 61 73 0a 58 20 |e as lar|ge as.X |
|00000ff0| 2a 20 6f 6e 65 20 6d 6f | 72 65 20 74 68 61 6e 20 |* one mo|re than |
|00001000| 74 68 65 20 73 69 7a 65 | 20 6f 66 20 74 68 65 20 |the size| of the |
|00001010| 6e 75 6d 62 65 72 20 61 | 74 20 74 68 65 20 74 6f |number a|t the to|
|00001020| 70 20 6c 65 76 65 6c 20 | 72 65 63 75 72 73 69 76 |p level |recursiv|
|00001030| 65 20 63 61 6c 6c 2e 0a | 58 20 2a 2f 0a 58 73 74 |e call..|X */.Xst|
|00001040| 61 74 69 63 20 4c 45 4e | 0a 58 64 6f 6d 75 6c 28 |atic LEN|.Xdomul(|
|00001050| 76 31 2c 20 73 69 7a 65 | 31 2c 20 76 32 2c 20 73 |v1, size|1, v2, s|
|00001060| 69 7a 65 32 2c 20 61 6e | 73 29 0a 58 09 48 41 4c |ize2, an|s).X.HAL|
|00001070| 46 20 2a 76 31 3b 09 09 | 2f 2a 20 66 69 72 73 74 |F *v1;..|/* first|
|00001080| 20 6e 75 6d 62 65 72 20 | 2a 2f 0a 58 09 4c 45 4e | number |*/.X.LEN|
|00001090| 20 73 69 7a 65 31 3b 09 | 09 2f 2a 20 73 69 7a 65 | size1;.|./* size|
|000010a0| 20 6f 66 20 66 69 72 73 | 74 20 6e 75 6d 62 65 72 | of firs|t number|
|000010b0| 20 2a 2f 0a 58 09 48 41 | 4c 46 20 2a 76 32 3b 09 | */.X.HA|LF *v2;.|
|000010c0| 09 2f 2a 20 73 65 63 6f | 6e 64 20 6e 75 6d 62 65 |./* seco|nd numbe|
|000010d0| 72 20 2a 2f 0a 58 09 4c | 45 4e 20 73 69 7a 65 32 |r */.X.L|EN size2|
|000010e0| 3b 09 09 2f 2a 20 73 69 | 7a 65 20 6f 66 20 73 65 |;../* si|ze of se|
|000010f0| 63 6f 6e 64 20 6e 75 6d | 62 65 72 20 2a 2f 0a 58 |cond num|ber */.X|
|00001100| 09 48 41 4c 46 20 2a 61 | 6e 73 3b 09 09 2f 2a 20 |.HALF *a|ns;../* |
|00001110| 6c 6f 63 61 74 69 6f 6e | 20 66 6f 72 20 72 65 73 |location| for res|
|00001120| 75 6c 74 20 2a 2f 0a 58 | 7b 0a 58 09 4c 45 4e 20 |ult */.X|{.X.LEN |
|00001130| 73 68 69 66 74 3b 09 09 | 2f 2a 20 61 6d 6f 75 6e |shift;..|/* amoun|
|00001140| 74 20 6e 75 6d 62 65 72 | 73 20 61 72 65 20 73 68 |t number|s are sh|
|00001150| 69 66 74 65 64 20 62 79 | 20 2a 2f 0a 58 09 4c 45 |ifted by| */.X.LE|
|00001160| 4e 20 73 69 7a 65 41 3b | 09 09 2f 2a 20 73 69 7a |N sizeA;|../* siz|
|00001170| 65 20 6f 66 20 6c 65 66 | 74 20 68 61 6c 66 20 6f |e of lef|t half o|
|00001180| 66 20 66 69 72 73 74 20 | 6e 75 6d 62 65 72 20 2a |f first |number *|
|00001190| 2f 0a 58 09 4c 45 4e 20 | 73 69 7a 65 42 3b 09 09 |/.X.LEN |sizeB;..|
|000011a0| 2f 2a 20 73 69 7a 65 20 | 6f 66 20 72 69 67 68 74 |/* size |of right|
|000011b0| 20 68 61 6c 66 20 6f 66 | 20 66 69 72 73 74 20 6e | half of| first n|
|000011c0| 75 6d 62 65 72 20 2a 2f | 0a 58 09 4c 45 4e 20 73 |umber */|.X.LEN s|
|000011d0| 69 7a 65 43 3b 09 09 2f | 2a 20 73 69 7a 65 20 6f |izeC;../|* size o|
|000011e0| 66 20 6c 65 66 74 20 68 | 61 6c 66 20 6f 66 20 73 |f left h|alf of s|
|000011f0| 65 63 6f 6e 64 20 6e 75 | 6d 62 65 72 20 2a 2f 0a |econd nu|mber */.|
|00001200| 58 09 4c 45 4e 20 73 69 | 7a 65 44 3b 09 09 2f 2a |X.LEN si|zeD;../*|
|00001210| 20 73 69 7a 65 20 6f 66 | 20 72 69 67 68 74 20 68 | size of| right h|
|00001220| 61 6c 66 20 6f 66 20 73 | 65 63 6f 6e 64 20 6e 75 |alf of s|econd nu|
|00001230| 6d 62 65 72 20 2a 2f 0a | 58 09 4c 45 4e 20 73 69 |mber */.|X.LEN si|
|00001240| 7a 65 41 42 3b 09 09 2f | 2a 20 73 69 7a 65 20 6f |zeAB;../|* size o|
|00001250| 66 20 73 75 62 74 72 61 | 63 74 69 6f 6e 20 6f 66 |f subtra|ction of|
|00001260| 20 41 20 61 6e 64 20 42 | 20 2a 2f 0a 58 09 4c 45 | A and B| */.X.LE|
|00001270| 4e 20 73 69 7a 65 44 43 | 3b 09 09 2f 2a 20 73 69 |N sizeDC|;../* si|
|00001280| 7a 65 20 6f 66 20 73 75 | 62 74 72 61 63 74 69 6f |ze of su|btractio|
|00001290| 6e 20 6f 66 20 44 20 61 | 6e 64 20 43 20 2a 2f 0a |n of D a|nd C */.|
|000012a0| 58 09 4c 45 4e 20 73 69 | 7a 65 41 42 44 43 3b 09 |X.LEN si|zeABDC;.|
|000012b0| 09 2f 2a 20 73 69 7a 65 | 20 6f 66 20 70 72 6f 64 |./* size| of prod|
|000012c0| 75 63 74 20 6f 66 20 61 | 62 6f 76 65 20 74 77 6f |uct of a|bove two|
|000012d0| 20 72 65 73 75 6c 74 73 | 20 2a 2f 0a 58 09 4c 45 | results| */.X.LE|
|000012e0| 4e 20 73 75 62 73 69 7a | 65 3b 09 09 2f 2a 20 73 |N subsiz|e;../* s|
|000012f0| 69 7a 65 20 6f 66 20 64 | 69 66 66 65 72 65 6e 63 |ize of d|ifferenc|
|00001300| 65 20 6f 66 20 68 61 6c | 66 73 20 2a 2f 0a 58 09 |e of hal|fs */.X.|
|00001310| 4c 45 4e 20 63 6f 70 79 | 73 69 7a 65 3b 09 09 2f |LEN copy|size;../|
|00001320| 2a 20 73 69 7a 65 20 6f | 66 20 6e 75 6d 62 65 72 |* size o|f number|
|00001330| 20 6c 65 66 74 20 74 6f | 20 63 6f 70 79 20 2a 2f | left to| copy */|
|00001340| 0a 58 09 4c 45 4e 20 73 | 69 7a 65 74 6f 74 61 6c |.X.LEN s|izetotal|
|00001350| 3b 09 09 2f 2a 20 74 6f | 74 61 6c 20 73 69 7a 65 |;../* to|tal size|
|00001360| 20 6f 66 20 70 72 6f 64 | 75 63 74 20 2a 2f 0a 58 | of prod|uct */.X|
|00001370| 09 4c 45 4e 20 6c 65 6e | 3b 09 09 2f 2a 20 74 65 |.LEN len|;../* te|
|00001380| 6d 70 6f 72 61 72 79 20 | 6c 65 6e 67 74 68 20 2a |mporary |length *|
|00001390| 2f 0a 58 09 48 41 4c 46 | 20 2a 62 61 73 65 41 3b |/.X.HALF| *baseA;|
|000013a0| 09 09 2f 2a 20 62 61 73 | 65 20 6f 66 20 6c 65 66 |../* bas|e of lef|
|000013b0| 74 20 68 61 6c 66 20 6f | 66 20 66 69 72 73 74 20 |t half o|f first |
|000013c0| 6e 75 6d 62 65 72 20 2a | 2f 0a 58 09 48 41 4c 46 |number *|/.X.HALF|
|000013d0| 20 2a 62 61 73 65 42 3b | 09 09 2f 2a 20 62 61 73 | *baseB;|../* bas|
|000013e0| 65 20 6f 66 20 72 69 67 | 68 74 20 68 61 6c 66 20 |e of rig|ht half |
|000013f0| 6f 66 20 66 69 72 73 74 | 20 6e 75 6d 62 65 72 20 |of first| number |
|00001400| 2a 2f 0a 58 09 48 41 4c | 46 20 2a 62 61 73 65 43 |*/.X.HAL|F *baseC|
|00001410| 3b 09 09 2f 2a 20 62 61 | 73 65 20 6f 66 20 6c 65 |;../* ba|se of le|
|00001420| 66 74 20 68 61 6c 66 20 | 6f 66 20 73 65 63 6f 6e |ft half |of secon|
|00001430| 64 20 6e 75 6d 62 65 72 | 20 2a 2f 0a 58 09 48 41 |d number| */.X.HA|
|00001440| 4c 46 20 2a 62 61 73 65 | 44 3b 09 09 2f 2a 20 62 |LF *base|D;../* b|
|00001450| 61 73 65 20 6f 66 20 72 | 69 67 68 74 20 68 61 6c |ase of r|ight hal|
|00001460| 66 20 6f 66 20 73 65 63 | 6f 6e 64 20 6e 75 6d 62 |f of sec|ond numb|
|00001470| 65 72 20 2a 2f 0a 58 09 | 48 41 4c 46 20 2a 62 61 |er */.X.|HALF *ba|
|00001480| 73 65 41 42 3b 09 09 2f | 2a 20 62 61 73 65 20 6f |seAB;../|* base o|
|00001490| 66 20 72 65 73 75 6c 74 | 20 6f 66 20 73 75 62 74 |f result| of subt|
|000014a0| 72 61 63 74 69 6f 6e 20 | 6f 66 20 41 20 61 6e 64 |raction |of A and|
|000014b0| 20 42 20 2a 2f 0a 58 09 | 48 41 4c 46 20 2a 62 61 | B */.X.|HALF *ba|
|000014c0| 73 65 44 43 3b 09 09 2f | 2a 20 62 61 73 65 20 6f |seDC;../|* base o|
|000014d0| 66 20 72 65 73 75 6c 74 | 20 6f 66 20 73 75 62 74 |f result| of subt|
|000014e0| 72 61 63 74 69 6f 6e 20 | 6f 66 20 44 20 61 6e 64 |raction |of D and|
|000014f0| 20 43 20 2a 2f 0a 58 09 | 48 41 4c 46 20 2a 62 61 | C */.X.|HALF *ba|
|00001500| 73 65 41 42 44 43 3b 09 | 09 2f 2a 20 62 61 73 65 |seABDC;.|./* base|
|00001510| 20 6f 66 20 70 72 6f 64 | 75 63 74 20 6f 66 20 61 | of prod|uct of a|
|00001520| 62 6f 76 65 20 74 77 6f | 20 72 65 73 75 6c 74 73 |bove two| results|
|00001530| 20 2a 2f 0a 58 09 48 41 | 4c 46 20 2a 62 61 73 65 | */.X.HA|LF *base|
|00001540| 41 43 3b 09 09 2f 2a 20 | 62 61 73 65 20 6f 66 20 |AC;../* |base of |
|00001550| 70 72 6f 64 75 63 74 20 | 6f 66 20 41 20 61 6e 64 |product |of A and|
|00001560| 20 43 20 2a 2f 0a 58 09 | 48 41 4c 46 20 2a 62 61 | C */.X.|HALF *ba|
|00001570| 73 65 42 44 3b 09 09 2f | 2a 20 62 61 73 65 20 6f |seBD;../|* base o|
|00001580| 66 20 70 72 6f 64 75 63 | 74 20 6f 66 20 42 20 61 |f produc|t of B a|
|00001590| 6e 64 20 44 20 2a 2f 0a | 58 09 46 55 4c 4c 20 63 |nd D */.|X.FULL c|
|000015a0| 61 72 72 79 3b 09 09 2f | 2a 20 63 61 72 72 79 20 |arry;../|* carry |
|000015b0| 64 69 67 69 74 20 66 6f | 72 20 73 6d 61 6c 6c 20 |digit fo|r small |
|000015c0| 6d 75 6c 74 69 70 6c 69 | 63 61 74 69 6f 6e 73 20 |multipli|cations |
|000015d0| 2a 2f 0a 58 09 46 55 4c | 4c 20 63 61 72 72 79 41 |*/.X.FUL|L carryA|
|000015e0| 43 42 44 3b 09 09 2f 2a | 20 63 61 72 72 79 20 66 |CBD;../*| carry f|
|000015f0| 72 6f 6d 20 61 64 64 69 | 74 69 6f 6e 20 6f 66 20 |rom addi|tion of |
|00001600| 41 2a 43 20 61 6e 64 20 | 42 2a 44 20 2a 2f 0a 58 |A*C and |B*D */.X|
|00001610| 09 46 55 4c 4c 20 64 69 | 67 69 74 3b 09 09 2f 2a |.FULL di|git;../*|
|00001620| 20 73 69 6e 67 6c 65 20 | 64 69 67 69 74 20 6d 75 | single |digit mu|
|00001630| 6c 74 69 70 6c 79 69 6e | 67 20 62 79 20 2a 2f 0a |ltiplyin|g by */.|
|00001640| 58 09 48 41 4c 46 20 2a | 74 65 6d 70 3b 09 09 2f |X.HALF *|temp;../|
|00001650| 2a 20 62 61 73 65 20 66 | 6f 72 20 74 65 6d 70 6f |* base f|or tempo|
|00001660| 72 61 72 79 20 63 61 6c | 63 75 6c 61 74 69 6f 6e |rary cal|culation|
|00001670| 73 20 2a 2f 0a 58 09 42 | 4f 4f 4c 20 6e 65 67 3b |s */.X.B|OOL neg;|
|00001680| 09 09 2f 2a 20 77 68 65 | 74 68 65 72 20 69 6d 74 |../* whe|ther imt|
|00001690| 65 72 6d 65 64 69 61 74 | 65 20 74 65 72 6d 20 69 |ermediat|e term i|
|000016a0| 73 20 6e 65 67 61 74 69 | 76 65 20 2a 2f 0a 58 09 |s negati|ve */.X.|
|000016b0| 72 65 67 69 73 74 65 72 | 20 48 41 4c 46 20 2a 68 |register| HALF *h|
|000016c0| 64 2c 20 2a 68 31 2c 20 | 2a 68 32 3b 09 2f 2a 20 |d, *h1, |*h2;./* |
|000016d0| 66 6f 72 20 69 6e 6e 65 | 72 20 6c 6f 6f 70 73 20 |for inne|r loops |
|000016e0| 2a 2f 0a 58 09 53 49 55 | 4e 49 4f 4e 20 73 69 76 |*/.X.SIU|NION siv|
|000016f0| 61 6c 3b 09 09 2f 2a 20 | 66 6f 72 20 61 64 64 69 |al;../* |for addi|
|00001700| 74 69 6f 6e 20 6f 66 20 | 64 69 67 69 74 73 20 2a |tion of |digits *|
|00001710| 2f 0a 58 0a 58 09 2f 2a | 0a 58 09 20 2a 20 54 72 |/.X.X./*|.X. * Tr|
|00001720| 69 6d 20 74 68 65 20 6e | 75 6d 62 65 72 73 20 6f |im the n|umbers o|
|00001730| 66 20 6c 65 61 64 69 6e | 67 20 7a 65 72 6f 65 73 |f leadin|g zeroes|
|00001740| 20 61 6e 64 20 69 6e 69 | 74 69 61 6c 69 7a 65 20 | and ini|tialize |
|00001750| 74 68 65 0a 58 09 20 2a | 20 65 73 74 69 6d 61 74 |the.X. *| estimat|
|00001760| 65 64 20 73 69 7a 65 20 | 6f 66 20 74 68 65 20 72 |ed size |of the r|
|00001770| 65 73 75 6c 74 2e 0a 58 | 09 20 2a 2f 0a 58 09 68 |esult..X|. */.X.h|
|00001780| 64 20 3d 20 26 76 31 5b | 73 69 7a 65 31 20 2d 20 |d = &v1[|size1 - |
|00001790| 31 5d 3b 0a 58 09 77 68 | 69 6c 65 20 28 28 2a 68 |1];.X.wh|ile ((*h|
|000017a0| 64 20 3d 3d 20 30 29 20 | 26 26 20 28 73 69 7a 65 |d == 0) |&& (size|
|000017b0| 31 20 3e 20 31 29 29 20 | 7b 0a 58 09 09 68 64 2d |1 > 1)) |{.X..hd-|
|000017c0| 2d 3b 0a 58 09 09 73 69 | 7a 65 31 2d 2d 3b 0a 58 |-;.X..si|ze1--;.X|
|000017d0| 09 7d 0a 58 09 68 64 20 | 3d 20 26 76 32 5b 73 69 |.}.X.hd |= &v2[si|
|000017e0| 7a 65 32 20 2d 20 31 5d | 3b 0a 58 09 77 68 69 6c |ze2 - 1]|;.X.whil|
|000017f0| 65 20 28 28 2a 68 64 20 | 3d 3d 20 30 29 20 26 26 |e ((*hd |== 0) &&|
|00001800| 20 28 73 69 7a 65 32 20 | 3e 20 31 29 29 20 7b 0a | (size2 |> 1)) {.|
|00001810| 58 09 09 68 64 2d 2d 3b | 0a 58 09 09 73 69 7a 65 |X..hd--;|.X..size|
|00001820| 32 2d 2d 3b 0a 58 09 7d | 0a 58 09 73 69 7a 65 74 |2--;.X.}|.X.sizet|
|00001830| 6f 74 61 6c 20 3d 20 73 | 69 7a 65 31 20 2b 20 73 |otal = s|ize1 + s|
|00001840| 69 7a 65 32 3b 0a 58 0a | 58 09 2f 2a 0a 58 09 20 |ize2;.X.|X./*.X. |
|00001850| 2a 20 46 69 72 73 74 20 | 63 68 65 63 6b 20 66 6f |* First |check fo|
|00001860| 72 20 7a 65 72 6f 20 61 | 6e 73 77 65 72 2e 0a 58 |r zero a|nswer..X|
|00001870| 09 20 2a 2f 0a 58 09 69 | 66 20 28 28 28 73 69 7a |. */.X.i|f (((siz|
|00001880| 65 31 20 3d 3d 20 31 29 | 20 26 26 20 28 2a 76 31 |e1 == 1)| && (*v1|
|00001890| 20 3d 3d 20 30 29 29 20 | 7c 7c 20 28 28 73 69 7a | == 0)) ||| ((siz|
|000018a0| 65 32 20 3d 3d 20 31 29 | 20 26 26 20 28 2a 76 32 |e2 == 1)| && (*v2|
|000018b0| 20 3d 3d 20 30 29 29 29 | 20 7b 0a 58 09 09 2a 61 | == 0)))| {.X..*a|
|000018c0| 6e 73 20 3d 20 30 3b 0a | 58 09 09 72 65 74 75 72 |ns = 0;.|X..retur|
|000018d0| 6e 20 31 3b 0a 58 09 7d | 0a 58 0a 58 09 2f 2a 0a |n 1;.X.}|.X.X./*.|
|000018e0| 58 09 20 2a 20 45 78 63 | 68 61 6e 67 65 20 74 68 |X. * Exc|hange th|
|000018f0| 65 20 74 77 6f 20 6e 75 | 6d 62 65 72 73 20 69 66 |e two nu|mbers if|
|00001900| 20 6e 65 63 65 73 73 61 | 72 79 20 74 6f 20 6d 61 | necessa|ry to ma|
|00001910| 6b 65 20 74 68 65 20 6e | 75 6d 62 65 72 20 6f 66 |ke the n|umber of|
|00001920| 0a 58 09 20 2a 20 64 69 | 67 69 74 73 20 6f 66 20 |.X. * di|gits of |
|00001930| 74 68 65 20 66 69 72 73 | 74 20 6e 75 6d 62 65 72 |the firs|t number|
|00001940| 20 62 65 20 67 72 65 61 | 74 65 72 20 74 68 61 6e | be grea|ter than|
|00001950| 20 6f 72 20 65 71 75 61 | 6c 20 74 6f 20 74 68 65 | or equa|l to the|
|00001960| 0a 58 09 20 2a 20 73 65 | 63 6f 6e 64 20 6e 75 6d |.X. * se|cond num|
|00001970| 62 65 72 2e 0a 58 09 20 | 2a 2f 0a 58 09 69 66 20 |ber..X. |*/.X.if |
|00001980| 28 73 69 7a 65 31 20 3c | 20 73 69 7a 65 32 29 20 |(size1 <| size2) |
|00001990| 7b 0a 58 09 09 6c 65 6e | 20 3d 20 73 69 7a 65 31 |{.X..len| = size1|
|000019a0| 3b 20 73 69 7a 65 31 20 | 3d 20 73 69 7a 65 32 3b |; size1 |= size2;|
|000019b0| 20 73 69 7a 65 32 20 3d | 20 6c 65 6e 3b 0a 58 09 | size2 =| len;.X.|
|000019c0| 09 68 64 20 3d 20 76 31 | 3b 20 76 31 20 3d 20 76 |.hd = v1|; v1 = v|
|000019d0| 32 3b 20 76 32 20 3d 20 | 68 64 3b 0a 58 09 7d 0a |2; v2 = |hd;.X.}.|
|000019e0| 58 0a 58 09 2f 2a 0a 58 | 09 20 2a 20 49 66 20 74 |X.X./*.X|. * If t|
|000019f0| 68 65 20 73 6d 61 6c 6c | 65 72 20 6e 75 6d 62 65 |he small|er numbe|
|00001a00| 72 20 68 61 73 20 6f 6e | 6c 79 20 61 20 66 65 77 |r has on|ly a few|
|00001a10| 20 64 69 67 69 74 73 2c | 20 74 68 65 6e 20 63 61 | digits,| then ca|
|00001a20| 6c 63 75 6c 61 74 65 0a | 58 09 20 2a 20 74 68 65 |lculate.|X. * the|
|00001a30| 20 72 65 73 75 6c 74 20 | 69 6e 20 74 68 65 20 6e | result |in the n|
|00001a40| 6f 72 6d 61 6c 20 6d 61 | 6e 6e 65 72 20 69 6e 20 |ormal ma|nner in |
|00001a50| 6f 72 64 65 72 20 74 6f | 20 61 76 6f 69 64 20 74 |order to| avoid t|
|00001a60| 68 65 20 6f 76 65 72 68 | 65 61 64 0a 58 09 20 2a |he overh|ead.X. *|
|00001a70| 20 6f 66 20 74 68 65 20 | 72 65 63 75 72 73 69 6f | of the |recursio|
|00001a80| 6e 20 66 6f 72 20 73 6d | 61 6c 6c 20 6e 75 6d 62 |n for sm|all numb|
|00001a90| 65 72 73 2e 20 20 54 68 | 65 20 6e 75 6d 62 65 72 |ers.  Th|e number|
|00001aa0| 20 6f 66 20 64 69 67 69 | 74 73 20 77 68 65 72 65 | of digi|ts where|
|00001ab0| 0a 58 09 20 2a 20 74 68 | 65 20 61 6c 67 6f 72 69 |.X. * th|e algori|
|00001ac0| 74 68 6d 20 63 68 61 6e | 67 65 73 20 69 73 20 73 |thm chan|ges is s|
|00001ad0| 65 74 74 61 62 6c 65 20 | 66 72 6f 6d 20 32 20 74 |ettable |from 2 t|
|00001ae0| 6f 20 6d 61 78 69 6e 74 | 2e 0a 58 09 20 2a 2f 0a |o maxint|..X. */.|
|00001af0| 58 09 69 66 20 28 73 69 | 7a 65 32 20 3c 20 5f 6d |X.if (si|ze2 < _m|
|00001b00| 75 6c 32 5f 29 20 7b 0a | 58 09 09 2f 2a 0a 58 09 |ul2_) {.|X../*.X.|
|00001b10| 09 20 2a 20 46 69 72 73 | 74 20 63 6c 65 61 72 20 |. * Firs|t clear |
|00001b20| 74 68 65 20 74 6f 70 20 | 70 61 72 74 20 6f 66 20 |the top |part of |
|00001b30| 74 68 65 20 72 65 73 75 | 6c 74 2c 20 61 6e 64 20 |the resu|lt, and |
|00001b40| 74 68 65 6e 20 6d 75 6c | 74 69 70 6c 79 0a 58 09 |then mul|tiply.X.|
|00001b50| 09 20 2a 20 62 79 20 74 | 68 65 20 6c 6f 77 65 73 |. * by t|he lowes|
|00001b60| 74 20 64 69 67 69 74 20 | 74 6f 20 67 65 74 20 74 |t digit |to get t|
|00001b70| 68 65 20 66 69 72 73 74 | 20 70 61 72 74 69 61 6c |he first| partial|
|00001b80| 20 73 75 6d 2e 20 20 4c | 61 74 65 72 0a 58 09 09 | sum.  L|ater.X..|
|00001b90| 20 2a 20 70 72 6f 64 75 | 63 74 73 20 77 69 6c 6c | * produ|cts will|
|00001ba0| 20 74 68 65 6e 20 61 64 | 64 20 69 6e 74 6f 20 74 | then ad|d into t|
|00001bb0| 68 69 73 20 72 65 73 75 | 6c 74 2e 0a 58 09 09 20 |his resu|lt..X.. |
|00001bc0| 2a 2f 0a 58 09 09 68 64 | 20 3d 20 26 61 6e 73 5b |*/.X..hd| = &ans[|
|00001bd0| 73 69 7a 65 31 5d 3b 0a | 58 09 09 6c 65 6e 20 3d |size1];.|X..len =|
|00001be0| 20 73 69 7a 65 32 3b 0a | 58 09 09 77 68 69 6c 65 | size2;.|X..while|
|00001bf0| 20 28 6c 65 6e 2d 2d 29 | 0a 58 09 09 09 2a 68 64 | (len--)|.X...*hd|
|00001c00| 2b 2b 20 3d 20 30 3b 0a | 58 0a 58 09 09 64 69 67 |++ = 0;.|X.X..dig|
|00001c10| 69 74 20 3d 20 2a 76 32 | 2b 2b 3b 0a 58 09 09 68 |it = *v2|++;.X..h|
|00001c20| 31 20 3d 20 76 31 3b 0a | 58 09 09 68 64 20 3d 20 |1 = v1;.|X..hd = |
|00001c30| 61 6e 73 3b 0a 58 09 09 | 63 61 72 72 79 20 3d 20 |ans;.X..|carry = |
|00001c40| 30 3b 0a 58 09 09 6c 65 | 6e 20 3d 20 73 69 7a 65 |0;.X..le|n = size|
|00001c50| 31 3b 0a 58 09 09 77 68 | 69 6c 65 20 28 6c 65 6e |1;.X..wh|ile (len|
|00001c60| 20 3e 3d 20 34 29 20 7b | 09 2f 2a 20 65 78 70 61 | >= 4) {|./* expa|
|00001c70| 6e 64 20 74 68 65 20 6c | 6f 6f 70 20 73 6f 6d 65 |nd the l|oop some|
|00001c80| 20 2a 2f 0a 58 09 09 09 | 6c 65 6e 20 2d 3d 20 34 | */.X...|len -= 4|
|00001c90| 3b 0a 58 09 09 09 73 69 | 76 61 6c 2e 69 76 61 6c |;.X...si|val.ival|
|00001ca0| 75 65 20 3d 20 28 28 46 | 55 4c 4c 29 20 2a 68 31 |ue = ((F|ULL) *h1|
|00001cb0| 2b 2b 29 20 2a 20 64 69 | 67 69 74 20 2b 20 63 61 |++) * di|git + ca|
|00001cc0| 72 72 79 3b 0a 58 09 09 | 09 2a 68 64 2b 2b 20 3d |rry;.X..|.*hd++ =|
|00001cd0| 20 73 69 76 61 6c 2e 73 | 69 6c 6f 77 3b 0a 58 09 | sival.s|ilow;.X.|
|00001ce0| 09 09 63 61 72 72 79 20 | 3d 20 73 69 76 61 6c 2e |..carry |= sival.|
|00001cf0| 73 69 68 69 67 68 3b 0a | 58 09 09 09 73 69 76 61 |sihigh;.|X...siva|
|00001d00| 6c 2e 69 76 61 6c 75 65 | 20 3d 20 28 28 46 55 4c |l.ivalue| = ((FUL|
|00001d10| 4c 29 20 2a 68 31 2b 2b | 29 20 2a 20 64 69 67 69 |L) *h1++|) * digi|
|00001d20| 74 20 2b 20 63 61 72 72 | 79 3b 0a 58 09 09 09 2a |t + carr|y;.X...*|
|00001d30| 68 64 2b 2b 20 3d 20 73 | 69 76 61 6c 2e 73 69 6c |hd++ = s|ival.sil|
|00001d40| 6f 77 3b 0a 58 09 09 09 | 63 61 72 72 79 20 3d 20 |ow;.X...|carry = |
|00001d50| 73 69 76 61 6c 2e 73 69 | 68 69 67 68 3b 0a 58 09 |sival.si|high;.X.|
|00001d60| 09 09 73 69 76 61 6c 2e | 69 76 61 6c 75 65 20 3d |..sival.|ivalue =|
|00001d70| 20 28 28 46 55 4c 4c 29 | 20 2a 68 31 2b 2b 29 20 | ((FULL)| *h1++) |
|00001d80| 2a 20 64 69 67 69 74 20 | 2b 20 63 61 72 72 79 3b |* digit |+ carry;|
|00001d90| 0a 58 09 09 09 2a 68 64 | 2b 2b 20 3d 20 73 69 76 |.X...*hd|++ = siv|
|00001da0| 61 6c 2e 73 69 6c 6f 77 | 3b 0a 58 09 09 09 63 61 |al.silow|;.X...ca|
|00001db0| 72 72 79 20 3d 20 73 69 | 76 61 6c 2e 73 69 68 69 |rry = si|val.sihi|
|00001dc0| 67 68 3b 0a 58 09 09 09 | 73 69 76 61 6c 2e 69 76 |gh;.X...|sival.iv|
|00001dd0| 61 6c 75 65 20 3d 20 28 | 28 46 55 4c 4c 29 20 2a |alue = (|(FULL) *|
|00001de0| 68 31 2b 2b 29 20 2a 20 | 64 69 67 69 74 20 2b 20 |h1++) * |digit + |
|00001df0| 63 61 72 72 79 3b 0a 58 | 09 09 09 2a 68 64 2b 2b |carry;.X|...*hd++|
|00001e00| 20 3d 20 73 69 76 61 6c | 2e 73 69 6c 6f 77 3b 0a | = sival|.silow;.|
|00001e10| 58 09 09 09 63 61 72 72 | 79 20 3d 20 73 69 76 61 |X...carr|y = siva|
|00001e20| 6c 2e 73 69 68 69 67 68 | 3b 0a 58 09 09 7d 0a 58 |l.sihigh|;.X..}.X|
|00001e30| 09 09 77 68 69 6c 65 20 | 28 6c 65 6e 2d 2d 29 20 |..while |(len--) |
|00001e40| 7b 0a 58 09 09 09 73 69 | 76 61 6c 2e 69 76 61 6c |{.X...si|val.ival|
|00001e50| 75 65 20 3d 20 28 28 46 | 55 4c 4c 29 20 2a 68 31 |ue = ((F|ULL) *h1|
|00001e60| 2b 2b 29 20 2a 20 64 69 | 67 69 74 20 2b 20 63 61 |++) * di|git + ca|
|00001e70| 72 72 79 3b 0a 58 09 09 | 09 2a 68 64 2b 2b 20 3d |rry;.X..|.*hd++ =|
|00001e80| 20 73 69 76 61 6c 2e 73 | 69 6c 6f 77 3b 0a 58 09 | sival.s|ilow;.X.|
|00001e90| 09 09 63 61 72 72 79 20 | 3d 20 73 69 76 61 6c 2e |..carry |= sival.|
|00001ea0| 73 69 68 69 67 68 3b 0a | 58 09 09 7d 0a 58 09 09 |sihigh;.|X..}.X..|
|00001eb0| 2a 68 64 20 3d 20 28 48 | 41 4c 46 29 63 61 72 72 |*hd = (H|ALF)carr|
|00001ec0| 79 3b 0a 58 0a 58 09 09 | 2f 2a 0a 58 09 09 20 2a |y;.X.X..|/*.X.. *|
|00001ed0| 20 4e 6f 77 20 6d 75 6c | 74 69 70 6c 79 20 62 79 | Now mul|tiply by|
|00001ee0| 20 74 68 65 20 72 65 6d | 61 69 6e 69 6e 67 20 64 | the rem|aining d|
|00001ef0| 69 67 69 74 73 20 6f 66 | 20 74 68 65 20 73 65 63 |igits of| the sec|
|00001f00| 6f 6e 64 20 6e 75 6d 62 | 65 72 2c 0a 58 09 09 20 |ond numb|er,.X.. |
|00001f10| 2a 20 61 64 64 69 6e 67 | 20 65 61 63 68 20 70 72 |* adding| each pr|
|00001f20| 6f 64 75 63 74 20 69 6e | 74 6f 20 74 68 65 20 66 |oduct in|to the f|
|00001f30| 69 6e 61 6c 20 72 65 73 | 75 6c 74 2e 0a 58 09 09 |inal res|ult..X..|
|00001f40| 20 2a 2f 0a 58 09 09 68 | 32 20 3d 20 61 6e 73 3b | */.X..h|2 = ans;|
|00001f50| 0a 58 09 09 77 68 69 6c | 65 20 28 2d 2d 73 69 7a |.X..whil|e (--siz|
|00001f60| 65 32 20 3e 20 30 29 20 | 7b 0a 58 09 09 09 64 69 |e2 > 0) |{.X...di|
|00001f70| 67 69 74 20 3d 20 2a 76 | 32 2b 2b 3b 0a 58 09 09 |git = *v|2++;.X..|
|00001f80| 09 68 31 20 3d 20 76 31 | 3b 0a 58 09 09 09 68 64 |.h1 = v1|;.X...hd|
|00001f90| 20 3d 20 2b 2b 68 32 3b | 0a 58 09 09 09 69 66 20 | = ++h2;|.X...if |
|00001fa0| 28 64 69 67 69 74 20 3d | 3d 20 30 29 0a 58 09 09 |(digit =|= 0).X..|
|00001fb0| 09 09 63 6f 6e 74 69 6e | 75 65 3b 0a 58 09 09 09 |..contin|ue;.X...|
|00001fc0| 63 61 72 72 79 20 3d 20 | 30 3b 0a 58 09 09 09 6c |carry = |0;.X...l|
|00001fd0| 65 6e 20 3d 20 73 69 7a | 65 31 3b 0a 58 09 09 09 |en = siz|e1;.X...|
|00001fe0| 77 68 69 6c 65 20 28 6c | 65 6e 20 3e 3d 20 34 29 |while (l|en >= 4)|
|00001ff0| 20 7b 09 2f 2a 20 65 78 | 70 61 6e 64 20 74 68 65 | {./* ex|pand the|
|00002000| 20 6c 6f 6f 70 20 73 6f | 6d 65 20 2a 2f 0a 58 09 | loop so|me */.X.|
|00002010| 09 09 09 6c 65 6e 20 2d | 3d 20 34 3b 0a 58 09 09 |...len -|= 4;.X..|
|00002020| 09 09 73 69 76 61 6c 2e | 69 76 61 6c 75 65 20 3d |..sival.|ivalue =|
|00002030| 20 28 28 46 55 4c 4c 29 | 20 2a 68 31 2b 2b 29 20 | ((FULL)| *h1++) |
|00002040| 2a 20 64 69 67 69 74 0a | 58 09 09 09 09 09 2b 20 |* digit.|X.....+ |
|00002050| 28 28 46 55 4c 4c 29 20 | 2a 68 64 29 20 2b 20 63 |((FULL) |*hd) + c|
|00002060| 61 72 72 79 3b 0a 58 09 | 09 09 09 2a 68 64 2b 2b |arry;.X.|...*hd++|
|00002070| 20 3d 20 73 69 76 61 6c | 2e 73 69 6c 6f 77 3b 0a | = sival|.silow;.|
|00002080| 58 09 09 09 09 63 61 72 | 72 79 20 3d 20 73 69 76 |X....car|ry = siv|
|00002090| 61 6c 2e 73 69 68 69 67 | 68 3b 0a 58 09 09 09 09 |al.sihig|h;.X....|
|000020a0| 73 69 76 61 6c 2e 69 76 | 61 6c 75 65 20 3d 20 28 |sival.iv|alue = (|
|000020b0| 28 46 55 4c 4c 29 20 2a | 68 31 2b 2b 29 20 2a 20 |(FULL) *|h1++) * |
|000020c0| 64 69 67 69 74 0a 58 09 | 09 09 09 09 2b 20 28 28 |digit.X.|....+ ((|
|000020d0| 46 55 4c 4c 29 20 2a 68 | 64 29 20 2b 20 63 61 72 |FULL) *h|d) + car|
|000020e0| 72 79 3b 0a 58 09 09 09 | 09 2a 68 64 2b 2b 20 3d |ry;.X...|.*hd++ =|
|000020f0| 20 73 69 76 61 6c 2e 73 | 69 6c 6f 77 3b 0a 58 09 | sival.s|ilow;.X.|
|00002100| 09 09 09 63 61 72 72 79 | 20 3d 20 73 69 76 61 6c |...carry| = sival|
|00002110| 2e 73 69 68 69 67 68 3b | 0a 58 09 09 09 09 73 69 |.sihigh;|.X....si|
|00002120| 76 61 6c 2e 69 76 61 6c | 75 65 20 3d 20 28 28 46 |val.ival|ue = ((F|
|00002130| 55 4c 4c 29 20 2a 68 31 | 2b 2b 29 20 2a 20 64 69 |ULL) *h1|++) * di|
|00002140| 67 69 74 0a 58 09 09 09 | 09 09 2b 20 28 28 46 55 |git.X...|..+ ((FU|
|00002150| 4c 4c 29 20 2a 68 64 29 | 20 2b 20 63 61 72 72 79 |LL) *hd)| + carry|
|00002160| 3b 0a 58 09 09 09 09 2a | 68 64 2b 2b 20 3d 20 73 |;.X....*|hd++ = s|
|00002170| 69 76 61 6c 2e 73 69 6c | 6f 77 3b 0a 58 09 09 09 |ival.sil|ow;.X...|
|00002180| 09 63 61 72 72 79 20 3d | 20 73 69 76 61 6c 2e 73 |.carry =| sival.s|
|00002190| 69 68 69 67 68 3b 0a 58 | 09 09 09 09 73 69 76 61 |ihigh;.X|....siva|
|000021a0| 6c 2e 69 76 61 6c 75 65 | 20 3d 20 28 28 46 55 4c |l.ivalue| = ((FUL|
|000021b0| 4c 29 20 2a 68 31 2b 2b | 29 20 2a 20 64 69 67 69 |L) *h1++|) * digi|
|000021c0| 74 0a 58 09 09 09 09 09 | 2b 20 28 28 46 55 4c 4c |t.X.....|+ ((FULL|
|000021d0| 29 20 2a 68 64 29 20 2b | 20 63 61 72 72 79 3b 0a |) *hd) +| carry;.|
|000021e0| 58 09 09 09 09 2a 68 64 | 2b 2b 20 3d 20 73 69 76 |X....*hd|++ = siv|
|000021f0| 61 6c 2e 73 69 6c 6f 77 | 3b 0a 58 09 09 09 09 63 |al.silow|;.X....c|
|00002200| 61 72 72 79 20 3d 20 73 | 69 76 61 6c 2e 73 69 68 |arry = s|ival.sih|
|00002210| 69 67 68 3b 0a 58 09 09 | 09 7d 0a 58 09 09 09 77 |igh;.X..|.}.X...w|
|00002220| 68 69 6c 65 20 28 6c 65 | 6e 2d 2d 29 20 7b 0a 58 |hile (le|n--) {.X|
|00002230| 09 09 09 09 73 69 76 61 | 6c 2e 69 76 61 6c 75 65 |....siva|l.ivalue|
|00002240| 20 3d 20 28 28 46 55 4c | 4c 29 20 2a 68 31 2b 2b | = ((FUL|L) *h1++|
|00002250| 29 20 2a 20 64 69 67 69 | 74 0a 58 09 09 09 09 09 |) * digi|t.X.....|
|00002260| 2b 20 28 28 46 55 4c 4c | 29 20 2a 68 64 29 20 2b |+ ((FULL|) *hd) +|
|00002270| 20 63 61 72 72 79 3b 0a | 58 09 09 09 09 2a 68 64 | carry;.|X....*hd|
|00002280| 2b 2b 20 3d 20 73 69 76 | 61 6c 2e 73 69 6c 6f 77 |++ = siv|al.silow|
|00002290| 3b 0a 58 09 09 09 09 63 | 61 72 72 79 20 3d 20 73 |;.X....c|arry = s|
|000022a0| 69 76 61 6c 2e 73 69 68 | 69 67 68 3b 0a 58 09 09 |ival.sih|igh;.X..|
|000022b0| 09 7d 0a 58 09 09 09 77 | 68 69 6c 65 20 28 63 61 |.}.X...w|hile (ca|
|000022c0| 72 72 79 29 20 7b 0a 58 | 09 09 09 09 73 69 76 61 |rry) {.X|....siva|
|000022d0| 6c 2e 69 76 61 6c 75 65 | 20 3d 20 28 28 46 55 4c |l.ivalue| = ((FUL|
|000022e0| 4c 29 20 2a 68 64 29 20 | 2b 20 63 61 72 72 79 3b |L) *hd) |+ carry;|
|000022f0| 0a 58 09 09 09 09 2a 68 | 64 2b 2b 20 3d 20 73 69 |.X....*h|d++ = si|
|00002300| 76 61 6c 2e 73 69 6c 6f | 77 3b 0a 58 09 09 09 09 |val.silo|w;.X....|
|00002310| 63 61 72 72 79 20 3d 20 | 73 69 76 61 6c 2e 73 69 |carry = |sival.si|
|00002320| 68 69 67 68 3b 0a 58 09 | 09 09 7d 0a 58 09 09 7d |high;.X.|..}.X..}|
|00002330| 0a 58 0a 58 09 09 2f 2a | 0a 58 09 09 20 2a 20 4e |.X.X../*|.X.. * N|
|00002340| 6f 77 20 72 65 74 75 72 | 6e 20 74 68 65 20 74 72 |ow retur|n the tr|
|00002350| 75 65 20 73 69 7a 65 20 | 6f 66 20 74 68 65 20 6e |ue size |of the n|
|00002360| 75 6d 62 65 72 2e 0a 58 | 09 09 20 2a 2f 0a 58 09 |umber..X|.. */.X.|
|00002370| 09 6c 65 6e 20 3d 20 73 | 69 7a 65 74 6f 74 61 6c |.len = s|izetotal|
|00002380| 3b 0a 58 09 09 68 64 20 | 3d 20 26 61 6e 73 5b 6c |;.X..hd |= &ans[l|
|00002390| 65 6e 20 2d 20 31 5d 3b | 0a 58 09 09 77 68 69 6c |en - 1];|.X..whil|
|000023a0| 65 20 28 28 2a 68 64 20 | 3d 3d 20 30 29 20 26 26 |e ((*hd |== 0) &&|
|000023b0| 20 28 6c 65 6e 20 3e 20 | 31 29 29 20 7b 0a 58 09 | (len > |1)) {.X.|
|000023c0| 09 09 68 64 2d 2d 3b 0a | 58 09 09 09 6c 65 6e 2d |..hd--;.|X...len-|
|000023d0| 2d 3b 0a 58 09 09 7d 0a | 58 09 09 72 65 74 75 72 |-;.X..}.|X..retur|
|000023e0| 6e 20 6c 65 6e 3b 0a 58 | 09 7d 0a 58 0a 58 09 2f |n len;.X|.}.X.X./|
|000023f0| 2a 0a 58 09 20 2a 20 4e | 65 65 64 20 74 6f 20 6d |*.X. * N|eed to m|
|00002400| 75 6c 74 69 70 6c 79 20 | 62 79 20 61 20 6c 61 72 |ultiply |by a lar|
|00002410| 67 65 20 6e 75 6d 62 65 | 72 2e 0a 58 09 20 2a 20 |ge numbe|r..X. * |
|00002420| 41 6c 6c 6f 63 61 74 65 | 20 74 65 6d 70 6f 72 61 |Allocate| tempora|
|00002430| 72 79 20 73 70 61 63 65 | 20 66 6f 72 20 63 61 6c |ry space| for cal|
|00002440| 63 75 6c 61 74 69 6f 6e | 73 2c 20 61 6e 64 20 63 |culation|s, and c|
|00002450| 61 6c 63 75 6c 61 74 65 | 20 74 68 65 0a 58 09 20 |alculate| the.X. |
|00002460| 2a 20 76 61 6c 75 65 20 | 66 6f 72 20 74 68 65 20 |* value |for the |
|00002470| 73 68 69 66 74 2e 20 20 | 54 68 65 20 73 68 69 66 |shift.  |The shif|
|00002480| 74 20 76 61 6c 75 65 20 | 69 73 20 31 2f 32 20 74 |t value |is 1/2 t|
|00002490| 68 65 20 73 69 7a 65 20 | 6f 66 20 74 68 65 0a 58 |he size |of the.X|
|000024a0| 09 20 2a 20 6c 61 72 67 | 65 72 20 28 66 69 72 73 |. * larg|er (firs|
|000024b0| 74 29 20 6e 75 6d 62 65 | 72 20 28 72 6f 75 6e 64 |t) numbe|r (round|
|000024c0| 65 64 20 75 70 29 2e 20 | 20 54 68 65 20 61 6d 6f |ed up). | The amo|
|000024d0| 75 6e 74 20 6f 66 20 74 | 65 6d 70 6f 72 61 72 79 |unt of t|emporary|
|000024e0| 0a 58 09 20 2a 20 73 70 | 61 63 65 20 6e 65 65 64 |.X. * sp|ace need|
|000024f0| 65 64 20 69 73 20 74 77 | 69 63 65 20 74 68 65 20 |ed is tw|ice the |
|00002500| 73 69 7a 65 20 6f 66 20 | 74 68 65 20 73 68 69 66 |size of |the shif|
|00002510| 74 2c 20 70 6c 75 73 20 | 6f 6e 65 20 6d 6f 72 65 |t, plus |one more|
|00002520| 20 77 6f 72 64 0a 58 09 | 20 2a 20 66 6f 72 20 74 | word.X.| * for t|
|00002530| 68 65 20 6d 75 6c 74 69 | 70 6c 79 20 74 6f 20 75 |he multi|ply to u|
|00002540| 73 65 2e 0a 58 09 20 2a | 2f 0a 58 09 73 68 69 66 |se..X. *|/.X.shif|
|00002550| 74 20 3d 20 28 73 69 7a | 65 31 20 2b 20 31 29 20 |t = (siz|e1 + 1) |
|00002560| 2f 20 32 3b 0a 58 09 74 | 65 6d 70 20 3d 20 74 65 |/ 2;.X.t|emp = te|
|00002570| 6d 70 62 75 66 3b 0a 58 | 09 74 65 6d 70 62 75 66 |mpbuf;.X|.tempbuf|
|00002580| 20 2b 3d 20 28 32 20 2a | 20 73 68 69 66 74 29 20 | += (2 *| shift) |
|00002590| 2b 20 31 3b 0a 58 0a 58 | 09 2f 2a 0a 58 09 20 2a |+ 1;.X.X|./*.X. *|
|000025a0| 20 44 65 74 65 72 6d 69 | 6e 65 20 74 68 65 20 73 | Determi|ne the s|
|000025b0| 69 7a 65 73 20 61 6e 64 | 20 6c 6f 63 61 74 69 6f |izes and| locatio|
|000025c0| 6e 73 20 6f 66 20 61 6c | 6c 20 74 68 65 20 6e 75 |ns of al|l the nu|
|000025d0| 6d 62 65 72 73 2e 0a 58 | 09 20 2a 20 54 68 65 20 |mbers..X|. * The |
|000025e0| 76 61 6c 75 65 20 6f 66 | 20 73 69 7a 65 43 20 63 |value of| sizeC c|
|000025f0| 61 6e 20 62 65 20 6e 65 | 67 61 74 69 76 65 2c 20 |an be ne|gative, |
|00002600| 61 6e 64 20 74 68 69 73 | 20 69 73 20 63 68 65 63 |and this| is chec|
|00002610| 6b 65 64 20 6c 61 74 65 | 72 2e 0a 58 09 20 2a 20 |ked late|r..X. * |
|00002620| 54 68 65 20 76 61 6c 75 | 65 20 6f 66 20 73 69 7a |The valu|e of siz|
|00002630| 65 44 20 69 73 20 6c 69 | 6d 69 74 65 64 20 62 79 |eD is li|mited by|
|00002640| 20 74 68 65 20 66 75 6c | 6c 20 73 69 7a 65 20 6f | the ful|l size o|
|00002650| 66 20 74 68 65 20 6e 75 | 6d 62 65 72 2e 0a 58 09 |f the nu|mber..X.|
|00002660| 20 2a 2f 0a 58 09 62 61 | 73 65 41 20 3d 20 76 31 | */.X.ba|seA = v1|
|00002670| 20 2b 20 73 68 69 66 74 | 3b 0a 58 09 62 61 73 65 | + shift|;.X.base|
|00002680| 42 20 3d 20 76 31 3b 0a | 58 09 62 61 73 65 43 20 |B = v1;.|X.baseC |
|00002690| 3d 20 76 32 20 2b 20 73 | 68 69 66 74 3b 0a 58 09 |= v2 + s|hift;.X.|
|000026a0| 62 61 73 65 44 20 3d 20 | 76 32 3b 0a 58 09 62 61 |baseD = |v2;.X.ba|
|000026b0| 73 65 41 42 20 3d 20 61 | 6e 73 3b 0a 58 09 62 61 |seAB = a|ns;.X.ba|
|000026c0| 73 65 44 43 20 3d 20 61 | 6e 73 20 2b 20 73 68 69 |seDC = a|ns + shi|
|000026d0| 66 74 3b 0a 58 09 62 61 | 73 65 41 43 20 3d 20 61 |ft;.X.ba|seAC = a|
|000026e0| 6e 73 20 2b 20 73 68 69 | 66 74 20 2a 20 32 3b 0a |ns + shi|ft * 2;.|
|000026f0| 58 09 62 61 73 65 42 44 | 20 3d 20 61 6e 73 3b 0a |X.baseBD| = ans;.|
|00002700| 58 0a 58 09 73 69 7a 65 | 41 20 3d 20 73 69 7a 65 |X.X.size|A = size|
|00002710| 31 20 2d 20 73 68 69 66 | 74 3b 0a 58 09 73 69 7a |1 - shif|t;.X.siz|
|00002720| 65 43 20 3d 20 73 69 7a | 65 32 20 2d 20 73 68 69 |eC = siz|e2 - shi|
|00002730| 66 74 3b 0a 58 0a 58 09 | 73 69 7a 65 42 20 3d 20 |ft;.X.X.|sizeB = |
|00002740| 73 68 69 66 74 3b 0a 58 | 09 68 64 20 3d 20 26 62 |shift;.X|.hd = &b|
|00002750| 61 73 65 42 5b 73 68 69 | 66 74 20 2d 20 31 5d 3b |aseB[shi|ft - 1];|
|00002760| 0a 58 09 77 68 69 6c 65 | 20 28 28 2a 68 64 20 3d |.X.while| ((*hd =|
|00002770| 3d 20 30 29 20 26 26 20 | 28 73 69 7a 65 42 20 3e |= 0) && |(sizeB >|
|00002780| 20 31 29 29 20 7b 0a 58 | 09 09 68 64 2d 2d 3b 0a | 1)) {.X|..hd--;.|
|00002790| 58 09 09 73 69 7a 65 42 | 2d 2d 3b 0a 58 09 7d 0a |X..sizeB|--;.X.}.|
|000027a0| 58 0a 58 09 73 69 7a 65 | 44 20 3d 20 73 68 69 66 |X.X.size|D = shif|
|000027b0| 74 3b 0a 58 09 69 66 20 | 28 73 69 7a 65 44 20 3e |t;.X.if |(sizeD >|
|000027c0| 20 73 69 7a 65 32 29 0a | 58 09 09 73 69 7a 65 44 | size2).|X..sizeD|
|000027d0| 20 3d 20 73 69 7a 65 32 | 3b 0a 58 09 68 64 20 3d | = size2|;.X.hd =|
|000027e0| 20 26 62 61 73 65 44 5b | 73 69 7a 65 44 20 2d 20 | &baseD[|sizeD - |
|000027f0| 31 5d 3b 0a 58 09 77 68 | 69 6c 65 20 28 28 2a 68 |1];.X.wh|ile ((*h|
|00002800| 64 20 3d 3d 20 30 29 20 | 26 26 20 28 73 69 7a 65 |d == 0) |&& (size|
|00002810| 44 20 3e 20 31 29 29 20 | 7b 0a 58 09 09 68 64 2d |D > 1)) |{.X..hd-|
|00002820| 2d 3b 0a 58 09 09 73 69 | 7a 65 44 2d 2d 3b 0a 58 |-;.X..si|zeD--;.X|
|00002830| 09 7d 0a 58 0a 58 09 2f | 2a 0a 58 09 20 2a 20 49 |.}.X.X./|*.X. * I|
|00002840| 66 20 74 68 65 20 73 6d | 61 6c 6c 65 72 20 6e 75 |f the sm|aller nu|
|00002850| 6d 62 65 72 20 68 61 73 | 20 61 20 68 69 67 68 20 |mber has| a high |
|00002860| 68 61 6c 66 20 6f 66 20 | 7a 65 72 6f 2c 20 74 68 |half of |zero, th|
|00002870| 65 6e 20 63 61 6c 63 75 | 6c 61 74 65 0a 58 09 20 |en calcu|late.X. |
|00002880| 2a 20 74 68 65 20 72 65 | 73 75 6c 74 20 62 79 20 |* the re|sult by |
|00002890| 62 72 65 61 6b 69 6e 67 | 20 75 70 20 74 68 65 20 |breaking| up the |
|000028a0| 66 69 72 73 74 20 6e 75 | 6d 62 65 72 20 69 6e 74 |first nu|mber int|
|000028b0| 6f 20 74 77 6f 20 6e 75 | 6d 62 65 72 73 0a 58 09 |o two nu|mbers.X.|
|000028c0| 20 2a 20 61 6e 64 20 63 | 6f 6d 62 69 6e 69 6e 67 | * and c|ombining|
|000028d0| 20 74 68 65 20 72 65 73 | 75 6c 74 73 20 75 73 69 | the res|ults usi|
|000028e0| 6e 67 20 74 68 65 20 6f | 62 76 69 6f 75 73 20 66 |ng the o|bvious f|
|000028f0| 6f 72 6d 75 6c 61 3a 0a | 58 09 20 2a 09 28 41 2a |ormula:.|X. *.(A*|
|00002900| 53 2b 42 29 20 2a 20 44 | 20 3d 20 28 41 2a 44 29 |S+B) * D| = (A*D)|
|00002910| 2a 53 20 2b 20 42 2a 44 | 0a 58 09 20 2a 2f 0a 58 |*S + B*D|.X. */.X|
|00002920| 09 69 66 20 28 73 69 7a | 65 43 20 3c 3d 20 30 29 |.if (siz|eC <= 0)|
|00002930| 20 7b 0a 58 09 09 6c 65 | 6e 20 3d 20 64 6f 6d 75 | {.X..le|n = domu|
|00002940| 6c 28 62 61 73 65 42 2c | 20 73 69 7a 65 42 2c 20 |l(baseB,| sizeB, |
|00002950| 62 61 73 65 44 2c 20 73 | 69 7a 65 44 2c 20 61 6e |baseD, s|izeD, an|
|00002960| 73 29 3b 0a 58 09 09 68 | 64 20 3d 20 26 61 6e 73 |s);.X..h|d = &ans|
|00002970| 5b 6c 65 6e 5d 3b 0a 58 | 09 09 6c 65 6e 20 3d 20 |[len];.X|..len = |
|00002980| 73 69 7a 65 74 6f 74 61 | 6c 20 2d 20 6c 65 6e 3b |sizetota|l - len;|
|00002990| 0a 58 09 09 77 68 69 6c | 65 20 28 6c 65 6e 2d 2d |.X..whil|e (len--|
|000029a0| 29 0a 58 09 09 09 2a 68 | 64 2b 2b 20 3d 20 30 3b |).X...*h|d++ = 0;|
|000029b0| 0a 58 0a 58 09 09 2f 2a | 0a 58 09 09 20 2a 20 41 |.X.X../*|.X.. * A|
|000029c0| 64 64 20 74 68 65 20 73 | 65 63 6f 6e 64 20 6e 75 |dd the s|econd nu|
|000029d0| 6d 62 65 72 20 69 6e 74 | 6f 20 74 68 65 20 66 69 |mber int|o the fi|
|000029e0| 72 73 74 20 6e 75 6d 62 | 65 72 2c 20 73 68 69 66 |rst numb|er, shif|
|000029f0| 74 65 64 0a 58 09 09 20 | 2a 20 6f 76 65 72 20 61 |ted.X.. |* over a|
|00002a00| 74 20 74 68 65 20 63 6f | 72 72 65 63 74 20 70 6f |t the co|rrect po|
|00002a10| 73 69 74 69 6f 6e 2e 0a | 58 09 09 20 2a 2f 0a 58 |sition..|X.. */.X|
|00002a20| 09 09 6c 65 6e 20 3d 20 | 64 6f 6d 75 6c 28 62 61 |..len = |domul(ba|
|00002a30| 73 65 41 2c 20 73 69 7a | 65 41 2c 20 62 61 73 65 |seA, siz|eA, base|
|00002a40| 44 2c 20 73 69 7a 65 44 | 2c 20 74 65 6d 70 29 3b |D, sizeD|, temp);|
|00002a50| 0a 58 09 09 68 31 20 3d | 20 74 65 6d 70 3b 0a 58 |.X..h1 =| temp;.X|
|00002a60| 09 09 68 64 20 3d 20 61 | 6e 73 20 2b 20 73 68 69 |..hd = a|ns + shi|
|00002a70| 66 74 3b 0a 58 09 09 63 | 61 72 72 79 20 3d 20 30 |ft;.X..c|arry = 0|
|00002a80| 3b 0a 58 09 09 77 68 69 | 6c 65 20 28 6c 65 6e 2d |;.X..whi|le (len-|
|00002a90| 2d 29 20 7b 0a 58 09 09 | 09 73 69 76 61 6c 2e 69 |-) {.X..|.sival.i|
|00002aa0| 76 61 6c 75 65 20 3d 20 | 28 28 46 55 4c 4c 29 20 |value = |((FULL) |
|00002ab0| 2a 68 31 2b 2b 29 20 2b | 20 28 28 46 55 4c 4c 29 |*h1++) +| ((FULL)|
|00002ac0| 20 2a 68 64 29 20 2b 20 | 63 61 72 72 79 3b 0a 58 | *hd) + |carry;.X|
|00002ad0| 09 09 09 2a 68 64 2b 2b | 20 3d 20 73 69 76 61 6c |...*hd++| = sival|
|00002ae0| 2e 73 69 6c 6f 77 3b 0a | 58 09 09 09 63 61 72 72 |.silow;.|X...carr|
|00002af0| 79 20 3d 20 73 69 76 61 | 6c 2e 73 69 68 69 67 68 |y = siva|l.sihigh|
|00002b00| 3b 0a 58 09 09 7d 0a 58 | 09 09 77 68 69 6c 65 20 |;.X..}.X|..while |
|00002b10| 28 63 61 72 72 79 29 20 | 7b 0a 58 09 09 09 73 69 |(carry) |{.X...si|
|00002b20| 76 61 6c 2e 69 76 61 6c | 75 65 20 3d 20 28 28 46 |val.ival|ue = ((F|
|00002b30| 55 4c 4c 29 20 2a 68 64 | 29 20 2b 20 63 61 72 72 |ULL) *hd|) + carr|
|00002b40| 79 3b 0a 58 09 09 09 2a | 68 64 2b 2b 20 3d 20 73 |y;.X...*|hd++ = s|
|00002b50| 69 76 61 6c 2e 73 69 6c | 6f 77 3b 0a 58 09 09 09 |ival.sil|ow;.X...|
|00002b60| 63 61 72 72 79 20 3d 20 | 73 69 76 61 6c 2e 73 69 |carry = |sival.si|
|00002b70| 68 69 67 68 3b 0a 58 09 | 09 7d 0a 58 0a 58 09 09 |high;.X.|.}.X.X..|
|00002b80| 2f 2a 0a 58 09 09 20 2a | 20 44 65 74 65 72 6d 69 |/*.X.. *| Determi|
|00002b90| 6e 65 20 74 68 65 20 66 | 69 6e 61 6c 20 73 69 7a |ne the f|inal siz|
|00002ba0| 65 20 6f 66 20 74 68 65 | 20 6e 75 6d 62 65 72 20 |e of the| number |
|00002bb0| 61 6e 64 20 72 65 74 75 | 72 6e 20 69 74 2e 0a 58 |and retu|rn it..X|
|00002bc0| 09 09 20 2a 2f 0a 58 09 | 09 6c 65 6e 20 3d 20 73 |.. */.X.|.len = s|
|00002bd0| 69 7a 65 74 6f 74 61 6c | 3b 0a 58 09 09 68 64 20 |izetotal|;.X..hd |
|00002be0| 3d 20 26 61 6e 73 5b 6c | 65 6e 20 2d 20 31 5d 3b |= &ans[l|en - 1];|
|00002bf0| 0a 58 09 09 77 68 69 6c | 65 20 28 28 2a 68 64 20 |.X..whil|e ((*hd |
|00002c00| 3d 3d 20 30 29 20 26 26 | 20 28 6c 65 6e 20 3e 20 |== 0) &&| (len > |
|00002c10| 31 29 29 20 7b 0a 58 09 | 09 09 68 64 2d 2d 3b 0a |1)) {.X.|..hd--;.|
|00002c20| 58 09 09 09 6c 65 6e 2d | 2d 3b 0a 58 09 09 7d 0a |X...len-|-;.X..}.|
|00002c30| 58 09 09 74 65 6d 70 62 | 75 66 20 3d 20 74 65 6d |X..tempb|uf = tem|
|00002c40| 70 3b 0a 58 09 09 72 65 | 74 75 72 6e 20 6c 65 6e |p;.X..re|turn len|
|00002c50| 3b 0a 58 09 7d 0a 58 0a | 58 09 2f 2a 0a 58 09 20 |;.X.}.X.|X./*.X. |
|00002c60| 2a 20 4e 6f 77 20 77 65 | 20 6b 6e 6f 77 20 74 68 |* Now we| know th|
|00002c70| 61 74 20 74 68 65 20 68 | 69 67 68 20 68 61 6c 66 |at the h|igh half|
|00002c80| 73 20 6f 66 20 74 68 65 | 20 6e 75 6d 62 65 72 73 |s of the| numbers|
|00002c90| 20 61 72 65 20 6e 6f 6e | 7a 65 72 6f 2c 0a 58 09 | are non|zero,.X.|
|00002ca0| 20 2a 20 73 6f 20 77 65 | 20 63 61 6e 20 75 73 65 | * so we| can use|
|00002cb0| 20 74 68 65 20 63 6f 6d | 70 6c 65 74 65 20 66 6f | the com|plete fo|
|00002cc0| 72 6d 75 6c 61 2e 0a 58 | 09 20 2a 09 28 41 2a 53 |rmula..X|. *.(A*S|
|00002cd0| 2b 42 29 2a 28 43 2a 53 | 2b 44 29 20 3d 20 28 53 |+B)*(C*S|+D) = (S|
|00002ce0| 5e 32 2b 53 29 2a 41 2a | 43 20 2b 20 53 2a 28 41 |^2+S)*A*|C + S*(A|
|00002cf0| 2d 42 29 2a 28 44 2d 43 | 29 20 2b 20 28 53 2b 31 |-B)*(D-C|) + (S+1|
|00002d00| 29 2a 42 2a 44 2e 0a 58 | 09 20 2a 20 54 68 65 20 |)*B*D..X|. * The |
|00002d10| 73 74 65 70 73 20 61 72 | 65 20 64 6f 6e 65 20 69 |steps ar|e done i|
|00002d20| 6e 20 74 68 65 20 66 6f | 6c 6c 6f 77 69 6e 67 20 |n the fo|llowing |
|00002d30| 6f 72 64 65 72 3a 0a 58 | 09 20 2a 09 41 2d 42 0a |order:.X|. *.A-B.|
|00002d40| 58 09 20 2a 09 44 2d 43 | 0a 58 09 20 2a 09 28 41 |X. *.D-C|.X. *.(A|
|00002d50| 2d 42 29 2a 28 44 2d 43 | 29 0a 58 09 20 2a 09 53 |-B)*(D-C|).X. *.S|
|00002d60| 5e 32 2a 41 2a 43 20 2b | 20 42 2a 44 0a 58 09 20 |^2*A*C +| B*D.X. |
|00002d70| 2a 09 28 53 5e 32 2b 53 | 29 2a 41 2a 43 20 2b 20 |*.(S^2+S|)*A*C + |
|00002d80| 28 53 2b 31 29 2a 42 2a | 44 09 09 09 09 28 2a 29 |(S+1)*B*|D....(*)|
|00002d90| 0a 58 09 20 2a 09 28 53 | 5e 32 2b 53 29 2a 41 2a |.X. *.(S|^2+S)*A*|
|00002da0| 43 20 2b 20 53 2a 28 41 | 2d 42 29 2a 28 44 2d 43 |C + S*(A|-B)*(D-C|
|00002db0| 29 20 2b 20 28 53 2b 31 | 29 2a 42 2a 44 0a 58 09 |) + (S+1|)*B*D.X.|
|00002dc0| 20 2a 0a 58 09 20 2a 20 | 4e 6f 74 65 3a 20 73 74 | *.X. * |Note: st|
|00002dd0| 65 70 20 28 2a 29 20 61 | 62 6f 76 65 20 63 61 6e |ep (*) a|bove can|
|00002de0| 20 70 72 6f 64 75 63 65 | 20 61 20 72 65 73 75 6c | produce| a resul|
|00002df0| 74 20 77 68 69 63 68 20 | 69 73 20 6c 61 72 67 65 |t which |is large|
|00002e00| 72 20 74 68 61 6e 0a 58 | 09 20 2a 20 74 68 65 20 |r than.X|. * the |
|00002e10| 66 69 6e 61 6c 20 70 72 | 6f 64 75 63 74 20 77 69 |final pr|oduct wi|
|00002e20| 6c 6c 20 62 65 2c 20 61 | 6e 64 20 74 68 69 73 20 |ll be, a|nd this |
|00002e30| 69 73 20 77 68 65 72 65 | 20 74 68 65 20 65 78 74 |is where| the ext|
|00002e40| 72 61 20 77 6f 72 64 0a | 58 09 20 2a 20 6e 65 65 |ra word.|X. * nee|
|00002e50| 64 65 64 20 69 6e 20 74 | 68 65 20 70 72 6f 64 75 |ded in t|he produ|
|00002e60| 63 74 20 63 6f 6d 65 73 | 20 66 72 6f 6d 2e 20 20 |ct comes| from.  |
|00002e70| 41 66 74 65 72 20 74 68 | 65 20 66 69 6e 61 6c 20 |After th|e final |
|00002e80| 73 75 62 74 72 61 63 74 | 69 6f 6e 20 69 73 0a 58 |subtract|ion is.X|
|00002e90| 09 20 2a 20 64 6f 6e 65 | 2c 20 74 68 65 20 72 65 |. * done|, the re|
|00002ea0| 73 75 6c 74 20 66 69 74 | 73 20 69 6e 20 74 68 65 |sult fit|s in the|
|00002eb0| 20 65 78 70 65 63 74 65 | 64 20 73 69 7a 65 2e 20 | expecte|d size. |
|00002ec0| 20 55 73 69 6e 67 20 74 | 68 65 20 65 78 74 72 61 | Using t|he extra|
|00002ed0| 20 77 6f 72 64 0a 58 09 | 20 2a 20 69 73 20 65 61 | word.X.| * is ea|
|00002ee0| 73 69 65 72 20 74 68 61 | 6e 20 73 75 70 70 72 65 |sier tha|n suppre|
|00002ef0| 73 73 69 6e 67 20 74 68 | 65 20 63 61 72 72 69 65 |ssing th|e carrie|
|00002f00| 73 20 61 6e 64 20 62 6f | 72 72 6f 77 73 20 65 76 |s and bo|rrows ev|
|00002f10| 65 72 79 77 68 65 72 65 | 2e 0a 58 09 20 2a 0a 58 |erywhere|..X. *.X|
|00002f20| 09 20 2a 20 42 65 67 69 | 6e 20 62 79 20 66 6f 72 |. * Begi|n by for|
|00002f30| 6d 69 6e 67 20 74 68 65 | 20 70 72 6f 64 75 63 74 |ming the| product|
|00002f40| 20 28 41 2d 42 29 2a 28 | 44 2d 43 29 20 69 6e 74 | (A-B)*(|D-C) int|
|00002f50| 6f 20 61 20 74 65 6d 70 | 6f 72 61 72 79 0a 58 09 |o a temp|orary.X.|
|00002f60| 20 2a 20 6c 6f 63 61 74 | 69 6f 6e 20 74 68 61 74 | * locat|ion that|
|00002f70| 20 77 65 20 73 61 76 65 | 20 75 6e 74 69 6c 20 74 | we save| until t|
|00002f80| 68 65 20 66 69 6e 61 6c | 20 73 74 65 70 2e 20 20 |he final| step.  |
|00002f90| 44 6f 20 65 61 63 68 20 | 73 75 62 74 72 61 63 74 |Do each |subtract|
|00002fa0| 69 6f 6e 0a 58 09 20 2a | 20 61 74 20 70 6f 73 69 |ion.X. *| at posi|
|00002fb0| 74 69 6f 6e 73 20 30 20 | 61 6e 64 20 53 2e 20 20 |tions 0 |and S.  |
|00002fc0| 42 65 20 76 65 72 79 20 | 63 61 72 65 66 75 6c 20 |Be very |careful |
|00002fd0| 61 62 6f 75 74 20 74 68 | 65 20 72 65 6c 61 74 69 |about th|e relati|
|00002fe0| 76 65 20 73 69 7a 65 73 | 0a 58 09 20 2a 20 6f 66 |ve sizes|.X. * of|
|00002ff0| 20 74 68 65 20 6e 75 6d | 62 65 72 73 20 73 69 6e | the num|bers sin|
|00003000| 63 65 20 74 68 69 73 20 | 72 65 73 75 6c 74 20 63 |ce this |result c|
|00003010| 61 6e 20 62 65 20 6e 65 | 67 61 74 69 76 65 2e 20 |an be ne|gative. |
|00003020| 20 46 6f 72 20 74 68 65 | 20 66 69 72 73 74 0a 58 | For the| first.X|
|00003030| 09 20 2a 20 73 74 65 70 | 20 63 61 6c 63 75 6c 61 |. * step| calcula|
|00003040| 74 65 20 74 68 65 20 61 | 62 73 6f 6c 75 74 65 20 |te the a|bsolute |
|00003050| 64 69 66 66 65 72 65 6e | 63 65 20 6f 66 20 41 20 |differen|ce of A |
|00003060| 61 6e 64 20 42 20 69 6e | 74 6f 20 61 20 74 65 6d |and B in|to a tem|
|00003070| 70 6f 72 61 72 79 0a 58 | 09 20 2a 20 6c 6f 63 61 |porary.X|. * loca|
|00003080| 74 69 6f 6e 20 61 74 20 | 70 6f 73 69 74 69 6f 6e |tion at |position|
|00003090| 20 30 20 6f 66 20 74 68 | 65 20 72 65 73 75 6c 74 | 0 of th|e result|
|000030a0| 2e 20 20 4e 65 67 61 74 | 65 20 74 68 65 20 73 69 |.  Negat|e the si|
|000030b0| 67 6e 20 69 66 20 41 20 | 69 73 0a 58 09 20 2a 20 |gn if A |is.X. * |
|000030c0| 73 6d 61 6c 6c 65 72 20 | 74 68 61 6e 20 42 2e 0a |smaller |than B..|
|000030d0| 58 09 20 2a 2f 0a 58 09 | 6e 65 67 20 3d 20 46 41 |X. */.X.|neg = FA|
|000030e0| 4c 53 45 3b 0a 58 09 69 | 66 20 28 73 69 7a 65 41 |LSE;.X.i|f (sizeA|
|000030f0| 20 3d 3d 20 73 69 7a 65 | 42 29 20 7b 0a 58 09 09 | == size|B) {.X..|
|00003100| 6c 65 6e 20 3d 20 73 69 | 7a 65 41 3b 0a 58 09 09 |len = si|zeA;.X..|
|00003110| 68 31 20 3d 20 26 62 61 | 73 65 41 5b 6c 65 6e 20 |h1 = &ba|seA[len |
|00003120| 2d 20 31 5d 3b 0a 58 09 | 09 68 32 20 3d 20 26 62 |- 1];.X.|.h2 = &b|
|00003130| 61 73 65 42 5b 6c 65 6e | 20 2d 20 31 5d 3b 0a 58 |aseB[len| - 1];.X|
|00003140| 09 09 77 68 69 6c 65 20 | 28 28 6c 65 6e 20 3e 20 |..while |((len > |
|00003150| 31 29 20 26 26 20 28 2a | 68 31 20 3d 3d 20 2a 68 |1) && (*|h1 == *h|
|00003160| 32 29 29 20 7b 0a 58 09 | 09 09 6c 65 6e 2d 2d 3b |2)) {.X.|..len--;|
|00003170| 0a 58 09 09 09 68 31 2d | 2d 3b 0a 58 09 09 09 68 |.X...h1-|-;.X...h|
|00003180| 32 2d 2d 3b 0a 58 09 09 | 7d 0a 58 09 7d 0a 58 09 |2--;.X..|}.X.}.X.|
|00003190| 69 66 20 28 28 73 69 7a | 65 41 20 3e 20 73 69 7a |if ((siz|eA > siz|
|000031a0| 65 42 29 20 7c 7c 20 28 | 28 73 69 7a 65 41 20 3d |eB) || (|(sizeA =|
|000031b0| 3d 20 73 69 7a 65 42 29 | 20 26 26 20 28 2a 68 31 |= sizeB)| && (*h1|
|000031c0| 20 3e 20 2a 68 32 29 29 | 29 0a 58 09 7b 0a 58 09 | > *h2))|).X.{.X.|
|000031d0| 09 68 31 20 3d 20 62 61 | 73 65 41 3b 0a 58 09 09 |.h1 = ba|seA;.X..|
|000031e0| 68 32 20 3d 20 62 61 73 | 65 42 3b 0a 58 09 09 73 |h2 = bas|eB;.X..s|
|000031f0| 69 7a 65 41 42 20 3d 20 | 73 69 7a 65 41 3b 0a 58 |izeAB = |sizeA;.X|
|00003200| 09 09 73 75 62 73 69 7a | 65 20 3d 20 73 69 7a 65 |..subsiz|e = size|
|00003210| 42 3b 0a 58 09 7d 20 65 | 6c 73 65 20 7b 0a 58 09 |B;.X.} e|lse {.X.|
|00003220| 09 6e 65 67 20 3d 20 21 | 6e 65 67 3b 0a 58 09 09 |.neg = !|neg;.X..|
|00003230| 68 31 20 3d 20 62 61 73 | 65 42 3b 0a 58 09 09 68 |h1 = bas|eB;.X..h|
|00003240| 32 20 3d 20 62 61 73 65 | 41 3b 0a 58 09 09 73 69 |2 = base|A;.X..si|
|00003250| 7a 65 41 42 20 3d 20 73 | 69 7a 65 42 3b 0a 58 09 |zeAB = s|izeB;.X.|
|00003260| 09 73 75 62 73 69 7a 65 | 20 3d 20 73 69 7a 65 41 |.subsize| = sizeA|
|00003270| 3b 0a 58 09 7d 0a 58 09 | 63 6f 70 79 73 69 7a 65 |;.X.}.X.|copysize|
|00003280| 20 3d 20 73 69 7a 65 41 | 42 20 2d 20 73 75 62 73 | = sizeA|B - subs|
|00003290| 69 7a 65 3b 0a 58 0a 58 | 09 68 64 20 3d 20 62 61 |ize;.X.X|.hd = ba|
|000032a0| 73 65 41 42 3b 0a 58 09 | 63 61 72 72 79 20 3d 20 |seAB;.X.|carry = |
|000032b0| 30 3b 0a 58 09 77 68 69 | 6c 65 20 28 73 75 62 73 |0;.X.whi|le (subs|
|000032c0| 69 7a 65 2d 2d 29 20 7b | 0a 58 09 09 73 69 76 61 |ize--) {|.X..siva|
|000032d0| 6c 2e 69 76 61 6c 75 65 | 20 3d 20 42 41 53 45 31 |l.ivalue| = BASE1|
|000032e0| 20 2d 20 28 28 46 55 4c | 4c 29 20 2a 68 31 2b 2b | - ((FUL|L) *h1++|
|000032f0| 29 20 2b 20 28 28 46 55 | 4c 4c 29 20 2a 68 32 2b |) + ((FU|LL) *h2+|
|00003300| 2b 29 20 2b 20 63 61 72 | 72 79 3b 0a 58 09 09 2a |+) + car|ry;.X..*|
|00003310| 68 64 2b 2b 20 3d 20 42 | 41 53 45 31 20 2d 20 73 |hd++ = B|ASE1 - s|
|00003320| 69 76 61 6c 2e 73 69 6c | 6f 77 3b 0a 58 09 09 63 |ival.sil|ow;.X..c|
|00003330| 61 72 72 79 20 3d 20 73 | 69 76 61 6c 2e 73 69 68 |arry = s|ival.sih|
|00003340| 69 67 68 3b 0a 58 09 7d | 0a 58 09 77 68 69 6c 65 |igh;.X.}|.X.while|
|00003350| 20 28 63 6f 70 79 73 69 | 7a 65 2d 2d 29 20 7b 0a | (copysi|ze--) {.|
|00003360| 58 09 09 73 69 76 61 6c | 2e 69 76 61 6c 75 65 20 |X..sival|.ivalue |
|00003370| 3d 20 28 42 41 53 45 31 | 20 2d 20 28 28 46 55 4c |= (BASE1| - ((FUL|
|00003380| 4c 29 20 2a 68 31 2b 2b | 29 29 20 2b 20 63 61 72 |L) *h1++|)) + car|
|00003390| 72 79 3b 0a 58 09 09 2a | 68 64 2b 2b 20 3d 20 42 |ry;.X..*|hd++ = B|
|000033a0| 41 53 45 31 20 2d 20 73 | 69 76 61 6c 2e 73 69 6c |ASE1 - s|ival.sil|
|000033b0| 6f 77 3b 0a 58 09 09 63 | 61 72 72 79 20 3d 20 73 |ow;.X..c|arry = s|
|000033c0| 69 76 61 6c 2e 73 69 68 | 69 67 68 3b 0a 58 09 7d |ival.sih|igh;.X.}|
|000033d0| 0a 58 0a 58 09 68 64 20 | 3d 20 26 62 61 73 65 41 |.X.X.hd |= &baseA|
|000033e0| 42 5b 73 69 7a 65 41 42 | 20 2d 20 31 5d 3b 0a 58 |B[sizeAB| - 1];.X|
|000033f0| 09 77 68 69 6c 65 20 28 | 28 2a 68 64 20 3d 3d 20 |.while (|(*hd == |
|00003400| 30 29 20 26 26 20 28 73 | 69 7a 65 41 42 20 3e 20 |0) && (s|izeAB > |
|00003410| 31 29 29 20 7b 0a 58 09 | 09 68 64 2d 2d 3b 0a 58 |1)) {.X.|.hd--;.X|
|00003420| 09 09 73 69 7a 65 41 42 | 2d 2d 3b 0a 58 09 7d 0a |..sizeAB|--;.X.}.|
|00003430| 58 0a 58 09 2f 2a 0a 58 | 09 20 2a 20 54 68 69 73 |X.X./*.X|. * This|
|00003440| 20 63 6f 6d 70 6c 65 74 | 65 73 20 74 68 65 20 63 | complet|es the c|
|00003450| 61 6c 63 75 6c 61 74 69 | 6f 6e 20 6f 66 20 61 62 |alculati|on of ab|
|00003460| 73 28 41 2d 42 29 2e 20 | 20 46 6f 72 20 74 68 65 |s(A-B). | For the|
|00003470| 20 6e 65 78 74 20 73 74 | 65 70 0a 58 09 20 2a 20 | next st|ep.X. * |
|00003480| 63 61 6c 63 75 6c 61 74 | 65 20 74 68 65 20 61 62 |calculat|e the ab|
|00003490| 73 6f 6c 75 74 65 20 64 | 69 66 66 65 72 65 6e 63 |solute d|ifferenc|
|000034a0| 65 20 6f 66 20 44 20 61 | 6e 64 20 43 20 69 6e 74 |e of D a|nd C int|
|000034b0| 6f 20 61 20 74 65 6d 70 | 6f 72 61 72 79 0a 58 09 |o a temp|orary.X.|
|000034c0| 20 2a 20 6c 6f 63 61 74 | 69 6f 6e 20 61 74 20 70 | * locat|ion at p|
|000034d0| 6f 73 69 74 69 6f 6e 20 | 53 20 6f 66 20 74 68 65 |osition |S of the|
|000034e0| 20 72 65 73 75 6c 74 2e | 20 20 4e 65 67 61 74 65 | result.|  Negate|
|000034f0| 20 74 68 65 20 73 69 67 | 6e 20 69 66 20 43 20 69 | the sig|n if C i|
|00003500| 73 0a 58 09 20 2a 20 6c | 61 72 67 65 72 20 74 68 |s.X. * l|arger th|
|00003510| 61 6e 20 44 2e 0a 58 09 | 20 2a 2f 0a 58 09 69 66 |an D..X.| */.X.if|
|00003520| 20 28 73 69 7a 65 43 20 | 3d 3d 20 73 69 7a 65 44 | (sizeC |== sizeD|
|00003530| 29 20 7b 0a 58 09 09 6c | 65 6e 20 3d 20 73 69 7a |) {.X..l|en = siz|
|00003540| 65 43 3b 0a 58 09 09 68 | 31 20 3d 20 26 62 61 73 |eC;.X..h|1 = &bas|
|00003550| 65 43 5b 6c 65 6e 20 2d | 20 31 5d 3b 0a 58 09 09 |eC[len -| 1];.X..|
|00003560| 68 32 20 3d 20 26 62 61 | 73 65 44 5b 6c 65 6e 20 |h2 = &ba|seD[len |
|00003570| 2d 20 31 5d 3b 0a 58 09 | 09 77 68 69 6c 65 20 28 |- 1];.X.|.while (|
|00003580| 28 6c 65 6e 20 3e 20 31 | 29 20 26 26 20 28 2a 68 |(len > 1|) && (*h|
|00003590| 31 20 3d 3d 20 2a 68 32 | 29 29 20 7b 0a 58 09 09 |1 == *h2|)) {.X..|
|000035a0| 09 6c 65 6e 2d 2d 3b 0a | 58 09 09 09 68 31 2d 2d |.len--;.|X...h1--|
|000035b0| 3b 0a 58 09 09 09 68 32 | 2d 2d 3b 0a 58 09 09 7d |;.X...h2|--;.X..}|
|000035c0| 0a 58 09 7d 0a 58 09 69 | 66 20 28 28 73 69 7a 65 |.X.}.X.i|f ((size|
|000035d0| 43 20 3e 20 73 69 7a 65 | 44 29 20 7c 7c 20 28 28 |C > size|D) || ((|
|000035e0| 73 69 7a 65 43 20 3d 3d | 20 73 69 7a 65 44 29 20 |sizeC ==| sizeD) |
|000035f0| 26 26 20 28 2a 68 31 20 | 3e 20 2a 68 32 29 29 29 |&& (*h1 |> *h2)))|
|00003600| 0a 58 09 7b 0a 58 09 09 | 6e 65 67 20 3d 20 21 6e |.X.{.X..|neg = !n|
|00003610| 65 67 3b 0a 58 09 09 68 | 31 20 3d 20 62 61 73 65 |eg;.X..h|1 = base|
|00003620| 43 3b 0a 58 09 09 68 32 | 20 3d 20 62 61 73 65 44 |C;.X..h2| = baseD|
|00003630| 3b 0a 58 09 09 73 69 7a | 65 44 43 20 3d 20 73 69 |;.X..siz|eDC = si|
|00003640| 7a 65 43 3b 0a 58 09 09 | 73 75 62 73 69 7a 65 20 |zeC;.X..|subsize |
|00003650| 3d 20 73 69 7a 65 44 3b | 0a 58 09 7d 20 65 6c 73 |= sizeD;|.X.} els|
|00003660| 65 20 7b 0a 58 09 09 68 | 31 20 3d 20 62 61 73 65 |e {.X..h|1 = base|
|00003670| 44 3b 0a 58 09 09 68 32 | 20 3d 20 62 61 73 65 43 |D;.X..h2| = baseC|
|00003680| 3b 0a 58 09 09 73 69 7a | 65 44 43 20 3d 20 73 69 |;.X..siz|eDC = si|
|00003690| 7a 65 44 3b 0a 58 09 09 | 73 75 62 73 69 7a 65 20 |zeD;.X..|subsize |
|000036a0| 3d 20 73 69 7a 65 43 3b | 0a 58 09 7d 0a 58 09 63 |= sizeC;|.X.}.X.c|
|000036b0| 6f 70 79 73 69 7a 65 20 | 3d 20 73 69 7a 65 44 43 |opysize |= sizeDC|
|000036c0| 20 2d 20 73 75 62 73 69 | 7a 65 3b 0a 58 0a 58 09 | - subsi|ze;.X.X.|
|000036d0| 68 64 20 3d 20 62 61 73 | 65 44 43 3b 0a 58 09 63 |hd = bas|eDC;.X.c|
|000036e0| 61 72 72 79 20 3d 20 30 | 3b 0a 58 09 77 68 69 6c |arry = 0|;.X.whil|
|000036f0| 65 20 28 73 75 62 73 69 | 7a 65 2d 2d 29 20 7b 0a |e (subsi|ze--) {.|
|00003700| 58 09 09 73 69 76 61 6c | 2e 69 76 61 6c 75 65 20 |X..sival|.ivalue |
|00003710| 3d 20 42 41 53 45 31 20 | 2d 20 28 28 46 55 4c 4c |= BASE1 |- ((FULL|
|00003720| 29 20 2a 68 31 2b 2b 29 | 20 2b 20 28 28 46 55 4c |) *h1++)| + ((FUL|
|00003730| 4c 29 20 2a 68 32 2b 2b | 29 20 2b 20 63 61 72 72 |L) *h2++|) + carr|
|00003740| 79 3b 0a 58 09 09 2a 68 | 64 2b 2b 20 3d 20 42 41 |y;.X..*h|d++ = BA|
|00003750| 53 45 31 20 2d 20 73 69 | 76 61 6c 2e 73 69 6c 6f |SE1 - si|val.silo|
|00003760| 77 3b 0a 58 09 09 63 61 | 72 72 79 20 3d 20 73 69 |w;.X..ca|rry = si|
|00003770| 76 61 6c 2e 73 69 68 69 | 67 68 3b 0a 58 09 7d 0a |val.sihi|gh;.X.}.|
|00003780| 58 09 77 68 69 6c 65 20 | 28 63 6f 70 79 73 69 7a |X.while |(copysiz|
|00003790| 65 2d 2d 29 20 7b 0a 58 | 09 09 73 69 76 61 6c 2e |e--) {.X|..sival.|
|000037a0| 69 76 61 6c 75 65 20 3d | 20 28 42 41 53 45 31 20 |ivalue =| (BASE1 |
|000037b0| 2d 20 28 28 46 55 4c 4c | 29 20 2a 68 31 2b 2b 29 |- ((FULL|) *h1++)|
|000037c0| 29 20 2b 20 63 61 72 72 | 79 3b 0a 58 09 09 2a 68 |) + carr|y;.X..*h|
|000037d0| 64 2b 2b 20 3d 20 42 41 | 53 45 31 20 2d 20 73 69 |d++ = BA|SE1 - si|
|000037e0| 76 61 6c 2e 73 69 6c 6f | 77 3b 0a 58 09 09 63 61 |val.silo|w;.X..ca|
|000037f0| 72 72 79 20 3d 20 73 69 | 76 61 6c 2e 73 69 68 69 |rry = si|val.sihi|
|00003800| 67 68 3b 0a 58 09 7d 0a | 58 09 68 64 20 3d 20 26 |gh;.X.}.|X.hd = &|
|00003810| 62 61 73 65 44 43 5b 73 | 69 7a 65 44 43 20 2d 20 |baseDC[s|izeDC - |
|00003820| 31 5d 3b 0a 58 09 77 68 | 69 6c 65 20 28 28 2a 68 |1];.X.wh|ile ((*h|
|00003830| 64 20 3d 3d 20 30 29 20 | 26 26 20 28 73 69 7a 65 |d == 0) |&& (size|
|00003840| 44 43 20 3e 20 31 29 29 | 20 7b 0a 58 09 09 68 64 |DC > 1))| {.X..hd|
|00003850| 2d 2d 3b 0a 58 09 09 73 | 69 7a 65 44 43 2d 2d 3b |--;.X..s|izeDC--;|
|00003860| 0a 58 09 7d 0a 58 0a 58 | 09 2f 2a 0a 58 09 20 2a |.X.}.X.X|./*.X. *|
|00003870| 20 54 68 69 73 20 63 6f | 6d 70 6c 65 74 65 73 20 | This co|mpletes |
|00003880| 74 68 65 20 63 61 6c 63 | 75 6c 61 74 69 6f 6e 20 |the calc|ulation |
|00003890| 6f 66 20 61 62 73 28 44 | 2d 43 29 2e 20 20 4e 6f |of abs(D|-C).  No|
|000038a0| 77 20 6d 75 6c 74 69 70 | 6c 79 0a 58 09 20 2a 20 |w multip|ly.X. * |
|000038b0| 74 6f 67 65 74 68 65 72 | 20 61 62 73 28 41 2d 42 |together| abs(A-B|
|000038c0| 29 20 61 6e 64 20 61 62 | 73 28 44 2d 43 29 20 69 |) and ab|s(D-C) i|
|000038d0| 6e 74 6f 20 61 20 74 65 | 6d 70 6f 72 61 72 79 20 |nto a te|mporary |
|000038e0| 6c 6f 63 61 74 69 6f 6e | 2c 0a 58 09 20 2a 20 77 |location|,.X. * w|
|000038f0| 68 69 63 68 20 69 73 20 | 70 72 65 73 65 72 76 65 |hich is |preserve|
|00003900| 64 20 75 6e 74 69 6c 20 | 74 68 65 20 66 69 6e 61 |d until |the fina|
|00003910| 6c 20 73 74 65 70 73 2e | 0a 58 09 20 2a 2f 0a 58 |l steps.|.X. */.X|
|00003920| 09 62 61 73 65 41 42 44 | 43 20 3d 20 74 65 6d 70 |.baseABD|C = temp|
|00003930| 3b 0a 58 09 73 69 7a 65 | 41 42 44 43 20 3d 20 64 |;.X.size|ABDC = d|
|00003940| 6f 6d 75 6c 28 62 61 73 | 65 41 42 2c 20 73 69 7a |omul(bas|eAB, siz|
|00003950| 65 41 42 2c 20 62 61 73 | 65 44 43 2c 20 73 69 7a |eAB, bas|eDC, siz|
|00003960| 65 44 43 2c 20 62 61 73 | 65 41 42 44 43 29 3b 0a |eDC, bas|eABDC);.|
|00003970| 58 0a 58 09 2f 2a 0a 58 | 09 20 2a 20 4e 6f 77 20 |X.X./*.X|. * Now |
|00003980| 63 61 6c 63 75 6c 61 74 | 65 20 42 2a 44 20 61 6e |calculat|e B*D an|
|00003990| 64 20 41 2a 43 20 69 6e | 74 6f 20 6f 6e 65 20 6f |d A*C in|to one o|
|000039a0| 66 20 74 68 65 69 72 20 | 74 77 6f 20 66 69 6e 61 |f their |two fina|
|000039b0| 6c 20 6c 6f 63 61 74 69 | 6f 6e 73 2e 0a 58 09 20 |l locati|ons..X. |
|000039c0| 2a 20 4d 61 6b 65 20 73 | 75 72 65 20 74 68 65 20 |* Make s|ure the |
|000039d0| 68 69 67 68 20 6f 72 64 | 65 72 20 64 69 67 69 74 |high ord|er digit|
|000039e0| 73 20 6f 66 20 74 68 65 | 20 70 72 6f 64 75 63 74 |s of the| product|
|000039f0| 73 20 61 72 65 20 7a 65 | 72 6f 65 64 20 73 69 6e |s are ze|roed sin|
|00003a00| 63 65 0a 58 09 20 2a 20 | 74 68 69 73 20 69 6e 69 |ce.X. * |this ini|
|00003a10| 74 69 61 6c 69 7a 65 73 | 20 74 68 65 20 66 69 6e |tializes| the fin|
|00003a20| 61 6c 20 72 65 73 75 6c | 74 2e 20 20 42 65 20 63 |al resul|t.  Be c|
|00003a30| 61 72 65 66 75 6c 20 61 | 62 6f 75 74 20 74 68 69 |areful a|bout thi|
|00003a40| 73 20 7a 65 72 6f 69 6e | 67 0a 58 09 20 2a 20 73 |s zeroin|g.X. * s|
|00003a50| 69 6e 63 65 20 74 68 65 | 20 73 69 7a 65 20 6f 66 |ince the| size of|
|00003a60| 20 74 68 65 20 68 69 67 | 68 20 6f 72 64 65 72 20 | the hig|h order |
|00003a70| 77 6f 72 64 73 20 6d 69 | 67 68 74 20 62 65 20 73 |words mi|ght be s|
|00003a80| 6d 61 6c 6c 65 72 20 74 | 68 61 6e 0a 58 09 20 2a |maller t|han.X. *|
|00003a90| 20 74 68 65 20 73 68 69 | 66 74 20 73 69 7a 65 2e | the shi|ft size.|
|00003aa0| 20 20 44 6f 20 42 2a 44 | 20 66 69 72 73 74 20 73 |  Do B*D| first s|
|00003ab0| 69 6e 63 65 20 74 68 65 | 20 6d 75 6c 74 69 70 6c |ince the| multipl|
|00003ac0| 69 65 73 20 75 73 65 20 | 6f 6e 65 20 6d 6f 72 65 |ies use |one more|
|00003ad0| 0a 58 09 20 2a 20 77 6f | 72 64 20 74 68 61 6e 20 |.X. * wo|rd than |
|00003ae0| 74 68 65 20 73 69 7a 65 | 20 6f 66 20 74 68 65 20 |the size| of the |
|00003af0| 70 72 6f 64 75 63 74 2e | 20 20 41 6c 73 6f 20 7a |product.|  Also z|
|00003b00| 65 72 6f 20 74 68 65 20 | 66 69 6e 61 6c 20 65 78 |ero the |final ex|
|00003b10| 74 72 61 0a 58 09 20 2a | 20 77 6f 72 64 20 69 6e |tra.X. *| word in|
|00003b20| 20 74 68 65 20 72 65 73 | 75 6c 74 20 66 6f 72 20 | the res|ult for |
|00003b30| 70 6f 73 73 69 62 6c 65 | 20 63 61 72 72 69 65 73 |possible| carries|
|00003b40| 20 74 6f 20 75 73 65 2e | 0a 58 09 20 2a 2f 0a 58 | to use.|.X. */.X|
|00003b50| 09 6c 65 6e 20 3d 20 64 | 6f 6d 75 6c 28 62 61 73 |.len = d|omul(bas|
|00003b60| 65 42 2c 20 73 69 7a 65 | 42 2c 20 62 61 73 65 44 |eB, size|B, baseD|
|00003b70| 2c 20 73 69 7a 65 44 2c | 20 62 61 73 65 42 44 29 |, sizeD,| baseBD)|
|00003b80| 3b 0a 58 09 68 64 20 3d | 20 26 62 61 73 65 42 44 |;.X.hd =| &baseBD|
|00003b90| 5b 6c 65 6e 5d 3b 0a 58 | 09 6c 65 6e 20 3d 20 73 |[len];.X|.len = s|
|00003ba0| 68 69 66 74 20 2a 20 32 | 20 2d 20 6c 65 6e 3b 0a |hift * 2| - len;.|
|00003bb0| 58 09 77 68 69 6c 65 20 | 28 6c 65 6e 2d 2d 29 0a |X.while |(len--).|
|00003bc0| 58 09 09 2a 68 64 2b 2b | 20 3d 20 30 3b 0a 58 0a |X..*hd++| = 0;.X.|
|00003bd0| 58 09 6c 65 6e 20 3d 20 | 64 6f 6d 75 6c 28 62 61 |X.len = |domul(ba|
|00003be0| 73 65 41 2c 20 73 69 7a | 65 41 2c 20 62 61 73 65 |seA, siz|eA, base|
|00003bf0| 43 2c 20 73 69 7a 65 43 | 2c 20 62 61 73 65 41 43 |C, sizeC|, baseAC|
|00003c00| 29 3b 0a 58 09 68 64 20 | 3d 20 26 62 61 73 65 41 |);.X.hd |= &baseA|
|00003c10| 43 5b 6c 65 6e 5d 3b 0a | 58 09 6c 65 6e 20 3d 20 |C[len];.|X.len = |
|00003c20| 73 69 7a 65 74 6f 74 61 | 6c 20 2d 20 73 68 69 66 |sizetota|l - shif|
|00003c30| 74 20 2a 20 32 20 2d 20 | 6c 65 6e 20 2b 20 31 3b |t * 2 - |len + 1;|
|00003c40| 0a 58 09 77 68 69 6c 65 | 20 28 6c 65 6e 2d 2d 29 |.X.while| (len--)|
|00003c50| 0a 58 09 09 2a 68 64 2b | 2b 20 3d 20 30 3b 0a 58 |.X..*hd+|+ = 0;.X|
|00003c60| 0a 58 09 2f 2a 0a 58 09 | 20 2a 20 4e 6f 77 20 61 |.X./*.X.| * Now a|
|00003c70| 64 64 20 69 6e 20 41 2a | 43 20 61 6e 64 20 42 2a |dd in A*|C and B*|
|00003c80| 44 20 69 6e 74 6f 20 74 | 68 65 6d 73 65 6c 76 65 |D into t|hemselve|
|00003c90| 73 20 61 74 20 74 68 65 | 20 6f 74 68 65 72 20 73 |s at the| other s|
|00003ca0| 68 69 66 74 65 64 0a 58 | 09 20 2a 20 70 6f 73 69 |hifted.X|. * posi|
|00003cb0| 74 69 6f 6e 20 74 68 61 | 74 20 74 68 65 79 20 6e |tion tha|t they n|
|00003cc0| 65 65 64 2e 20 20 54 68 | 69 73 20 61 64 64 69 74 |eed.  Th|is addit|
|00003cd0| 69 6f 6e 20 69 73 20 74 | 72 69 63 6b 79 20 69 6e |ion is t|ricky in|
|00003ce0| 20 6f 72 64 65 72 20 74 | 6f 0a 58 09 20 2a 20 6d | order t|o.X. * m|
|00003cf0| 61 6b 65 20 73 75 72 65 | 20 74 68 61 74 20 74 68 |ake sure| that th|
|00003d00| 65 20 74 77 6f 20 61 64 | 64 69 74 69 6f 6e 73 20 |e two ad|ditions |
|00003d10| 63 61 6e 6e 6f 74 20 69 | 6e 74 65 72 66 65 72 65 |cannot i|nterfere|
|00003d20| 20 77 69 74 68 20 65 61 | 63 68 20 6f 74 68 65 72 | with ea|ch other|
|00003d30| 2e 0a 58 09 20 2a 20 54 | 68 65 72 65 66 6f 72 65 |..X. * T|herefore|
|00003d40| 20 77 65 20 66 69 72 73 | 74 20 61 64 64 20 69 6e | we firs|t add in|
|00003d50| 20 74 68 65 20 74 6f 70 | 20 68 61 6c 66 20 6f 66 | the top| half of|
|00003d60| 20 42 2a 44 20 61 6e 64 | 20 74 68 65 20 6c 6f 77 | B*D and| the low|
|00003d70| 65 72 20 68 61 6c 66 0a | 58 09 20 2a 20 6f 66 20 |er half.|X. * of |
|00003d80| 41 2a 43 2e 20 20 54 68 | 65 20 73 6f 75 72 63 65 |A*C.  Th|e source|
|00003d90| 73 20 61 6e 64 20 64 65 | 73 74 69 6e 61 74 69 6f |s and de|stinatio|
|00003da0| 6e 73 20 6f 66 20 74 68 | 65 73 65 20 74 77 6f 20 |ns of th|ese two |
|00003db0| 61 64 64 69 74 69 6f 6e | 73 0a 58 09 20 2a 20 6f |addition|s.X. * o|
|00003dc0| 76 65 72 6c 61 70 2c 20 | 61 6e 64 20 73 6f 20 74 |verlap, |and so t|
|00003dd0| 68 65 20 73 61 6d 65 20 | 61 6e 73 77 65 72 20 72 |he same |answer r|
|00003de0| 65 73 75 6c 74 73 20 66 | 72 6f 6d 20 74 68 65 20 |esults f|rom the |
|00003df0| 74 77 6f 20 61 64 64 69 | 74 69 6f 6e 73 2c 0a 58 |two addi|tions,.X|
|00003e00| 09 20 2a 20 74 68 75 73 | 20 6f 6e 6c 79 20 74 77 |. * thus| only tw|
|00003e10| 6f 20 70 6f 69 6e 74 65 | 72 73 20 73 75 66 66 69 |o pointe|rs suffi|
|00003e20| 63 65 20 66 6f 72 20 62 | 6f 74 68 20 61 64 64 69 |ce for b|oth addi|
|00003e30| 74 69 6f 6e 73 2e 20 20 | 4b 65 65 70 20 74 68 65 |tions.  |Keep the|
|00003e40| 0a 58 09 20 2a 20 66 69 | 6e 61 6c 20 63 61 72 72 |.X. * fi|nal carr|
|00003e50| 79 20 66 72 6f 6d 20 74 | 68 65 73 65 20 61 64 64 |y from t|hese add|
|00003e60| 69 74 69 6f 6e 73 20 66 | 6f 72 20 6c 61 74 65 72 |itions f|or later|
|00003e70| 20 75 73 65 20 73 69 6e | 63 65 20 77 65 20 63 61 | use sin|ce we ca|
|00003e80| 6e 6e 6f 74 0a 58 09 20 | 2a 20 61 66 66 6f 72 64 |nnot.X. |* afford|
|00003e90| 20 74 6f 20 63 68 61 6e | 67 65 20 74 68 65 20 74 | to chan|ge the t|
|00003ea0| 6f 70 20 68 61 6c 66 20 | 6f 66 20 41 2a 43 20 79 |op half |of A*C y|
|00003eb0| 65 74 2e 0a 58 09 20 2a | 2f 0a 58 09 68 31 20 3d |et..X. *|/.X.h1 =|
|00003ec0| 20 62 61 73 65 42 44 20 | 2b 20 73 68 69 66 74 3b | baseBD |+ shift;|
|00003ed0| 0a 58 09 68 32 20 3d 20 | 62 61 73 65 41 43 3b 0a |.X.h2 = |baseAC;.|
|00003ee0| 58 09 63 61 72 72 79 41 | 43 42 44 20 3d 20 30 3b |X.carryA|CBD = 0;|
|00003ef0| 0a 58 09 6c 65 6e 20 3d | 20 73 68 69 66 74 3b 0a |.X.len =| shift;.|
|00003f00| 58 09 77 68 69 6c 65 20 | 28 6c 65 6e 2d 2d 29 20 |X.while |(len--) |
|00003f10| 7b 0a 58 09 09 73 69 76 | 61 6c 2e 69 76 61 6c 75 |{.X..siv|al.ivalu|
|00003f20| 65 20 3d 20 28 28 46 55 | 4c 4c 29 20 2a 68 31 29 |e = ((FU|LL) *h1)|
|00003f30| 20 2b 20 28 28 46 55 4c | 4c 29 20 2a 68 32 29 20 | + ((FUL|L) *h2) |
|00003f40| 2b 20 63 61 72 72 79 41 | 43 42 44 3b 0a 58 09 09 |+ carryA|CBD;.X..|
|00003f50| 2a 68 31 2b 2b 20 3d 20 | 73 69 76 61 6c 2e 73 69 |*h1++ = |sival.si|
|00003f60| 6c 6f 77 3b 0a 58 09 09 | 2a 68 32 2b 2b 20 3d 20 |low;.X..|*h2++ = |
|00003f70| 73 69 76 61 6c 2e 73 69 | 6c 6f 77 3b 0a 58 09 09 |sival.si|low;.X..|
|00003f80| 63 61 72 72 79 41 43 42 | 44 20 3d 20 73 69 76 61 |carryACB|D = siva|
|00003f90| 6c 2e 73 69 68 69 67 68 | 3b 0a 58 09 7d 0a 58 0a |l.sihigh|;.X.}.X.|
|00003fa0| 58 09 2f 2a 0a 58 09 20 | 2a 20 4e 6f 77 20 61 64 |X./*.X. |* Now ad|
|00003fb0| 64 20 69 6e 20 74 68 65 | 20 62 6f 74 74 6f 6d 20 |d in the| bottom |
|00003fc0| 68 61 6c 66 20 6f 66 20 | 42 2a 44 20 61 6e 64 20 |half of |B*D and |
|00003fd0| 74 68 65 20 74 6f 70 20 | 68 61 6c 66 20 6f 66 20 |the top |half of |
|00003fe0| 41 2a 43 2e 0a 58 09 20 | 2a 20 54 68 65 73 65 20 |A*C..X. |* These |
|00003ff0| 61 64 64 69 74 69 6f 6e | 73 20 61 72 65 20 73 74 |addition|s are st|
|00004000| 72 61 69 67 68 74 66 6f | 72 77 61 72 64 2c 20 65 |raightfo|rward, e|
|00004010| 78 63 65 70 74 20 74 68 | 61 74 20 41 2a 43 20 73 |xcept th|at A*C s|
|00004020| 68 6f 75 6c 64 0a 58 09 | 20 2a 20 62 65 20 64 6f |hould.X.| * be do|
|00004030| 6e 65 20 66 69 72 73 74 | 20 62 65 63 61 75 73 65 |ne first| because|
|00004040| 20 6f 66 20 70 6f 73 73 | 69 62 6c 65 20 63 61 72 | of poss|ible car|
|00004050| 72 69 65 73 20 66 72 6f | 6d 20 42 2a 44 2c 20 61 |ries fro|m B*D, a|
|00004060| 6e 64 20 74 68 65 0a 58 | 09 20 2a 20 74 6f 70 20 |nd the.X|. * top |
|00004070| 68 61 6c 66 20 6f 66 20 | 41 2a 43 20 6d 69 67 68 |half of |A*C migh|
|00004080| 74 20 6e 6f 74 20 65 78 | 69 73 74 2e 20 20 41 64 |t not ex|ist.  Ad|
|00004090| 64 20 69 6e 20 6f 6e 65 | 20 6f 66 20 74 68 65 20 |d in one| of the |
|000040a0| 63 61 72 72 69 65 73 0a | 58 09 20 2a 20 66 72 6f |carries.|X. * fro|
|000040b0| 6d 20 74 68 65 20 70 72 | 65 76 69 6f 75 73 20 61 |m the pr|evious a|
|000040c0| 64 64 69 74 69 6f 6e 20 | 77 68 69 6c 65 20 77 65 |ddition |while we|
|000040d0| 20 61 72 65 20 61 74 20 | 69 74 2e 0a 58 09 20 2a | are at |it..X. *|
|000040e0| 2f 0a 58 09 68 31 20 3d | 20 62 61 73 65 41 43 20 |/.X.h1 =| baseAC |
|000040f0| 2b 20 73 68 69 66 74 3b | 0a 58 09 68 64 20 3d 20 |+ shift;|.X.hd = |
|00004100| 62 61 73 65 41 43 3b 0a | 58 09 63 61 72 72 79 20 |baseAC;.|X.carry |
|00004110| 3d 20 63 61 72 72 79 41 | 43 42 44 3b 0a 58 09 6c |= carryA|CBD;.X.l|
|00004120| 65 6e 20 3d 20 73 69 7a | 65 74 6f 74 61 6c 20 2d |en = siz|etotal -|
|00004130| 20 33 20 2a 20 73 68 69 | 66 74 3b 0a 58 09 77 68 | 3 * shi|ft;.X.wh|
|00004140| 69 6c 65 20 28 6c 65 6e | 2d 2d 29 20 7b 0a 58 09 |ile (len|--) {.X.|
|00004150| 09 73 69 76 61 6c 2e 69 | 76 61 6c 75 65 20 3d 20 |.sival.i|value = |
|00004160| 28 28 46 55 4c 4c 29 20 | 2a 68 31 2b 2b 29 20 2b |((FULL) |*h1++) +|
|00004170| 20 28 28 46 55 4c 4c 29 | 20 2a 68 64 29 20 2b 20 | ((FULL)| *hd) + |
|00004180| 63 61 72 72 79 3b 0a 58 | 09 09 2a 68 64 2b 2b 20 |carry;.X|..*hd++ |
|00004190| 3d 20 73 69 76 61 6c 2e | 73 69 6c 6f 77 3b 0a 58 |= sival.|silow;.X|
|000041a0| 09 09 63 61 72 72 79 20 | 3d 20 73 69 76 61 6c 2e |..carry |= sival.|
|000041b0| 73 69 68 69 67 68 3b 0a | 58 09 7d 0a 58 09 77 68 |sihigh;.|X.}.X.wh|
|000041c0| 69 6c 65 20 28 63 61 72 | 72 79 29 20 7b 0a 58 09 |ile (car|ry) {.X.|
|000041d0| 09 73 69 76 61 6c 2e 69 | 76 61 6c 75 65 20 3d 20 |.sival.i|value = |
|000041e0| 28 28 46 55 4c 4c 29 20 | 2a 68 64 29 20 2b 20 63 |((FULL) |*hd) + c|
|000041f0| 61 72 72 79 3b 0a 58 09 | 09 2a 68 64 2b 2b 20 3d |arry;.X.|.*hd++ =|
|00004200| 20 73 69 76 61 6c 2e 73 | 69 6c 6f 77 3b 0a 58 09 | sival.s|ilow;.X.|
|00004210| 09 63 61 72 72 79 20 3d | 20 73 69 76 61 6c 2e 73 |.carry =| sival.s|
|00004220| 69 68 69 67 68 3b 0a 58 | 09 7d 0a 58 0a 58 09 68 |ihigh;.X|.}.X.X.h|
|00004230| 31 20 3d 20 62 61 73 65 | 42 44 3b 0a 58 09 68 64 |1 = base|BD;.X.hd|
|00004240| 20 3d 20 62 61 73 65 42 | 44 20 2b 20 73 68 69 66 | = baseB|D + shif|
|00004250| 74 3b 0a 58 09 63 61 72 | 72 79 20 3d 20 30 3b 0a |t;.X.car|ry = 0;.|
|00004260| 58 09 6c 65 6e 20 3d 20 | 73 68 69 66 74 3b 0a 58 |X.len = |shift;.X|
|00004270| 09 77 68 69 6c 65 20 28 | 6c 65 6e 2d 2d 29 20 7b |.while (|len--) {|
|00004280| 0a 58 09 09 73 69 76 61 | 6c 2e 69 76 61 6c 75 65 |.X..siva|l.ivalue|
|00004290| 20 3d 20 28 28 46 55 4c | 4c 29 20 2a 68 31 2b 2b | = ((FUL|L) *h1++|
|000042a0| 29 20 2b 20 28 28 46 55 | 4c 4c 29 20 2a 68 64 29 |) + ((FU|LL) *hd)|
|000042b0| 20 2b 20 63 61 72 72 79 | 3b 0a 58 09 09 2a 68 64 | + carry|;.X..*hd|
|000042c0| 2b 2b 20 3d 20 73 69 76 | 61 6c 2e 73 69 6c 6f 77 |++ = siv|al.silow|
|000042d0| 3b 0a 58 09 09 63 61 72 | 72 79 20 3d 20 73 69 76 |;.X..car|ry = siv|
|000042e0| 61 6c 2e 73 69 68 69 67 | 68 3b 0a 58 09 7d 0a 58 |al.sihig|h;.X.}.X|
|000042f0| 09 77 68 69 6c 65 20 28 | 63 61 72 72 79 29 20 7b |.while (|carry) {|
|00004300| 0a 58 09 09 73 69 76 61 | 6c 2e 69 76 61 6c 75 65 |.X..siva|l.ivalue|
|00004310| 20 3d 20 28 28 46 55 4c | 4c 29 20 2a 68 64 29 20 | = ((FUL|L) *hd) |
|00004320| 2b 20 63 61 72 72 79 3b | 0a 58 09 09 2a 68 64 2b |+ carry;|.X..*hd+|
|00004330| 2b 20 3d 20 73 69 76 61 | 6c 2e 73 69 6c 6f 77 3b |+ = siva|l.silow;|
|00004340| 0a 58 09 09 63 61 72 72 | 79 20 3d 20 73 69 76 61 |.X..carr|y = siva|
|00004350| 6c 2e 73 69 68 69 67 68 | 3b 0a 58 09 7d 0a 58 0a |l.sihigh|;.X.}.X.|
|00004360| 58 09 2f 2a 0a 58 09 20 | 2a 20 4e 6f 77 20 66 69 |X./*.X. |* Now fi|
|00004370| 6e 61 6c 6c 79 20 61 64 | 64 20 69 6e 20 74 68 65 |nally ad|d in the|
|00004380| 20 6f 74 68 65 72 20 64 | 65 6c 61 79 65 64 20 63 | other d|elayed c|
|00004390| 61 72 72 79 20 66 72 6f | 6d 20 74 68 65 0a 58 09 |arry fro|m the.X.|
|000043a0| 20 2a 20 61 62 6f 76 65 | 20 61 64 64 69 74 69 6f | * above| additio|
|000043b0| 6e 2e 0a 58 09 20 2a 2f | 0a 58 09 68 64 20 3d 20 |n..X. */|.X.hd = |
|000043c0| 62 61 73 65 41 43 20 2b | 20 73 68 69 66 74 3b 0a |baseAC +| shift;.|
|000043d0| 58 09 77 68 69 6c 65 20 | 28 63 61 72 72 79 41 43 |X.while |(carryAC|
|000043e0| 42 44 29 20 7b 0a 58 09 | 09 73 69 76 61 6c 2e 69 |BD) {.X.|.sival.i|
|000043f0| 76 61 6c 75 65 20 3d 20 | 28 28 46 55 4c 4c 29 20 |value = |((FULL) |
|00004400| 2a 68 64 29 20 2b 20 63 | 61 72 72 79 41 43 42 44 |*hd) + c|arryACBD|
|00004410| 3b 0a 58 09 09 2a 68 64 | 2b 2b 20 3d 20 73 69 76 |;.X..*hd|++ = siv|
|00004420| 61 6c 2e 73 69 6c 6f 77 | 3b 0a 58 09 09 63 61 72 |al.silow|;.X..car|
|00004430| 72 79 41 43 42 44 20 3d | 20 73 69 76 61 6c 2e 73 |ryACBD =| sival.s|
|00004440| 69 68 69 67 68 3b 0a 58 | 09 7d 0a 58 0a 58 09 2f |ihigh;.X|.}.X.X./|
|00004450| 2a 0a 58 09 20 2a 20 4e | 6f 77 20 66 69 6e 61 6c |*.X. * N|ow final|
|00004460| 6c 79 20 61 64 64 20 6f | 72 20 73 75 62 74 72 61 |ly add o|r subtra|
|00004470| 63 74 20 28 41 2d 42 29 | 2a 28 44 2d 43 29 20 69 |ct (A-B)|*(D-C) i|
|00004480| 6e 74 6f 20 74 68 65 20 | 66 69 6e 61 6c 20 72 65 |nto the |final re|
|00004490| 73 75 6c 74 20 61 74 0a | 58 09 20 2a 20 74 68 65 |sult at.|X. * the|
|000044a0| 20 63 6f 72 72 65 63 74 | 20 70 6f 73 69 74 69 6f | correct| positio|
|000044b0| 6e 20 28 53 29 2c 20 61 | 63 63 6f 72 64 69 6e 67 |n (S), a|ccording|
|000044c0| 20 74 6f 20 77 68 65 74 | 68 65 72 20 69 74 20 69 | to whet|her it i|
|000044d0| 73 20 70 6f 73 69 74 69 | 76 65 20 6f 72 0a 58 09 |s positi|ve or.X.|
|000044e0| 20 2a 20 6e 65 67 61 74 | 69 76 65 2e 20 20 57 68 | * negat|ive.  Wh|
|000044f0| 65 6e 20 73 75 62 74 72 | 61 63 74 69 6e 67 2c 20 |en subtr|acting, |
|00004500| 74 68 65 20 61 6e 73 77 | 65 72 20 63 61 6e 6e 6f |the answ|er canno|
|00004510| 74 20 67 6f 20 6e 65 67 | 61 74 69 76 65 2e 0a 58 |t go neg|ative..X|
|00004520| 09 20 2a 2f 0a 58 09 68 | 31 20 3d 20 62 61 73 65 |. */.X.h|1 = base|
|00004530| 41 42 44 43 3b 0a 58 09 | 68 64 20 3d 20 61 6e 73 |ABDC;.X.|hd = ans|
|00004540| 20 2b 20 73 68 69 66 74 | 3b 0a 58 09 63 61 72 72 | + shift|;.X.carr|
|00004550| 79 20 3d 20 30 3b 0a 58 | 09 6c 65 6e 20 3d 20 73 |y = 0;.X|.len = s|
|00004560| 69 7a 65 41 42 44 43 3b | 0a 58 09 69 66 20 28 6e |izeABDC;|.X.if (n|
|00004570| 65 67 29 20 7b 0a 58 09 | 09 77 68 69 6c 65 20 28 |eg) {.X.|.while (|
|00004580| 6c 65 6e 2d 2d 29 20 7b | 0a 58 09 09 09 73 69 76 |len--) {|.X...siv|
|00004590| 61 6c 2e 69 76 61 6c 75 | 65 20 3d 20 42 41 53 45 |al.ivalu|e = BASE|
|000045a0| 31 20 2d 20 28 28 46 55 | 4c 4c 29 20 2a 68 64 29 |1 - ((FU|LL) *hd)|
|000045b0| 20 2b 0a 58 09 09 09 09 | 28 28 46 55 4c 4c 29 20 | +.X....|((FULL) |
|000045c0| 2a 68 31 2b 2b 29 20 2b | 20 63 61 72 72 79 3b 0a |*h1++) +| carry;.|
|000045d0| 58 09 09 09 2a 68 64 2b | 2b 20 3d 20 42 41 53 45 |X...*hd+|+ = BASE|
|000045e0| 31 20 2d 20 73 69 76 61 | 6c 2e 73 69 6c 6f 77 3b |1 - siva|l.silow;|
|000045f0| 0a 58 09 09 09 63 61 72 | 72 79 20 3d 20 73 69 76 |.X...car|ry = siv|
|00004600| 61 6c 2e 73 69 68 69 67 | 68 3b 0a 58 09 09 7d 0a |al.sihig|h;.X..}.|
|00004610| 58 09 09 77 68 69 6c 65 | 20 28 63 61 72 72 79 29 |X..while| (carry)|
|00004620| 20 7b 0a 58 09 09 09 73 | 69 76 61 6c 2e 69 76 61 | {.X...s|ival.iva|
|00004630| 6c 75 65 20 3d 20 42 41 | 53 45 31 20 2d 20 28 28 |lue = BA|SE1 - ((|
|00004640| 46 55 4c 4c 29 20 2a 68 | 64 29 20 2b 20 63 61 72 |FULL) *h|d) + car|
|00004650| 72 79 3b 0a 58 09 09 09 | 2a 68 64 2b 2b 20 3d 20 |ry;.X...|*hd++ = |
|00004660| 42 41 53 45 31 20 2d 20 | 73 69 76 61 6c 2e 73 69 |BASE1 - |sival.si|
|00004670| 6c 6f 77 3b 0a 58 09 09 | 09 63 61 72 72 79 20 3d |low;.X..|.carry =|
|00004680| 20 73 69 76 61 6c 2e 73 | 69 68 69 67 68 3b 0a 58 | sival.s|ihigh;.X|
|00004690| 09 09 7d 0a 58 09 7d 20 | 65 6c 73 65 20 7b 0a 58 |..}.X.} |else {.X|
|000046a0| 09 09 77 68 69 6c 65 20 | 28 6c 65 6e 2d 2d 29 20 |..while |(len--) |
|000046b0| 7b 0a 58 09 09 09 73 69 | 76 61 6c 2e 69 76 61 6c |{.X...si|val.ival|
|000046c0| 75 65 20 3d 20 28 28 46 | 55 4c 4c 29 20 2a 68 31 |ue = ((F|ULL) *h1|
|000046d0| 2b 2b 29 20 2b 20 28 28 | 46 55 4c 4c 29 20 2a 68 |++) + ((|FULL) *h|
|000046e0| 64 29 20 2b 20 63 61 72 | 72 79 3b 0a 58 09 09 09 |d) + car|ry;.X...|
|000046f0| 2a 68 64 2b 2b 20 3d 20 | 73 69 76 61 6c 2e 73 69 |*hd++ = |sival.si|
|00004700| 6c 6f 77 3b 0a 58 09 09 | 09 63 61 72 72 79 20 3d |low;.X..|.carry =|
|00004710| 20 73 69 76 61 6c 2e 73 | 69 68 69 67 68 3b 0a 58 | sival.s|ihigh;.X|
|00004720| 09 09 7d 0a 58 09 09 77 | 68 69 6c 65 20 28 63 61 |..}.X..w|hile (ca|
|00004730| 72 72 79 29 20 7b 0a 58 | 09 09 09 73 69 76 61 6c |rry) {.X|...sival|
|00004740| 2e 69 76 61 6c 75 65 20 | 3d 20 28 28 46 55 4c 4c |.ivalue |= ((FULL|
|00004750| 29 20 2a 68 64 29 20 2b | 20 63 61 72 72 79 3b 0a |) *hd) +| carry;.|
|00004760| 58 09 09 09 2a 68 64 2b | 2b 20 3d 20 73 69 76 61 |X...*hd+|+ = siva|
|00004770| 6c 2e 73 69 6c 6f 77 3b | 0a 58 09 09 09 63 61 72 |l.silow;|.X...car|
|00004780| 72 79 20 3d 20 73 69 76 | 61 6c 2e 73 69 68 69 67 |ry = siv|al.sihig|
|00004790| 68 3b 0a 58 09 09 7d 0a | 58 09 7d 0a 58 0a 58 09 |h;.X..}.|X.}.X.X.|
|000047a0| 2f 2a 0a 58 09 20 2a 20 | 46 69 6e 61 6c 6c 79 20 |/*.X. * |Finally |
|000047b0| 64 65 74 65 72 6d 69 6e | 65 20 74 68 65 20 73 69 |determin|e the si|
|000047c0| 7a 65 20 6f 66 20 74 68 | 65 20 66 69 6e 61 6c 20 |ze of th|e final |
|000047d0| 72 65 73 75 6c 74 20 61 | 6e 64 20 72 65 74 75 72 |result a|nd retur|
|000047e0| 6e 20 74 68 61 74 2e 0a | 58 09 20 2a 2f 0a 58 09 |n that..|X. */.X.|
|000047f0| 6c 65 6e 20 3d 20 73 69 | 7a 65 74 6f 74 61 6c 3b |len = si|zetotal;|
|00004800| 0a 58 09 68 64 20 3d 20 | 26 61 6e 73 5b 6c 65 6e |.X.hd = |&ans[len|
|00004810| 20 2d 20 31 5d 3b 0a 58 | 09 77 68 69 6c 65 20 28 | - 1];.X|.while (|
|00004820| 28 2a 68 64 20 3d 3d 20 | 30 29 20 26 26 20 28 6c |(*hd == |0) && (l|
|00004830| 65 6e 20 3e 20 31 29 29 | 20 7b 0a 58 09 09 68 64 |en > 1))| {.X..hd|
|00004840| 2d 2d 3b 0a 58 09 09 6c | 65 6e 2d 2d 3b 0a 58 09 |--;.X..l|en--;.X.|
|00004850| 7d 0a 58 09 74 65 6d 70 | 62 75 66 20 3d 20 74 65 |}.X.temp|buf = te|
|00004860| 6d 70 3b 0a 58 09 72 65 | 74 75 72 6e 20 6c 65 6e |mp;.X.re|turn len|
|00004870| 3b 0a 58 7d 0a 58 0a 58 | 0a 58 2f 2a 0a 58 20 2a |;.X}.X.X|.X/*.X *|
|00004880| 20 53 71 75 61 72 65 20 | 61 20 6e 75 6d 62 65 72 | Square |a number|
|00004890| 20 62 79 20 75 73 69 6e | 67 20 74 68 65 20 66 6f | by usin|g the fo|
|000048a0| 6c 6c 6f 77 69 6e 67 20 | 66 6f 72 6d 75 6c 61 20 |llowing |formula |
|000048b0| 72 65 63 75 72 73 69 76 | 65 6c 79 3a 0a 58 20 2a |recursiv|ely:.X *|
|000048c0| 09 28 41 2a 53 2b 42 29 | 5e 32 20 3d 20 28 53 5e |.(A*S+B)|^2 = (S^|
|000048d0| 32 2b 53 29 2a 41 5e 32 | 20 2b 20 28 53 2b 31 29 |2+S)*A^2| + (S+1)|
|000048e0| 2a 42 5e 32 20 2d 20 53 | 2a 28 41 2d 42 29 5e 32 |*B^2 - S|*(A-B)^2|
|000048f0| 0a 58 20 2a 20 77 68 65 | 72 65 20 53 20 69 73 20 |.X * whe|re S is |
|00004900| 61 20 70 6f 77 65 72 20 | 6f 66 20 32 5e 31 36 2c |a power |of 2^16,|
|00004910| 20 61 6e 64 20 73 6f 20 | 6d 75 6c 74 69 70 6c 69 | and so |multipli|
|00004920| 65 73 20 62 79 20 69 74 | 20 61 72 65 20 73 68 69 |es by it| are shi|
|00004930| 66 74 73 2c 0a 58 20 2a | 20 61 6e 64 20 41 20 61 |fts,.X *| and A a|
|00004940| 6e 64 20 42 20 61 72 65 | 20 74 68 65 20 6c 65 66 |nd B are| the lef|
|00004950| 74 20 61 6e 64 20 72 69 | 67 68 74 20 68 61 6c 66 |t and ri|ght half|
|00004960| 73 20 6f 66 20 74 68 65 | 20 6e 75 6d 62 65 72 20 |s of the| number |
|00004970| 74 6f 20 73 71 75 61 72 | 65 2e 0a 58 20 2a 2f 0a |to squar|e..X */.|
|00004980| 58 76 6f 69 64 0a 58 7a | 73 71 75 61 72 65 28 7a |Xvoid.Xz|square(z|
|00004990| 2c 20 72 65 73 29 0a 58 | 09 5a 56 41 4c 55 45 20 |, res).X|.ZVALUE |
|000049a0| 7a 2c 20 2a 72 65 73 3b | 0a 58 7b 0a 58 09 4c 45 |z, *res;|.X{.X.LE|
|000049b0| 4e 20 6c 65 6e 3b 0a 58 | 0a 58 09 69 66 20 28 69 |N len;.X|.X.if (i|
|000049c0| 73 7a 65 72 6f 28 7a 29 | 29 20 7b 0a 58 09 09 2a |szero(z)|) {.X..*|
|000049d0| 72 65 73 20 3d 20 5f 7a | 65 72 6f 5f 3b 0a 58 09 |res = _z|ero_;.X.|
|000049e0| 09 72 65 74 75 72 6e 3b | 0a 58 09 7d 0a 58 09 69 |.return;|.X.}.X.i|
|000049f0| 66 20 28 69 73 75 6e 69 | 74 28 7a 29 29 20 7b 0a |f (isuni|t(z)) {.|
|00004a00| 58 09 09 2a 72 65 73 20 | 3d 20 5f 6f 6e 65 5f 3b |X..*res |= _one_;|
|00004a10| 0a 58 09 09 72 65 74 75 | 72 6e 3b 0a 58 09 7d 0a |.X..retu|rn;.X.}.|
|00004a20| 58 0a 58 09 2f 2a 0a 58 | 09 20 2a 20 41 6c 6c 6f |X.X./*.X|. * Allo|
|00004a30| 63 61 74 65 20 61 20 74 | 65 6d 70 6f 72 61 72 79 |cate a t|emporary|
|00004a40| 20 61 72 72 61 79 20 69 | 66 20 6e 65 63 65 73 73 | array i|f necess|
|00004a50| 61 72 79 20 66 6f 72 20 | 74 68 65 20 72 65 63 75 |ary for |the recu|
|00004a60| 72 73 69 6f 6e 20 74 6f | 20 75 73 65 2e 0a 58 09 |rsion to| use..X.|
|00004a70| 20 2a 20 54 68 65 20 61 | 72 72 61 79 20 6e 65 65 | * The a|rray nee|
|00004a80| 64 73 20 74 6f 20 62 65 | 20 61 6c 6c 6f 63 61 74 |ds to be| allocat|
|00004a90| 65 64 20 6c 61 72 67 65 | 20 65 6e 6f 75 67 68 20 |ed large| enough |
|00004aa0| 66 6f 72 20 61 6c 6c 20 | 6f 66 20 74 68 65 0a 58 |for all |of the.X|
|00004ab0| 09 20 2a 20 74 65 6d 70 | 6f 72 61 72 79 20 72 65 |. * temp|orary re|
|00004ac0| 73 75 6c 74 73 20 74 6f | 20 66 69 74 20 69 6e 2e |sults to| fit in.|
|00004ad0| 20 20 54 68 69 73 20 73 | 69 7a 65 20 69 73 20 61 |  This s|ize is a|
|00004ae0| 62 6f 75 74 20 33 20 74 | 69 6d 65 73 20 74 68 65 |bout 3 t|imes the|
|00004af0| 0a 58 09 20 2a 20 73 69 | 7a 65 20 6f 66 20 74 68 |.X. * si|ze of th|
|00004b00| 65 20 6f 72 69 67 69 6e | 61 6c 20 6e 75 6d 62 65 |e origin|al numbe|
|00004b10| 72 2c 20 73 69 6e 63 65 | 20 65 61 63 68 20 72 65 |r, since| each re|
|00004b20| 63 75 72 73 69 6f 6e 20 | 6c 65 76 65 6c 20 75 73 |cursion |level us|
|00004b30| 65 73 20 33 2f 32 0a 58 | 09 20 2a 20 6f 66 20 74 |es 3/2.X|. * of t|
|00004b40| 68 65 20 73 69 7a 65 20 | 6f 66 20 69 74 73 20 67 |he size |of its g|
|00004b50| 69 76 65 6e 20 6e 75 6d | 62 65 72 2c 20 61 6e 64 |iven num|ber, and|
|00004b60| 20 77 68 6f 73 65 20 73 | 69 7a 65 20 69 73 20 31 | whose s|ize is 1|
|00004b70| 2f 32 20 74 68 65 20 73 | 69 7a 65 0a 58 09 20 2a |/2 the s|ize.X. *|
|00004b80| 20 6f 66 20 74 68 65 20 | 70 72 65 76 69 6f 75 73 | of the |previous|
|00004b90| 20 6c 65 76 65 6c 2e 20 | 20 54 68 65 20 73 75 6d | level. | The sum|
|00004ba0| 20 6f 66 20 74 68 65 20 | 69 6e 66 69 6e 69 74 65 | of the |infinite|
|00004bb0| 20 73 65 72 69 65 73 20 | 69 73 20 33 2e 0a 58 09 | series |is 3..X.|
|00004bc0| 20 2a 20 41 6c 6c 6f 63 | 61 74 65 20 73 6f 6d 65 | * Alloc|ate some|
|00004bd0| 20 65 78 74 72 61 20 77 | 6f 72 64 73 20 66 6f 72 | extra w|ords for|
|00004be0| 20 72 6f 75 6e 64 69 6e | 67 20 75 70 20 74 68 65 | roundin|g up the|
|00004bf0| 20 73 69 7a 65 73 2e 0a | 58 09 20 2a 2f 0a 58 09 | sizes..|X. */.X.|
|00004c00| 6c 65 6e 20 3d 20 33 20 | 2a 20 7a 2e 6c 65 6e 20 |len = 3 |* z.len |
|00004c10| 2b 20 33 32 3b 0a 58 09 | 74 65 6d 70 62 75 66 20 |+ 32;.X.|tempbuf |
|00004c20| 3d 20 7a 61 6c 6c 6f 63 | 74 65 6d 70 28 6c 65 6e |= zalloc|temp(len|
|00004c30| 29 3b 0a 58 0a 58 09 72 | 65 73 2d 3e 73 69 67 6e |);.X.X.r|es->sign|
|00004c40| 20 3d 20 30 3b 0a 58 09 | 72 65 73 2d 3e 76 20 3d | = 0;.X.|res->v =|
|00004c50| 20 61 6c 6c 6f 63 28 7a | 2e 6c 65 6e 20 2a 20 32 | alloc(z|.len * 2|
|00004c60| 29 3b 0a 58 09 72 65 73 | 2d 3e 6c 65 6e 20 3d 20 |);.X.res|->len = |
|00004c70| 64 6f 73 71 75 61 72 65 | 28 7a 2e 76 2c 20 7a 2e |dosquare|(z.v, z.|
|00004c80| 6c 65 6e 2c 20 72 65 73 | 2d 3e 76 29 3b 0a 58 7d |len, res|->v);.X}|
|00004c90| 0a 58 0a 58 0a 58 2f 2a | 0a 58 20 2a 20 52 65 63 |.X.X.X/*|.X * Rec|
|00004ca0| 75 72 73 69 76 65 20 72 | 6f 75 74 69 6e 65 20 74 |ursive r|outine t|
|00004cb0| 6f 20 73 71 75 61 72 65 | 20 61 20 6e 75 6d 62 65 |o square| a numbe|
|00004cc0| 72 20 62 79 20 73 70 6c | 69 74 74 69 6e 67 20 69 |r by spl|itting i|
|00004cd0| 74 20 75 70 20 69 6e 74 | 6f 20 74 77 6f 20 6e 75 |t up int|o two nu|
|00004ce0| 6d 62 65 72 73 0a 58 20 | 2a 20 6f 66 20 68 61 6c |mbers.X |* of hal|
|00004cf0| 66 20 74 68 65 20 73 69 | 7a 65 2c 20 61 6e 64 20 |f the si|ze, and |
|00004d00| 75 73 69 6e 67 20 74 68 | 65 20 72 65 73 75 6c 74 |using th|e result|
|00004d10| 73 20 6f 66 20 73 71 75 | 61 72 69 6e 67 20 74 68 |s of squ|aring th|
|00004d20| 65 20 73 75 62 70 69 65 | 63 65 73 2e 0a 58 20 2a |e subpie|ces..X *|
|00004d30| 20 54 68 65 20 72 65 73 | 75 6c 74 20 69 73 20 70 | The res|ult is p|
|00004d40| 6c 61 63 65 64 20 69 6e | 20 74 68 65 20 69 6e 64 |laced in| the ind|
|00004d50| 69 63 61 74 65 64 20 61 | 72 72 61 79 2c 20 77 68 |icated a|rray, wh|
|00004d60| 69 63 68 20 6d 75 73 74 | 20 62 65 20 6c 61 72 67 |ich must| be larg|
|00004d70| 65 0a 58 20 2a 20 65 6e | 6f 75 67 68 20 66 6f 72 |e.X * en|ough for|
|00004d80| 20 74 68 65 20 72 65 73 | 75 6c 74 20 28 73 69 7a | the res|ult (siz|
|00004d90| 65 20 2a 20 32 29 2e 20 | 20 52 65 74 75 72 6e 73 |e * 2). | Returns|
|00004da0| 20 74 68 65 20 73 69 7a | 65 20 6f 66 20 74 68 65 | the siz|e of the|
|00004db0| 20 72 65 73 75 6c 74 2e | 0a 58 20 2a 20 54 68 69 | result.|.X * Thi|
|00004dc0| 73 20 75 73 65 73 20 61 | 20 74 65 6d 70 6f 72 61 |s uses a| tempora|
|00004dd0| 72 79 20 61 72 72 61 79 | 20 77 68 69 63 68 20 6d |ry array| which m|
|00004de0| 75 73 74 20 62 65 20 33 | 20 74 69 6d 65 73 20 61 |ust be 3| times a|
|00004df0| 73 20 6c 61 72 67 65 20 | 61 73 20 74 68 65 0a 58 |s large |as the.X|
|00004e00| 20 2a 20 73 69 7a 65 20 | 6f 66 20 74 68 65 20 6e | * size |of the n|
|00004e10| 75 6d 62 65 72 20 61 74 | 20 74 68 65 20 74 6f 70 |umber at| the top|
|00004e20| 20 6c 65 76 65 6c 20 72 | 65 63 75 72 73 69 76 65 | level r|ecursive|
|00004e30| 20 63 61 6c 6c 2e 0a 58 | 20 2a 2f 0a 58 73 74 61 | call..X| */.Xsta|
|00004e40| 74 69 63 20 4c 45 4e 0a | 58 64 6f 73 71 75 61 72 |tic LEN.|Xdosquar|
|00004e50| 65 28 76 70 2c 20 73 69 | 7a 65 2c 20 61 6e 73 29 |e(vp, si|ze, ans)|
|00004e60| 0a 58 09 48 41 4c 46 20 | 2a 76 70 3b 09 09 2f 2a |.X.HALF |*vp;../*|
|00004e70| 20 76 61 6c 75 65 20 74 | 6f 20 62 65 20 73 71 75 | value t|o be squ|
|00004e80| 61 72 65 64 20 2a 2f 0a | 58 09 4c 45 4e 20 73 69 |ared */.|X.LEN si|
|00004e90| 7a 65 3b 09 09 2f 2a 20 | 6c 65 6e 67 74 68 20 6f |ze;../* |length o|
|00004ea0| 66 20 76 61 6c 75 65 20 | 74 6f 20 73 71 75 61 72 |f value |to squar|
|00004eb0| 65 20 2a 2f 0a 58 09 48 | 41 4c 46 20 2a 61 6e 73 |e */.X.H|ALF *ans|
|00004ec0| 3b 09 09 2f 2a 20 6c 6f | 63 61 74 69 6f 6e 20 66 |;../* lo|cation f|
|00004ed0| 6f 72 20 72 65 73 75 6c | 74 20 2a 2f 0a 58 7b 0a |or resul|t */.X{.|
|00004ee0| 58 09 4c 45 4e 20 73 68 | 69 66 74 3b 09 09 2f 2a |X.LEN sh|ift;../*|
|00004ef0| 20 61 6d 6f 75 6e 74 20 | 6e 75 6d 62 65 72 73 20 | amount |numbers |
|00004f00| 61 72 65 20 73 68 69 66 | 74 65 64 20 62 79 20 2a |are shif|ted by *|
|00004f10| 2f 0a 58 09 4c 45 4e 20 | 73 69 7a 65 41 3b 09 09 |/.X.LEN |sizeA;..|
|00004f20| 2f 2a 20 73 69 7a 65 20 | 6f 66 20 6c 65 66 74 20 |/* size |of left |
|00004f30| 68 61 6c 66 20 6f 66 20 | 6e 75 6d 62 65 72 20 74 |half of |number t|
|00004f40| 6f 20 73 71 75 61 72 65 | 20 2a 2f 0a 58 09 4c 45 |o square| */.X.LE|
|00004f50| 4e 20 73 69 7a 65 42 3b | 09 09 2f 2a 20 73 69 7a |N sizeB;|../* siz|
|00004f60| 65 20 6f 66 20 72 69 67 | 68 74 20 68 61 6c 66 20 |e of rig|ht half |
|00004f70| 6f 66 20 6e 75 6d 62 65 | 72 20 74 6f 20 73 71 75 |of numbe|r to squ|
|00004f80| 61 72 65 20 2a 2f 0a 58 | 09 4c 45 4e 20 73 69 7a |are */.X|.LEN siz|
|00004f90| 65 41 41 3b 09 09 2f 2a | 20 73 69 7a 65 20 6f 66 |eAA;../*| size of|
|00004fa0| 20 73 71 75 61 72 65 20 | 6f 66 20 6c 65 66 74 20 | square |of left |
|00004fb0| 68 61 6c 66 20 2a 2f 0a | 58 09 4c 45 4e 20 73 69 |half */.|X.LEN si|
|00004fc0| 7a 65 42 42 3b 09 09 2f | 2a 20 73 69 7a 65 20 6f |zeBB;../|* size o|
|00004fd0| 66 20 73 71 75 61 72 65 | 20 6f 66 20 72 69 67 68 |f square| of righ|
|00004fe0| 74 20 68 61 6c 66 20 2a | 2f 0a 58 09 4c 45 4e 20 |t half *|/.X.LEN |
|00004ff0| 73 69 7a 65 41 41 42 42 | 3b 09 09 2f 2a 20 73 69 |sizeAABB|;../* si|
|00005000| 7a 65 20 6f 66 20 73 75 | 6d 20 6f 66 20 73 71 75 |ze of su|m of squ|
|00005010| 61 72 65 73 20 6f 66 20 | 41 20 61 6e 64 20 42 20 |ares of |A and B |
|00005020| 2a 2f 0a 58 09 4c 45 4e | 20 73 69 7a 65 41 42 3b |*/.X.LEN| sizeAB;|
|00005030| 09 09 2f 2a 20 73 69 7a | 65 20 6f 66 20 64 69 66 |../* siz|e of dif|
|00005040| 66 65 72 65 6e 63 65 20 | 6f 66 20 41 20 61 6e 64 |ference |of A and|
|00005050| 20 42 20 2a 2f 0a 58 09 | 4c 45 4e 20 73 69 7a 65 | B */.X.|LEN size|
|00005060| 41 42 41 42 3b 09 09 2f | 2a 20 73 69 7a 65 20 6f |ABAB;../|* size o|
|00005070| 66 20 73 71 75 61 72 65 | 20 6f 66 20 64 69 66 66 |f square| of diff|
|00005080| 65 72 65 6e 63 65 20 6f | 66 20 41 20 61 6e 64 20 |erence o|f A and |
|00005090| 42 20 2a 2f 0a 58 09 4c | 45 4e 20 73 75 62 73 69 |B */.X.L|EN subsi|
|000050a0| 7a 65 3b 09 09 2f 2a 20 | 73 69 7a 65 20 6f 66 20 |ze;../* |size of |
|000050b0| 64 69 66 66 65 72 65 6e | 63 65 20 6f 66 20 68 61 |differen|ce of ha|
|000050c0| 6c 66 73 20 2a 2f 0a 58 | 09 4c 45 4e 20 63 6f 70 |lfs */.X|.LEN cop|
|000050d0| 79 73 69 7a 65 3b 09 09 | 2f 2a 20 73 69 7a 65 20 |ysize;..|/* size |
|000050e0| 6f 66 20 6e 75 6d 62 65 | 72 20 6c 65 66 74 20 74 |of numbe|r left t|
|000050f0| 6f 20 63 6f 70 79 20 2a | 2f 0a 58 09 4c 45 4e 20 |o copy *|/.X.LEN |
|00005100| 73 75 6d 73 69 7a 65 3b | 09 09 2f 2a 20 73 69 7a |sumsize;|../* siz|
|00005110| 65 20 6f 66 20 73 75 6d | 20 2a 2f 0a 58 09 4c 45 |e of sum| */.X.LE|
|00005120| 4e 20 73 69 7a 65 74 6f | 74 61 6c 3b 09 09 2f 2a |N sizeto|tal;../*|
|00005130| 20 74 6f 74 61 6c 20 73 | 69 7a 65 20 6f 66 20 73 | total s|ize of s|
|00005140| 71 75 61 72 65 20 2a 2f | 0a 58 09 4c 45 4e 20 6c |quare */|.X.LEN l|
|00005150| 65 6e 3b 09 09 2f 2a 20 | 74 65 6d 70 6f 72 61 72 |en;../* |temporar|
|00005160| 79 20 6c 65 6e 67 74 68 | 20 2a 2f 0a 58 09 4c 45 |y length| */.X.LE|
|00005170| 4e 20 6c 65 6e 31 3b 09 | 09 2f 2a 20 61 6e 6f 74 |N len1;.|./* anot|
|00005180| 68 65 72 20 74 65 6d 70 | 6f 72 61 72 79 20 6c 65 |her temp|orary le|
|00005190| 6e 67 74 68 20 2a 2f 0a | 58 09 46 55 4c 4c 20 63 |ngth */.|X.FULL c|
|000051a0| 61 72 72 79 3b 09 09 2f | 2a 20 63 61 72 72 79 20 |arry;../|* carry |
|000051b0| 64 69 67 69 74 20 66 6f | 72 20 73 6d 61 6c 6c 20 |digit fo|r small |
|000051c0| 6d 75 6c 74 69 70 6c 69 | 63 61 74 69 6f 6e 73 20 |multipli|cations |
|000051d0| 2a 2f 0a 58 09 46 55 4c | 4c 20 64 69 67 69 74 3b |*/.X.FUL|L digit;|
|000051e0| 09 09 2f 2a 20 73 69 6e | 67 6c 65 20 64 69 67 69 |../* sin|gle digi|
|000051f0| 74 20 6d 75 6c 74 69 70 | 6c 79 69 6e 67 20 62 79 |t multip|lying by|
|00005200| 20 2a 2f 0a 58 09 48 41 | 4c 46 20 2a 74 65 6d 70 | */.X.HA|LF *temp|
|00005210| 3b 09 09 2f 2a 20 62 61 | 73 65 20 66 6f 72 20 74 |;../* ba|se for t|
|00005220| 65 6d 70 6f 72 61 72 79 | 20 63 61 6c 63 75 6c 61 |emporary| calcula|
|00005230| 74 69 6f 6e 73 20 2a 2f | 0a 58 09 48 41 4c 46 20 |tions */|.X.HALF |
|00005240| 2a 62 61 73 65 41 3b 09 | 09 2f 2a 20 62 61 73 65 |*baseA;.|./* base|
|00005250| 20 6f 66 20 6c 65 66 74 | 20 68 61 6c 66 20 6f 66 | of left| half of|
|00005260| 20 6e 75 6d 62 65 72 20 | 2a 2f 0a 58 09 48 41 4c | number |*/.X.HAL|
|00005270| 46 20 2a 62 61 73 65 42 | 3b 09 09 2f 2a 20 62 61 |F *baseB|;../* ba|
|00005280| 73 65 20 6f 66 20 72 69 | 67 68 74 20 68 61 6c 66 |se of ri|ght half|
|00005290| 20 6f 66 20 6e 75 6d 62 | 65 72 20 2a 2f 0a 58 09 | of numb|er */.X.|
|000052a0| 48 41 4c 46 20 2a 62 61 | 73 65 41 41 3b 09 09 2f |HALF *ba|seAA;../|
|000052b0| 2a 20 62 61 73 65 20 6f | 66 20 73 71 75 61 72 65 |* base o|f square|
|000052c0| 20 6f 66 20 6c 65 66 74 | 20 68 61 6c 66 20 6f 66 | of left| half of|
|000052d0| 20 6e 75 6d 62 65 72 20 | 2a 2f 0a 58 09 48 41 4c | number |*/.X.HAL|
|000052e0| 46 20 2a 62 61 73 65 42 | 42 3b 09 09 2f 2a 20 62 |F *baseB|B;../* b|
|000052f0| 61 73 65 20 6f 66 20 73 | 71 75 61 72 65 20 6f 66 |ase of s|quare of|
|00005300| 20 72 69 67 68 74 20 68 | 61 6c 66 20 6f 66 20 6e | right h|alf of n|
|00005310| 75 6d 62 65 72 20 2a 2f | 0a 58 09 48 41 4c 46 20 |umber */|.X.HALF |
|00005320| 2a 62 61 73 65 41 41 42 | 42 3b 09 09 2f 2a 20 62 |*baseAAB|B;../* b|
|00005330| 61 73 65 20 6f 66 20 73 | 75 6d 20 6f 66 20 73 71 |ase of s|um of sq|
|00005340| 75 61 72 65 73 20 6f 66 | 20 41 20 61 6e 64 20 42 |uares of| A and B|
|00005350| 20 2a 2f 0a 58 09 48 41 | 4c 46 20 2a 62 61 73 65 | */.X.HA|LF *base|
|00005360| 41 42 3b 09 09 2f 2a 20 | 62 61 73 65 20 6f 66 20 |AB;../* |base of |
|00005370| 64 69 66 66 65 72 65 6e | 63 65 20 6f 66 20 41 20 |differen|ce of A |
|00005380| 61 6e 64 20 42 20 2a 2f | 0a 58 09 48 41 4c 46 20 |and B */|.X.HALF |
|00005390| 2a 62 61 73 65 41 42 41 | 42 3b 09 09 2f 2a 20 62 |*baseABA|B;../* b|
|000053a0| 61 73 65 20 6f 66 20 73 | 71 75 61 72 65 20 6f 66 |ase of s|quare of|
|000053b0| 20 64 69 66 66 65 72 65 | 6e 63 65 20 6f 66 20 41 | differe|nce of A|
|000053c0| 20 61 6e 64 20 42 20 2a | 2f 0a 58 09 72 65 67 69 | and B *|/.X.regi|
|000053d0| 73 74 65 72 20 48 41 4c | 46 20 2a 68 64 2c 20 2a |ster HAL|F *hd, *|
|000053e0| 68 31 2c 20 2a 68 32 2c | 20 2a 68 33 3b 09 2f 2a |h1, *h2,| *h3;./*|
|000053f0| 20 66 6f 72 20 69 6e 6e | 65 72 20 6c 6f 6f 70 73 | for inn|er loops|
|00005400| 20 2a 2f 0a 58 09 53 49 | 55 4e 49 4f 4e 20 73 69 | */.X.SI|UNION si|
|00005410| 76 61 6c 3b 09 09 2f 2a | 20 66 6f 72 20 61 64 64 |val;../*| for add|
|00005420| 69 74 69 6f 6e 20 6f 66 | 20 64 69 67 69 74 73 20 |ition of| digits |
|00005430| 2a 2f 0a 58 0a 58 09 2f | 2a 0a 58 09 20 2a 20 46 |*/.X.X./|*.X. * F|
|00005440| 69 72 73 74 20 74 72 69 | 6d 20 74 68 65 20 6e 75 |irst tri|m the nu|
|00005450| 6d 62 65 72 20 6f 66 20 | 6c 65 61 64 69 6e 67 20 |mber of |leading |
|00005460| 7a 65 72 6f 65 73 2e 0a | 58 09 20 2a 2f 0a 58 09 |zeroes..|X. */.X.|
|00005470| 68 64 20 3d 20 26 76 70 | 5b 73 69 7a 65 20 2d 20 |hd = &vp|[size - |
|00005480| 31 5d 3b 0a 58 09 77 68 | 69 6c 65 20 28 28 2a 68 |1];.X.wh|ile ((*h|
|00005490| 64 20 3d 3d 20 30 29 20 | 26 26 20 28 73 69 7a 65 |d == 0) |&& (size|
|000054a0| 20 3e 20 31 29 29 20 7b | 0a 58 09 09 73 69 7a 65 | > 1)) {|.X..size|
|000054b0| 2d 2d 3b 0a 58 09 09 68 | 64 2d 2d 3b 0a 58 09 7d |--;.X..h|d--;.X.}|
|000054c0| 0a 58 09 73 69 7a 65 74 | 6f 74 61 6c 20 3d 20 73 |.X.sizet|otal = s|
|000054d0| 69 7a 65 20 2b 20 73 69 | 7a 65 3b 0a 58 0a 58 09 |ize + si|ze;.X.X.|
|000054e0| 2f 2a 0a 58 09 20 2a 20 | 49 66 20 74 68 65 20 6e |/*.X. * |If the n|
|000054f0| 75 6d 62 65 72 20 68 61 | 73 20 6f 6e 6c 79 20 61 |umber ha|s only a|
|00005500| 20 73 6d 61 6c 6c 20 6e | 75 6d 62 65 72 20 6f 66 | small n|umber of|
|00005510| 20 64 69 67 69 74 73 2c | 20 74 68 65 6e 20 75 73 | digits,| then us|
|00005520| 65 20 74 68 65 0a 58 09 | 20 2a 20 28 61 6c 6d 6f |e the.X.| * (almo|
|00005530| 73 74 29 20 6e 6f 72 6d | 61 6c 20 6d 75 6c 74 69 |st) norm|al multi|
|00005540| 70 6c 69 63 61 74 69 6f | 6e 20 6d 65 74 68 6f 64 |plicatio|n method|
|00005550| 2e 20 20 4d 75 6c 74 69 | 70 6c 79 20 65 61 63 68 |.  Multi|ply each|
|00005560| 20 68 61 6c 66 77 6f 72 | 64 0a 58 09 20 2a 20 6f | halfwor|d.X. * o|
|00005570| 6e 6c 79 20 62 79 20 74 | 68 6f 73 65 20 68 61 6c |nly by t|hose hal|
|00005580| 66 77 61 72 64 73 20 66 | 75 72 74 68 65 72 20 6f |fwards f|urther o|
|00005590| 6e 20 69 6e 20 74 68 65 | 20 6e 75 6d 62 65 72 2e |n in the| number.|
|000055a0| 20 20 4d 69 73 73 65 64 | 20 74 65 72 6d 73 0a 58 |  Missed| terms.X|
|000055b0| 09 20 2a 20 77 69 6c 6c | 20 74 68 65 6e 20 62 65 |. * will| then be|
|000055c0| 20 74 68 65 20 73 61 6d | 65 20 70 61 69 72 73 20 | the sam|e pairs |
|000055d0| 6f 66 20 70 72 6f 64 75 | 63 74 73 20 72 65 70 65 |of produ|cts repe|
|000055e0| 61 74 65 64 2c 20 61 6e | 64 20 74 68 65 20 73 71 |ated, an|d the sq|
|000055f0| 75 61 72 65 73 0a 58 09 | 20 2a 20 6f 66 20 65 61 |uares.X.| * of ea|
|00005600| 63 68 20 68 61 6c 66 77 | 6f 72 64 2e 20 20 54 68 |ch halfw|ord.  Th|
|00005610| 65 20 66 69 72 73 74 20 | 63 61 73 65 20 69 73 20 |e first |case is |
|00005620| 68 61 6e 64 6c 65 64 20 | 62 79 20 64 6f 75 62 6c |handled |by doubl|
|00005630| 69 6e 67 20 74 68 65 0a | 58 09 20 2a 20 72 65 73 |ing the.|X. * res|
|00005640| 75 6c 74 2e 20 20 54 68 | 65 20 73 65 63 6f 6e 64 |ult.  Th|e second|
|00005650| 20 63 61 73 65 20 69 73 | 20 68 61 6e 64 6c 65 64 | case is| handled|
|00005660| 20 65 78 70 6c 69 63 69 | 74 6c 79 2e 20 20 54 68 | explici|tly.  Th|
|00005670| 65 20 6e 75 6d 62 65 72 | 20 6f 66 0a 58 09 20 2a |e number| of.X. *|
|00005680| 20 64 69 67 69 74 73 20 | 77 68 65 72 65 20 74 68 | digits |where th|
|00005690| 65 20 61 6c 67 6f 72 69 | 74 68 6d 20 63 68 61 6e |e algori|thm chan|
|000056a0| 67 65 73 20 69 73 20 73 | 65 74 74 61 62 6c 65 20 |ges is s|ettable |
|000056b0| 66 72 6f 6d 20 32 20 74 | 6f 20 6d 61 78 69 6e 74 |from 2 t|o maxint|
|000056c0| 2e 0a 58 09 20 2a 2f 0a | 58 09 69 66 20 28 73 69 |..X. */.|X.if (si|
|000056d0| 7a 65 20 3c 20 5f 73 71 | 32 5f 29 20 7b 0a 58 09 |ze < _sq|2_) {.X.|
|000056e0| 09 68 64 20 3d 20 61 6e | 73 3b 0a 58 09 09 6c 65 |.hd = an|s;.X..le|
|000056f0| 6e 20 3d 20 73 69 7a 65 | 74 6f 74 61 6c 3b 0a 58 |n = size|total;.X|
|00005700| 09 09 77 68 69 6c 65 20 | 28 6c 65 6e 2d 2d 29 0a |..while |(len--).|
|00005710| 58 09 09 09 2a 68 64 2b | 2b 20 3d 20 30 3b 0a 58 |X...*hd+|+ = 0;.X|
|00005720| 0a 58 09 09 68 32 20 3d | 20 76 70 3b 0a 58 09 09 |.X..h2 =| vp;.X..|
|00005730| 68 64 20 3d 20 61 6e 73 | 20 2b 20 31 3b 0a 58 09 |hd = ans| + 1;.X.|
|00005740| 09 66 6f 72 20 28 6c 65 | 6e 20 3d 20 73 69 7a 65 |.for (le|n = size|
|00005750| 3b 20 6c 65 6e 2d 2d 3b | 20 68 64 20 2b 3d 20 32 |; len--;| hd += 2|
|00005760| 29 20 7b 0a 58 09 09 09 | 64 69 67 69 74 20 3d 20 |) {.X...|digit = |
|00005770| 28 46 55 4c 4c 29 20 2a | 68 32 2b 2b 3b 0a 58 09 |(FULL) *|h2++;.X.|
|00005780| 09 09 69 66 20 28 64 69 | 67 69 74 20 3d 3d 20 30 |..if (di|git == 0|
|00005790| 29 0a 58 09 09 09 09 63 | 6f 6e 74 69 6e 75 65 3b |).X....c|ontinue;|
|000057a0| 0a 58 09 09 09 68 33 20 | 3d 20 68 32 3b 0a 58 09 |.X...h3 |= h2;.X.|
|000057b0| 09 09 68 31 20 3d 20 68 | 64 3b 0a 58 09 09 09 63 |..h1 = h|d;.X...c|
|000057c0| 61 72 72 79 20 3d 20 30 | 3b 0a 58 09 09 09 6c 65 |arry = 0|;.X...le|
|000057d0| 6e 31 20 3d 20 6c 65 6e | 3b 0a 58 09 09 09 77 68 |n1 = len|;.X...wh|
|000057e0| 69 6c 65 20 28 6c 65 6e | 31 20 3e 3d 20 34 29 20 |ile (len|1 >= 4) |
|000057f0| 7b 09 2f 2a 20 65 78 70 | 61 6e 64 20 74 68 65 20 |{./* exp|and the |
|00005800| 6c 6f 6f 70 20 73 6f 6d | 65 20 2a 2f 0a 58 09 09 |loop som|e */.X..|
|00005810| 09 09 6c 65 6e 31 20 2d | 3d 20 34 3b 0a 58 09 09 |..len1 -|= 4;.X..|
|00005820| 09 09 73 69 76 61 6c 2e | 69 76 61 6c 75 65 20 3d |..sival.|ivalue =|
|00005830| 20 28 64 69 67 69 74 20 | 2a 20 28 28 46 55 4c 4c | (digit |* ((FULL|
|00005840| 29 20 2a 68 33 2b 2b 29 | 29 0a 58 09 09 09 09 09 |) *h3++)|).X.....|
|00005850| 2b 20 28 28 46 55 4c 4c | 29 20 2a 68 31 29 20 2b |+ ((FULL|) *h1) +|
|00005860| 20 63 61 72 72 79 3b 0a | 58 09 09 09 09 2a 68 31 | carry;.|X....*h1|
|00005870| 2b 2b 20 3d 20 73 69 76 | 61 6c 2e 73 69 6c 6f 77 |++ = siv|al.silow|
|00005880| 3b 0a 58 09 09 09 09 73 | 69 76 61 6c 2e 69 76 61 |;.X....s|ival.iva|
|00005890| 6c 75 65 20 3d 20 28 64 | 69 67 69 74 20 2a 20 28 |lue = (d|igit * (|
|000058a0| 28 46 55 4c 4c 29 20 2a | 68 33 2b 2b 29 29 0a 58 |(FULL) *|h3++)).X|
|000058b0| 09 09 09 09 09 2b 20 28 | 28 46 55 4c 4c 29 20 2a |.....+ (|(FULL) *|
|000058c0| 68 31 29 20 2b 20 28 28 | 46 55 4c 4c 29 20 73 69 |h1) + ((|FULL) si|
|000058d0| 76 61 6c 2e 73 69 68 69 | 67 68 29 3b 0a 58 09 09 |val.sihi|gh);.X..|
|000058e0| 09 09 2a 68 31 2b 2b 20 | 3d 20 73 69 76 61 6c 2e |..*h1++ |= sival.|
|000058f0| 73 69 6c 6f 77 3b 0a 58 | 09 09 09 09 73 69 76 61 |silow;.X|....siva|
|00005900| 6c 2e 69 76 61 6c 75 65 | 20 3d 20 28 64 69 67 69 |l.ivalue| = (digi|
|00005910| 74 20 2a 20 28 28 46 55 | 4c 4c 29 20 2a 68 33 2b |t * ((FU|LL) *h3+|
|00005920| 2b 29 29 0a 58 09 09 09 | 09 09 2b 20 28 28 46 55 |+)).X...|..+ ((FU|
|00005930| 4c 4c 29 20 2a 68 31 29 | 20 2b 20 28 28 46 55 4c |LL) *h1)| + ((FUL|
|00005940| 4c 29 20 73 69 76 61 6c | 2e 73 69 68 69 67 68 29 |L) sival|.sihigh)|
|00005950| 3b 0a 58 09 09 09 09 2a | 68 31 2b 2b 20 3d 20 73 |;.X....*|h1++ = s|
|00005960| 69 76 61 6c 2e 73 69 6c | 6f 77 3b 0a 58 09 09 09 |ival.sil|ow;.X...|
|00005970| 09 73 69 76 61 6c 2e 69 | 76 61 6c 75 65 20 3d 20 |.sival.i|value = |
|00005980| 28 64 69 67 69 74 20 2a | 20 28 28 46 55 4c 4c 29 |(digit *| ((FULL)|
|00005990| 20 2a 68 33 2b 2b 29 29 | 0a 58 09 09 09 09 09 2b | *h3++))|.X.....+|
|000059a0| 20 28 28 46 55 4c 4c 29 | 20 2a 68 31 29 20 2b 20 | ((FULL)| *h1) + |
|000059b0| 28 28 46 55 4c 4c 29 20 | 73 69 76 61 6c 2e 73 69 |((FULL) |sival.si|
|000059c0| 68 69 67 68 29 3b 0a 58 | 09 09 09 09 2a 68 31 2b |high);.X|....*h1+|
|000059d0| 2b 20 3d 20 73 69 76 61 | 6c 2e 73 69 6c 6f 77 3b |+ = siva|l.silow;|
|000059e0| 0a 58 09 09 09 09 63 61 | 72 72 79 20 3d 20 73 69 |.X....ca|rry = si|
|000059f0| 76 61 6c 2e 73 69 68 69 | 67 68 3b 0a 58 09 09 09 |val.sihi|gh;.X...|
|00005a00| 7d 0a 58 09 09 09 77 68 | 69 6c 65 20 28 6c 65 6e |}.X...wh|ile (len|
|00005a10| 31 2d 2d 29 20 7b 0a 58 | 09 09 09 09 73 69 76 61 |1--) {.X|....siva|
|00005a20| 6c 2e 69 76 61 6c 75 65 | 20 3d 20 28 64 69 67 69 |l.ivalue| = (digi|
|00005a30| 74 20 2a 20 28 28 46 55 | 4c 4c 29 20 2a 68 33 2b |t * ((FU|LL) *h3+|
|00005a40| 2b 29 29 0a 58 09 09 09 | 09 09 2b 20 28 28 46 55 |+)).X...|..+ ((FU|
|00005a50| 4c 4c 29 20 2a 68 31 29 | 20 2b 20 63 61 72 72 79 |LL) *h1)| + carry|
|00005a60| 3b 0a 58 09 09 09 09 2a | 68 31 2b 2b 20 3d 20 73 |;.X....*|h1++ = s|
|00005a70| 69 76 61 6c 2e 73 69 6c | 6f 77 3b 0a 58 09 09 09 |ival.sil|ow;.X...|
|00005a80| 09 63 61 72 72 79 20 3d | 20 73 69 76 61 6c 2e 73 |.carry =| sival.s|
|00005a90| 69 68 69 67 68 3b 0a 58 | 09 09 09 7d 0a 58 09 09 |ihigh;.X|...}.X..|
|00005aa0| 09 77 68 69 6c 65 20 28 | 63 61 72 72 79 29 20 7b |.while (|carry) {|
|00005ab0| 0a 58 09 09 09 09 73 69 | 76 61 6c 2e 69 76 61 6c |.X....si|val.ival|
|00005ac0| 75 65 20 3d 20 28 28 46 | 55 4c 4c 29 20 2a 68 31 |ue = ((F|ULL) *h1|
|00005ad0| 29 20 2b 20 63 61 72 72 | 79 3b 0a 58 09 09 09 09 |) + carr|y;.X....|
|00005ae0| 2a 68 31 2b 2b 20 3d 20 | 73 69 76 61 6c 2e 73 69 |*h1++ = |sival.si|
|00005af0| 6c 6f 77 3b 0a 58 09 09 | 09 09 63 61 72 72 79 20 |low;.X..|..carry |
|00005b00| 3d 20 73 69 76 61 6c 2e | 73 69 68 69 67 68 3b 0a |= sival.|sihigh;.|
|00005b10| 58 09 09 09 7d 0a 58 09 | 09 7d 0a 58 0a 58 09 09 |X...}.X.|.}.X.X..|
|00005b20| 2f 2a 0a 58 09 09 20 2a | 20 4e 6f 77 20 64 6f 75 |/*.X.. *| Now dou|
|00005b30| 62 6c 65 20 74 68 65 20 | 72 65 73 75 6c 74 2e 0a |ble the |result..|
|00005b40| 58 09 09 20 2a 20 54 68 | 65 72 65 20 69 73 20 6e |X.. * Th|ere is n|
|00005b50| 6f 20 66 69 6e 61 6c 20 | 63 61 72 72 79 20 74 6f |o final |carry to|
|00005b60| 20 77 6f 72 72 79 20 61 | 62 6f 75 74 20 62 65 63 | worry a|bout bec|
|00005b70| 61 75 73 65 20 77 65 0a | 58 09 09 20 2a 20 68 61 |ause we.|X.. * ha|
|00005b80| 6e 64 6c 65 20 61 6c 6c | 20 64 69 67 69 74 73 20 |ndle all| digits |
|00005b90| 6f 66 20 74 68 65 20 72 | 65 73 75 6c 74 20 77 68 |of the r|esult wh|
|00005ba0| 69 63 68 20 6d 75 73 74 | 20 66 69 74 2e 0a 58 09 |ich must| fit..X.|
|00005bb0| 09 20 2a 2f 0a 58 09 09 | 63 61 72 72 79 20 3d 20 |. */.X..|carry = |
|00005bc0| 30 3b 0a 58 09 09 68 64 | 20 3d 20 61 6e 73 3b 0a |0;.X..hd| = ans;.|
|00005bd0| 58 09 09 6c 65 6e 20 3d | 20 73 69 7a 65 74 6f 74 |X..len =| sizetot|
|00005be0| 61 6c 3b 0a 58 09 09 77 | 68 69 6c 65 20 28 6c 65 |al;.X..w|hile (le|
|00005bf0| 6e 2d 2d 29 20 7b 0a 58 | 09 09 09 64 69 67 69 74 |n--) {.X|...digit|
|00005c00| 20 3d 20 28 28 46 55 4c | 4c 29 20 2a 68 64 29 3b | = ((FUL|L) *hd);|
|00005c10| 0a 58 09 09 09 73 69 76 | 61 6c 2e 69 76 61 6c 75 |.X...siv|al.ivalu|
|00005c20| 65 20 3d 20 64 69 67 69 | 74 20 2b 20 64 69 67 69 |e = digi|t + digi|
|00005c30| 74 20 2b 20 63 61 72 72 | 79 3b 0a 58 09 09 09 2a |t + carr|y;.X...*|
|00005c40| 68 64 2b 2b 20 3d 20 73 | 69 76 61 6c 2e 73 69 6c |hd++ = s|ival.sil|
|00005c50| 6f 77 3b 0a 58 09 09 09 | 63 61 72 72 79 20 3d 20 |ow;.X...|carry = |
|00005c60| 73 69 76 61 6c 2e 73 69 | 68 69 67 68 3b 0a 58 09 |sival.si|high;.X.|
|00005c70| 09 7d 0a 58 0a 58 09 09 | 2f 2a 0a 58 09 09 20 2a |.}.X.X..|/*.X.. *|
|00005c80| 20 4e 6f 77 20 61 64 64 | 20 69 6e 20 74 68 65 20 | Now add| in the |
|00005c90| 73 71 75 61 72 65 73 20 | 6f 66 20 65 61 63 68 20 |squares |of each |
|00005ca0| 68 61 6c 66 77 6f 72 64 | 2e 0a 58 09 09 20 2a 2f |halfword|..X.. */|
|00005cb0| 0a 58 09 09 63 61 72 72 | 79 20 3d 20 30 3b 0a 58 |.X..carr|y = 0;.X|
|00005cc0| 09 09 68 64 20 3d 20 61 | 6e 73 3b 0a 58 09 09 68 |..hd = a|ns;.X..h|
|00005cd0| 33 20 3d 20 76 70 3b 0a | 58 09 09 6c 65 6e 20 3d |3 = vp;.|X..len =|
|00005ce0| 20 73 69 7a 65 3b 0a 58 | 09 09 77 68 69 6c 65 20 | size;.X|..while |
|00005cf0| 28 6c 65 6e 2d 2d 29 20 | 7b 0a 58 09 09 09 64 69 |(len--) |{.X...di|
|00005d00| 67 69 74 20 3d 20 28 28 | 46 55 4c 4c 29 20 2a 68 |git = ((|FULL) *h|
|00005d10| 33 2b 2b 29 3b 0a 58 09 | 09 09 73 69 76 61 6c 2e |3++);.X.|..sival.|
|00005d20| 69 76 61 6c 75 65 20 3d | 20 64 69 67 69 74 20 2a |ivalue =| digit *|
|00005d30| 20 64 69 67 69 74 20 2b | 20 28 28 46 55 4c 4c 29 | digit +| ((FULL)|
|00005d40| 20 2a 68 64 29 20 2b 20 | 63 61 72 72 79 3b 0a 58 | *hd) + |carry;.X|
|00005d50| 09 09 09 2a 68 64 2b 2b | 20 3d 20 73 69 76 61 6c |...*hd++| = sival|
|00005d60| 2e 73 69 6c 6f 77 3b 0a | 58 09 09 09 63 61 72 72 |.silow;.|X...carr|
|00005d70| 79 20 3d 20 73 69 76 61 | 6c 2e 73 69 68 69 67 68 |y = siva|l.sihigh|
|00005d80| 3b 0a 58 09 09 09 73 69 | 76 61 6c 2e 69 76 61 6c |;.X...si|val.ival|
|00005d90| 75 65 20 3d 20 28 28 46 | 55 4c 4c 29 20 2a 68 64 |ue = ((F|ULL) *hd|
|00005da0| 29 20 2b 20 63 61 72 72 | 79 3b 0a 58 09 09 09 2a |) + carr|y;.X...*|
|00005db0| 68 64 2b 2b 20 3d 20 73 | 69 76 61 6c 2e 73 69 6c |hd++ = s|ival.sil|
|00005dc0| 6f 77 3b 0a 58 09 09 09 | 63 61 72 72 79 20 3d 20 |ow;.X...|carry = |
|00005dd0| 73 69 76 61 6c 2e 73 69 | 68 69 67 68 3b 0a 58 09 |sival.si|high;.X.|
|00005de0| 09 7d 0a 58 09 09 77 68 | 69 6c 65 20 28 63 61 72 |.}.X..wh|ile (car|
|00005df0| 72 79 29 20 7b 0a 58 09 | 09 09 73 69 76 61 6c 2e |ry) {.X.|..sival.|
|00005e00| 69 76 61 6c 75 65 20 3d | 20 28 28 46 55 4c 4c 29 |ivalue =| ((FULL)|
|00005e10| 20 2a 68 64 29 20 2b 20 | 63 61 72 72 79 3b 0a 58 | *hd) + |carry;.X|
|00005e20| 09 09 09 2a 68 64 2b 2b | 20 3d 20 73 69 76 61 6c |...*hd++| = sival|
|00005e30| 2e 73 69 6c 6f 77 3b 0a | 58 09 09 09 63 61 72 72 |.silow;.|X...carr|
|00005e40| 79 20 3d 20 73 69 76 61 | 6c 2e 73 69 68 69 67 68 |y = siva|l.sihigh|
|00005e50| 3b 0a 58 09 09 7d 0a 58 | 0a 58 09 09 2f 2a 0a 58 |;.X..}.X|.X../*.X|
|00005e60| 09 09 20 2a 20 46 69 6e | 61 6c 6c 79 20 72 65 74 |.. * Fin|ally ret|
|00005e70| 75 72 6e 20 74 68 65 20 | 73 69 7a 65 20 6f 66 20 |urn the |size of |
|00005e80| 74 68 65 20 72 65 73 75 | 6c 74 2e 0a 58 09 09 20 |the resu|lt..X.. |
|00005e90| 2a 2f 0a 58 09 09 6c 65 | 6e 20 3d 20 73 69 7a 65 |*/.X..le|n = size|
|00005ea0| 74 6f 74 61 6c 3b 0a 58 | 09 09 68 64 20 3d 20 26 |total;.X|..hd = &|
|00005eb0| 61 6e 73 5b 6c 65 6e 20 | 2d 20 31 5d 3b 0a 58 09 |ans[len |- 1];.X.|
|00005ec0| 09 77 68 69 6c 65 20 28 | 28 2a 68 64 20 3d 3d 20 |.while (|(*hd == |
|00005ed0| 30 29 20 26 26 20 28 6c | 65 6e 20 3e 20 31 29 29 |0) && (l|en > 1))|
|00005ee0| 20 7b 0a 58 09 09 09 6c | 65 6e 2d 2d 3b 0a 58 09 | {.X...l|en--;.X.|
|00005ef0| 09 09 68 64 2d 2d 3b 0a | 58 09 09 7d 0a 58 09 09 |..hd--;.|X..}.X..|
|00005f00| 72 65 74 75 72 6e 20 6c | 65 6e 3b 0a 58 09 7d 0a |return l|en;.X.}.|
|00005f10| 58 0a 58 09 2f 2a 0a 58 | 09 20 2a 20 54 68 65 20 |X.X./*.X|. * The |
|00005f20| 6e 75 6d 62 65 72 20 74 | 6f 20 62 65 20 73 71 75 |number t|o be squ|
|00005f30| 61 72 65 64 20 69 73 20 | 6c 61 72 67 65 2e 0a 58 |ared is |large..X|
|00005f40| 09 20 2a 20 41 6c 6c 6f | 63 61 74 65 20 74 65 6d |. * Allo|cate tem|
|00005f50| 70 6f 72 61 72 79 20 73 | 70 61 63 65 20 61 6e 64 |porary s|pace and|
|00005f60| 20 64 65 74 65 72 6d 69 | 6e 65 20 74 68 65 20 73 | determi|ne the s|
|00005f70| 69 7a 65 73 20 61 6e 64 | 0a 58 09 20 2a 20 70 6f |izes and|.X. * po|
|00005f80| 73 69 74 69 6f 6e 73 20 | 6f 66 20 74 68 65 20 76 |sitions |of the v|
|00005f90| 61 6c 75 65 73 20 74 6f | 20 62 65 20 63 61 6c 63 |alues to| be calc|
|00005fa0| 75 6c 61 74 65 64 2e 0a | 58 09 20 2a 2f 0a 58 09 |ulated..|X. */.X.|
|00005fb0| 74 65 6d 70 20 3d 20 74 | 65 6d 70 62 75 66 3b 0a |temp = t|empbuf;.|
|00005fc0| 58 09 74 65 6d 70 62 75 | 66 20 2b 3d 20 28 33 20 |X.tempbu|f += (3 |
|00005fd0| 2a 20 28 73 69 7a 65 20 | 2b 20 31 29 20 2f 20 32 |* (size |+ 1) / 2|
|00005fe0| 29 3b 0a 58 0a 58 09 73 | 69 7a 65 41 20 3d 20 73 |);.X.X.s|izeA = s|
|00005ff0| 69 7a 65 20 2f 20 32 3b | 0a 58 09 73 69 7a 65 42 |ize / 2;|.X.sizeB|
|00006000| 20 3d 20 73 69 7a 65 20 | 2d 20 73 69 7a 65 41 3b | = size |- sizeA;|
|00006010| 0a 58 09 73 68 69 66 74 | 20 3d 20 73 69 7a 65 42 |.X.shift| = sizeB|
|00006020| 3b 0a 58 09 62 61 73 65 | 41 20 3d 20 76 70 20 2b |;.X.base|A = vp +|
|00006030| 20 73 69 7a 65 42 3b 0a | 58 09 62 61 73 65 42 20 | sizeB;.|X.baseB |
|00006040| 3d 20 76 70 3b 0a 58 09 | 62 61 73 65 41 41 20 3d |= vp;.X.|baseAA =|
|00006050| 20 26 61 6e 73 5b 73 68 | 69 66 74 20 2a 20 32 5d | &ans[sh|ift * 2]|
|00006060| 3b 0a 58 09 62 61 73 65 | 42 42 20 3d 20 61 6e 73 |;.X.base|BB = ans|
|00006070| 3b 0a 58 09 62 61 73 65 | 41 41 42 42 20 3d 20 74 |;.X.base|AABB = t|
|00006080| 65 6d 70 3b 0a 58 09 62 | 61 73 65 41 42 20 3d 20 |emp;.X.b|aseAB = |
|00006090| 74 65 6d 70 3b 0a 58 09 | 62 61 73 65 41 42 41 42 |temp;.X.|baseABAB|
|000060a0| 20 3d 20 26 74 65 6d 70 | 5b 73 68 69 66 74 5d 3b | = &temp|[shift];|
|000060b0| 0a 58 0a 58 09 2f 2a 0a | 58 09 20 2a 20 54 72 69 |.X.X./*.|X. * Tri|
|000060c0| 6d 20 74 68 65 20 73 65 | 63 6f 6e 64 20 6e 75 6d |m the se|cond num|
|000060d0| 62 65 72 20 6f 66 20 6c | 65 61 64 69 6e 67 20 7a |ber of l|eading z|
|000060e0| 65 72 6f 65 73 2e 0a 58 | 09 20 2a 2f 0a 58 09 68 |eroes..X|. */.X.h|
|000060f0| 64 20 3d 20 26 62 61 73 | 65 42 5b 73 69 7a 65 42 |d = &bas|eB[sizeB|
|00006100| 20 2d 20 31 5d 3b 0a 58 | 09 77 68 69 6c 65 20 28 | - 1];.X|.while (|
|00006110| 28 2a 68 64 20 3d 3d 20 | 30 29 20 26 26 20 28 73 |(*hd == |0) && (s|
|00006120| 69 7a 65 42 20 3e 20 31 | 29 29 20 7b 0a 58 09 09 |izeB > 1|)) {.X..|
|00006130| 73 69 7a 65 42 2d 2d 3b | 0a 58 09 09 68 64 2d 2d |sizeB--;|.X..hd--|
|00006140| 3b 0a 58 09 7d 0a 58 0a | 58 09 2f 2a 0a 58 09 20 |;.X.}.X.|X./*.X. |
|00006150| 2a 20 4e 6f 77 20 74 6f | 20 70 72 6f 63 65 65 64 |* Now to| proceed|
|00006160| 20 74 6f 20 63 61 6c 63 | 75 6c 61 74 65 20 74 68 | to calc|ulate th|
|00006170| 65 20 72 65 73 75 6c 74 | 20 75 73 69 6e 67 20 74 |e result| using t|
|00006180| 68 65 20 66 6f 72 6d 75 | 6c 61 2e 0a 58 09 20 2a |he formu|la..X. *|
|00006190| 09 28 41 2a 53 2b 42 29 | 5e 32 20 3d 20 28 53 5e |.(A*S+B)|^2 = (S^|
|000061a0| 32 2b 53 29 2a 41 5e 32 | 20 2b 20 28 53 2b 31 29 |2+S)*A^2| + (S+1)|
|000061b0| 2a 42 5e 32 20 2d 20 53 | 2a 28 41 2d 42 29 5e 32 |*B^2 - S|*(A-B)^2|
|000061c0| 2e 0a 58 09 20 2a 20 54 | 68 65 20 73 74 65 70 73 |..X. * T|he steps|
|000061d0| 20 61 72 65 20 64 6f 6e | 65 20 69 6e 20 74 68 65 | are don|e in the|
|000061e0| 20 66 6f 6c 6c 6f 77 69 | 6e 67 20 6f 72 64 65 72 | followi|ng order|
|000061f0| 3a 0a 58 09 20 2a 09 53 | 5e 32 2a 41 5e 32 20 2b |:.X. *.S|^2*A^2 +|
|00006200| 20 42 5e 32 0a 58 09 20 | 2a 09 41 5e 32 20 2b 20 | B^2.X. |*.A^2 + |
|00006210| 42 5e 32 0a 58 09 20 2a | 09 28 53 5e 32 2b 53 29 |B^2.X. *|.(S^2+S)|
|00006220| 2a 41 5e 32 20 2b 20 28 | 53 2b 31 29 2a 42 5e 32 |*A^2 + (|S+1)*B^2|
|00006230| 0a 58 09 20 2a 09 28 41 | 2d 42 29 5e 32 0a 58 09 |.X. *.(A|-B)^2.X.|
|00006240| 20 2a 09 28 53 5e 32 2b | 53 29 2a 41 5e 32 20 2b | *.(S^2+|S)*A^2 +|
|00006250| 20 28 53 2b 31 29 2a 42 | 5e 32 20 2d 20 53 2a 28 | (S+1)*B|^2 - S*(|
|00006260| 41 2d 42 29 5e 32 2e 0a | 58 09 20 2a 0a 58 09 20 |A-B)^2..|X. *.X. |
|00006270| 2a 20 42 65 67 69 6e 20 | 62 79 20 66 6f 72 6d 69 |* Begin |by formi|
|00006280| 6e 67 20 74 68 65 20 73 | 71 75 61 72 65 73 20 6f |ng the s|quares o|
|00006290| 66 20 74 77 6f 20 74 68 | 65 20 68 61 6c 66 73 20 |f two th|e halfs |
|000062a0| 63 6f 6e 63 61 74 65 6e | 61 74 65 64 0a 58 09 20 |concaten|ated.X. |
|000062b0| 2a 20 74 6f 67 65 74 68 | 65 72 20 69 6e 20 74 68 |* togeth|er in th|
|000062c0| 65 20 66 69 6e 61 6c 20 | 72 65 73 75 6c 74 20 6c |e final |result l|
|000062d0| 6f 63 61 74 69 6f 6e 2e | 20 20 4d 61 6b 65 20 73 |ocation.|  Make s|
|000062e0| 75 72 65 20 74 68 61 74 | 20 74 68 65 0a 58 09 20 |ure that| the.X. |
|000062f0| 2a 20 68 69 67 68 65 73 | 74 20 77 6f 72 64 73 20 |* highes|t words |
|00006300| 6f 66 20 74 68 65 20 72 | 65 73 75 6c 74 73 20 61 |of the r|esults a|
|00006310| 72 65 20 7a 65 72 6f 2e | 0a 58 09 20 2a 2f 0a 58 |re zero.|.X. */.X|
|00006320| 09 73 69 7a 65 42 42 20 | 3d 20 64 6f 73 71 75 61 |.sizeBB |= dosqua|
|00006330| 72 65 28 62 61 73 65 42 | 2c 20 73 69 7a 65 42 2c |re(baseB|, sizeB,|
|00006340| 20 62 61 73 65 42 42 29 | 3b 0a 58 09 68 64 20 3d | baseBB)|;.X.hd =|
|00006350| 20 26 62 61 73 65 42 42 | 5b 73 69 7a 65 42 42 5d | &baseBB|[sizeBB]|
|00006360| 3b 0a 58 09 6c 65 6e 20 | 3d 20 73 68 69 66 74 20 |;.X.len |= shift |
|00006370| 2a 20 32 20 2d 20 73 69 | 7a 65 42 42 3b 0a 58 09 |* 2 - si|zeBB;.X.|
|00006380| 77 68 69 6c 65 20 28 6c | 65 6e 2d 2d 29 0a 58 09 |while (l|en--).X.|
|00006390| 09 2a 68 64 2b 2b 20 3d | 20 30 3b 0a 58 0a 58 09 |.*hd++ =| 0;.X.X.|
|000063a0| 73 69 7a 65 41 41 20 3d | 20 64 6f 73 71 75 61 72 |sizeAA =| dosquar|
|000063b0| 65 28 62 61 73 65 41 2c | 20 73 69 7a 65 41 2c 20 |e(baseA,| sizeA, |
|000063c0| 62 61 73 65 41 41 29 3b | 0a 58 09 68 64 20 3d 20 |baseAA);|.X.hd = |
|000063d0| 26 62 61 73 65 41 41 5b | 73 69 7a 65 41 41 5d 3b |&baseAA[|sizeAA];|
|000063e0| 0a 58 09 6c 65 6e 20 3d | 20 73 69 7a 65 74 6f 74 |.X.len =| sizetot|
|000063f0| 61 6c 20 2d 20 73 68 69 | 66 74 20 2a 20 32 20 2d |al - shi|ft * 2 -|
+--------+-------------------------+-------------------------+--------+--------+