/ Usenet 1994 January / usenetsourcesnewsgroupsinfomagicjanuary1994.iso / sources / unix / volume10 / cbw / part08

< prev

next >

MacOS 8.1

Win98

DOS

browse contents

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 50 61 74 68 3a 20 73 65 | 69 73 6d 6f 21 75 75 6e |Path: se|ismo!uun|
|00000010| 65 74 21 72 73 0a 46 72 | 6f 6d 3a 20 72 73 40 75 |et!rs.Fr|om: rs@u|
|00000020| 75 6e 65 74 2e 55 55 2e | 4e 45 54 20 28 52 69 63 |unet.UU.|NET (Ric|
|00000030| 68 20 53 61 6c 7a 29 0a | 4e 65 77 73 67 72 6f 75 |h Salz).|Newsgrou|
|00000040| 70 73 3a 20 63 6f 6d 70 | 2e 73 6f 75 72 63 65 73 |ps: comp|.sources|
|00000050| 2e 75 6e 69 78 0a 53 75 | 62 6a 65 63 74 3a 20 76 |.unix.Su|bject: v|
|00000060| 31 30 69 30 30 38 3a 20 | 20 43 72 79 70 74 20 42 |10i008: | Crypt B|
|00000070| 72 65 61 6b 65 72 27 73 | 20 57 6f 72 6b 62 65 6e |reaker's| Workben|
|00000080| 63 68 2c 20 50 61 72 74 | 30 38 2f 31 31 0a 4d 65 |ch, Part|08/11.Me|
|00000090| 73 73 61 67 65 2d 49 44 | 3a 20 3c 33 39 30 40 75 |ssage-ID|: <390@u|
|000000a0| 75 6e 65 74 2e 55 55 2e | 4e 45 54 3e 0a 44 61 74 |unet.UU.|NET>.Dat|
|000000b0| 65 3a 20 31 39 20 4a 75 | 6e 20 38 37 20 30 33 3a |e: 19 Ju|n 87 03:|
|000000c0| 31 37 3a 31 37 20 47 4d | 54 0a 4f 72 67 61 6e 69 |17:17 GM|T.Organi|
|000000d0| 7a 61 74 69 6f 6e 3a 20 | 55 55 4e 45 54 20 43 6f |zation: |UUNET Co|
|000000e0| 6d 6d 75 6e 69 63 61 74 | 69 6f 6e 73 20 53 65 72 |mmunicat|ions Ser|
|000000f0| 76 69 63 65 73 2c 20 41 | 72 6c 69 6e 67 74 6f 6e |vices, A|rlington|
|00000100| 2c 20 56 41 0a 4c 69 6e | 65 73 3a 20 31 37 30 35 |, VA.Lin|es: 1705|
|00000110| 0a 41 70 70 72 6f 76 65 | 64 3a 20 72 73 40 75 75 |.Approve|d: rs@uu|
|00000120| 6e 65 74 2e 75 75 2e 6e | 65 74 0a 0a 53 75 62 6d |net.uu.n|et..Subm|
|00000130| 69 74 74 65 64 20 62 79 | 3a 20 52 6f 62 65 72 74 |itted by|: Robert|
|00000140| 20 57 2e 20 42 61 6c 64 | 77 69 6e 20 3c 42 41 4c | W. Bald|win <BAL|
|00000150| 44 57 49 4e 40 58 58 2e | 4c 43 53 2e 4d 49 54 2e |DWIN@XX.|LCS.MIT.|
|00000160| 45 44 55 3e 0a 4d 6f 64 | 2e 73 6f 75 72 63 65 73 |EDU>.Mod|.sources|
|00000170| 3a 20 56 6f 6c 75 6d 65 | 20 31 30 2c 20 49 73 73 |: Volume| 10, Iss|
|00000180| 75 65 20 38 0a 41 72 63 | 68 69 76 65 2d 6e 61 6d |ue 8.Arc|hive-nam|
|00000190| 65 3a 20 63 62 77 2f 50 | 61 72 74 30 38 0a 0a 0a |e: cbw/P|art08...|
|000001a0| 23 21 20 2f 62 69 6e 2f | 73 68 0a 23 20 54 68 69 |#! /bin/|sh.# Thi|
|000001b0| 73 20 69 73 20 61 20 73 | 68 65 6c 6c 20 61 72 63 |s is a s|hell arc|
|000001c0| 68 69 76 65 2e 20 20 52 | 65 6d 6f 76 65 20 61 6e |hive.  R|emove an|
|000001d0| 79 74 68 69 6e 67 20 62 | 65 66 6f 72 65 20 74 68 |ything b|efore th|
|000001e0| 69 73 20 6c 69 6e 65 2c | 20 74 68 65 6e 20 75 6e |is line,| then un|
|000001f0| 70 61 63 6b 0a 23 20 69 | 74 20 62 79 20 73 61 76 |pack.# i|t by sav|
|00000200| 69 6e 67 20 69 74 20 69 | 6e 74 6f 20 61 20 66 69 |ing it i|nto a fi|
|00000210| 6c 65 20 61 6e 64 20 74 | 79 70 69 6e 67 20 22 73 |le and t|yping "s|
|00000220| 68 20 66 69 6c 65 22 2e | 20 20 54 6f 20 6f 76 65 |h file".|  To ove|
|00000230| 72 77 72 69 74 65 20 65 | 78 69 73 74 69 6e 67 0a |rwrite e|xisting.|
|00000240| 23 20 66 69 6c 65 73 2c | 20 74 79 70 65 20 22 73 |# files,| type "s|
|00000250| 68 20 66 69 6c 65 20 2d | 63 22 2e 20 20 59 6f 75 |h file -|c".  You|
|00000260| 20 63 61 6e 20 61 6c 73 | 6f 20 66 65 65 64 20 74 | can als|o feed t|
|00000270| 68 69 73 20 61 73 20 73 | 74 61 6e 64 61 72 64 20 |his as s|tandard |
|00000280| 69 6e 70 75 74 20 76 69 | 61 0a 23 20 75 6e 73 68 |input vi|a.# unsh|
|00000290| 61 72 2c 20 6f 72 20 62 | 79 20 74 79 70 69 6e 67 |ar, or b|y typing|
|000002a0| 20 22 73 68 20 3c 66 69 | 6c 65 22 2c 20 65 2e 67 | "sh <fi|le", e.g|
|000002b0| 2e 2e 20 20 49 66 20 74 | 68 69 73 20 61 72 63 68 |..  If t|his arch|
|000002c0| 69 76 65 20 69 73 20 63 | 6f 6d 70 6c 65 74 65 2c |ive is c|omplete,|
|000002d0| 20 79 6f 75 0a 23 20 77 | 69 6c 6c 20 73 65 65 20 | you.# w|ill see |
|000002e0| 74 68 65 20 66 6f 6c 6c | 6f 77 69 6e 67 20 6d 65 |the foll|owing me|
|000002f0| 73 73 61 67 65 20 61 74 | 20 74 68 65 20 65 6e 64 |ssage at| the end|
|00000300| 3a 0a 23 09 09 22 45 6e | 64 20 6f 66 20 61 72 63 |:.#.."En|d of arc|
|00000310| 68 69 76 65 20 38 20 28 | 6f 66 20 31 31 29 2e 22 |hive 8 (|of 11)."|
|00000320| 0a 23 20 43 6f 6e 74 65 | 6e 74 73 3a 20 20 73 74 |.# Conte|nts:  st|
|00000330| 61 74 73 2e 63 20 74 65 | 72 6d 69 6e 61 6c 2e 63 |ats.c te|rminal.c|
|00000340| 0a 23 20 57 72 61 70 70 | 65 64 20 62 79 20 72 73 |.# Wrapp|ed by rs|
|00000350| 40 75 75 6e 65 74 20 6f | 6e 20 57 65 64 20 4a 75 |@uunet o|n Wed Ju|
|00000360| 6e 20 31 37 20 31 38 3a | 31 37 3a 32 36 20 31 39 |n 17 18:|17:26 19|
|00000370| 38 37 0a 50 41 54 48 3d | 2f 62 69 6e 3a 2f 75 73 |87.PATH=|/bin:/us|
|00000380| 72 2f 62 69 6e 3a 2f 75 | 73 72 2f 75 63 62 20 3b |r/bin:/u|sr/ucb ;|
|00000390| 20 65 78 70 6f 72 74 20 | 50 41 54 48 0a 69 66 20 | export |PATH.if |
|000003a0| 74 65 73 74 20 2d 66 20 | 73 74 61 74 73 2e 63 20 |test -f |stats.c |
|000003b0| 2d 61 20 22 24 7b 31 7d | 22 20 21 3d 20 22 2d 63 |-a "${1}|" != "-c|
|000003c0| 22 20 3b 20 74 68 65 6e | 20 0a 20 20 65 63 68 6f |" ; then| .  echo|
|000003d0| 20 73 68 61 72 3a 20 57 | 69 6c 6c 20 6e 6f 74 20 | shar: W|ill not |
|000003e0| 6f 76 65 72 2d 77 72 69 | 74 65 20 65 78 69 73 74 |over-wri|te exist|
|000003f0| 69 6e 67 20 66 69 6c 65 | 20 5c 22 73 74 61 74 73 |ing file| \"stats|
|00000400| 2e 63 5c 22 0a 65 6c 73 | 65 0a 65 63 68 6f 20 73 |.c\".els|e.echo s|
|00000410| 68 61 72 3a 20 45 78 74 | 72 61 63 74 69 6e 67 20 |har: Ext|racting |
|00000420| 5c 22 73 74 61 74 73 2e | 63 5c 22 20 5c 28 31 38 |\"stats.|c\" \(18|
|00000430| 37 39 36 20 63 68 61 72 | 61 63 74 65 72 73 5c 29 |796 char|acters\)|
|00000440| 0a 73 65 64 20 22 73 2f | 5e 58 2f 2f 22 20 3e 73 |.sed "s/|^X//" >s|
|00000450| 74 61 74 73 2e 63 20 3c | 3c 27 45 4e 44 5f 4f 46 |tats.c <|<'END_OF|
|00000460| 5f 73 74 61 74 73 2e 63 | 27 0a 58 2f 2a 0a 58 20 |_stats.c|'.X/*.X |
|00000470| 2a 20 53 74 61 74 69 73 | 74 69 63 73 20 72 6f 75 |* Statis|tics rou|
|00000480| 74 69 6e 65 73 2e 0a 58 | 20 2a 0a 58 20 2a 20 52 |tines..X| *.X * R|
|00000490| 6f 62 65 72 74 20 57 2e | 20 42 61 6c 64 77 69 6e |obert W.| Baldwin|
|000004a0| 2c 20 4a 61 6e 75 61 72 | 79 20 31 39 38 35 2e 0a |, Januar|y 1985..|
|000004b0| 58 20 2a 20 53 63 6f 72 | 69 6e 67 20 62 61 73 65 |X * Scor|ing base|
|000004c0| 64 20 6f 6e 20 6c 65 74 | 74 65 72 20 70 61 69 72 |d on let|ter pair|
|000004d0| 73 20 61 64 64 65 64 20 | 62 79 20 42 6f 62 20 42 |s added |by Bob B|
|000004e0| 61 6c 64 77 69 6e 20 4d | 61 79 20 31 39 38 35 2e |aldwin M|ay 1985.|
|000004f0| 0a 58 20 2a 2f 0a 58 0a | 58 0a 58 23 69 6e 63 6c |.X */.X.|X.X#incl|
|00000500| 75 64 65 09 3c 73 74 64 | 69 6f 2e 68 3e 0a 58 23 |ude.<std|io.h>.X#|
|00000510| 69 6e 63 6c 75 64 65 09 | 3c 6d 61 74 68 2e 68 3e |include.|<math.h>|
|00000520| 0a 58 23 69 6e 63 6c 75 | 64 65 09 22 77 69 6e 64 |.X#inclu|de."wind|
|00000530| 6f 77 2e 68 22 0a 58 23 | 69 6e 63 6c 75 64 65 09 |ow.h".X#|include.|
|00000540| 22 73 70 65 63 73 2e 68 | 22 0a 58 23 69 6e 63 6c |"specs.h|".X#incl|
|00000550| 75 64 65 09 22 63 69 70 | 68 65 72 2e 68 22 0a 58 |ude."cip|her.h".X|
|00000560| 0a 58 0a 58 23 64 65 66 | 69 6e 65 20 53 54 41 4e |.X.X#def|ine STAN|
|00000570| 44 41 4c 4f 4e 45 09 46 | 41 4c 53 45 0a 58 0a 58 |DALONE.F|ALSE.X.X|
|00000580| 0a 58 0a 58 2f 2a 20 47 | 6c 6f 62 61 6c 73 20 2a |.X.X/* G|lobals *|
|00000590| 2f 0a 58 69 6e 74 09 09 | 73 74 61 74 73 31 6c 6f |/.Xint..|stats1lo|
|000005a0| 61 64 65 64 20 3d 20 46 | 41 4c 53 45 3b 09 2f 2a |aded = F|ALSE;./*|
|000005b0| 20 54 72 75 65 20 69 66 | 20 6c 65 74 74 65 72 20 | True if| letter |
|000005c0| 73 74 61 74 73 20 6c 6f | 61 64 65 64 2e 20 2a 2f |stats lo|aded. */|
|000005d0| 0a 58 63 68 61 72 09 2a | 6c 65 74 74 65 72 73 74 |.Xchar.*|letterst|
|000005e0| 61 74 73 3b 09 09 09 2f | 2a 20 46 69 6c 65 6e 61 |ats;.../|* Filena|
|000005f0| 6d 65 20 74 6f 20 66 69 | 6e 64 20 73 69 6e 67 6c |me to fi|nd singl|
|00000600| 65 20 6c 65 74 74 65 72 | 20 63 6f 75 6e 74 73 2e |e letter| counts.|
|00000610| 20 2a 2f 0a 58 69 6e 74 | 09 09 73 74 61 74 73 32 | */.Xint|..stats2|
|00000620| 6c 6f 61 64 65 64 20 3d | 20 46 41 4c 53 45 3b 09 |loaded =| FALSE;.|
|00000630| 2f 2a 20 54 72 75 65 20 | 69 66 20 6c 65 74 74 65 |/* True |if lette|
|00000640| 72 20 70 61 69 72 20 73 | 74 61 74 73 20 6c 6f 61 |r pair s|tats loa|
|00000650| 64 65 64 2e 20 2a 2f 0a | 58 63 68 61 72 09 2a 62 |ded. */.|Xchar.*b|
|00000660| 69 67 72 61 6d 73 74 61 | 74 73 3b 09 09 09 2f 2a |igramsta|ts;.../*|
|00000670| 20 46 69 6c 65 6e 61 6d | 65 20 74 6f 20 66 69 6e | Filenam|e to fin|
|00000680| 64 20 6c 65 74 74 65 72 | 20 70 61 69 72 20 63 6f |d letter| pair co|
|00000690| 75 6e 74 73 2e 20 2a 2f | 0a 58 0a 58 0a 58 2f 2a |unts. */|.X.X.X/*|
|000006a0| 20 54 68 69 73 20 61 72 | 72 61 79 20 63 6f 6e 74 | This ar|ray cont|
|000006b0| 61 69 6e 73 20 74 68 65 | 20 73 69 6e 67 6c 65 20 |ains the| single |
|000006c0| 6c 65 74 74 65 72 20 66 | 72 65 71 75 65 6e 63 69 |letter f|requenci|
|000006d0| 65 73 2e 0a 58 20 2a 20 | 54 68 61 74 20 69 73 2c |es..X * |That is,|
|000006e0| 20 70 72 6f 62 5b 69 5d | 20 69 73 20 74 68 65 20 | prob[i]| is the |
|000006f0| 70 72 6f 62 61 62 69 6c | 69 74 79 20 74 68 61 74 |probabil|ity that|
|00000700| 20 61 20 72 61 6e 64 6f | 6d 6c 79 20 73 65 6c 65 | a rando|mly sele|
|00000710| 63 74 65 64 0a 58 20 2a | 20 70 6c 61 69 6e 74 65 |cted.X *| plainte|
|00000720| 78 74 20 63 68 61 72 61 | 63 74 65 72 20 68 61 73 |xt chara|cter has|
|00000730| 20 61 73 63 69 69 20 76 | 61 6c 75 65 20 69 2e 0a | ascii v|alue i..|
|00000740| 58 20 2a 2f 0a 58 66 6c | 6f 61 74 09 70 72 6f 62 |X */.Xfl|oat.prob|
|00000750| 5b 4d 41 58 43 48 41 52 | 2b 31 5d 3b 0a 58 66 6c |[MAXCHAR|+1];.Xfl|
|00000760| 6f 61 74 09 70 6d 65 61 | 6e 2c 20 70 76 61 72 3b |oat.pmea|n, pvar;|
|00000770| 09 09 09 2f 2a 20 4d 65 | 61 6e 20 61 6e 64 20 76 |.../* Me|an and v|
|00000780| 61 72 69 61 6e 63 65 20 | 6f 66 20 61 62 6f 76 65 |ariance |of above|
|00000790| 2e 20 2a 2f 0a 58 0a 58 | 0a 58 2f 2a 20 54 68 69 |. */.X.X|.X/* Thi|
|000007a0| 73 20 61 72 72 61 79 20 | 63 6f 6e 74 61 69 6e 73 |s array |contains|
|000007b0| 20 74 68 65 20 62 61 73 | 65 20 74 65 6e 20 6c 6f | the bas|e ten lo|
|000007c0| 67 61 72 69 74 68 6d 73 | 20 6f 66 20 74 68 65 20 |garithms| of the |
|000007d0| 73 69 6e 67 6c 65 0a 58 | 20 2a 20 6c 65 74 74 65 |single.X| * lette|
|000007e0| 72 20 66 72 65 71 75 65 | 6e 63 69 65 73 20 28 70 |r freque|ncies (p|
|000007f0| 72 6f 62 61 62 69 6c 69 | 74 69 65 73 29 20 6f 66 |robabili|ties) of|
|00000800| 20 41 53 43 49 49 20 63 | 68 61 72 61 63 74 65 72 | ASCII c|haracter|
|00000810| 73 2e 0a 58 20 2a 20 54 | 68 65 20 66 72 65 71 75 |s..X * T|he frequ|
|00000820| 65 6e 63 69 65 73 20 61 | 72 65 20 62 65 74 77 65 |encies a|re betwe|
|00000830| 65 6e 20 30 20 61 6e 64 | 20 31 2c 20 73 6f 20 74 |en 0 and| 1, so t|
|00000840| 68 65 20 65 6e 74 65 72 | 69 65 73 20 69 6e 20 74 |he enter|ies in t|
|00000850| 68 65 0a 58 20 2a 20 74 | 61 62 6c 65 20 61 72 65 |he.X * t|able are|
|00000860| 20 62 65 74 77 65 65 6e | 20 6d 69 6e 75 73 20 69 | between| minus i|
|00000870| 6e 66 69 6e 69 74 79 20 | 61 6e 64 20 30 2e 0a 58 |nfinity |and 0..X|
|00000880| 20 2a 20 54 68 65 20 6c | 6f 67 20 6f 66 20 61 20 | * The l|og of a |
|00000890| 76 61 6e 69 73 68 69 6e | 67 6c 79 20 73 6d 61 6c |vanishin|gly smal|
|000008a0| 6c 20 66 72 65 71 75 65 | 6e 63 79 20 69 73 20 72 |l freque|ncy is r|
|000008b0| 65 70 72 65 73 65 6e 74 | 65 64 20 61 73 0a 58 20 |epresent|ed as.X |
|000008c0| 2a 20 7a 65 72 6f 20 72 | 61 74 68 65 72 20 74 68 |* zero r|ather th|
|000008d0| 61 6e 20 73 6f 6d 65 20 | 6c 61 72 67 65 20 6e 65 |an some |large ne|
|000008e0| 67 61 74 69 76 65 20 6e | 75 6d 62 65 72 2e 0a 58 |gative n|umber..X|
|000008f0| 20 2a 20 54 68 65 20 65 | 78 70 65 63 74 65 64 20 | * The e|xpected |
|00000900| 76 61 6c 75 65 20 6f 66 | 20 6c 6f 67 70 72 6f 62 |value of| logprob|
|00000910| 5b 63 5d 20 66 6f 72 20 | 61 20 72 61 6e 64 6f 6d |[c] for |a random|
|00000920| 6c 79 20 73 65 6c 65 63 | 74 65 64 0a 58 20 2a 20 |ly selec|ted.X * |
|00000930| 63 68 61 72 61 63 74 65 | 72 2c 20 63 2c 20 69 73 |characte|r, c, is|
|00000940| 20 67 69 76 65 6e 20 62 | 79 20 6c 6f 67 6d 65 61 | given b|y logmea|
|00000950| 6e 2e 20 20 54 68 65 20 | 76 61 72 69 61 6e 63 65 |n.  The |variance|
|00000960| 20 6f 66 20 6c 6f 67 70 | 72 6f 62 5b 63 5d 0a 58 | of logp|rob[c].X|
|00000970| 20 2a 20 69 73 20 67 69 | 76 65 6e 20 62 79 20 6c | * is gi|ven by l|
|00000980| 6f 67 76 61 72 2e 20 20 | 54 68 65 20 73 74 61 6e |ogvar.  |The stan|
|00000990| 64 61 72 64 20 64 65 76 | 69 61 74 69 6f 6e 20 69 |dard dev|iation i|
|000009a0| 73 20 69 6e 20 6c 6f 67 | 73 64 2e 0a 58 20 2a 2f |s in log|sd..X */|
|000009b0| 0a 58 66 6c 6f 61 74 09 | 6c 6f 67 70 72 6f 62 5b |.Xfloat.|logprob[|
|000009c0| 4d 41 58 43 48 41 52 2b | 31 5d 3b 0a 58 66 6c 6f |MAXCHAR+|1];.Xflo|
|000009d0| 61 74 09 6c 6f 67 6d 65 | 61 6e 2c 20 6c 6f 67 76 |at.logme|an, logv|
|000009e0| 61 72 2c 20 6c 6f 67 73 | 64 3b 0a 58 0a 58 0a 58 |ar, logs|d;.X.X.X|
|000009f0| 2f 2a 20 54 68 69 73 20 | 61 72 72 61 79 20 63 6f |/* This |array co|
|00000a00| 6e 74 61 69 6e 73 20 74 | 68 65 20 62 69 67 72 61 |ntains t|he bigra|
|00000a10| 6d 20 28 6c 65 74 74 65 | 72 20 70 61 69 72 29 20 |m (lette|r pair) |
|00000a20| 66 72 65 71 75 65 6e 63 | 69 65 73 2e 0a 58 20 2a |frequenc|ies..X *|
|00000a30| 20 54 68 61 74 20 69 73 | 2c 20 62 69 70 72 6f 62 | That is|, biprob|
|00000a40| 5b 69 5d 5b 6a 5d 20 69 | 73 20 74 68 65 20 70 72 |[i][j] i|s the pr|
|00000a50| 6f 62 61 62 69 6c 69 74 | 79 20 74 68 61 74 20 61 |obabilit|y that a|
|00000a60| 20 72 61 6e 64 6f 6d 6c | 79 20 73 65 6c 65 63 74 | randoml|y select|
|00000a70| 65 64 0a 58 20 2a 20 70 | 61 69 72 20 6f 66 20 61 |ed.X * p|air of a|
|00000a80| 64 6a 61 63 65 6e 74 20 | 70 6c 61 69 6e 74 65 78 |djacent |plaintex|
|00000a90| 74 20 63 68 61 72 61 63 | 74 65 72 73 20 69 73 20 |t charac|ters is |
|00000aa0| 41 69 2c 20 41 6a 2e 20 | 20 57 68 65 72 65 20 41 |Ai, Aj. | Where A|
|00000ab0| 69 0a 58 20 2a 20 69 73 | 20 74 68 65 20 69 74 68 |i.X * is| the ith|
|00000ac0| 20 6c 65 74 74 65 72 20 | 6f 66 20 74 68 65 20 61 | letter |of the a|
|00000ad0| 6c 70 68 61 62 65 74 20 | 28 30 20 3d 20 27 41 27 |lphabet |(0 = 'A'|
|00000ae0| 20 6f 72 20 27 61 27 2c | 20 61 6e 64 0a 58 20 2a | or 'a',| and.X *|
|00000af0| 20 32 36 20 3d 20 73 70 | 61 63 65 20 6f 72 20 6f | 26 = sp|ace or o|
|00000b00| 74 68 65 72 20 6e 6f 6e | 2d 61 6c 70 68 61 62 65 |ther non|-alphabe|
|00000b10| 74 69 63 29 2e 0a 58 20 | 2a 20 45 76 65 6e 74 75 |tic)..X |* Eventu|
|00000b20| 61 6c 6c 79 20 74 68 69 | 73 20 77 69 6c 6c 20 62 |ally thi|s will b|
|00000b30| 65 20 67 65 6e 65 72 61 | 6c 69 7a 65 64 20 74 6f |e genera|lized to|
|00000b40| 20 69 6e 63 6c 75 64 65 | 20 61 6e 20 61 72 62 69 | include| an arbi|
|00000b50| 74 72 61 72 79 0a 58 20 | 2a 20 63 68 61 72 61 63 |trary.X |* charac|
|00000b60| 74 65 72 20 74 72 61 6e | 73 6c 61 74 69 6f 6e 20 |ter tran|slation |
|00000b70| 74 61 62 6c 65 20 74 6f | 20 68 61 6e 64 6c 65 20 |table to| handle |
|00000b80| 70 75 6e 63 74 75 61 74 | 69 6f 6e 20 61 6e 64 20 |punctuat|ion and |
|00000b90| 74 6f 20 61 6c 6c 6f 77 | 0a 58 20 2a 20 67 72 6f |to allow|.X * gro|
|00000ba0| 75 70 73 20 6f 66 20 63 | 68 61 72 61 63 74 65 72 |ups of c|haracter|
|00000bb0| 73 20 73 75 63 68 20 61 | 73 20 28 2c 20 3c 2c 20 |s such a|s (, <, |
|00000bc0| 7b 2c 20 61 6e 64 20 5b | 20 74 6f 20 62 65 20 74 |{, and [| to be t|
|00000bd0| 72 65 61 74 65 64 20 74 | 68 65 20 73 61 6d 65 2e |reated t|he same.|
|00000be0| 0a 58 20 2a 20 0a 58 20 | 2a 20 54 68 65 20 61 72 |.X * .X |* The ar|
|00000bf0| 72 61 79 20 73 6c 62 69 | 70 72 6f 62 20 69 73 20 |ray slbi|prob is |
|00000c00| 74 68 65 20 73 69 6e 67 | 6c 65 20 6c 65 74 74 65 |the sing|le lette|
|00000c10| 72 20 66 72 65 71 75 65 | 6e 63 69 65 73 20 74 61 |r freque|ncies ta|
|00000c20| 6b 65 6e 20 66 72 6f 6d | 20 74 68 65 0a 58 20 2a |ken from| the.X *|
|00000c30| 20 73 61 6d 65 20 73 6f | 75 72 63 65 73 20 61 73 | same so|urces as|
|00000c40| 20 62 69 70 72 6f 62 5b | 5d 5b 5d 2e 0a 58 20 2a | biprob[|][]..X *|
|00000c50| 2f 0a 58 66 6c 6f 61 74 | 09 62 69 70 72 6f 62 5b |/.Xfloat|.biprob[|
|00000c60| 4d 58 42 49 49 4e 44 45 | 58 5d 5b 4d 58 42 49 49 |MXBIINDE|X][MXBII|
|00000c70| 4e 44 45 58 5d 3b 0a 58 | 66 6c 6f 61 74 09 73 6c |NDEX];.X|float.sl|
|00000c80| 62 69 70 72 6f 62 5b 4d | 58 42 49 49 4e 44 45 58 |biprob[M|XBIINDEX|
|00000c90| 5d 3b 0a 58 0a 58 0a 58 | 2f 2a 20 54 68 65 20 61 |];.X.X.X|/* The a|
|00000ca0| 72 72 61 79 20 69 73 20 | 75 73 65 64 20 74 6f 20 |rray is |used to |
|00000cb0| 6d 61 70 20 66 72 6f 6d | 20 37 2d 62 69 74 20 61 |map from| 7-bit a|
|00000cc0| 73 63 69 69 20 74 6f 20 | 69 6e 64 69 63 65 73 20 |scii to |indices |
|00000cd0| 69 6e 20 74 68 65 20 62 | 69 70 72 6f 62 0a 58 20 |in the b|iprob.X |
|00000ce0| 2a 20 61 6e 64 20 72 65 | 6c 61 74 65 64 20 61 72 |* and re|lated ar|
|00000cf0| 72 61 79 73 2e 20 20 54 | 68 65 20 76 61 72 69 61 |rays.  T|he varia|
|00000d00| 62 6c 65 20 6e 62 69 63 | 68 61 72 73 20 69 73 20 |ble nbic|hars is |
|00000d10| 73 65 74 20 74 6f 20 74 | 68 65 20 6e 65 78 74 20 |set to t|he next |
|00000d20| 69 6e 64 65 78 0a 58 20 | 2a 20 74 6f 20 75 73 65 |index.X |* to use|
|00000d30| 20 69 6e 20 74 68 65 20 | 62 69 70 72 6f 62 20 61 | in the |biprob a|
|00000d40| 72 72 61 79 2e 0a 58 20 | 2a 2f 0a 58 69 6e 74 09 |rray..X |*/.Xint.|
|00000d50| 09 63 68 61 72 5f 62 69 | 6d 61 70 5b 4d 41 58 43 |.char_bi|map[MAXC|
|00000d60| 48 41 52 2b 31 5d 3b 0a | 58 69 6e 74 09 09 6e 62 |HAR+1];.|Xint..nb|
|00000d70| 69 63 68 61 72 73 3b 0a | 58 0a 58 0a 58 2f 2a 20 |ichars;.|X.X.X/* |
|00000d80| 54 68 69 73 20 61 72 72 | 61 79 20 63 6f 6e 74 61 |This arr|ay conta|
|00000d90| 69 6e 73 20 74 68 65 20 | 62 61 73 65 20 74 65 6e |ins the |base ten|
|00000da0| 20 6c 6f 67 61 72 69 74 | 68 6d 73 20 6f 66 20 74 | logarit|hms of t|
|00000db0| 68 65 20 6c 65 74 74 65 | 72 20 70 61 69 72 0a 58 |he lette|r pair.X|
|00000dc0| 20 2a 20 66 72 65 71 75 | 65 6e 63 69 65 73 20 28 | * frequ|encies (|
|00000dd0| 62 69 70 72 6f 62 5b 5d | 5b 5d 29 2e 0a 58 20 2a |biprob[]|[])..X *|
|00000de0| 20 54 68 65 20 66 72 65 | 71 75 65 6e 63 69 65 73 | The fre|quencies|
|00000df0| 20 61 72 65 20 62 65 74 | 77 65 65 6e 20 30 20 61 | are bet|ween 0 a|
|00000e00| 6e 64 20 31 2c 20 73 6f | 20 74 68 65 20 65 6e 74 |nd 1, so| the ent|
|00000e10| 65 72 69 65 73 20 69 6e | 20 74 68 65 0a 58 20 2a |eries in| the.X *|
|00000e20| 20 74 61 62 6c 65 20 61 | 72 65 20 62 65 74 77 65 | table a|re betwe|
|00000e30| 65 6e 20 6d 69 6e 75 73 | 20 69 6e 66 69 6e 69 74 |en minus| infinit|
|00000e40| 79 20 61 6e 64 20 30 2e | 0a 58 20 2a 20 54 68 65 |y and 0.|.X * The|
|00000e50| 20 6c 6f 67 20 6f 66 20 | 61 20 76 61 6e 69 73 68 | log of |a vanish|
|00000e60| 69 6e 67 6c 79 20 73 6d | 61 6c 6c 20 66 72 65 71 |ingly sm|all freq|
|00000e70| 75 65 6e 63 79 20 69 73 | 20 72 65 70 72 65 73 65 |uency is| represe|
|00000e80| 6e 74 65 64 20 61 73 0a | 58 20 2a 20 7a 65 72 6f |nted as.|X * zero|
|00000e90| 20 72 61 74 68 65 72 20 | 74 68 61 6e 20 73 6f 6d | rather |than som|
|00000ea0| 65 20 6c 61 72 67 65 20 | 6e 65 67 61 74 69 76 65 |e large |negative|
|00000eb0| 20 6e 75 6d 62 65 72 2e | 0a 58 20 2a 20 54 68 65 | number.|.X * The|
|00000ec0| 20 65 78 70 65 63 74 65 | 64 20 76 61 6c 75 65 20 | expecte|d value |
|00000ed0| 6f 66 20 62 69 6c 6f 67 | 70 72 6f 62 5b 63 5d 20 |of bilog|prob[c] |
|00000ee0| 66 6f 72 20 61 20 72 61 | 6e 64 6f 6d 6c 79 20 73 |for a ra|ndomly s|
|00000ef0| 65 6c 65 63 74 65 64 0a | 58 20 2a 20 63 68 61 72 |elected.|X * char|
|00000f00| 61 63 74 65 72 2c 20 63 | 2c 20 69 73 20 67 69 76 |acter, c|, is giv|
|00000f10| 65 6e 20 62 79 20 62 69 | 6c 6f 67 6d 65 61 6e 2e |en by bi|logmean.|
|00000f20| 20 20 54 68 65 20 76 61 | 72 69 61 6e 63 65 20 6f |  The va|riance o|
|00000f30| 66 20 62 69 6c 6f 67 70 | 72 6f 62 5b 63 5d 0a 58 |f bilogp|rob[c].X|
|00000f40| 20 2a 20 69 73 20 67 69 | 76 65 6e 20 62 79 20 62 | * is gi|ven by b|
|00000f50| 69 6c 6f 67 76 61 72 2e | 20 20 54 68 65 20 73 74 |ilogvar.|  The st|
|00000f60| 61 6e 64 61 72 64 20 64 | 65 76 69 61 74 69 6f 6e |andard d|eviation|
|00000f70| 20 69 73 20 69 6e 20 62 | 69 6c 6f 67 73 64 2e 0a | is in b|ilogsd..|
|00000f80| 58 20 2a 2f 0a 58 66 6c | 6f 61 74 09 62 69 6c 6f |X */.Xfl|oat.bilo|
|00000f90| 67 70 72 6f 62 5b 4d 58 | 42 49 49 4e 44 45 58 5d |gprob[MX|BIINDEX]|
|00000fa0| 5b 4d 58 42 49 49 4e 44 | 45 58 5d 3b 0a 58 0a 58 |[MXBIIND|EX];.X.X|
|00000fb0| 0a 58 2f 2a 20 54 68 69 | 73 20 76 65 63 74 6f 72 |.X/* Thi|s vector|
|00000fc0| 20 63 6f 6e 74 61 69 6e | 73 20 74 68 65 20 62 61 | contain|s the ba|
|00000fd0| 73 65 20 74 65 6e 20 6c | 6f 67 61 72 69 74 68 6d |se ten l|ogarithm|
|00000fe0| 73 20 6f 66 20 74 68 65 | 20 73 69 6e 67 6c 65 20 |s of the| single |
|00000ff0| 6c 65 74 74 65 72 0a 58 | 20 2a 20 66 72 65 71 75 |letter.X| * frequ|
|00001000| 65 6e 63 69 65 73 20 77 | 68 69 63 68 20 61 72 65 |encies w|hich are|
|00001010| 20 64 65 72 69 76 65 64 | 20 66 72 6f 6d 20 62 69 | derived| from bi|
|00001020| 70 72 6f 62 5b 5d 5b 5d | 2e 0a 58 20 2a 20 54 68 |prob[][]|..X * Th|
|00001030| 65 79 20 61 72 65 20 75 | 73 65 64 20 74 6f 20 63 |ey are u|sed to c|
|00001040| 6f 6d 70 75 74 65 20 74 | 68 65 20 6c 6f 67 20 6f |ompute t|he log o|
|00001050| 66 20 74 68 65 20 63 6f | 6e 64 69 74 69 6f 6e 61 |f the co|nditiona|
|00001060| 6c 20 70 72 6f 62 61 62 | 69 6c 69 74 69 65 73 0a |l probab|ilities.|
|00001070| 58 20 2a 20 6f 66 20 6c | 65 74 74 65 72 20 70 61 |X * of l|etter pa|
|00001080| 69 72 73 2c 20 67 69 76 | 65 6e 20 61 20 6b 6e 6f |irs, giv|en a kno|
|00001090| 77 6e 20 76 61 6c 75 65 | 20 66 6f 72 20 65 69 74 |wn value| for eit|
|000010a0| 68 65 72 20 74 68 65 20 | 66 69 72 73 74 20 6f 72 |her the |first or|
|000010b0| 0a 58 20 2a 20 73 65 63 | 6f 6e 64 20 63 68 61 72 |.X * sec|ond char|
|000010c0| 61 63 74 65 72 2e 0a 58 | 20 2a 20 53 70 65 63 69 |acter..X| * Speci|
|000010d0| 66 69 63 61 6c 6c 79 3a | 20 6c 6f 67 28 20 70 72 |fically:| log( pr|
|000010e0| 6f 62 28 58 59 20 67 69 | 76 65 6e 20 58 3d 41 69 |ob(XY gi|ven X=Ai|
|000010f0| 29 20 29 20 65 71 75 61 | 6c 73 0a 58 20 2a 20 20 |) ) equa|ls.X *  |
|00001100| 20 20 6c 6f 67 28 20 70 | 72 6f 62 28 58 59 29 20 |  log( p|rob(XY) |
|00001110| 2f 20 70 72 6f 62 28 41 | 69 29 20 29 20 77 68 69 |/ prob(A|i) ) whi|
|00001120| 63 68 20 65 71 75 61 6c | 73 20 0a 58 20 2a 20 20 |ch equal|s .X *  |
|00001130| 20 20 62 69 6c 6f 67 70 | 72 6f 62 5b 58 5d 5b 59 |  bilogp|rob[X][Y|
|00001140| 5d 20 2d 20 73 6c 6c 6f | 67 70 72 6f 62 5b 41 69 |] - sllo|gprob[Ai|
|00001150| 5d 2e 0a 58 20 2a 2f 0a | 58 66 6c 6f 61 74 09 73 |]..X */.|Xfloat.s|
|00001160| 6c 6c 6f 67 70 72 6f 62 | 5b 4d 58 42 49 49 4e 44 |llogprob|[MXBIIND|
|00001170| 45 58 5d 3b 0a 58 0a 58 | 0a 58 2f 2a 20 54 68 65 |EX];.X.X|.X/* The|
|00001180| 20 73 63 6f 72 69 6e 67 | 20 66 75 6e 63 74 69 6f | scoring| functio|
|00001190| 6e 20 74 68 61 74 20 75 | 73 65 73 20 6c 65 74 74 |n that u|ses lett|
|000011a0| 65 72 20 70 61 69 72 20 | 66 72 65 71 75 65 6e 63 |er pair |frequenc|
|000011b0| 69 65 73 20 69 73 20 62 | 61 73 65 64 0a 58 20 2a |ies is b|ased.X *|
|000011c0| 20 6f 6e 20 61 20 73 74 | 61 74 69 73 74 69 63 20 | on a st|atistic |
|000011d0| 74 68 61 74 20 68 61 73 | 20 61 20 63 6f 6d 70 75 |that has| a compu|
|000011e0| 74 61 62 6c 65 20 6d 65 | 61 6e 20 61 6e 64 20 76 |table me|an and v|
|000011f0| 61 72 69 61 6e 63 65 20 | 28 61 6e 64 0a 58 20 2a |ariance |(and.X *|
|00001200| 20 73 74 61 6e 64 61 72 | 64 20 64 65 76 69 61 74 | standar|d deviat|
|00001210| 69 6f 6e 29 2e 20 20 54 | 68 65 20 61 72 65 20 73 |ion).  T|he are s|
|00001220| 74 6f 72 65 64 20 69 6e | 20 74 68 65 20 66 6f 6c |tored in| the fol|
|00001230| 6c 6f 77 69 6e 67 20 76 | 61 72 69 61 62 6c 65 73 |lowing v|ariables|
|00001240| 2e 0a 58 20 2a 2f 0a 58 | 66 6c 6f 61 74 09 73 63 |..X */.X|float.sc|
|00001250| 6f 72 65 32 5f 6d 65 61 | 6e 2c 20 73 63 6f 72 65 |ore2_mea|n, score|
|00001260| 32 5f 76 61 72 2c 20 73 | 63 6f 72 65 32 5f 73 64 |2_var, s|core2_sd|
|00001270| 2c 20 73 63 6f 72 65 32 | 5f 73 63 61 6c 65 3b 0a |, score2|_scale;.|
|00001280| 58 66 6c 6f 61 74 09 73 | 63 6f 72 65 31 5f 6d 65 |Xfloat.s|core1_me|
|00001290| 61 6e 2c 20 73 63 6f 72 | 65 31 5f 76 61 72 2c 20 |an, scor|e1_var, |
|000012a0| 73 63 6f 72 65 31 5f 73 | 64 2c 20 73 63 6f 72 65 |score1_s|d, score|
|000012b0| 31 5f 73 63 61 6c 65 3b | 0a 58 0a 58 0a 58 23 69 |1_scale;|.X.X.X#i|
|000012c0| 66 20 53 54 41 4e 44 41 | 4c 4f 4e 45 0a 58 23 64 |f STANDA|LONE.X#d|
|000012d0| 65 66 69 6e 65 09 66 69 | 6c 65 6e 61 6d 65 09 09 |efine.fi|lename..|
|000012e0| 22 2f 75 73 72 2f 62 61 | 6c 64 77 69 6e 2f 45 63 |"/usr/ba|ldwin/Ec|
|000012f0| 72 79 70 74 2f 6d 73 73 | 2d 62 69 67 72 61 6d 2e |rypt/mss|-bigram.|
|00001300| 73 74 61 74 73 22 0a 58 | 6d 61 69 6e 28 29 0a 58 |stats".X|main().X|
|00001310| 7b 0a 58 09 46 49 4c 45 | 09 2a 69 6e 70 3b 0a 58 |{.X.FILE|.*inp;.X|
|00001320| 0a 58 09 6c 6f 61 64 5f | 32 73 74 61 74 73 5f 66 |.X.load_|2stats_f|
|00001330| 72 6f 6d 28 66 69 6c 65 | 6e 61 6d 65 29 3b 0a 58 |rom(file|name);.X|
|00001340| 09 70 72 69 6e 74 5f 32 | 73 74 61 74 73 28 73 74 |.print_2|stats(st|
|00001350| 64 6f 75 74 29 3b 0a 58 | 7d 0a 58 23 65 6e 64 69 |dout);.X|}.X#endi|
|00001360| 66 0a 58 0a 58 0a 58 0a | 58 2f 2a 20 53 63 6f 72 |f.X.X.X.|X/* Scor|
|00001370| 65 20 74 68 65 20 67 69 | 76 65 6e 20 70 6c 61 69 |e the gi|ven plai|
|00001380| 6e 74 65 78 74 20 62 6c | 6f 63 6b 2e 20 20 52 65 |ntext bl|ock.  Re|
|00001390| 74 75 72 6e 73 20 61 20 | 66 6c 6f 61 74 69 6e 67 |turns a |floating|
|000013a0| 20 70 6f 69 6e 74 20 6e | 75 6d 62 65 72 2e 0a 58 | point n|umber..X|
|000013b0| 20 2a 20 46 6f 72 20 6e | 6f 77 20 61 20 73 74 75 | * For n|ow a stu|
|000013c0| 64 2e 0a 58 20 2a 2f 0a | 58 66 6c 6f 61 74 09 73 |d..X */.|Xfloat.s|
|000013d0| 63 6f 72 65 28 70 62 6c | 6f 63 6b 29 0a 58 69 6e |core(pbl|ock).Xin|
|000013e0| 74 09 09 70 62 6c 6f 63 | 6b 5b 5d 3b 0a 58 7b 0a |t..pbloc|k[];.X{.|
|000013f0| 58 09 69 6e 74 09 09 70 | 63 68 61 72 3b 0a 58 09 |X.int..p|char;.X.|
|00001400| 69 6e 74 09 09 69 3b 0a | 58 09 66 6c 6f 61 74 09 |int..i;.|X.float.|
|00001410| 73 63 6f 72 65 3b 0a 58 | 0a 58 09 73 63 6f 72 65 |score;.X|.X.score|
|00001420| 20 3d 20 30 2e 30 3b 0a | 58 09 66 6f 72 20 28 69 | = 0.0;.|X.for (i|
|00001430| 20 3d 20 30 20 3b 20 69 | 20 3c 20 42 4c 4f 43 4b | = 0 ; i| < BLOCK|
|00001440| 53 49 5a 45 20 3b 20 69 | 2b 2b 29 20 20 7b 0a 58 |SIZE ; i|++)  {.X|
|00001450| 09 09 70 63 68 61 72 20 | 3d 20 70 62 6c 6f 63 6b |..pchar |= pblock|
|00001460| 5b 69 5d 3b 0a 58 09 09 | 69 66 20 28 70 63 68 61 |[i];.X..|if (pcha|
|00001470| 72 20 3d 3d 20 2d 31 29 | 20 20 63 6f 6e 74 69 6e |r == -1)|  contin|
|00001480| 75 65 3b 0a 58 09 09 69 | 66 20 28 70 63 68 61 72 |ue;.X..i|f (pchar|
|00001490| 20 3d 3d 20 27 20 27 29 | 09 09 09 7b 73 63 6f 72 | == ' ')|...{scor|
|000014a0| 65 20 2b 3d 20 30 2e 31 | 3b 7d 0a 58 09 09 65 6c |e += 0.1|;}.X..el|
|000014b0| 73 65 20 69 66 20 28 6c | 6c 65 74 74 65 72 28 70 |se if (l|letter(p|
|000014c0| 63 68 61 72 29 29 09 7b | 73 63 6f 72 65 20 2b 3d |char)).{|score +=|
|000014d0| 20 30 2e 31 3b 7d 0a 58 | 09 09 65 6c 73 65 20 69 | 0.1;}.X|..else i|
|000014e0| 66 20 28 75 6c 65 74 74 | 65 72 28 70 63 68 61 72 |f (ulett|er(pchar|
|000014f0| 29 29 09 7b 73 63 6f 72 | 65 20 2b 3d 20 30 2e 30 |)).{scor|e += 0.0|
|00001500| 35 3b 7d 0a 58 09 09 65 | 6c 73 65 20 69 66 20 28 |5;}.X..e|lse if (|
|00001510| 70 72 69 6e 74 61 62 6c | 65 28 70 63 68 61 72 29 |printabl|e(pchar)|
|00001520| 29 09 7b 73 63 6f 72 65 | 20 2b 3d 20 30 2e 30 32 |).{score| += 0.02|
|00001530| 3b 7d 0a 58 09 09 65 6c | 73 65 20 69 66 20 28 70 |;}.X..el|se if (p|
|00001540| 63 68 61 72 20 3d 3d 20 | 27 5c 74 27 20 7c 7c 20 |char == |'\t' || |
|00001550| 70 63 68 61 72 20 3d 3d | 20 27 5c 6e 27 20 7c 7c |pchar ==| '\n' |||
|00001560| 20 70 63 68 61 72 20 3d | 3d 20 27 5c 66 27 29 0a | pchar =|= '\f').|
|00001570| 58 09 09 09 09 09 09 09 | 09 09 7b 73 63 6f 72 65 |X.......|..{score|
|00001580| 20 2b 3d 20 30 2e 30 35 | 3b 7d 0a 58 09 09 65 6c | += 0.05|;}.X..el|
|00001590| 73 65 20 69 66 20 28 27 | 30 27 20 3c 3d 20 70 63 |se if ('|0' <= pc|
|000015a0| 68 61 72 20 26 26 20 70 | 63 68 61 72 20 3c 3d 20 |har && p|char <= |
|000015b0| 27 39 27 29 0a 58 09 09 | 09 09 09 09 09 09 09 7b |'9').X..|.......{|
|000015c0| 73 63 6f 72 65 20 2b 3d | 20 30 2e 30 35 3b 7d 0a |score +=| 0.05;}.|
|000015d0| 58 09 09 65 6c 73 65 20 | 09 09 09 09 09 09 7b 73 |X..else |......{s|
|000015e0| 63 6f 72 65 20 2d 3d 20 | 30 2e 34 3b 7d 0a 58 09 |core -= |0.4;}.X.|
|000015f0| 09 7d 0a 58 09 72 65 74 | 75 72 6e 28 73 63 6f 72 |.}.X.ret|urn(scor|
|00001600| 65 29 3b 0a 58 7d 0a 58 | 0a 58 0a 58 2f 2a 20 53 |e);.X}.X|.X.X/* S|
|00001610| 63 6f 72 65 20 61 20 76 | 65 63 74 6f 72 20 6f 66 |core a v|ector of|
|00001620| 20 69 6e 74 65 67 65 72 | 73 20 74 68 61 74 20 72 | integer|s that r|
|00001630| 65 70 72 65 73 65 6e 74 | 20 63 68 61 72 61 63 74 |epresent| charact|
|00001640| 65 72 73 2e 0a 58 20 2a | 20 54 68 65 20 76 65 63 |ers..X *| The vec|
|00001650| 74 6f 72 20 69 73 20 74 | 65 72 6d 69 6e 61 74 65 |tor is t|erminate|
|00001660| 64 20 62 79 20 61 20 76 | 61 6c 75 65 20 6f 66 20 |d by a v|alue of |
|00001670| 4e 4f 4e 45 2e 0a 58 20 | 2a 20 54 68 65 20 72 65 |NONE..X |* The re|
|00001680| 74 75 72 6e 65 64 20 73 | 63 6f 72 65 20 69 73 20 |turned s|core is |
|00001690| 74 68 65 20 6e 75 6d 62 | 65 72 20 6f 66 20 73 74 |the numb|er of st|
|000016a0| 61 6e 64 61 72 64 20 64 | 65 76 69 61 74 69 6f 6e |andard d|eviation|
|000016b0| 73 0a 58 20 2a 20 74 68 | 61 74 20 74 68 65 20 6f |s.X * th|at the o|
|000016c0| 62 73 65 72 76 65 64 20 | 73 74 61 74 69 73 74 69 |bserved |statisti|
|000016d0| 63 73 20 64 69 66 66 65 | 72 73 20 66 72 6f 6d 20 |cs diffe|rs from |
|000016e0| 69 74 73 20 65 78 70 65 | 63 74 65 64 20 76 61 6c |its expe|cted val|
|000016f0| 75 65 2e 0a 58 20 2a 20 | 53 63 6f 72 65 73 20 61 |ue..X * |Scores a|
|00001700| 72 65 20 70 6f 73 69 74 | 69 76 65 20 77 69 74 68 |re posit|ive with|
|00001710| 20 6c 6f 77 20 73 63 6f | 72 65 73 20 62 65 69 6e | low sco|res bein|
|00001720| 67 20 62 65 74 74 65 72 | 2e 0a 58 20 2a 20 41 20 |g better|..X * A |
|00001730| 6e 65 67 61 74 69 76 65 | 20 73 63 6f 72 65 20 69 |negative| score i|
|00001740| 6e 64 69 63 61 74 65 73 | 20 61 6e 20 69 6d 70 6f |ndicates| an impo|
|00001750| 73 73 69 62 6c 65 20 70 | 6c 61 69 6e 74 65 78 74 |ssible p|laintext|
|00001760| 20 76 61 6c 75 65 2e 0a | 58 20 2a 2f 0a 58 66 6c | value..|X */.Xfl|
|00001770| 6f 61 74 09 70 76 65 63 | 5f 31 73 63 6f 72 65 28 |oat.pvec|_1score(|
|00001780| 70 76 65 63 29 0a 58 69 | 6e 74 09 2a 70 76 65 63 |pvec).Xi|nt.*pvec|
|00001790| 3b 0a 58 7b 0a 58 09 69 | 6e 74 09 09 69 3b 0a 58 |;.X{.X.i|nt..i;.X|
|000017a0| 09 69 6e 74 09 09 63 3b | 0a 58 09 66 6c 6f 61 74 |.int..c;|.X.float|
|000017b0| 09 74 6d 70 2c 20 73 75 | 6d 2c 20 63 6f 75 6e 74 |.tmp, su|m, count|
|000017c0| 2c 20 73 63 6f 72 65 3b | 0a 58 0a 58 09 69 66 20 |, score;|.X.X.if |
|000017d0| 28 21 73 74 61 74 73 31 | 6c 6f 61 64 65 64 29 20 |(!stats1|loaded) |
|000017e0| 20 7b 0a 58 09 09 6c 6f | 61 64 5f 31 73 74 61 74 | {.X..lo|ad_1stat|
|000017f0| 73 5f 66 72 6f 6d 28 6c | 65 74 74 65 72 73 74 61 |s_from(l|ettersta|
|00001800| 74 73 29 3b 0a 58 09 09 | 7d 0a 58 0a 58 09 63 6f |ts);.X..|}.X.X.co|
|00001810| 75 6e 74 20 3d 20 30 2e | 30 3b 0a 58 09 73 75 6d |unt = 0.|0;.X.sum|
|00001820| 20 3d 20 30 2e 30 3b 0a | 58 09 77 68 69 6c 65 20 | = 0.0;.|X.while |
|00001830| 28 2a 70 76 65 63 20 21 | 3d 20 4e 4f 4e 45 29 20 |(*pvec !|= NONE) |
|00001840| 20 7b 0a 58 09 09 63 6f | 75 6e 74 20 2b 3d 20 31 | {.X..co|unt += 1|
|00001850| 2e 30 3b 0a 58 09 09 63 | 20 3d 20 2a 70 76 65 63 |.0;.X..c| = *pvec|
|00001860| 2b 2b 3b 0a 58 09 09 69 | 66 20 28 63 20 21 3d 20 |++;.X..i|f (c != |
|00001870| 63 20 26 20 43 48 41 52 | 4d 41 53 4b 29 20 20 72 |c & CHAR|MASK)  r|
|00001880| 65 74 75 72 6e 28 2d 31 | 2e 30 29 3b 0a 58 09 09 |eturn(-1|.0);.X..|
|00001890| 74 6d 70 20 3d 20 6c 6f | 67 70 72 6f 62 5b 63 20 |tmp = lo|gprob[c |
|000018a0| 26 20 43 48 41 52 4d 41 | 53 4b 5d 3b 0a 58 09 09 |& CHARMA|SK];.X..|
|000018b0| 69 66 20 28 74 6d 70 20 | 3d 3d 20 30 2e 30 29 20 |if (tmp |== 0.0) |
|000018c0| 20 72 65 74 75 72 6e 28 | 2d 31 2e 30 29 3b 0a 58 | return(|-1.0);.X|
|000018d0| 09 09 73 75 6d 20 2b 3d | 20 74 6d 70 3b 0a 58 09 |..sum +=| tmp;.X.|
|000018e0| 09 7d 0a 58 0a 58 09 69 | 66 20 28 63 6f 75 6e 74 |.}.X.X.i|f (count|
|000018f0| 20 3d 3d 20 30 2e 30 29 | 20 20 72 65 74 75 72 6e | == 0.0)|  return|
|00001900| 28 2d 31 2e 30 29 3b 0a | 58 09 74 6d 70 20 3d 20 |(-1.0);.|X.tmp = |
|00001910| 28 73 75 6d 20 2f 20 63 | 6f 75 6e 74 29 20 2d 20 |(sum / c|ount) - |
|00001920| 6c 6f 67 6d 65 61 6e 3b | 0a 58 09 74 6d 70 20 3d |logmean;|.X.tmp =|
|00001930| 20 74 6d 70 20 3e 20 30 | 20 3f 20 74 6d 70 20 3a | tmp > 0| ? tmp :|
|00001940| 20 30 2e 30 20 2d 20 74 | 6d 70 3b 0a 58 09 73 63 | 0.0 - t|mp;.X.sc|
|00001950| 6f 72 65 20 3d 20 74 6d | 70 20 2f 20 28 6c 6f 67 |ore = tm|p / (log|
|00001960| 73 64 20 2f 20 73 71 72 | 74 28 63 6f 75 6e 74 29 |sd / sqr|t(count)|
|00001970| 29 3b 0a 58 2f 2a 09 70 | 72 69 6e 74 66 28 22 20 |);.X/*.p|rintf(" |
|00001980| 20 64 69 76 69 64 69 6e | 67 20 62 79 20 6c 6f 67 | dividin|g by log|
|00001990| 73 64 20 79 69 65 6c 64 | 73 20 25 67 22 2c 20 73 |sd yield|s %g", s|
|000019a0| 63 6f 72 65 29 3b 0a 58 | 09 74 6d 70 20 3d 20 74 |core);.X|.tmp = t|
|000019b0| 6d 70 20 2a 20 74 6d 70 | 3b 0a 58 09 74 6d 70 20 |mp * tmp|;.X.tmp |
|000019c0| 3d 20 28 74 6d 70 20 2a | 20 63 6f 75 6e 74 29 20 |= (tmp *| count) |
|000019d0| 2f 20 6c 6f 67 76 61 72 | 3b 0a 58 09 74 6d 70 20 |/ logvar|;.X.tmp |
|000019e0| 3d 20 65 78 70 28 2d 30 | 2e 35 20 2a 20 74 6d 70 |= exp(-0|.5 * tmp|
|000019f0| 29 3b 0a 58 09 70 72 69 | 6e 74 66 28 22 5c 6e 54 |);.X.pri|ntf("\nT|
|00001a00| 68 65 20 65 78 70 6f 6e | 65 6e 74 69 61 6c 20 79 |he expon|ential y|
|00001a10| 69 65 6c 64 73 20 25 67 | 22 2c 20 74 6d 70 29 3b |ields %g|", tmp);|
|00001a20| 0a 58 09 70 72 69 6e 74 | 66 28 22 5c 6e 22 29 3b |.X.print|f("\n");|
|00001a30| 0a 58 09 73 63 6f 72 65 | 20 3d 20 73 71 72 74 28 |.X.score| = sqrt(|
|00001a40| 63 6f 75 6e 74 29 20 2a | 20 65 78 70 28 28 30 20 |count) *| exp((0 |
|00001a50| 2d 20 74 6d 70 29 2f 32 | 2e 30 29 3b 0a 58 2a 2f |- tmp)/2|.0);.X*/|
|00001a60| 0a 58 09 72 65 74 75 72 | 6e 28 73 63 6f 72 65 29 |.X.retur|n(score)|
|00001a70| 3b 0a 58 7d 0a 58 0a 58 | 0a 58 2f 2a 20 53 63 6f |;.X}.X.X|.X/* Sco|
|00001a80| 72 65 20 61 20 76 65 63 | 74 6f 72 20 6f 66 20 69 |re a vec|tor of i|
|00001a90| 6e 74 65 67 65 72 73 20 | 74 68 61 74 20 72 65 70 |ntegers |that rep|
|00001aa0| 72 65 73 65 6e 74 20 63 | 68 61 72 61 63 74 65 72 |resent c|haracter|
|00001ab0| 73 2e 0a 58 20 2a 20 54 | 68 65 20 76 65 63 74 6f |s..X * T|he vecto|
|00001ac0| 72 20 69 73 20 74 65 72 | 6d 69 6e 61 74 65 64 20 |r is ter|minated |
|00001ad0| 62 79 20 61 20 76 61 6c | 75 65 20 6f 66 20 4e 4f |by a val|ue of NO|
|00001ae0| 4e 45 2e 0a 58 20 2a 20 | 53 63 6f 72 69 6e 67 20 |NE..X * |Scoring |
|00001af0| 69 73 20 62 61 73 65 64 | 20 6f 6e 20 72 61 74 69 |is based| on rati|
|00001b00| 6f 20 6f 66 20 6f 62 73 | 65 72 76 65 64 20 61 6e |o of obs|erved an|
|00001b10| 64 20 65 78 70 65 63 74 | 65 64 20 76 61 72 69 61 |d expect|ed varia|
|00001b20| 6e 63 65 2e 0a 58 20 2a | 2f 0a 58 66 6c 6f 61 74 |nce..X *|/.Xfloat|
|00001b30| 09 76 61 72 5f 31 73 63 | 6f 72 65 28 70 76 65 63 |.var_1sc|ore(pvec|
|00001b40| 29 0a 58 69 6e 74 09 2a | 70 76 65 63 3b 0a 58 7b |).Xint.*|pvec;.X{|
|00001b50| 0a 58 09 69 6e 74 09 09 | 69 3b 0a 58 09 69 6e 74 |.X.int..|i;.X.int|
|00001b60| 09 09 63 3b 0a 58 09 66 | 6c 6f 61 74 09 74 6d 70 |..c;.X.f|loat.tmp|
|00001b70| 2c 20 73 75 6d 2c 20 63 | 6f 75 6e 74 2c 20 73 63 |, sum, c|ount, sc|
|00001b80| 6f 72 65 3b 0a 58 0a 58 | 09 69 66 20 28 21 73 74 |ore;.X.X|.if (!st|
|00001b90| 61 74 73 31 6c 6f 61 64 | 65 64 29 20 20 7b 0a 58 |ats1load|ed)  {.X|
|00001ba0| 09 09 6c 6f 61 64 5f 31 | 73 74 61 74 73 5f 66 72 |..load_1|stats_fr|
|00001bb0| 6f 6d 28 6c 65 74 74 65 | 72 73 74 61 74 73 29 3b |om(lette|rstats);|
|00001bc0| 0a 58 09 09 7d 0a 58 0a | 58 09 63 6f 75 6e 74 20 |.X..}.X.|X.count |
|00001bd0| 3d 20 30 2e 30 3b 0a 58 | 09 73 75 6d 20 3d 20 30 |= 0.0;.X|.sum = 0|
|00001be0| 2e 30 3b 0a 58 09 77 68 | 69 6c 65 20 28 2a 70 76 |.0;.X.wh|ile (*pv|
|00001bf0| 65 63 20 21 3d 20 4e 4f | 4e 45 29 20 20 7b 0a 58 |ec != NO|NE)  {.X|
|00001c00| 09 09 63 6f 75 6e 74 20 | 2b 3d 20 31 2e 30 3b 0a |..count |+= 1.0;.|
|00001c10| 58 09 09 63 20 3d 20 2a | 70 76 65 63 2b 2b 3b 0a |X..c = *|pvec++;.|
|00001c20| 58 09 09 69 66 20 28 63 | 20 21 3d 20 63 20 26 20 |X..if (c| != c & |
|00001c30| 43 48 41 52 4d 41 53 4b | 29 20 20 72 65 74 75 72 |CHARMASK|)  retur|
|00001c40| 6e 28 30 2e 30 29 3b 0a | 58 09 09 74 6d 70 20 3d |n(0.0);.|X..tmp =|
|00001c50| 20 6c 6f 67 70 72 6f 62 | 5b 63 20 26 20 43 48 41 | logprob|[c & CHA|
|00001c60| 52 4d 41 53 4b 5d 3b 0a | 58 09 09 69 66 20 28 74 |RMASK];.|X..if (t|
|00001c70| 6d 70 20 3d 3d 20 30 2e | 30 29 20 20 72 65 74 75 |mp == 0.|0)  retu|
|00001c80| 72 6e 28 30 2e 30 29 3b | 0a 58 09 09 74 6d 70 20 |rn(0.0);|.X..tmp |
|00001c90| 3d 20 74 6d 70 20 2d 20 | 6c 6f 67 6d 65 61 6e 3b |= tmp - |logmean;|
|00001ca0| 0a 58 09 09 74 6d 70 20 | 3d 20 74 6d 70 20 2a 20 |.X..tmp |= tmp * |
|00001cb0| 74 6d 70 3b 0a 58 09 09 | 73 75 6d 20 2b 3d 20 74 |tmp;.X..|sum += t|
|00001cc0| 6d 70 3b 0a 58 09 09 7d | 0a 58 0a 58 09 69 66 20 |mp;.X..}|.X.X.if |
|00001cd0| 28 63 6f 75 6e 74 20 3d | 3d 20 30 2e 30 29 20 20 |(count =|= 0.0)  |
|00001ce0| 72 65 74 75 72 6e 28 30 | 2e 30 29 3b 0a 58 09 73 |return(0|.0);.X.s|
|00001cf0| 63 6f 72 65 20 3d 20 73 | 75 6d 20 2f 20 28 63 6f |core = s|um / (co|
|00001d00| 75 6e 74 20 2a 20 6c 6f | 67 76 61 72 29 3b 0a 58 |unt * lo|gvar);.X|
|00001d10| 09 72 65 74 75 72 6e 28 | 73 63 6f 72 65 29 3b 0a |.return(|score);.|
|00001d20| 58 7d 0a 58 0a 58 0a 58 | 2f 2a 20 53 63 6f 72 65 |X}.X.X.X|/* Score|
|00001d30| 20 61 20 76 65 63 74 6f | 72 20 6f 66 20 69 6e 74 | a vecto|r of int|
|00001d40| 65 67 65 72 73 20 74 68 | 61 74 20 72 65 70 72 65 |egers th|at repre|
|00001d50| 73 65 6e 74 20 63 68 61 | 72 61 63 74 65 72 73 2e |sent cha|racters.|
|00001d60| 0a 58 20 2a 20 54 68 65 | 20 76 65 63 74 6f 72 20 |.X * The| vector |
|00001d70| 69 73 20 74 65 72 6d 69 | 6e 61 74 65 64 20 62 79 |is termi|nated by|
|00001d80| 20 61 20 76 61 6c 75 65 | 20 6f 66 20 4e 4f 4e 45 | a value| of NONE|
|00001d90| 2e 0a 58 20 2a 20 53 63 | 6f 72 65 20 69 73 20 74 |..X * Sc|ore is t|
|00001da0| 68 65 20 70 72 6f 62 61 | 62 69 6c 69 74 79 20 74 |he proba|bility t|
|00001db0| 68 61 74 20 74 68 65 20 | 67 69 76 65 6e 20 63 68 |hat the |given ch|
|00001dc0| 61 72 61 63 74 65 72 73 | 0a 58 20 2a 20 77 65 72 |aracters|.X * wer|
|00001dd0| 65 20 64 72 61 77 6e 20 | 66 72 6f 6d 20 65 6e 67 |e drawn |from eng|
|00001de0| 6c 69 73 68 2e 0a 58 20 | 2a 20 4e 4f 54 45 3a 20 |lish..X |* NOTE: |
|00001df0| 64 6f 65 73 6e 27 74 20 | 63 6f 72 72 65 63 74 6c |doesn't |correctl|
|00001e00| 79 20 68 61 6e 64 6c 65 | 20 72 65 70 65 61 74 65 |y handle| repeate|
|00001e10| 64 20 6c 65 74 74 65 72 | 73 2e 0a 58 20 2a 2f 0a |d letter|s..X */.|
|00001e20| 58 66 6c 6f 61 74 09 70 | 72 6f 62 5f 31 73 63 6f |Xfloat.p|rob_1sco|
|00001e30| 72 65 28 70 76 65 63 29 | 0a 58 69 6e 74 09 2a 70 |re(pvec)|.Xint.*p|
|00001e40| 76 65 63 3b 0a 58 7b 0a | 58 09 69 6e 74 09 09 69 |vec;.X{.|X.int..i|
|00001e50| 3b 0a 58 09 69 6e 74 09 | 09 63 3b 0a 58 09 66 6c |;.X.int.|.c;.X.fl|
|00001e60| 6f 61 74 09 74 6d 70 2c | 20 70 72 6f 64 75 63 74 |oat.tmp,| product|
|00001e70| 2c 20 63 6f 75 6e 74 2c | 20 73 63 6f 72 65 3b 0a |, count,| score;.|
|00001e80| 58 0a 58 09 69 66 20 28 | 21 73 74 61 74 73 31 6c |X.X.if (|!stats1l|
|00001e90| 6f 61 64 65 64 29 20 20 | 7b 0a 58 09 09 6c 6f 61 |oaded)  |{.X..loa|
|00001ea0| 64 5f 31 73 74 61 74 73 | 5f 66 72 6f 6d 28 6c 65 |d_1stats|_from(le|
|00001eb0| 74 74 65 72 73 74 61 74 | 73 29 3b 0a 58 09 09 7d |tterstat|s);.X..}|
|00001ec0| 0a 58 0a 58 09 63 6f 75 | 6e 74 20 3d 20 30 2e 30 |.X.X.cou|nt = 0.0|
|00001ed0| 3b 0a 58 09 70 72 6f 64 | 75 63 74 20 3d 20 31 2e |;.X.prod|uct = 1.|
|00001ee0| 30 3b 0a 58 09 77 68 69 | 6c 65 20 28 2a 70 76 65 |0;.X.whi|le (*pve|
|00001ef0| 63 20 21 3d 20 4e 4f 4e | 45 29 20 20 7b 0a 58 09 |c != NON|E)  {.X.|
|00001f00| 09 63 6f 75 6e 74 20 2b | 3d 20 31 2e 30 3b 0a 58 |.count +|= 1.0;.X|
|00001f10| 09 09 63 20 3d 20 2a 70 | 76 65 63 2b 2b 3b 0a 58 |..c = *p|vec++;.X|
|00001f20| 09 09 69 66 20 28 63 20 | 21 3d 20 63 20 26 20 43 |..if (c |!= c & C|
|00001f30| 48 41 52 4d 41 53 4b 29 | 20 20 72 65 74 75 72 6e |HARMASK)|  return|
|00001f40| 28 30 2e 30 29 3b 0a 58 | 09 09 70 72 6f 64 75 63 |(0.0);.X|..produc|
|00001f50| 74 20 2a 3d 20 70 72 6f | 62 5b 63 5d 20 2a 20 63 |t *= pro|b[c] * c|
|00001f60| 6f 75 6e 74 3b 0a 58 09 | 09 7d 0a 58 0a 58 09 69 |ount;.X.|.}.X.X.i|
|00001f70| 66 20 28 63 6f 75 6e 74 | 20 3d 3d 20 30 2e 30 29 |f (count| == 0.0)|
|00001f80| 20 20 72 65 74 75 72 6e | 28 30 2e 30 29 3b 0a 58 |  return|(0.0);.X|
|00001f90| 09 73 63 6f 72 65 20 3d | 20 70 72 6f 64 75 63 74 |.score =| product|
|00001fa0| 3b 0a 58 09 72 65 74 75 | 72 6e 28 73 63 6f 72 65 |;.X.retu|rn(score|
|00001fb0| 29 3b 0a 58 7d 0a 58 0a | 58 0a 58 2f 2a 20 53 63 |);.X}.X.|X.X/* Sc|
|00001fc0| 6f 72 65 20 61 20 67 75 | 65 73 73 20 62 61 73 65 |ore a gu|ess base|
|00001fd0| 64 20 6f 6e 20 6c 65 74 | 74 65 72 20 70 61 69 72 |d on let|ter pair|
|00001fe0| 20 66 72 65 71 75 65 6e | 63 69 65 73 2e 0a 58 20 | frequen|cies..X |
|00001ff0| 2a 20 54 68 65 20 72 65 | 74 75 72 6e 65 64 20 73 |* The re|turned s|
|00002000| 63 6f 72 65 20 69 73 20 | 74 68 65 20 6e 75 6d 62 |core is |the numb|
|00002010| 65 72 20 6f 66 20 73 74 | 61 6e 64 61 72 64 20 64 |er of st|andard d|
|00002020| 65 76 69 61 74 69 6f 6e | 73 0a 58 20 2a 20 74 68 |eviation|s.X * th|
|00002030| 61 74 20 74 68 65 20 6f | 62 73 65 72 76 65 64 20 |at the o|bserved |
|00002040| 73 74 61 74 69 73 74 69 | 63 73 20 64 69 66 66 65 |statisti|cs diffe|
|00002050| 72 73 20 66 72 6f 6d 20 | 69 74 73 20 65 78 70 65 |rs from |its expe|
|00002060| 63 74 65 64 20 76 61 6c | 75 65 2e 0a 58 20 2a 20 |cted val|ue..X * |
|00002070| 53 63 6f 72 65 73 20 61 | 72 65 20 70 6f 73 69 74 |Scores a|re posit|
|00002080| 69 76 65 20 77 69 74 68 | 20 6c 6f 77 20 73 63 6f |ive with| low sco|
|00002090| 72 65 73 20 62 65 69 6e | 67 20 62 65 74 74 65 72 |res bein|g better|
|000020a0| 2e 0a 58 20 2a 20 41 20 | 6e 65 67 61 74 69 76 65 |..X * A |negative|
|000020b0| 20 73 63 6f 72 65 20 69 | 6e 64 69 63 61 74 65 73 | score i|ndicates|
|000020c0| 20 61 6e 20 69 6d 70 6f | 73 73 69 62 6c 65 20 70 | an impo|ssible p|
|000020d0| 6c 61 69 6e 74 65 78 74 | 20 76 61 6c 75 65 2e 0a |laintext| value..|
|000020e0| 58 20 2a 2f 0a 58 66 6c | 6f 61 74 09 67 73 69 5f |X */.Xfl|oat.gsi_|
|000020f0| 32 73 63 6f 72 65 28 67 | 73 69 29 0a 58 72 65 67 |2score(g|si).Xreg|
|00002100| 09 09 67 73 69 6e 66 6f | 09 2a 67 73 69 3b 0a 58 |..gsinfo|.*gsi;.X|
|00002110| 7b 0a 58 09 66 6c 6f 61 | 74 09 73 63 6f 72 65 3b |{.X.floa|t.score;|
|00002120| 0a 58 09 66 6c 6f 61 74 | 09 74 6f 74 61 6c 3b 0a |.X.float|.total;.|
|00002130| 58 09 66 6c 6f 61 74 09 | 74 6d 70 3b 0a 58 09 69 |X.float.|tmp;.X.i|
|00002140| 6e 74 09 09 6e 63 68 61 | 72 73 3b 0a 58 09 69 6e |nt..ncha|rs;.X.in|
|00002150| 74 09 09 69 3b 0a 58 72 | 65 67 09 69 6e 74 09 09 |t..i;.Xr|eg.int..|
|00002160| 70 6f 73 3b 09 0a 58 72 | 65 67 09 69 6e 74 09 09 |pos;..Xr|eg.int..|
|00002170| 63 3b 0a 58 72 65 67 09 | 69 6e 74 09 09 63 65 6e |c;.Xreg.|int..cen|
|00002180| 74 65 72 5f 6c 65 74 74 | 65 72 3b 0a 58 09 69 6e |ter_lett|er;.X.in|
|00002190| 74 09 09 6c 65 66 74 5f | 6c 65 74 74 65 72 2c 20 |t..left_|letter, |
|000021a0| 72 69 67 68 74 5f 6c 65 | 74 74 65 72 3b 0a 58 09 |right_le|tter;.X.|
|000021b0| 66 6c 6f 61 74 09 70 61 | 69 72 5f 73 63 6f 72 65 |float.pa|ir_score|
|000021c0| 3b 0a 58 0a 58 09 69 66 | 20 28 21 73 74 61 74 73 |;.X.X.if| (!stats|
|000021d0| 32 6c 6f 61 64 65 64 29 | 20 20 7b 0a 58 09 09 6c |2loaded)|  {.X..l|
|000021e0| 6f 61 64 5f 32 73 74 61 | 74 73 5f 66 72 6f 6d 28 |oad_2sta|ts_from(|
|000021f0| 62 69 67 72 61 6d 73 74 | 61 74 73 29 3b 0a 58 09 |bigramst|ats);.X.|
|00002200| 09 7d 0a 58 0a 58 09 6e | 63 68 61 72 73 20 3d 20 |.}.X.X.n|chars = |
|00002210| 30 3b 0a 58 09 74 6f 74 | 61 6c 20 3d 20 30 2e 30 |0;.X.tot|al = 0.0|
|00002220| 3b 0a 58 09 66 6f 72 20 | 28 69 20 3d 20 30 20 3b |;.X.for |(i = 0 ;|
|00002230| 20 28 70 6f 73 20 3d 20 | 67 73 69 2d 3e 63 70 6f | (pos = |gsi->cpo|
|00002240| 73 5b 69 5d 29 20 21 3d | 20 4e 4f 4e 45 20 3b 20 |s[i]) !=| NONE ; |
|00002250| 69 2b 2b 29 20 20 7b 0a | 58 09 09 6e 63 68 61 72 |i++)  {.|X..nchar|
|00002260| 73 2b 2b 3b 0a 58 09 09 | 63 20 3d 20 28 67 73 69 |s++;.X..|c = (gsi|
|00002270| 2d 3e 63 67 75 65 73 73 | 65 64 29 5b 70 6f 73 5d |->cguess|ed)[pos]|
|00002280| 3b 0a 58 09 09 63 65 6e | 74 65 72 5f 6c 65 74 74 |;.X..cen|ter_lett|
|00002290| 65 72 20 3d 20 63 68 61 | 72 5f 62 69 6d 61 70 5b |er = cha|r_bimap[|
|000022a0| 63 20 26 20 43 48 41 52 | 4d 41 53 4b 5d 3b 0a 58 |c & CHAR|MASK];.X|
|000022b0| 09 09 69 66 20 28 73 6c | 6c 6f 67 70 72 6f 62 5b |..if (sl|logprob[|
|000022c0| 63 65 6e 74 65 72 5f 6c | 65 74 74 65 72 5d 20 3d |center_l|etter] =|
|000022d0| 3d 20 30 2e 30 29 0a 58 | 09 09 09 72 65 74 75 72 |= 0.0).X|...retur|
|000022e0| 6e 28 2d 31 2e 30 29 3b | 0a 58 0a 58 09 09 69 66 |n(-1.0);|.X.X..if|
|000022f0| 20 28 70 6f 73 20 3d 3d | 20 30 29 20 7b 0a 58 09 | (pos ==| 0) {.X.|
|00002300| 09 09 74 6f 74 61 6c 20 | 2b 3d 20 73 6c 6c 6f 67 |..total |+= sllog|
|00002310| 70 72 6f 62 5b 63 65 6e | 74 65 72 5f 6c 65 74 74 |prob[cen|ter_lett|
|00002320| 65 72 5d 3b 0a 58 09 09 | 09 7d 0a 58 09 09 65 6c |er];.X..|.}.X..el|
|00002330| 73 65 20 7b 0a 58 09 09 | 09 63 20 3d 20 28 67 73 |se {.X..|.c = (gs|
|00002340| 69 2d 3e 63 6b 6e 6f 77 | 6e 29 5b 70 6f 73 20 2d |i->cknow|n)[pos -|
|00002350| 20 31 5d 3b 0a 58 09 09 | 09 69 66 20 28 63 20 3d | 1];.X..|.if (c =|
|00002360| 3d 20 4e 4f 4e 45 29 20 | 20 7b 0a 58 09 09 09 09 |= NONE) | {.X....|
|00002370| 63 20 3d 20 28 67 73 69 | 2d 3e 63 67 75 65 73 73 |c = (gsi|->cguess|
|00002380| 65 64 29 5b 70 6f 73 20 | 2d 20 31 5d 3b 0a 58 09 |ed)[pos |- 1];.X.|
|00002390| 09 09 09 7d 0a 58 09 09 | 09 69 66 20 28 63 20 3d |...}.X..|.if (c =|
|000023a0| 3d 20 4e 4f 4e 45 29 20 | 20 7b 0a 58 09 09 09 09 |= NONE) | {.X....|
|000023b0| 74 6f 74 61 6c 20 2b 3d | 20 73 6c 6c 6f 67 70 72 |total +=| sllogpr|
|000023c0| 6f 62 5b 63 65 6e 74 65 | 72 5f 6c 65 74 74 65 72 |ob[cente|r_letter|
|000023d0| 5d 3b 0a 58 09 09 09 09 | 7d 0a 58 09 09 09 65 6c |];.X....|}.X...el|
|000023e0| 73 65 20 7b 0a 58 09 09 | 09 09 6c 65 66 74 5f 6c |se {.X..|..left_l|
|000023f0| 65 74 74 65 72 20 3d 20 | 63 68 61 72 5f 62 69 6d |etter = |char_bim|
|00002400| 61 70 5b 63 20 26 20 43 | 48 41 52 4d 41 53 4b 5d |ap[c & C|HARMASK]|
|00002410| 3b 0a 58 09 09 09 09 70 | 61 69 72 5f 73 63 6f 72 |;.X....p|air_scor|
|00002420| 65 20 3d 20 62 69 6c 6f | 67 70 72 6f 62 5b 6c 65 |e = bilo|gprob[le|
|00002430| 66 74 5f 6c 65 74 74 65 | 72 5d 5b 63 65 6e 74 65 |ft_lette|r][cente|
|00002440| 72 5f 6c 65 74 74 65 72 | 5d 3b 0a 58 09 09 09 09 |r_letter|];.X....|
|00002450| 69 66 20 28 70 61 69 72 | 5f 73 63 6f 72 65 20 3d |if (pair|_score =|
|00002460| 3d 20 30 2e 30 29 0a 58 | 09 09 09 09 09 72 65 74 |= 0.0).X|.....ret|
|00002470| 75 72 6e 28 2d 31 2e 30 | 29 3b 0a 58 09 09 09 09 |urn(-1.0|);.X....|
|00002480| 74 6f 74 61 6c 20 2b 3d | 20 70 61 69 72 5f 73 63 |total +=| pair_sc|
|00002490| 6f 72 65 20 2d 20 73 6c | 6c 6f 67 70 72 6f 62 5b |ore - sl|logprob[|
|000024a0| 63 65 6e 74 65 72 5f 6c | 65 74 74 65 72 5d 3b 0a |center_l|etter];.|
|000024b0| 58 09 09 09 09 7d 0a 58 | 09 09 09 7d 0a 58 0a 58 |X....}.X|...}.X.X|
|000024c0| 09 09 69 66 20 28 70 6f | 73 20 3d 3d 20 28 42 4c |..if (po|s == (BL|
|000024d0| 4f 43 4b 53 49 5a 45 20 | 2d 20 31 29 29 20 7b 0a |OCKSIZE |- 1)) {.|
|000024e0| 58 09 09 09 74 6f 74 61 | 6c 20 2b 3d 20 73 6c 6c |X...tota|l += sll|
|000024f0| 6f 67 70 72 6f 62 5b 63 | 65 6e 74 65 72 5f 6c 65 |ogprob[c|enter_le|
|00002500| 74 74 65 72 5d 3b 0a 58 | 09 09 09 7d 0a 58 09 09 |tter];.X|...}.X..|
|00002510| 65 6c 73 65 20 7b 0a 58 | 09 09 09 63 20 3d 20 28 |else {.X|...c = (|
|00002520| 67 73 69 2d 3e 63 6b 6e | 6f 77 6e 29 5b 70 6f 73 |gsi->ckn|own)[pos|
|00002530| 20 2b 20 31 5d 3b 0a 58 | 09 09 09 69 66 20 28 63 | + 1];.X|...if (c|
|00002540| 20 3d 3d 20 4e 4f 4e 45 | 29 20 20 7b 0a 58 09 09 | == NONE|)  {.X..|
|00002550| 09 09 63 20 3d 20 28 67 | 73 69 2d 3e 63 67 75 65 |..c = (g|si->cgue|
|00002560| 73 73 65 64 29 5b 70 6f | 73 20 2b 20 31 5d 3b 0a |ssed)[po|s + 1];.|
|00002570| 58 09 09 09 09 7d 0a 58 | 09 09 09 69 66 20 28 63 |X....}.X|...if (c|
|00002580| 20 3d 3d 20 4e 4f 4e 45 | 29 20 20 7b 0a 58 09 09 | == NONE|)  {.X..|
|00002590| 09 09 74 6f 74 61 6c 20 | 2b 3d 20 73 6c 6c 6f 67 |..total |+= sllog|
|000025a0| 70 72 6f 62 5b 63 65 6e | 74 65 72 5f 6c 65 74 74 |prob[cen|ter_lett|
|000025b0| 65 72 5d 3b 0a 58 09 09 | 09 09 7d 0a 58 09 09 09 |er];.X..|..}.X...|
|000025c0| 65 6c 73 65 20 7b 0a 58 | 09 09 09 09 72 69 67 68 |else {.X|....righ|
|000025d0| 74 5f 6c 65 74 74 65 72 | 20 3d 20 63 68 61 72 5f |t_letter| = char_|
|000025e0| 62 69 6d 61 70 5b 63 20 | 26 20 43 48 41 52 4d 41 |bimap[c |& CHARMA|
|000025f0| 53 4b 5d 3b 0a 58 09 09 | 09 09 70 61 69 72 5f 73 |SK];.X..|..pair_s|
|00002600| 63 6f 72 65 20 3d 20 62 | 69 6c 6f 67 70 72 6f 62 |core = b|ilogprob|
|00002610| 5b 63 65 6e 74 65 72 5f | 6c 65 74 74 65 72 5d 5b |[center_|letter][|
|00002620| 72 69 67 68 74 5f 6c 65 | 74 74 65 72 5d 3b 0a 58 |right_le|tter];.X|
|00002630| 09 09 09 09 69 66 20 28 | 70 61 69 72 5f 73 63 6f |....if (|pair_sco|
|00002640| 72 65 20 3d 3d 20 30 2e | 30 29 0a 58 09 09 09 09 |re == 0.|0).X....|
|00002650| 09 72 65 74 75 72 6e 28 | 2d 31 2e 30 29 3b 0a 58 |.return(|-1.0);.X|
|00002660| 09 09 09 09 74 6f 74 61 | 6c 20 2b 3d 20 70 61 69 |....tota|l += pai|
|00002670| 72 5f 73 63 6f 72 65 20 | 2d 20 73 6c 6c 6f 67 70 |r_score |- sllogp|
|00002680| 72 6f 62 5b 63 65 6e 74 | 65 72 5f 6c 65 74 74 65 |rob[cent|er_lette|
|00002690| 72 5d 3b 0a 58 09 09 09 | 09 7d 0a 58 09 09 09 7d |r];.X...|.}.X...}|
|000026a0| 0a 58 09 09 7d 0a 58 0a | 58 09 69 66 20 28 6e 63 |.X..}.X.|X.if (nc|
|000026b0| 68 61 72 73 20 3d 3d 20 | 30 29 0a 58 09 09 72 65 |hars == |0).X..re|
|000026c0| 74 75 72 6e 28 2d 31 2e | 30 29 3b 0a 58 09 74 6d |turn(-1.|0);.X.tm|
|000026d0| 70 20 3d 20 28 74 6f 74 | 61 6c 20 2f 20 6e 63 68 |p = (tot|al / nch|
|000026e0| 61 72 73 29 20 2d 20 73 | 63 6f 72 65 32 5f 6d 65 |ars) - s|core2_me|
|000026f0| 61 6e 3b 0a 58 09 74 6d | 70 20 3d 20 74 6d 70 20 |an;.X.tm|p = tmp |
|00002700| 3e 20 30 2e 30 20 3f 20 | 74 6d 70 20 3a 20 30 2e |> 0.0 ? |tmp : 0.|
|00002710| 30 20 2d 20 74 6d 70 3b | 0a 58 09 73 63 6f 72 65 |0 - tmp;|.X.score|
|00002720| 20 3d 20 74 6d 70 20 2f | 20 28 73 63 6f 72 65 32 | = tmp /| (score2|
|00002730| 5f 73 64 20 2f 20 69 73 | 71 72 74 5b 6e 63 68 61 |_sd / is|qrt[ncha|
|00002740| 72 73 5d 29 3b 0a 58 09 | 72 65 74 75 72 6e 28 73 |rs]);.X.|return(s|
|00002750| 63 6f 72 65 29 3b 0a 58 | 7d 0a 58 0a 58 0a 58 2f |core);.X|}.X.X.X/|
|00002760| 2a 20 53 63 6f 72 65 20 | 61 20 67 75 65 73 73 20 |* Score |a guess |
|00002770| 62 61 73 65 64 20 6f 6e | 20 73 69 6e 67 6c 65 20 |based on| single |
|00002780| 6c 65 74 74 65 72 20 66 | 72 65 71 75 65 6e 63 69 |letter f|requenci|
|00002790| 65 73 2e 0a 58 20 2a 20 | 54 68 65 20 72 65 74 75 |es..X * |The retu|
|000027a0| 72 6e 65 64 20 73 63 6f | 72 65 20 69 73 20 74 68 |rned sco|re is th|
|000027b0| 65 20 6e 75 6d 62 65 72 | 20 6f 66 20 73 74 61 6e |e number| of stan|
|000027c0| 64 61 72 64 20 64 65 76 | 69 61 74 69 6f 6e 73 0a |dard dev|iations.|
|000027d0| 58 20 2a 20 74 68 61 74 | 20 74 68 65 20 6f 62 73 |X * that| the obs|
|000027e0| 65 72 76 65 64 20 73 74 | 61 74 69 73 74 69 63 73 |erved st|atistics|
|000027f0| 20 64 69 66 66 65 72 73 | 20 66 72 6f 6d 20 69 74 | differs| from it|
|00002800| 73 20 65 78 70 65 63 74 | 65 64 20 76 61 6c 75 65 |s expect|ed value|
|00002810| 2e 0a 58 20 2a 20 53 63 | 6f 72 65 73 20 61 72 65 |..X * Sc|ores are|
|00002820| 20 70 6f 73 69 74 69 76 | 65 20 77 69 74 68 20 6c | positiv|e with l|
|00002830| 6f 77 20 73 63 6f 72 65 | 73 20 62 65 69 6e 67 20 |ow score|s being |
|00002840| 62 65 74 74 65 72 2e 0a | 58 20 2a 20 41 20 6e 65 |better..|X * A ne|
|00002850| 67 61 74 69 76 65 20 73 | 63 6f 72 65 20 69 6e 64 |gative s|core ind|
|00002860| 69 63 61 74 65 73 20 61 | 6e 20 69 6d 70 6f 73 73 |icates a|n imposs|
|00002870| 69 62 6c 65 20 70 6c 61 | 69 6e 74 65 78 74 20 76 |ible pla|intext v|
|00002880| 61 6c 75 65 2e 0a 58 20 | 2a 2f 0a 58 66 6c 6f 61 |alue..X |*/.Xfloa|
|00002890| 74 09 67 73 69 5f 31 73 | 63 6f 72 65 28 67 73 69 |t.gsi_1s|core(gsi|
|000028a0| 29 0a 58 72 65 67 09 09 | 67 73 69 6e 66 6f 09 2a |).Xreg..|gsinfo.*|
|000028b0| 67 73 69 3b 0a 58 7b 0a | 58 72 65 67 09 69 6e 74 |gsi;.X{.|Xreg.int|
|000028c0| 09 09 70 6f 73 3b 0a 58 | 09 69 6e 74 09 09 69 3b |..pos;.X|.int..i;|
|000028d0| 0a 58 09 69 6e 74 09 09 | 63 3b 0a 58 09 69 6e 74 |.X.int..|c;.X.int|
|000028e0| 09 09 6e 63 68 61 72 73 | 3b 0a 58 09 66 6c 6f 61 |..nchars|;.X.floa|
|000028f0| 74 09 73 75 6d 2c 20 73 | 63 6f 72 65 3b 0a 58 72 |t.sum, s|core;.Xr|
|00002900| 65 67 09 66 6c 6f 61 74 | 09 74 6d 70 3b 0a 58 0a |eg.float|.tmp;.X.|
|00002910| 58 09 69 66 20 28 21 73 | 74 61 74 73 31 6c 6f 61 |X.if (!s|tats1loa|
|00002920| 64 65 64 29 20 20 7b 0a | 58 09 09 6c 6f 61 64 5f |ded)  {.|X..load_|
|00002930| 31 73 74 61 74 73 5f 66 | 72 6f 6d 28 6c 65 74 74 |1stats_f|rom(lett|
|00002940| 65 72 73 74 61 74 73 29 | 3b 0a 58 09 09 7d 0a 58 |erstats)|;.X..}.X|
|00002950| 0a 58 09 6e 63 68 61 72 | 73 20 3d 20 30 3b 0a 58 |.X.nchar|s = 0;.X|
|00002960| 09 73 75 6d 20 3d 20 30 | 2e 30 3b 0a 58 09 66 6f |.sum = 0|.0;.X.fo|
|00002970| 72 20 28 69 20 3d 20 30 | 20 3b 20 28 70 6f 73 20 |r (i = 0| ; (pos |
|00002980| 3d 20 67 73 69 2d 3e 63 | 70 6f 73 5b 69 5d 29 20 |= gsi->c|pos[i]) |
|00002990| 21 3d 20 4e 4f 4e 45 20 | 3b 20 69 2b 2b 29 20 20 |!= NONE |; i++)  |
|000029a0| 7b 0a 58 09 09 6e 63 68 | 61 72 73 2b 2b 3b 0a 58 |{.X..nch|ars++;.X|
|000029b0| 09 09 63 20 3d 20 28 67 | 73 69 2d 3e 63 67 75 65 |..c = (g|si->cgue|
|000029c0| 73 73 65 64 29 5b 70 6f | 73 5d 3b 0a 58 09 09 74 |ssed)[po|s];.X..t|
|000029d0| 6d 70 20 3d 20 6c 6f 67 | 70 72 6f 62 5b 63 20 26 |mp = log|prob[c &|
|000029e0| 20 43 48 41 52 4d 41 53 | 4b 5d 3b 0a 58 09 09 69 | CHARMAS|K];.X..i|
|000029f0| 66 20 28 74 6d 70 20 3d | 3d 20 30 2e 30 29 0a 58 |f (tmp =|= 0.0).X|
|00002a00| 09 09 09 72 65 74 75 72 | 6e 28 2d 31 2e 30 29 3b |...retur|n(-1.0);|
|00002a10| 0a 58 09 09 73 75 6d 20 | 2b 3d 20 74 6d 70 3b 0a |.X..sum |+= tmp;.|
|00002a20| 58 09 09 7d 0a 58 0a 58 | 09 69 66 20 28 6e 63 68 |X..}.X.X|.if (nch|
|00002a30| 61 72 73 20 3d 3d 20 30 | 29 0a 58 09 09 72 65 74 |ars == 0|).X..ret|
|00002a40| 75 72 6e 28 2d 31 2e 30 | 29 3b 0a 58 09 74 6d 70 |urn(-1.0|);.X.tmp|
|00002a50| 20 3d 20 28 73 75 6d 20 | 2f 20 6e 63 68 61 72 73 | = (sum |/ nchars|
|00002a60| 29 20 2d 20 6c 6f 67 6d | 65 61 6e 3b 0a 58 09 74 |) - logm|ean;.X.t|
|00002a70| 6d 70 20 3d 20 74 6d 70 | 20 3e 20 30 20 3f 20 74 |mp = tmp| > 0 ? t|
|00002a80| 6d 70 20 3a 20 30 2e 30 | 20 2d 20 74 6d 70 3b 0a |mp : 0.0| - tmp;.|
|00002a90| 58 09 73 63 6f 72 65 20 | 3d 20 74 6d 70 20 2f 20 |X.score |= tmp / |
|00002aa0| 28 6c 6f 67 73 64 20 2f | 20 69 73 71 72 74 5b 6e |(logsd /| isqrt[n|
|00002ab0| 63 68 61 72 73 5d 29 3b | 0a 58 0a 58 09 72 65 74 |chars]);|.X.X.ret|
|00002ac0| 75 72 6e 28 73 63 6f 72 | 65 29 3b 0a 58 7d 0a 58 |urn(scor|e);.X}.X|
|00002ad0| 0a 58 0a 58 2f 2a 20 43 | 6f 6d 70 75 74 65 20 65 |.X.X/* C|ompute e|
|00002ae0| 78 70 65 63 74 65 64 20 | 76 61 6c 75 65 20 6f 66 |xpected |value of|
|00002af0| 20 61 20 73 63 6f 72 69 | 6e 67 20 66 75 6e 63 74 | a scori|ng funct|
|00002b00| 69 6f 6e 20 67 69 76 65 | 6e 0a 58 20 2a 20 61 20 |ion give|n.X * a |
|00002b10| 76 65 63 74 6f 72 20 6f | 66 20 70 72 6f 62 61 62 |vector o|f probab|
|00002b20| 69 6c 69 74 69 65 73 20 | 66 6f 72 20 76 61 6c 75 |ilities |for valu|
|00002b30| 65 73 20 61 6e 64 20 61 | 20 76 65 63 74 6f 72 20 |es and a| vector |
|00002b40| 6f 66 0a 58 20 2a 20 73 | 63 6f 72 65 73 20 66 6f |of.X * s|cores fo|
|00002b50| 72 20 65 61 63 68 20 76 | 61 6c 75 65 2e 0a 58 20 |r each v|alue..X |
|00002b60| 2a 2f 0a 58 66 6c 6f 61 | 74 09 76 65 63 5f 6d 65 |*/.Xfloa|t.vec_me|
|00002b70| 61 6e 28 70 72 6f 62 76 | 65 63 2c 20 73 63 6f 72 |an(probv|ec, scor|
|00002b80| 65 76 65 63 2c 20 6d 61 | 78 69 6e 64 65 78 29 0a |evec, ma|xindex).|
|00002b90| 58 66 6c 6f 61 74 09 2a | 70 72 6f 62 76 65 63 3b |Xfloat.*|probvec;|
|00002ba0| 0a 58 66 6c 6f 61 74 09 | 2a 73 63 6f 72 65 76 65 |.Xfloat.|*scoreve|
|00002bb0| 63 3b 0a 58 69 6e 74 09 | 09 6d 61 78 69 6e 64 65 |c;.Xint.|.maxinde|
|00002bc0| 78 3b 0a 58 7b 0a 58 09 | 69 6e 74 09 09 69 3b 0a |x;.X{.X.|int..i;.|
|00002bd0| 58 09 66 6c 6f 61 74 09 | 6d 65 61 6e 3b 0a 58 0a |X.float.|mean;.X.|
|00002be0| 58 09 6d 65 61 6e 20 3d | 20 30 2e 30 3b 0a 58 09 |X.mean =| 0.0;.X.|
|00002bf0| 66 6f 72 20 28 69 20 3d | 20 30 20 3b 20 69 20 3c |for (i =| 0 ; i <|
|00002c00| 3d 20 6d 61 78 69 6e 64 | 65 78 20 3b 20 69 2b 2b |= maxind|ex ; i++|
|00002c10| 29 20 7b 0a 58 09 09 6d | 65 61 6e 20 2b 3d 20 28 |) {.X..m|ean += (|
|00002c20| 2a 70 72 6f 62 76 65 63 | 2b 2b 29 20 2a 20 28 2a |*probvec|++) * (*|
|00002c30| 73 63 6f 72 65 76 65 63 | 2b 2b 29 3b 0a 58 09 09 |scorevec|++);.X..|
|00002c40| 7d 0a 58 09 72 65 74 75 | 72 6e 28 6d 65 61 6e 29 |}.X.retu|rn(mean)|
|00002c50| 3b 0a 58 7d 0a 58 0a 58 | 0a 58 2f 2a 20 43 6f 6d |;.X}.X.X|.X/* Com|
|00002c60| 70 75 74 65 20 76 61 72 | 69 61 6e 63 65 20 6f 66 |pute var|iance of|
|00002c70| 20 61 20 73 63 6f 72 69 | 6e 67 20 66 75 6e 63 74 | a scori|ng funct|
|00002c80| 69 6f 6e 20 67 69 76 65 | 6e 0a 58 20 2a 20 61 20 |ion give|n.X * a |
|00002c90| 76 65 63 74 6f 72 20 6f | 66 20 70 72 6f 62 61 62 |vector o|f probab|
|00002ca0| 69 6c 69 74 69 65 73 20 | 66 6f 72 20 76 61 6c 75 |ilities |for valu|
|00002cb0| 65 73 20 61 6e 64 20 61 | 20 76 65 63 74 6f 72 20 |es and a| vector |
|00002cc0| 6f 66 0a 58 20 2a 20 73 | 63 6f 72 65 73 20 66 6f |of.X * s|cores fo|
|00002cd0| 72 20 65 61 63 68 20 76 | 61 6c 75 65 2e 0a 58 20 |r each v|alue..X |
|00002ce0| 2a 2f 0a 58 66 6c 6f 61 | 74 09 76 65 63 5f 76 61 |*/.Xfloa|t.vec_va|
|00002cf0| 72 69 61 6e 63 65 28 70 | 72 6f 62 76 65 63 2c 20 |riance(p|robvec, |
|00002d00| 73 63 6f 72 65 76 65 63 | 2c 20 6d 61 78 69 6e 64 |scorevec|, maxind|
|00002d10| 65 78 29 0a 58 66 6c 6f | 61 74 09 2a 70 72 6f 62 |ex).Xflo|at.*prob|
|00002d20| 76 65 63 3b 0a 58 66 6c | 6f 61 74 09 2a 73 63 6f |vec;.Xfl|oat.*sco|
|00002d30| 72 65 76 65 63 3b 0a 58 | 69 6e 74 09 09 6d 61 78 |revec;.X|int..max|
|00002d40| 69 6e 64 65 78 3b 0a 58 | 7b 0a 58 09 69 6e 74 09 |index;.X|{.X.int.|
|00002d50| 09 69 3b 0a 58 09 66 6c | 6f 61 74 09 76 61 72 2c |.i;.X.fl|oat.var,|
|00002d60| 20 6d 65 61 6e 3b 0a 58 | 09 66 6c 6f 61 74 09 64 | mean;.X|.float.d|
|00002d70| 65 6c 74 61 3b 0a 58 0a | 58 09 6d 65 61 6e 20 3d |elta;.X.|X.mean =|
|00002d80| 20 76 65 63 5f 6d 65 61 | 6e 28 70 72 6f 62 76 65 | vec_mea|n(probve|
|00002d90| 63 2c 20 73 63 6f 72 65 | 76 65 63 2c 20 6d 61 78 |c, score|vec, max|
|00002da0| 69 6e 64 65 78 29 3b 0a | 58 09 76 61 72 20 3d 20 |index);.|X.var = |
|00002db0| 30 2e 30 3b 0a 58 09 66 | 6f 72 20 28 69 20 3d 20 |0.0;.X.f|or (i = |
|00002dc0| 30 20 3b 20 69 20 3c 3d | 20 6d 61 78 69 6e 64 65 |0 ; i <=| maxinde|
|00002dd0| 78 20 3b 20 69 2b 2b 29 | 20 7b 0a 58 09 09 64 65 |x ; i++)| {.X..de|
|00002de0| 6c 74 61 20 3d 20 28 2a | 73 63 6f 72 65 76 65 63 |lta = (*|scorevec|
|00002df0| 2b 2b 29 20 2d 20 6d 65 | 61 6e 3b 0a 58 09 09 76 |++) - me|an;.X..v|
|00002e00| 61 72 20 2b 3d 20 28 2a | 70 72 6f 62 76 65 63 2b |ar += (*|probvec+|
|00002e10| 2b 29 20 2a 20 28 64 65 | 6c 74 61 20 2a 20 64 65 |+) * (de|lta * de|
|00002e20| 6c 74 61 29 3b 0a 58 09 | 09 7d 0a 58 09 72 65 74 |lta);.X.|.}.X.ret|
|00002e30| 75 72 6e 28 76 61 72 29 | 3b 0a 58 7d 0a 58 0a 58 |urn(var)|;.X}.X.X|
|00002e40| 0a 58 2f 2a 20 52 65 61 | 64 20 66 72 6f 6d 20 67 |.X/* Rea|d from g|
|00002e50| 69 76 65 6e 20 73 74 72 | 65 61 6d 20 74 6f 20 73 |iven str|eam to s|
|00002e60| 65 74 20 75 70 20 6c 6f | 67 70 72 6f 62 20 74 61 |et up lo|gprob ta|
|00002e70| 62 6c 65 20 61 6e 64 20 | 63 6f 6e 73 74 61 6e 74 |ble and |constant|
|00002e80| 73 0a 58 20 2a 20 6c 6f | 67 6d 65 61 6e 20 61 6e |s.X * lo|gmean an|
|00002e90| 64 20 6c 6f 67 76 61 72 | 2e 0a 58 20 2a 0a 58 20 |d logvar|..X *.X |
|00002ea0| 2a 20 54 68 65 20 74 61 | 62 6c 65 20 66 6f 72 6d |* The ta|ble form|
|00002eb0| 61 74 20 69 73 3a 0a 58 | 20 2a 20 3c 54 6f 74 61 |at is:.X| * <Tota|
|00002ec0| 6c 20 63 6f 75 6e 74 3e | 0a 58 20 2a 20 3c 42 6c |l count>|.X * <Bl|
|00002ed0| 61 6e 6b 6c 69 6e 65 3e | 0a 58 20 2a 20 3c 43 6f |ankline>|.X * <Co|
|00002ee0| 75 6e 74 3e 3c 73 70 61 | 63 65 3e 3c 4f 6e 65 20 |unt><spa|ce><One |
|00002ef0| 6f 72 20 6d 6f 72 65 20 | 73 6c 61 73 68 69 66 69 |or more |slashifi|
|00002f00| 65 64 20 63 68 61 72 61 | 63 74 65 72 73 20 74 6f |ed chara|cters to|
|00002f10| 20 73 68 61 72 65 20 74 | 68 61 74 20 63 6f 75 6e | share t|hat coun|
|00002f20| 74 3e 0a 58 20 2a 20 20 | 2e 2e 2e 0a 58 20 2a 20 |t>.X *  |....X * |
|00002f30| 3c 43 6f 75 6e 74 3e 3c | 73 70 61 63 65 3e 3c 4f |<Count><|space><O|
|00002f40| 6e 65 20 6f 72 20 6d 6f | 72 65 20 73 6c 61 73 68 |ne or mo|re slash|
|00002f50| 69 66 69 65 64 20 63 68 | 61 72 61 63 74 65 72 73 |ified ch|aracters|
|00002f60| 20 74 6f 20 73 68 61 72 | 65 20 74 68 61 74 20 63 | to shar|e that c|
|00002f70| 6f 75 6e 74 3e 0a 58 20 | 2a 20 3c 42 6c 61 6e 6b |ount>.X |* <Blank|
|00002f80| 6c 69 6e 65 3e 0a 58 20 | 2a 20 3c 45 4f 46 3e 0a |line>.X |* <EOF>.|
|00002f90| 58 20 2a 2f 0a 58 6c 6f | 61 64 5f 31 73 74 61 74 |X */.Xlo|ad_1stat|
|00002fa0| 73 28 69 6e 70 29 0a 58 | 46 49 4c 45 09 2a 69 6e |s(inp).X|FILE.*in|
|00002fb0| 70 3b 0a 58 7b 0a 58 09 | 69 6e 74 09 09 69 2c 6e |p;.X{.X.|int..i,n|
|00002fc0| 3b 0a 58 09 69 6e 74 09 | 09 74 6d 70 3b 0a 58 09 |;.X.int.|.tmp;.X.|
|00002fd0| 69 6e 74 09 09 63 3b 0a | 58 09 66 6c 6f 61 74 09 |int..c;.|X.float.|
|00002fe0| 76 2c 20 6c 76 2c 20 66 | 76 3b 0a 58 09 66 6c 6f |v, lv, f|v;.X.flo|
|00002ff0| 61 74 09 65 74 6f 74 61 | 6c 2c 20 63 74 6f 74 61 |at.etota|l, ctota|
|00003000| 6c 3b 0a 58 0a 58 09 73 | 74 61 74 73 31 6c 6f 61 |l;.X.X.s|tats1loa|
|00003010| 64 65 64 20 3d 20 54 52 | 55 45 3b 0a 58 0a 58 09 |ded = TR|UE;.X.X.|
|00003020| 66 6f 72 20 28 69 20 3d | 20 30 20 3b 20 69 20 3c |for (i =| 0 ; i <|
|00003030| 3d 20 4d 41 58 43 48 41 | 52 20 3b 20 69 2b 2b 29 |= MAXCHA|R ; i++)|
|00003040| 20 20 6c 6f 67 70 72 6f | 62 5b 69 5d 20 3d 20 30 |  logpro|b[i] = 0|
|00003050| 2e 30 3b 0a 58 0a 58 09 | 69 66 20 28 66 73 63 61 |.0;.X.X.|if (fsca|
|00003060| 6e 66 28 69 6e 70 2c 20 | 22 25 64 22 2c 20 26 74 |nf(inp, |"%d", &t|
|00003070| 6d 70 29 20 21 3d 20 31 | 29 20 20 7b 0a 58 09 09 |mp) != 1|)  {.X..|
|00003080| 70 72 69 6e 74 66 28 22 | 45 72 72 6f 72 20 77 68 |printf("|Error wh|
|00003090| 69 6c 65 20 67 65 74 74 | 69 6e 67 20 74 6f 74 61 |ile gett|ing tota|
|000030a0| 6c 22 29 3b 0a 58 09 09 | 72 65 74 75 72 6e 3b 0a |l");.X..|return;.|
|000030b0| 58 09 09 7d 0a 58 09 65 | 74 6f 74 61 6c 20 3d 20 |X..}.X.e|total = |
|000030c0| 74 6d 70 3b 0a 58 09 63 | 74 6f 74 61 6c 20 3d 20 |tmp;.X.c|total = |
|000030d0| 30 2e 30 3b 0a 58 0a 58 | 09 69 66 20 28 66 73 63 |0.0;.X.X|.if (fsc|
|000030e0| 61 6e 66 28 69 6e 70 2c | 20 22 5c 6e 22 29 20 21 |anf(inp,| "\n") !|
|000030f0| 3d 20 30 29 20 20 7b 0a | 58 09 09 70 72 69 6e 74 |= 0)  {.|X..print|
|00003100| 66 28 22 45 72 72 6f 72 | 20 77 68 69 6c 65 20 73 |f("Error| while s|
|00003110| 6b 69 70 70 69 6e 67 20 | 62 6c 61 6e 6b 20 6c 69 |kipping |blank li|
|00003120| 6e 65 22 29 3b 0a 58 09 | 09 72 65 74 75 72 6e 3b |ne");.X.|.return;|
|00003130| 0a 58 09 09 7d 0a 58 0a | 58 09 77 68 69 6c 65 20 |.X..}.X.|X.while |
|00003140| 28 54 52 55 45 29 20 7b | 0a 58 09 09 69 66 20 28 |(TRUE) {|.X..if (|
|00003150| 28 6e 20 3d 20 66 73 63 | 61 6e 66 28 69 6e 70 2c |(n = fsc|anf(inp,|
|00003160| 20 22 25 64 22 2c 20 26 | 74 6d 70 29 29 20 21 3d | "%d", &|tmp)) !=|
|00003170| 20 31 29 20 20 7b 0a 58 | 09 09 09 69 66 20 28 6e | 1)  {.X|...if (n|
|00003180| 20 3d 3d 20 30 29 20 62 | 72 65 61 6b 3b 0a 58 09 | == 0) b|reak;.X.|
|00003190| 09 09 69 66 20 28 6e 20 | 3d 3d 20 45 4f 46 29 20 |..if (n |== EOF) |
|000031a0| 62 72 65 61 6b 3b 0a 58 | 09 09 09 70 72 69 6e 74 |break;.X|...print|
|000031b0| 66 28 22 45 72 72 6f 72 | 20 77 68 69 6c 65 20 67 |f("Error| while g|
|000031c0| 65 74 74 69 6e 67 20 63 | 68 61 72 61 63 74 65 72 |etting c|haracter|
|000031d0| 20 63 6f 75 6e 74 22 29 | 3b 0a 58 09 09 09 72 65 | count")|;.X...re|
|000031e0| 74 75 72 6e 3b 0a 58 09 | 09 09 7d 0a 58 09 09 76 |turn;.X.|..}.X..v|
|000031f0| 20 3d 20 74 6d 70 3b 0a | 58 09 09 63 74 6f 74 61 | = tmp;.|X..ctota|
|00003200| 6c 20 2b 3d 20 76 3b 0a | 58 09 09 66 76 20 3d 20 |l += v;.|X..fv = |
|00003210| 76 2f 65 74 6f 74 61 6c | 3b 0a 58 09 09 69 66 20 |v/etotal|;.X..if |
|00003220| 28 66 76 20 21 3d 20 30 | 2e 30 29 20 20 7b 6c 76 |(fv != 0|.0)  {lv|
|00003230| 20 3d 20 6c 6f 67 31 30 | 28 66 76 29 3b 7d 0a 58 | = log10|(fv);}.X|
|00003240| 09 09 65 6c 73 65 20 7b | 6c 76 20 3d 20 30 2e 30 |..else {|lv = 0.0|
|00003250| 3b 7d 0a 58 0a 58 09 09 | 63 20 3d 20 72 65 61 64 |;}.X.X..|c = read|
|00003260| 5f 63 68 61 72 28 69 6e | 70 29 3b 09 09 2f 2a 20 |_char(in|p);../* |
|00003270| 53 6b 69 70 20 74 68 65 | 20 73 70 61 63 65 2e 20 |Skip the| space. |
|00003280| 2a 2f 0a 58 09 09 77 68 | 69 6c 65 20 28 54 52 55 |*/.X..wh|ile (TRU|
|00003290| 45 29 20 7b 0a 58 09 09 | 09 63 20 3d 20 72 65 61 |E) {.X..|.c = rea|
|000032a0| 64 5f 63 68 61 72 28 69 | 6e 70 29 3b 0a 58 09 09 |d_char(i|np);.X..|
|000032b0| 09 69 66 20 28 63 20 3d | 3d 20 45 4f 4c 29 20 20 |.if (c =|= EOL)  |
|000032c0| 62 72 65 61 6b 3b 0a 58 | 09 09 09 70 72 6f 62 5b |break;.X|...prob[|
|000032d0| 63 26 43 48 41 52 4d 41 | 53 4b 5d 20 3d 20 66 76 |c&CHARMA|SK] = fv|
|000032e0| 3b 0a 58 09 09 09 6c 6f | 67 70 72 6f 62 5b 63 26 |;.X...lo|gprob[c&|
|000032f0| 43 48 41 52 4d 41 53 4b | 5d 20 3d 20 6c 76 3b 0a |CHARMASK|] = lv;.|
|00003300| 58 09 09 09 7d 0a 58 09 | 09 7d 0a 58 0a 58 09 69 |X...}.X.|.}.X.X.i|
|00003310| 66 20 28 65 74 6f 74 61 | 6c 20 21 3d 20 63 74 6f |f (etota|l != cto|
|00003320| 74 61 6c 29 20 7b 0a 58 | 09 09 70 72 69 6e 74 66 |tal) {.X|..printf|
|00003330| 28 22 45 78 70 65 63 74 | 65 64 20 74 6f 74 61 6c |("Expect|ed total|
|00003340| 20 69 73 20 25 66 2e 20 | 20 41 63 74 75 61 6c 20 | is %f. | Actual |
|00003350| 74 6f 74 61 6c 20 69 73 | 20 25 66 2e 5c 6e 22 2c |total is| %f.\n",|
|00003360| 65 74 6f 74 61 6c 2c 63 | 74 6f 74 61 6c 29 3b 0a |etotal,c|total);.|
|00003370| 58 09 09 7d 0a 58 0a 58 | 09 6c 6f 67 6d 65 61 6e |X..}.X.X|.logmean|
|00003380| 20 3d 20 76 65 63 5f 6d | 65 61 6e 28 70 72 6f 62 | = vec_m|ean(prob|
|00003390| 2c 20 6c 6f 67 70 72 6f | 62 2c 20 4d 41 58 43 48 |, logpro|b, MAXCH|
|000033a0| 41 52 29 3b 0a 58 09 6c | 6f 67 76 61 72 20 20 3d |AR);.X.l|ogvar  =|
|000033b0| 20 76 65 63 5f 76 61 72 | 69 61 6e 63 65 28 70 72 | vec_var|iance(pr|
|000033c0| 6f 62 2c 20 6c 6f 67 70 | 72 6f 62 2c 20 4d 41 58 |ob, logp|rob, MAX|
|000033d0| 43 48 41 52 29 3b 0a 58 | 09 6c 6f 67 73 64 20 3d |CHAR);.X|.logsd =|
|000033e0| 20 73 71 72 74 28 6c 6f | 67 76 61 72 29 3b 0a 58 | sqrt(lo|gvar);.X|
|000033f0| 09 73 63 6f 72 65 31 5f | 6d 65 61 6e 20 3d 20 6c |.score1_|mean = l|
|00003400| 6f 67 6d 65 61 6e 3b 0a | 58 09 73 63 6f 72 65 31 |ogmean;.|X.score1|
|00003410| 5f 76 61 72 20 3d 20 6c | 6f 67 76 61 72 3b 0a 58 |_var = l|ogvar;.X|
|00003420| 09 73 63 6f 72 65 31 5f | 73 64 20 3d 20 6c 6f 67 |.score1_|sd = log|
|00003430| 73 64 3b 0a 58 09 73 63 | 6f 72 65 31 5f 73 63 61 |sd;.X.sc|ore1_sca|
|00003440| 6c 65 20 3d 20 73 71 72 | 74 28 32 20 2a 20 50 49 |le = sqr|t(2 * PI|
|00003450| 20 2a 20 73 63 6f 72 65 | 31 5f 76 61 72 29 3b 0a | * score|1_var);.|
|00003460| 58 09 70 6d 65 61 6e 20 | 3d 20 76 65 63 5f 6d 65 |X.pmean |= vec_me|
|00003470| 61 6e 28 70 72 6f 62 2c | 20 70 72 6f 62 2c 20 4d |an(prob,| prob, M|
|00003480| 41 58 43 48 41 52 29 3b | 0a 58 09 70 76 61 72 20 |AXCHAR);|.X.pvar |
|00003490| 3d 20 76 65 63 5f 76 61 | 72 69 61 6e 63 65 28 70 |= vec_va|riance(p|
|000034a0| 72 6f 62 2c 20 70 72 6f | 62 2c 20 4d 41 58 43 48 |rob, pro|b, MAXCH|
|000034b0| 41 52 29 3b 0a 58 7d 0a | 58 0a 58 0a 58 2f 2a 20 |AR);.X}.|X.X.X/* |
|000034c0| 4c 6f 61 64 20 74 68 65 | 20 6c 65 74 74 65 72 20 |Load the| letter |
|000034d0| 70 61 69 72 20 73 74 61 | 74 69 73 74 69 63 73 20 |pair sta|tistics |
|000034e0| 66 72 6f 6d 20 74 68 65 | 20 67 69 76 65 6e 20 66 |from the| given f|
|000034f0| 69 6c 65 20 6e 61 6d 65 | 2e 0a 58 20 2a 2f 0a 58 |ile name|..X */.X|
|00003500| 6c 6f 61 64 5f 32 73 74 | 61 74 73 5f 66 72 6f 6d |load_2st|ats_from|
|00003510| 28 73 74 61 74 66 6e 61 | 6d 65 29 0a 58 63 68 61 |(statfna|me).Xcha|
|00003520| 72 09 2a 73 74 61 74 66 | 6e 61 6d 65 3b 09 09 2f |r.*statf|name;../|
|00003530| 2a 20 46 75 6c 6c 20 70 | 61 74 68 20 6e 61 6d 65 |* Full p|ath name|
|00003540| 20 6f 66 20 66 69 6c 65 | 20 77 69 74 68 20 73 74 | of file| with st|
|00003550| 61 74 69 73 74 69 63 73 | 2e 20 2a 2f 0a 58 7b 0a |atistics|. */.X{.|
|00003560| 58 09 46 49 4c 45 09 2a | 69 6e 70 3b 0a 58 0a 58 |X.FILE.*|inp;.X.X|
|00003570| 09 69 66 20 28 28 69 6e | 70 20 3d 20 66 6f 70 65 |.if ((in|p = fope|
|00003580| 6e 28 73 74 61 74 66 6e | 61 6d 65 2c 20 22 72 22 |n(statfn|ame, "r"|
|00003590| 29 29 20 3d 3d 20 4e 55 | 4c 4c 29 20 7b 0a 58 09 |)) == NU|LL) {.X.|
|000035a0| 09 70 72 69 6e 74 66 28 | 22 5c 6e 43 61 6e 27 74 |.printf(|"\nCan't|
|000035b0| 20 6f 70 65 6e 20 25 73 | 20 74 6f 20 72 65 61 64 | open %s| to read|
|000035c0| 20 6c 65 74 74 65 72 20 | 73 74 61 74 69 73 74 69 | letter |statisti|
|000035d0| 63 73 5c 6e 22 2c 20 73 | 74 61 74 66 6e 61 6d 65 |cs\n", s|tatfname|
|000035e0| 29 3b 0a 58 09 09 65 78 | 69 74 28 30 29 3b 0a 58 |);.X..ex|it(0);.X|
|000035f0| 09 09 7d 0a 58 09 6c 6f | 61 64 5f 32 73 74 61 74 |..}.X.lo|ad_2stat|
|00003600| 73 28 69 6e 70 29 3b 0a | 58 09 66 63 6c 6f 73 65 |s(inp);.|X.fclose|
|00003610| 28 69 6e 70 29 3b 0a 58 | 7d 0a 58 0a 58 0a 58 2f |(inp);.X|}.X.X.X/|
|00003620| 2a 20 52 65 61 64 20 66 | 72 6f 6d 20 67 69 76 65 |* Read f|rom give|
|00003630| 6e 20 73 74 72 65 61 6d | 20 74 6f 20 73 65 74 20 |n stream| to set |
|00003640| 75 70 20 62 69 6c 6f 67 | 70 72 6f 62 20 74 61 62 |up bilog|prob tab|
|00003650| 6c 65 20 61 6e 64 20 63 | 6f 6e 73 74 61 6e 74 73 |le and c|onstants|
|00003660| 0a 58 20 2a 20 62 69 6c | 6f 67 6d 65 61 6e 2c 20 |.X * bil|ogmean, |
|00003670| 62 69 6c 6f 67 73 64 2c | 20 61 6e 64 20 62 69 6c |bilogsd,| and bil|
|00003680| 6f 67 76 61 72 2e 0a 58 | 20 2a 0a 58 20 2a 20 54 |ogvar..X| *.X * T|
|00003690| 68 65 20 66 6f 72 6d 61 | 74 20 6f 66 20 74 68 65 |he forma|t of the|
|000036a0| 20 73 74 61 74 69 73 74 | 69 63 73 20 66 69 6c 65 | statist|ics file|
|000036b0| 20 69 73 3a 20 5b 54 68 | 69 73 20 73 68 6f 75 6c | is: [Th|is shoul|
|000036c0| 64 20 62 65 20 6d 6f 72 | 65 20 67 65 6e 65 72 61 |d be mor|e genera|
|000036d0| 6c 2e 5d 0a 58 20 2a 20 | 3c 54 6f 74 61 6c 20 63 |l.].X * |<Total c|
|000036e0| 6f 75 6e 74 3e 0a 58 20 | 2a 20 3c 42 6c 61 6e 6b |ount>.X |* <Blank|
|000036f0| 6c 69 6e 65 3e 0a 58 20 | 2a 20 3c 73 69 6e 67 6c |line>.X |* <singl|
|00003700| 65 20 6c 65 74 74 65 72 | 20 63 6f 75 6e 74 73 3e |e letter| counts>|
|00003710| 0a 58 20 2a 20 3c 6c 69 | 6e 65 20 77 69 74 68 20 |.X * <li|ne with |
|00003720| 74 68 65 20 63 68 61 72 | 73 20 27 2a 2a 2a 27 3e |the char|s '***'>|
|00003730| 0a 58 20 2a 20 3c 64 6f | 75 62 6c 65 20 6c 65 74 |.X * <do|uble let|
|00003740| 74 65 72 20 63 6f 75 6e | 74 73 3e 0a 58 20 2a 20 |ter coun|ts>.X * |
|00003750| 3c 6c 69 6e 65 20 77 69 | 74 68 20 74 68 65 20 63 |<line wi|th the c|
|00003760| 68 61 72 73 20 27 2a 2a | 2a 27 3e 0a 58 20 2a 20 |hars '**|*'>.X * |
|00003770| 3c 6d 65 61 6e 20 6f 66 | 20 6d 61 74 72 69 78 3e |<mean of| matrix>|
|00003780| 0a 58 20 2a 20 3c 76 61 | 72 69 61 6e 63 65 20 6f |.X * <va|riance o|
|00003790| 66 20 6d 61 74 72 69 78 | 3e 0a 58 20 2a 20 3c 73 |f matrix|>.X * <s|
|000037a0| 74 61 6e 64 61 72 64 20 | 64 65 76 69 61 74 69 6f |tandard |deviatio|
|000037b0| 6e 20 6f 66 20 6d 61 74 | 72 69 78 3e 0a 58 20 2a |n of mat|rix>.X *|
|000037c0| 20 3c 42 6c 61 6e 6b 6c | 69 6e 65 3e 0a 58 20 2a | <Blankl|ine>.X *|
|000037d0| 20 3c 45 4f 46 3e 0a 58 | 20 2a 0a 58 20 2a 20 57 | <EOF>.X| *.X * W|
|000037e0| 68 65 72 65 20 73 69 6e | 67 6c 65 20 6c 65 74 74 |here sin|gle lett|
|000037f0| 65 72 20 63 6f 75 6e 74 | 73 20 61 6c 73 6f 20 64 |er count|s also d|
|00003800| 65 66 69 6e 65 20 74 68 | 65 20 6d 61 70 70 69 6e |efine th|e mappin|
|00003810| 67 20 66 72 6f 6d 20 61 | 73 63 69 69 20 63 68 61 |g from a|scii cha|
|00003820| 72 73 20 74 6f 0a 58 20 | 2a 20 64 69 73 74 69 6e |rs to.X |* distin|
|00003830| 67 75 69 73 68 65 64 20 | 6c 65 74 74 65 72 73 20 |guished |letters |
|00003840| 28 69 2e 65 2e 2c 20 61 | 6c 6c 20 6f 70 65 6e 20 |(i.e., a|ll open |
|00003850| 62 72 61 63 6b 65 74 73 | 20 61 72 65 20 74 72 65 |brackets| are tre|
|00003860| 61 74 65 64 20 74 68 65 | 20 73 61 6d 65 29 2e 0a |ated the| same)..|
|00003870| 58 20 2a 20 54 68 65 20 | 73 69 6e 67 6c 65 20 6c |X * The |single l|
|00003880| 65 74 74 65 72 20 66 6f | 72 6d 61 74 20 69 73 3a |etter fo|rmat is:|
|00003890| 0a 58 20 2a 20 3c 43 6f | 75 6e 74 3e 3c 73 70 61 |.X * <Co|unt><spa|
|000038a0| 63 65 3e 3c 4f 6e 65 20 | 6f 72 20 6d 6f 72 65 20 |ce><One |or more |
|000038b0| 73 6c 61 73 68 69 66 69 | 65 64 20 63 68 61 72 61 |slashifi|ed chara|
|000038c0| 63 74 65 72 73 20 74 6f | 20 73 68 61 72 65 20 74 |cters to| share t|
|000038d0| 68 61 74 20 63 6f 75 6e | 74 3e 0a 58 20 2a 20 20 |hat coun|t>.X *  |
|000038e0| 2e 2e 2e 0a 58 20 2a 20 | 3c 43 6f 75 6e 74 3e 3c |....X * |<Count><|
|000038f0| 73 70 61 63 65 3e 3c 4f | 6e 65 20 6f 72 20 6d 6f |space><O|ne or mo|
|00003900| 72 65 20 73 6c 61 73 68 | 69 66 69 65 64 20 63 68 |re slash|ified ch|
|00003910| 61 72 61 63 74 65 72 73 | 20 74 6f 20 73 68 61 72 |aracters| to shar|
|00003920| 65 20 74 68 61 74 20 63 | 6f 75 6e 74 3e 0a 58 20 |e that c|ount>.X |
|00003930| 2a 20 4e 4f 54 45 3a 20 | 74 68 65 20 66 69 72 73 |* NOTE: |the firs|
|00003940| 74 20 65 6e 74 72 79 20 | 73 68 6f 75 6c 64 20 62 |t entry |should b|
|00003950| 65 20 66 6f 72 20 61 20 | 6c 6f 77 20 70 72 6f 62 |e for a |low prob|
|00003960| 61 62 69 6c 69 74 79 20 | 6c 65 74 74 65 72 20 62 |ability |letter b|
|00003970| 65 63 61 75 73 65 20 74 | 68 65 0a 58 20 2a 20 64 |ecause t|he.X * d|
|00003980| 65 66 61 75 6c 74 20 6d | 61 70 70 69 6e 67 20 66 |efault m|apping f|
|00003990| 6f 72 20 75 6e 6b 6e 6f | 77 6e 20 63 68 61 72 73 |or unkno|wn chars|
|000039a0| 20 69 73 20 7a 65 72 6f | 2e 20 20 53 65 65 20 63 | is zero|.  See c|
|000039b0| 6f 64 65 20 66 6f 72 20 | 64 65 74 61 69 6c 73 2e |ode for |details.|
|000039c0| 0a 58 20 2a 0a 58 20 2a | 20 54 68 65 20 64 6f 75 |.X *.X *| The dou|
|000039d0| 62 6c 65 20 6c 65 74 74 | 65 72 20 66 6f 72 6d 61 |ble lett|er forma|
|000039e0| 74 20 69 73 3a 0a 58 20 | 2a 20 3c 43 6f 75 6e 74 |t is:.X |* <Count|
|000039f0| 3e 3c 73 70 61 63 65 3e | 3c 52 65 70 72 65 73 65 |><space>|<Represe|
|00003a00| 6e 74 61 74 69 76 65 20 | 6f 66 20 66 69 72 73 74 |ntative |of first|
|00003a10| 20 6c 65 74 74 65 72 20 | 67 72 6f 75 70 3e 3c 52 | letter |group><R|
|00003a20| 65 70 20 6f 66 20 73 65 | 63 6f 6e 64 20 6c 65 74 |ep of se|cond let|
|00003a30| 74 65 72 3e 0a 58 20 2a | 20 20 2e 2e 2e 0a 58 20 |ter>.X *|  ....X |
|00003a40| 2a 20 3c 43 6f 75 6e 74 | 3e 3c 73 70 61 63 65 3e |* <Count|><space>|
|00003a50| 3c 52 65 70 72 65 73 65 | 6e 74 61 74 69 76 65 20 |<Represe|ntative |
|00003a60| 6f 66 20 66 69 72 73 74 | 20 6c 65 74 74 65 72 20 |of first| letter |
|00003a70| 67 72 6f 75 70 3e 3c 52 | 65 70 20 6f 66 20 73 65 |group><R|ep of se|
|00003a80| 63 6f 6e 64 20 6c 65 74 | 74 65 72 3e 0a 58 20 2a |cond let|ter>.X *|
|00003a90| 0a 58 20 2a 20 46 6f 72 | 20 65 78 61 6d 70 6c 65 |.X * For| example|
|00003aa0| 20 69 66 20 27 54 27 20 | 61 6e 64 20 27 74 27 20 | if 'T' |and 't' |
|00003ab0| 61 72 65 20 74 72 65 61 | 74 65 64 20 74 68 65 20 |are trea|ted the |
|00003ac0| 73 61 6d 65 2c 20 61 20 | 64 6f 75 62 6c 65 20 6c |same, a |double l|
|00003ad0| 65 74 74 65 72 20 65 6e | 74 72 79 0a 58 20 2a 20 |etter en|try.X * |
|00003ae0| 6d 69 67 68 74 20 6c 6f | 6f 6b 20 6c 69 6b 65 3a |might lo|ok like:|
|00003af0| 20 22 31 32 34 37 20 54 | 54 22 20 61 6e 64 20 63 | "1247 T|T" and c|
|00003b00| 6f 75 6e 74 20 66 6f 72 | 20 54 74 2c 20 74 54 2c |ount for| Tt, tT,|
|00003b10| 20 74 74 2c 20 61 6e 64 | 20 54 54 2e 0a 58 20 2a | tt, and| TT..X *|
|00003b20| 2f 0a 58 6c 6f 61 64 5f | 32 73 74 61 74 73 28 69 |/.Xload_|2stats(i|
|00003b30| 6e 70 29 0a 58 46 49 4c | 45 09 2a 69 6e 70 3b 0a |np).XFIL|E.*inp;.|
|00003b40| 58 7b 0a 58 72 65 67 69 | 73 74 65 72 09 69 6e 74 |X{.Xregi|ster.int|
|00003b50| 09 09 69 2c 6a 3b 0a 58 | 20 20 20 20 69 6e 74 09 |..i,j;.X|    int.|
|00003b60| 09 6e 3b 0a 58 09 69 6e | 74 09 09 74 6d 70 3b 0a |.n;.X.in|t..tmp;.|
|00003b70| 58 09 69 6e 74 09 09 63 | 3b 0a 58 09 69 6e 74 09 |X.int..c|;.X.int.|
|00003b80| 09 6c 65 66 74 5f 69 6e | 64 65 78 2c 20 72 69 67 |.left_in|dex, rig|
|00003b90| 68 74 5f 69 6e 64 65 78 | 3b 0a 58 09 66 6c 6f 61 |ht_index|;.X.floa|
|00003ba0| 74 09 76 2c 20 6c 76 2c | 20 66 76 3b 0a 58 09 66 |t.v, lv,| fv;.X.f|
|00003bb0| 6c 6f 61 74 09 65 74 6f | 74 61 6c 2c 20 63 74 6f |loat.eto|tal, cto|
|00003bc0| 74 61 6c 3b 0a 58 09 63 | 68 61 72 09 6c 69 6e 65 |tal;.X.c|har.line|
|00003bd0| 62 75 66 5b 33 30 30 5d | 3b 0a 58 0a 58 09 73 74 |buf[300]|;.X.X.st|
|00003be0| 61 74 73 32 6c 6f 61 64 | 65 64 20 3d 20 54 52 55 |ats2load|ed = TRU|
|00003bf0| 45 3b 0a 58 09 6e 62 69 | 63 68 61 72 73 20 3d 20 |E;.X.nbi|chars = |
|00003c00| 30 3b 0a 58 0a 58 09 66 | 6f 72 20 28 69 20 3d 20 |0;.X.X.f|or (i = |
|00003c10| 30 20 3b 20 69 20 3c 20 | 4d 58 42 49 49 4e 44 45 |0 ; i < |MXBIINDE|
|00003c20| 58 20 3b 20 69 2b 2b 29 | 20 20 7b 0a 58 09 09 73 |X ; i++)|  {.X..s|
|00003c30| 6c 6c 6f 67 70 72 6f 62 | 5b 69 5d 20 3d 20 30 2e |llogprob|[i] = 0.|
|00003c40| 30 3b 0a 58 09 09 73 6c | 62 69 70 72 6f 62 5b 69 |0;.X..sl|biprob[i|
|00003c50| 5d 20 3d 20 30 2e 30 3b | 0a 58 09 09 66 6f 72 20 |] = 0.0;|.X..for |
|00003c60| 28 6a 20 3d 20 30 20 3b | 20 6a 20 3c 20 4d 58 42 |(j = 0 ;| j < MXB|
|00003c70| 49 49 4e 44 45 58 20 3b | 20 6a 2b 2b 29 20 20 7b |IINDEX ;| j++)  {|
|00003c80| 0a 58 09 20 09 09 62 69 | 6c 6f 67 70 72 6f 62 5b |.X. ..bi|logprob[|
|00003c90| 69 5d 5b 6a 5d 20 3d 20 | 30 2e 30 3b 0a 58 09 09 |i][j] = |0.0;.X..|
|00003ca0| 09 62 69 70 72 6f 62 5b | 69 5d 5b 6a 5d 20 3d 20 |.biprob[|i][j] = |
|00003cb0| 30 2e 30 3b 0a 58 09 09 | 09 7d 0a 58 09 09 7d 0a |0.0;.X..|.}.X..}.|
|00003cc0| 58 0a 58 09 66 6f 72 20 | 28 69 20 3d 20 30 20 3b |X.X.for |(i = 0 ;|
|00003cd0| 20 69 20 3c 20 4d 41 58 | 43 48 41 52 2b 31 20 3b | i < MAX|CHAR+1 ;|
|00003ce0| 20 69 2b 2b 29 0a 58 09 | 09 63 68 61 72 5f 62 69 | i++).X.|.char_bi|
|00003cf0| 6d 61 70 5b 69 5d 20 3d | 20 30 3b 09 09 2f 2a 20 |map[i] =| 0;../* |
|00003d00| 44 65 66 61 75 6c 74 20 | 69 6e 64 65 78 20 69 66 |Default |index if|
|00003d10| 20 63 68 61 72 20 75 6e | 6b 6e 6f 77 6e 2e 20 2a | char un|known. *|
|00003d20| 2f 0a 58 0a 58 09 69 66 | 20 28 66 73 63 61 6e 66 |/.X.X.if| (fscanf|
|00003d30| 28 69 6e 70 2c 20 22 25 | 64 22 2c 20 26 74 6d 70 |(inp, "%|d", &tmp|
|00003d40| 29 20 21 3d 20 31 29 20 | 20 7b 0a 58 09 09 70 72 |) != 1) | {.X..pr|
|00003d50| 69 6e 74 66 28 22 45 72 | 72 6f 72 20 77 68 69 6c |intf("Er|ror whil|
|00003d60| 65 20 67 65 74 74 69 6e | 67 20 74 6f 74 61 6c 22 |e gettin|g total"|
|00003d70| 29 3b 0a 58 09 09 65 78 | 69 74 28 30 29 3b 0a 58 |);.X..ex|it(0);.X|
|00003d80| 09 09 7d 0a 58 09 65 74 | 6f 74 61 6c 20 3d 20 74 |..}.X.et|otal = t|
|00003d90| 6d 70 3b 0a 58 0a 58 09 | 69 66 20 28 66 73 63 61 |mp;.X.X.|if (fsca|
|00003da0| 6e 66 28 69 6e 70 2c 20 | 22 5c 6e 22 29 20 21 3d |nf(inp, |"\n") !=|
|00003db0| 20 30 29 20 20 7b 0a 58 | 09 09 70 72 69 6e 74 66 | 0)  {.X|..printf|
|00003dc0| 28 22 45 72 72 6f 72 20 | 77 68 69 6c 65 20 73 6b |("Error |while sk|
|00003dd0| 69 70 70 69 6e 67 20 62 | 6c 61 6e 6b 20 6c 69 6e |ipping b|lank lin|
|00003de0| 65 20 62 65 66 6f 72 65 | 20 73 69 6e 67 6c 65 20 |e before| single |
|00003df0| 6c 65 74 74 65 72 73 22 | 29 3b 0a 58 09 09 65 78 |letters"|);.X..ex|
|00003e00| 69 74 28 30 29 3b 0a 58 | 09 09 7d 0a 58 0a 58 09 |it(0);.X|..}.X.X.|
|00003e10| 63 74 6f 74 61 6c 20 3d | 20 30 2e 30 3b 0a 58 09 |ctotal =| 0.0;.X.|
|00003e20| 77 68 69 6c 65 20 28 54 | 52 55 45 29 20 7b 0a 58 |while (T|RUE) {.X|
|00003e30| 09 09 69 66 20 28 28 6e | 20 3d 20 66 73 63 61 6e |..if ((n| = fscan|
|00003e40| 66 28 69 6e 70 2c 20 22 | 25 64 22 2c 20 26 74 6d |f(inp, "|%d", &tm|
|00003e50| 70 29 29 20 21 3d 20 31 | 29 20 20 7b 0a 58 09 09 |p)) != 1|)  {.X..|
|00003e60| 09 69 66 20 28 6e 20 3d | 3d 20 30 29 20 62 72 65 |.if (n =|= 0) bre|
|00003e70| 61 6b 3b 0a 58 09 09 09 | 69 66 20 28 6e 20 3d 3d |ak;.X...|if (n ==|
|00003e80| 20 45 4f 46 29 20 62 72 | 65 61 6b 3b 0a 58 09 09 | EOF) br|eak;.X..|
|00003e90| 09 70 72 69 6e 74 66 28 | 22 45 72 72 6f 72 20 77 |.printf(|"Error w|
|00003ea0| 68 69 6c 65 20 67 65 74 | 74 69 6e 67 20 63 68 61 |hile get|ting cha|
|00003eb0| 72 61 63 74 65 72 20 63 | 6f 75 6e 74 20 28 73 69 |racter c|ount (si|
|00003ec0| 6e 67 6c 65 73 29 22 29 | 3b 0a 58 09 09 09 65 78 |ngles)")|;.X...ex|
|00003ed0| 69 74 28 30 29 3b 0a 58 | 09 09 09 7d 0a 58 09 09 |it(0);.X|...}.X..|
|00003ee0| 76 20 3d 20 74 6d 70 3b | 0a 58 09 09 63 74 6f 74 |v = tmp;|.X..ctot|
|00003ef0| 61 6c 20 2b 3d 20 76 3b | 0a 58 09 09 66 76 20 3d |al += v;|.X..fv =|
|00003f00| 20 76 2f 65 74 6f 74 61 | 6c 3b 0a 58 09 09 69 66 | v/etota|l;.X..if|
|00003f10| 20 28 66 76 20 3d 3d 20 | 30 2e 30 29 0a 58 09 09 | (fv == |0.0).X..|
|00003f20| 09 6c 76 20 3d 20 30 2e | 30 3b 0a 58 09 09 65 6c |.lv = 0.|0;.X..el|
|00003f30| 73 65 0a 58 09 09 09 6c | 76 20 3d 20 6c 6f 67 31 |se.X...l|v = log1|
|00003f40| 30 28 66 76 29 3b 0a 58 | 0a 58 09 09 63 20 3d 20 |0(fv);.X|.X..c = |
|00003f50| 72 65 61 64 5f 63 68 61 | 72 28 69 6e 70 29 3b 09 |read_cha|r(inp);.|
|00003f60| 09 2f 2a 20 53 6b 69 70 | 20 74 68 65 20 73 70 61 |./* Skip| the spa|
|00003f70| 63 65 2e 20 2a 2f 0a 58 | 09 09 77 68 69 6c 65 20 |ce. */.X|..while |
|00003f80| 28 54 52 55 45 29 20 7b | 0a 58 09 09 09 63 20 3d |(TRUE) {|.X...c =|
|00003f90| 20 72 65 61 64 5f 63 68 | 61 72 28 69 6e 70 29 3b | read_ch|ar(inp);|
|00003fa0| 0a 58 09 09 09 69 66 20 | 28 63 20 3d 3d 20 45 4f |.X...if |(c == EO|
|00003fb0| 4c 29 20 20 62 72 65 61 | 6b 3b 0a 58 09 09 09 63 |L)  brea|k;.X...c|
|00003fc0| 68 61 72 5f 62 69 6d 61 | 70 5b 63 20 26 20 43 48 |har_bima|p[c & CH|
|00003fd0| 41 52 4d 41 53 4b 5d 20 | 3d 20 6e 62 69 63 68 61 |ARMASK] |= nbicha|
|00003fe0| 72 73 3b 0a 58 09 09 09 | 73 6c 62 69 70 72 6f 62 |rs;.X...|slbiprob|
|00003ff0| 5b 6e 62 69 63 68 61 72 | 73 5d 20 3d 20 66 76 3b |[nbichar|s] = fv;|
|00004000| 0a 58 09 09 09 73 6c 6c | 6f 67 70 72 6f 62 5b 6e |.X...sll|ogprob[n|
|00004010| 62 69 63 68 61 72 73 5d | 20 3d 20 6c 76 3b 0a 58 |bichars]| = lv;.X|
|00004020| 09 09 09 7d 0a 58 09 09 | 6e 62 69 63 68 61 72 73 |...}.X..|nbichars|
|00004030| 2b 2b 3b 0a 58 09 09 7d | 0a 58 0a 58 09 69 66 20 |++;.X..}|.X.X.if |
|00004040| 28 65 74 6f 74 61 6c 20 | 21 3d 20 63 74 6f 74 61 |(etotal |!= ctota|
|00004050| 6c 29 20 7b 0a 58 09 09 | 70 72 69 6e 74 66 28 22 |l) {.X..|printf("|
|00004060| 45 78 70 65 63 74 65 64 | 20 74 6f 74 61 6c 20 69 |Expected| total i|
|00004070| 73 20 25 66 2e 20 20 41 | 63 74 75 61 6c 20 74 6f |s %f.  A|ctual to|
|00004080| 74 61 6c 20 69 73 20 25 | 66 20 66 6f 72 20 73 69 |tal is %|f for si|
|00004090| 6e 67 6c 65 73 2e 5c 6e | 22 2c 0a 58 09 09 09 09 |ngles.\n|",.X....|
|000040a0| 65 74 6f 74 61 6c 2c 20 | 63 74 6f 74 61 6c 29 3b |etotal, |ctotal);|
|000040b0| 0a 58 09 09 65 78 69 74 | 28 30 29 3b 0a 58 09 09 |.X..exit|(0);.X..|
|000040c0| 7d 0a 58 0a 58 0a 58 09 | 69 66 20 28 66 73 63 61 |}.X.X.X.|if (fsca|
|000040d0| 6e 66 28 69 6e 70 2c 20 | 22 2a 2a 2a 5c 6e 22 29 |nf(inp, |"***\n")|
|000040e0| 20 21 3d 20 30 29 20 20 | 7b 0a 58 09 09 70 72 69 | != 0)  |{.X..pri|
|000040f0| 6e 74 66 28 22 45 72 72 | 6f 72 20 6f 6e 20 64 65 |ntf("Err|or on de|
|00004100| 6c 69 6d 69 74 65 72 20 | 62 65 66 6f 72 65 20 6c |limiter |before l|
|00004110| 65 74 74 65 72 20 70 61 | 69 72 73 22 29 3b 0a 58 |etter pa|irs");.X|
|00004120| 09 09 65 78 69 74 28 30 | 29 3b 0a 58 09 09 7d 0a |..exit(0|);.X..}.|
|00004130| 58 0a 58 09 63 74 6f 74 | 61 6c 20 3d 20 30 2e 30 |X.X.ctot|al = 0.0|
|00004140| 3b 0a 58 09 77 68 69 6c | 65 20 28 54 52 55 45 29 |;.X.whil|e (TRUE)|
|00004150| 20 7b 0a 58 09 09 69 66 | 20 28 28 6e 20 3d 20 66 | {.X..if| ((n = f|
|00004160| 73 63 61 6e 66 28 69 6e | 70 2c 20 22 25 64 22 2c |scanf(in|p, "%d",|
|00004170| 20 26 74 6d 70 29 29 20 | 21 3d 20 31 29 20 20 7b | &tmp)) |!= 1)  {|
|00004180| 0a 58 09 09 09 69 66 20 | 28 6e 20 3d 3d 20 30 29 |.X...if |(n == 0)|
|00004190| 20 62 72 65 61 6b 3b 0a | 58 09 09 09 69 66 20 28 | break;.|X...if (|
|000041a0| 6e 20 3d 3d 20 45 4f 46 | 29 20 62 72 65 61 6b 3b |n == EOF|) break;|
|000041b0| 0a 58 09 09 09 70 72 69 | 6e 74 66 28 22 45 72 72 |.X...pri|ntf("Err|
|000041c0| 6f 72 20 77 68 69 6c 65 | 20 67 65 74 74 69 6e 67 |or while| getting|
|000041d0| 20 63 68 61 72 61 63 74 | 65 72 20 63 6f 75 6e 74 | charact|er count|
|000041e0| 20 28 70 61 69 72 73 29 | 22 29 3b 0a 58 09 09 09 | (pairs)|");.X...|
|000041f0| 65 78 69 74 28 30 29 3b | 0a 58 09 09 09 7d 0a 58 |exit(0);|.X...}.X|
|00004200| 09 09 76 20 3d 20 74 6d | 70 3b 0a 58 09 09 63 74 |..v = tm|p;.X..ct|
|00004210| 6f 74 61 6c 20 2b 3d 20 | 76 3b 0a 58 09 09 66 76 |otal += |v;.X..fv|
|00004220| 20 3d 20 76 2f 65 74 6f | 74 61 6c 3b 0a 58 09 09 | = v/eto|tal;.X..|
|00004230| 69 66 20 28 66 76 20 3d | 3d 20 30 2e 30 29 0a 58 |if (fv =|= 0.0).X|
|00004240| 09 09 09 6c 76 20 3d 20 | 30 2e 30 3b 0a 58 09 09 |...lv = |0.0;.X..|
|00004250| 65 6c 73 65 0a 58 09 09 | 09 6c 76 20 3d 20 6c 6f |else.X..|.lv = lo|
|00004260| 67 31 30 28 66 76 29 3b | 0a 58 0a 58 09 09 63 20 |g10(fv);|.X.X..c |
|00004270| 3d 20 72 65 61 64 5f 63 | 68 61 72 28 69 6e 70 29 |= read_c|har(inp)|
|00004280| 3b 09 09 2f 2a 20 53 6b | 69 70 20 74 68 65 20 73 |;../* Sk|ip the s|
|00004290| 70 61 63 65 2e 20 2a 2f | 0a 58 09 09 63 20 3d 20 |pace. */|.X..c = |
|000042a0| 72 65 61 64 5f 63 68 61 | 72 28 69 6e 70 29 3b 09 |read_cha|r(inp);.|
|000042b0| 09 2f 2a 20 46 69 72 73 | 74 20 6c 65 74 74 65 72 |./* Firs|t letter|
|000042c0| 2e 20 2a 2f 0a 58 09 09 | 69 66 20 28 63 20 3d 3d |. */.X..|if (c ==|
|000042d0| 20 45 4f 4c 29 20 20 7b | 0a 58 09 09 09 70 72 69 | EOL)  {|.X...pri|
|000042e0| 6e 74 66 28 22 4c 69 6e | 65 20 65 6e 64 73 20 62 |ntf("Lin|e ends b|
|000042f0| 65 66 6f 72 65 20 6c 65 | 74 74 65 72 20 70 61 69 |efore le|tter pai|
|00004300| 72 22 29 3b 0a 58 09 09 | 09 65 78 69 74 28 30 29 |r");.X..|.exit(0)|
|00004310| 3b 0a 58 09 09 09 7d 0a | 58 09 09 6c 65 66 74 5f |;.X...}.|X..left_|
|00004320| 69 6e 64 65 78 20 3d 20 | 63 68 61 72 5f 62 69 6d |index = |char_bim|
|00004330| 61 70 5b 63 20 26 20 43 | 48 41 52 4d 41 53 4b 5d |ap[c & C|HARMASK]|
|00004340| 3b 0a 58 09 09 63 20 3d | 20 72 65 61 64 5f 63 68 |;.X..c =| read_ch|
|00004350| 61 72 28 69 6e 70 29 3b | 09 09 2f 2a 20 53 65 63 |ar(inp);|../* Sec|
|00004360| 6f 6e 64 20 6c 65 74 74 | 65 72 2e 20 2a 2f 0a 58 |ond lett|er. */.X|
|00004370| 09 09 69 66 20 28 63 20 | 3d 3d 20 45 4f 4c 29 20 |..if (c |== EOL) |
|00004380| 20 7b 0a 58 09 09 09 70 | 72 69 6e 74 66 28 22 4c | {.X...p|rintf("L|
|00004390| 69 6e 65 20 65 6e 64 73 | 20 69 6e 20 6d 69 64 64 |ine ends| in midd|
|000043a0| 6c 65 20 6f 66 20 6c 65 | 74 74 65 72 20 70 61 69 |le of le|tter pai|
|000043b0| 72 22 29 3b 0a 58 09 09 | 09 65 78 69 74 28 30 29 |r");.X..|.exit(0)|
|000043c0| 3b 0a 58 09 09 09 7d 0a | 58 09 09 72 69 67 68 74 |;.X...}.|X..right|
|000043d0| 5f 69 6e 64 65 78 20 3d | 20 63 68 61 72 5f 62 69 |_index =| char_bi|
|000043e0| 6d 61 70 5b 63 20 26 20 | 43 48 41 52 4d 41 53 4b |map[c & |CHARMASK|
|000043f0| 5d 3b 0a 58 0a 58 09 09 | 62 69 70 72 6f 62 5b 6c |];.X.X..|biprob[l|
|00004400| 65 66 74 5f 69 6e 64 65 | 78 5d 5b 72 69 67 68 74 |eft_inde|x][right|
|00004410| 5f 69 6e 64 65 78 5d 20 | 3d 20 66 76 3b 0a 58 09 |_index] |= fv;.X.|
|00004420| 09 62 69 6c 6f 67 70 72 | 6f 62 5b 6c 65 66 74 5f |.bilogpr|ob[left_|
|00004430| 69 6e 64 65 78 5d 5b 72 | 69 67 68 74 5f 69 6e 64 |index][r|ight_ind|
|00004440| 65 78 5d 20 3d 20 6c 76 | 3b 0a 58 09 09 7d 0a 58 |ex] = lv|;.X..}.X|
|00004450| 0a 58 09 69 66 20 28 65 | 74 6f 74 61 6c 20 21 3d |.X.if (e|total !=|
|00004460| 20 63 74 6f 74 61 6c 29 | 20 7b 0a 58 09 09 70 72 | ctotal)| {.X..pr|
|00004470| 69 6e 74 66 28 22 45 78 | 70 65 63 74 65 64 20 74 |intf("Ex|pected t|
|00004480| 6f 74 61 6c 20 69 73 20 | 25 66 2e 20 20 41 63 74 |otal is |%f.  Act|
|00004490| 75 61 6c 20 74 6f 74 61 | 6c 20 69 73 20 25 66 20 |ual tota|l is %f |
|000044a0| 66 6f 72 20 70 61 69 72 | 73 2e 5c 6e 22 2c 0a 58 |for pair|s.\n",.X|
|000044b0| 09 09 09 09 65 74 6f 74 | 61 6c 2c 20 63 74 6f 74 |....etot|al, ctot|
|000044c0| 61 6c 29 3b 0a 58 09 09 | 65 78 69 74 28 30 29 3b |al);.X..|exit(0);|
|000044d0| 0a 58 09 09 7d 0a 58 0a | 58 09 69 66 20 28 66 73 |.X..}.X.|X.if (fs|
|000044e0| 63 61 6e 66 28 69 6e 70 | 2c 20 22 2a 2a 2a 5c 6e |canf(inp|, "***\n|
|000044f0| 22 29 20 3d 3d 20 30 29 | 20 20 7b 0a 58 09 09 69 |") == 0)|  {.X..i|
|00004500| 66 20 28 66 73 63 61 6e | 66 28 69 6e 70 2c 20 22 |f (fscan|f(inp, "|
|00004510| 25 66 22 2c 20 26 73 63 | 6f 72 65 32 5f 6d 65 61 |%f", &sc|ore2_mea|
|00004520| 6e 29 20 21 3d 20 31 29 | 20 20 7b 0a 58 09 09 09 |n) != 1)|  {.X...|
|00004530| 70 72 69 6e 74 66 28 22 | 45 72 72 6f 72 20 72 65 |printf("|Error re|
|00004540| 61 64 69 6e 67 20 6d 65 | 61 6e 2e 22 29 3b 0a 58 |ading me|an.");.X|
|00004550| 09 09 09 65 78 69 74 28 | 30 29 3b 0a 58 09 09 09 |...exit(|0);.X...|
|00004560| 7d 0a 58 09 09 69 66 20 | 28 66 73 63 61 6e 66 28 |}.X..if |(fscanf(|
|00004570| 69 6e 70 2c 20 22 25 66 | 22 2c 20 26 73 63 6f 72 |inp, "%f|", &scor|
|00004580| 65 32 5f 76 61 72 29 20 | 21 3d 20 31 29 20 20 7b |e2_var) |!= 1)  {|
|00004590| 0a 58 09 09 09 70 72 69 | 6e 74 66 28 22 45 72 72 |.X...pri|ntf("Err|
|000045a0| 6f 72 20 72 65 61 64 69 | 6e 67 20 76 61 72 69 61 |or readi|ng varia|
|000045b0| 6e 63 65 2e 22 29 3b 0a | 58 09 09 09 65 78 69 74 |nce.");.|X...exit|
|000045c0| 28 30 29 3b 0a 58 09 09 | 09 7d 0a 58 09 09 69 66 |(0);.X..|.}.X..if|
|000045d0| 20 28 66 73 63 61 6e 66 | 28 69 6e 70 2c 20 22 25 | (fscanf|(inp, "%|
|000045e0| 66 22 2c 20 26 73 63 6f | 72 65 32 5f 73 64 29 20 |f", &sco|re2_sd) |
|000045f0| 21 3d 20 31 29 20 20 7b | 0a 58 09 09 09 70 72 69 |!= 1)  {|.X...pri|
|00004600| 6e 74 66 28 22 45 72 72 | 6f 72 20 72 65 61 64 69 |ntf("Err|or readi|
|00004610| 6e 67 20 73 74 61 6e 64 | 61 72 64 20 64 65 76 69 |ng stand|ard devi|
|00004620| 61 74 69 6f 6e 73 2e 22 | 29 3b 0a 58 09 09 09 65 |ations."|);.X...e|
|00004630| 78 69 74 28 30 29 3b 0a | 58 09 09 09 7d 0a 58 09 |xit(0);.|X...}.X.|
|00004640| 09 73 63 6f 72 65 32 5f | 73 63 61 6c 65 20 3d 20 |.score2_|scale = |
|00004650| 73 71 72 74 28 32 20 2a | 20 50 49 20 2a 20 73 63 |sqrt(2 *| PI * sc|
|00004660| 6f 72 65 32 5f 76 61 72 | 29 3b 0a 58 09 09 61 70 |ore2_var|);.X..ap|
|00004670| 70 72 6f 78 5f 69 6e 69 | 74 28 29 3b 0a 58 09 09 |prox_ini|t();.X..|
|00004680| 72 65 74 75 72 6e 3b 0a | 58 09 09 7d 0a 58 0a 58 |return;.|X..}.X.X|
|00004690| 09 73 74 61 74 73 32 28 | 29 3b 0a 58 09 70 72 69 |.stats2(|);.X.pri|
|000046a0| 6e 74 66 28 22 4d 65 61 | 6e 3a 20 25 66 2c 20 56 |ntf("Mea|n: %f, V|
|000046b0| 61 72 3a 20 25 66 2c 20 | 53 44 3a 20 25 66 5c 6e |ar: %f, |SD: %f\n|
|000046c0| 22 2c 20 73 63 6f 72 65 | 32 5f 6d 65 61 6e 2c 20 |", score|2_mean, |
|000046d0| 73 63 6f 72 65 32 5f 76 | 61 72 2c 20 73 63 6f 72 |score2_v|ar, scor|
|000046e0| 65 32 5f 73 64 29 3b 0a | 58 7d 0a 58 0a 58 0a 58 |e2_sd);.|X}.X.X.X|
|000046f0| 2f 2a 20 43 6f 6d 70 75 | 74 65 20 73 63 6f 72 69 |/* Compu|te scori|
|00004700| 6e 67 20 73 74 61 74 69 | 73 74 69 63 73 20 66 6f |ng stati|stics fo|
|00004710| 72 20 74 68 65 20 6c 65 | 74 74 65 72 20 70 61 69 |r the le|tter pai|
|00004720| 72 20 66 72 65 71 75 65 | 6e 63 69 65 73 2e 0a 58 |r freque|ncies..X|
|00004730| 20 2a 20 55 73 65 73 20 | 74 68 65 20 67 6c 6f 62 | * Uses |the glob|
|00004740| 61 6c 73 3a 20 62 69 70 | 72 6f 62 5b 5d 5b 5d 2c |als: bip|rob[][],|
|00004750| 20 73 6c 6c 6f 67 62 69 | 70 72 6f 62 5b 5d 2c 20 | sllogbi|prob[], |
|00004760| 61 6e 64 20 62 69 6c 6f | 67 70 72 6f 62 5b 5d 5b |and bilo|gprob[][|
|00004770| 5d 2e 0a 58 20 2a 20 53 | 65 74 73 20 67 6f 62 61 |]..X * S|ets goba|
|00004780| 6c 73 3a 20 73 63 6f 72 | 65 32 5f 6d 65 61 6e 2c |ls: scor|e2_mean,|
|00004790| 20 73 63 6f 72 65 32 5f | 76 61 72 2c 20 73 63 6f | score2_|var, sco|
|000047a0| 72 65 32 5f 73 64 2e 0a | 58 20 2a 2f 0a 58 73 74 |re2_sd..|X */.Xst|
|000047b0| 61 74 73 32 28 29 0a 58 | 7b 0a 58 72 65 67 69 73 |ats2().X|{.Xregis|
|000047c0| 74 65 72 09 69 6e 74 09 | 69 2c 6a 2c 6b 3b 0a 58 |ter.int.|i,j,k;.X|
|000047d0| 09 66 6c 6f 61 74 09 6d | 65 61 6e 2c 20 76 61 72 |.float.m|ean, var|
|000047e0| 3b 0a 58 09 66 6c 6f 61 | 74 09 77 65 69 67 68 74 |;.X.floa|t.weight|
|000047f0| 2c 20 73 63 6f 72 65 3b | 0a 58 0a 58 09 6d 65 61 |, score;|.X.X.mea|
|00004800| 6e 20 3d 20 30 2e 30 3b | 0a 58 09 76 61 72 20 3d |n = 0.0;|.X.var =|
|00004810| 20 30 2e 30 3b 0a 58 09 | 66 6f 72 20 28 69 20 3d | 0.0;.X.|for (i =|
|00004820| 20 30 20 3b 20 69 20 3c | 20 6e 62 69 63 68 61 72 | 0 ; i <| nbichar|
|00004830| 73 20 3b 20 69 2b 2b 29 | 0a 58 09 09 66 6f 72 20 |s ; i++)|.X..for |
|00004840| 28 6a 20 3d 20 30 20 3b | 20 6a 20 3c 20 6e 62 69 |(j = 0 ;| j < nbi|
|00004850| 63 68 61 72 73 20 3b 20 | 6a 2b 2b 29 20 20 7b 0a |chars ; |j++)  {.|
|00004860| 58 09 09 09 69 66 20 28 | 73 6c 62 69 70 72 6f 62 |X...if (|slbiprob|
|00004870| 5b 6a 5d 20 3d 3d 20 30 | 2e 30 29 0a 58 09 09 09 |[j] == 0|.0).X...|
|00004880| 09 63 6f 6e 74 69 6e 75 | 65 3b 0a 58 09 09 09 66 |.continu|e;.X...f|
|00004890| 6f 72 20 28 6b 20 3d 20 | 30 20 3b 20 6b 20 3c 20 |or (k = |0 ; k < |
|000048a0| 6e 62 69 63 68 61 72 73 | 20 3b 20 6b 2b 2b 29 20 |nbichars| ; k++) |
|000048b0| 7b 0a 58 09 09 09 09 77 | 65 69 67 68 74 20 3d 20 |{.X....w|eight = |
|000048c0| 62 69 70 72 6f 62 5b 69 | 5d 5b 6a 5d 20 2a 20 62 |biprob[i|][j] * b|
|000048d0| 69 70 72 6f 62 5b 6a 5d | 5b 6b 5d 20 2f 20 73 6c |iprob[j]|[k] / sl|
|000048e0| 62 69 70 72 6f 62 5b 6a | 5d 3b 0a 58 09 09 09 09 |biprob[j|];.X....|
|000048f0| 73 63 6f 72 65 20 3d 20 | 62 69 6c 6f 67 70 72 6f |score = |bilogpro|
|00004900| 62 5b 69 5d 5b 6a 5d 20 | 2b 20 62 69 6c 6f 67 70 |b[i][j] |+ bilogp|
|00004910| 72 6f 62 5b 6a 5d 5b 6b | 5d 20 2d 20 32 20 2a 20 |rob[j][k|] - 2 * |
|00004920| 73 6c 6c 6f 67 70 72 6f | 62 5b 6a 5d 3b 0a 58 09 |sllogpro|b[j];.X.|
|00004930| 09 09 09 6d 65 61 6e 20 | 2b 3d 20 77 65 69 67 68 |...mean |+= weigh|
|00004940| 74 20 2a 20 73 63 6f 72 | 65 3b 0a 58 09 09 09 09 |t * scor|e;.X....|
|00004950| 76 61 72 20 2b 3d 20 77 | 65 69 67 68 74 20 2a 20 |var += w|eight * |
|00004960| 73 63 6f 72 65 20 2a 20 | 73 63 6f 72 65 3b 0a 58 |score * |score;.X|
|00004970| 09 09 09 09 7d 0a 58 09 | 09 09 7d 0a 58 09 76 61 |....}.X.|..}.X.va|
|00004980| 72 20 2d 3d 20 6d 65 61 | 6e 20 2a 20 6d 65 61 6e |r -= mea|n * mean|
|00004990| 3b 0a 58 0a 58 09 73 63 | 6f 72 65 32 5f 6d 65 61 |;.X.X.sc|ore2_mea|
|000049a0| 6e 20 3d 20 6d 65 61 6e | 3b 0a 58 09 73 63 6f 72 |n = mean|;.X.scor|
|000049b0| 65 32 5f 76 61 72 20 3d | 20 76 61 72 3b 0a 58 09 |e2_var =| var;.X.|
|000049c0| 73 63 6f 72 65 32 5f 73 | 64 20 3d 20 73 71 72 74 |score2_s|d = sqrt|
|000049d0| 28 73 63 6f 72 65 32 5f | 76 61 72 29 3b 0a 58 09 |(score2_|var);.X.|
|000049e0| 73 63 6f 72 65 32 5f 73 | 63 61 6c 65 20 3d 20 73 |score2_s|cale = s|
|000049f0| 71 72 74 28 32 20 2a 20 | 50 49 20 2a 20 73 63 6f |qrt(2 * |PI * sco|
|00004a00| 72 65 32 5f 76 61 72 29 | 3b 0a 58 0a 58 09 61 70 |re2_var)|;.X.X.ap|
|00004a10| 70 72 6f 78 5f 69 6e 69 | 74 28 29 3b 0a 58 7d 0a |prox_ini|t();.X}.|
|00004a20| 58 0a 58 0a 58 2f 2a 20 | 50 72 69 6e 74 20 74 68 |X.X.X/* |Print th|
|00004a30| 65 20 62 69 67 72 61 6d | 20 73 74 61 74 69 73 74 |e bigram| statist|
|00004a40| 69 63 73 2e 0a 58 20 2a | 2f 0a 58 70 72 69 6e 74 |ics..X *|/.Xprint|
|00004a50| 5f 32 73 74 61 74 73 28 | 6f 75 74 29 0a 58 46 49 |_2stats(|out).XFI|
|00004a60| 4c 45 09 2a 6f 75 74 3b | 0a 58 7b 0a 58 09 66 6c |LE.*out;|.X{.X.fl|
|00004a70| 6f 61 74 09 73 6c 6c 6f | 67 5f 6d 65 61 6e 3b 0a |oat.sllo|g_mean;.|
|00004a80| 58 09 66 6c 6f 61 74 09 | 73 6c 6c 6f 67 5f 76 61 |X.float.|sllog_va|
|00004a90| 72 3b 0a 58 09 66 6c 6f | 61 74 09 6c 65 76 5f 6d |r;.X.flo|at.lev_m|
|00004aa0| 65 61 6e 2c 20 6c 65 76 | 5f 76 61 72 3b 0a 58 09 |ean, lev|_var;.X.|
|00004ab0| 66 6c 6f 61 74 09 72 65 | 76 5f 6d 65 61 6e 2c 20 |float.re|v_mean, |
|00004ac0| 72 65 76 5f 76 61 72 3b | 0a 58 0a 58 09 66 70 72 |rev_var;|.X.X.fpr|
|00004ad0| 69 6e 74 66 28 6f 75 74 | 2c 20 22 5c 74 5c 74 42 |intf(out|, "\t\tB|
|00004ae0| 69 67 72 61 6d 20 53 74 | 61 74 69 73 74 69 63 73 |igram St|atistics|
|00004af0| 5c 6e 22 29 3b 0a 58 09 | 66 70 72 69 6e 74 66 28 |\n");.X.|fprintf(|
|00004b00| 6f 75 74 2c 20 22 53 63 | 6f 72 65 32 5f 6d 65 61 |out, "Sc|ore2_mea|
|00004b10| 6e 20 69 73 20 25 66 22 | 2c 20 73 63 6f 72 65 32 |n is %f"|, score2|
|00004b20| 5f 6d 65 61 6e 29 3b 0a | 58 09 66 70 72 69 6e 74 |_mean);.|X.fprint|
|00004b30| 66 28 6f 75 74 2c 20 22 | 2c 20 73 63 6f 72 65 32 |f(out, "|, score2|
|00004b40| 5f 76 61 72 20 69 73 20 | 25 66 22 2c 20 73 63 6f |_var is |%f", sco|
|00004b50| 72 65 32 5f 76 61 72 29 | 3b 0a 58 09 66 70 72 69 |re2_var)|;.X.fpri|
|00004b60| 6e 74 66 28 6f 75 74 2c | 20 22 2c 20 73 63 6f 72 |ntf(out,| ", scor|
|00004b70| 65 32 5f 73 64 20 69 73 | 20 25 66 22 2c 20 73 63 |e2_sd is| %f", sc|
|00004b80| 6f 72 65 32 5f 73 64 29 | 3b 0a 58 09 66 70 72 69 |ore2_sd)|;.X.fpri|
|00004b90| 6e 74 66 28 6f 75 74 2c | 20 22 5c 6e 6e 62 69 63 |ntf(out,| "\nnbic|
|00004ba0| 68 61 72 73 20 69 73 20 | 25 64 22 2c 20 6e 62 69 |hars is |%d", nbi|
|00004bb0| 63 68 61 72 73 29 3b 0a | 58 09 66 70 72 69 6e 74 |chars);.|X.fprint|
|00004bc0| 66 28 6f 75 74 2c 20 22 | 5c 6e 22 29 3b 0a 58 0a |f(out, "|\n");.X.|
|00004bd0| 58 09 73 6c 6c 6f 67 5f | 6d 65 61 6e 20 3d 20 76 |X.sllog_|mean = v|
|00004be0| 65 63 5f 6d 65 61 6e 28 | 73 6c 62 69 70 72 6f 62 |ec_mean(|slbiprob|
|00004bf0| 2c 20 73 6c 6c 6f 67 70 | 72 6f 62 2c 20 6e 62 69 |, sllogp|rob, nbi|
|00004c00| 63 68 61 72 73 29 3b 0a | 58 09 73 6c 6c 6f 67 5f |chars);.|X.sllog_|
|00004c10| 76 61 72 20 3d 20 76 65 | 63 5f 76 61 72 69 61 6e |var = ve|c_varian|
|00004c20| 63 65 28 73 6c 62 69 70 | 72 6f 62 2c 20 73 6c 6c |ce(slbip|rob, sll|
|00004c30| 6f 67 70 72 6f 62 2c 20 | 6e 62 69 63 68 61 72 73 |ogprob, |nbichars|
|00004c40| 29 3b 0a 58 09 66 70 72 | 69 6e 74 66 28 6f 75 74 |);.X.fpr|intf(out|
|00004c50| 2c 20 22 73 6c 6c 6f 67 | 5f 6d 65 61 6e 20 69 73 |, "sllog|_mean is|
|00004c60| 20 25 66 22 2c 20 73 6c | 6c 6f 67 5f 6d 65 61 6e | %f", sl|log_mean|
|00004c70| 29 3b 0a 58 09 66 70 72 | 69 6e 74 66 28 6f 75 74 |);.X.fpr|intf(out|
|00004c80| 2c 20 22 2c 20 73 6c 6c | 6f 67 5f 76 61 72 20 69 |, ", sll|og_var i|
|00004c90| 73 20 25 66 22 2c 20 73 | 6c 6c 6f 67 5f 76 61 72 |s %f", s|llog_var|
|00004ca0| 29 3b 0a 58 09 66 70 72 | 69 6e 74 66 28 6f 75 74 |);.X.fpr|intf(out|
|00004cb0| 2c 20 22 5c 6e 22 29 3b | 0a 58 7d 0a 58 0a 58 0a |, "\n");|.X}.X.X.|
|00004cc0| 58 0a 58 2f 2a 20 50 72 | 69 6e 74 20 74 68 65 20 |X.X/* Pr|int the |
|00004cd0| 66 69 72 73 74 20 6f 72 | 64 65 72 20 6c 6f 67 20 |first or|der log |
|00004ce0| 73 74 61 74 69 73 74 69 | 63 73 20 6f 6e 20 61 20 |statisti|cs on a |
|00004cf0| 73 74 72 65 61 6d 2e 0a | 58 20 2a 2f 0a 58 70 72 |stream..|X */.Xpr|
|00004d00| 69 6e 74 5f 31 73 74 61 | 74 73 28 6f 75 74 29 0a |int_1sta|ts(out).|
|00004d10| 58 46 49 4c 45 09 2a 6f | 75 74 3b 0a 58 7b 0a 58 |XFILE.*o|ut;.X{.X|
|00004d20| 0a 58 09 66 70 72 69 6e | 74 66 28 6f 75 74 2c 20 |.X.fprin|tf(out, |
|00004d30| 22 53 69 6e 67 6c 65 20 | 6c 65 74 74 65 72 20 66 |"Single |letter f|
|00004d40| 72 65 71 75 65 6e 63 69 | 65 73 5c 6e 22 29 3b 0a |requenci|es\n");.|
|00004d50| 58 09 66 70 72 69 6e 74 | 66 28 6f 75 74 2c 20 22 |X.fprint|f(out, "|
|00004d60| 5c 6e 45 78 70 65 63 74 | 65 64 20 76 61 6c 75 65 |\nExpect|ed value|
|00004d70| 20 6f 66 20 70 72 6f 62 | 20 69 73 20 25 66 2e 20 | of prob| is %f. |
|00004d80| 20 56 61 72 69 61 6e 63 | 65 20 69 73 20 25 66 2e | Varianc|e is %f.|
|00004d90| 5c 6e 22 2c 0a 58 09 20 | 20 20 20 20 20 20 70 6d |\n",.X. |      pm|
|00004da0| 65 61 6e 2c 20 70 76 61 | 72 29 3b 0a 58 09 70 72 |ean, pva|r);.X.pr|
|00004db0| 69 6e 74 5f 73 74 61 74 | 5f 74 61 62 28 6f 75 74 |int_stat|_tab(out|
|00004dc0| 2c 20 70 72 6f 62 2c 20 | 4d 41 58 43 48 41 52 29 |, prob, |MAXCHAR)|
|00004dd0| 3b 0a 58 0a 58 09 66 70 | 72 69 6e 74 66 28 6f 75 |;.X.X.fp|rintf(ou|
|00004de0| 74 2c 20 22 5c 6e 45 78 | 70 65 63 74 65 64 20 76 |t, "\nEx|pected v|
|00004df0| 61 6c 75 65 20 6f 66 20 | 6c 6f 67 70 72 6f 62 20 |alue of |logprob |
|00004e00| 69 73 20 25 66 2e 20 20 | 56 61 72 69 61 6e 63 65 |is %f.  |Variance|
|00004e10| 20 69 73 20 25 66 2e 5c | 6e 22 2c 0a 58 09 20 20 | is %f.\|n",.X.  |
|00004e20| 20 20 20 20 20 6c 6f 67 | 6d 65 61 6e 2c 20 6c 6f |     log|mean, lo|
|00004e30| 67 76 61 72 29 3b 0a 58 | 09 66 70 72 69 6e 74 66 |gvar);.X|.fprintf|
|00004e40| 28 6f 75 74 2c 20 22 4c | 6f 67 20 6f 66 20 73 69 |(out, "L|og of si|
|00004e50| 6e 67 6c 65 20 6c 65 74 | 74 65 72 20 66 72 65 71 |ngle let|ter freq|
|00004e60| 75 65 6e 63 69 65 73 5c | 6e 22 29 3b 0a 58 09 70 |uencies\|n");.X.p|
|00004e70| 72 69 6e 74 5f 73 74 61 | 74 5f 74 61 62 28 6f 75 |rint_sta|t_tab(ou|
|00004e80| 74 2c 20 6c 6f 67 70 72 | 6f 62 2c 20 4d 41 58 43 |t, logpr|ob, MAXC|
|00004e90| 48 41 52 29 3b 0a 58 7d | 0a 58 0a 58 0a 58 2f 2a |HAR);.X}|.X.X.X/*|
|00004ea0| 20 44 75 6d 70 61 20 73 | 74 61 74 69 73 74 69 63 | Dumpa s|tatistic|
|00004eb0| 73 20 74 61 62 6c 65 20 | 6f 6e 20 74 6f 20 61 20 |s table |on to a |
|00004ec0| 73 74 72 65 61 6d 2e 0a | 58 20 2a 2f 0a 58 70 72 |stream..|X */.Xpr|
|00004ed0| 69 6e 74 5f 73 74 61 74 | 5f 74 61 62 28 6f 75 74 |int_stat|_tab(out|
|00004ee0| 2c 20 74 61 62 6c 65 2c | 20 6d 61 78 69 6e 64 65 |, table,| maxinde|
|00004ef0| 78 29 0a 58 46 49 4c 45 | 09 2a 6f 75 74 3b 0a 58 |x).XFILE|.*out;.X|
|00004f00| 66 6c 6f 61 74 09 74 61 | 62 6c 65 5b 5d 3b 0a 58 |float.ta|ble[];.X|
|00004f10| 69 6e 74 09 09 6d 61 78 | 69 6e 64 65 78 3b 0a 58 |int..max|index;.X|
|00004f20| 7b 0a 58 09 69 6e 74 09 | 09 69 3b 0a 58 0a 58 09 |{.X.int.|.i;.X.X.|
|00004f30| 66 6f 72 20 28 69 20 3d | 20 30 20 3b 20 69 20 3c |for (i =| 0 ; i <|
|00004f40| 3d 20 6d 61 78 69 6e 64 | 65 78 20 3b 20 69 2b 2b |= maxind|ex ; i++|
|00004f50| 29 20 7b 0a 58 09 09 69 | 66 20 28 69 20 25 20 38 |) {.X..i|f (i % 8|
|00004f60| 20 3d 3d 20 30 29 20 20 | 66 70 72 69 6e 74 66 28 | == 0)  |fprintf(|
|00004f70| 6f 75 74 2c 20 22 5c 6e | 22 29 3b 0a 58 09 09 66 |out, "\n|");.X..f|
|00004f80| 70 72 69 6e 74 66 28 6f | 75 74 2c 20 22 25 37 2e |printf(o|ut, "%7.|
|00004f90| 34 66 20 22 2c 20 74 61 | 62 6c 65 5b 69 5d 29 3b |4f ", ta|ble[i]);|
|00004fa0| 0a 58 09 09 7d 0a 58 09 | 66 70 72 69 6e 74 66 28 |.X..}.X.|fprintf(|
|00004fb0| 6f 75 74 2c 20 22 5c 6e | 22 29 3b 0a 58 7d 0a 58 |out, "\n|");.X}.X|
|00004fc0| 0a 58 0a 58 2f 2a 20 4c | 6f 61 64 20 74 68 65 20 |.X.X/* L|oad the |
|00004fd0| 66 69 72 73 74 20 6f 72 | 64 65 72 20 73 74 61 74 |first or|der stat|
|00004fe0| 69 73 74 69 63 73 20 66 | 72 6f 6d 20 74 68 65 20 |istics f|rom the |
|00004ff0| 67 69 76 65 6e 20 66 69 | 6c 65 20 6e 61 6d 65 2e |given fi|le name.|
|00005000| 0a 58 20 2a 2f 0a 58 6c | 6f 61 64 5f 31 73 74 61 |.X */.Xl|oad_1sta|
|00005010| 74 73 5f 66 72 6f 6d 28 | 73 74 61 74 66 6e 61 6d |ts_from(|statfnam|
|00005020| 65 29 0a 58 63 68 61 72 | 09 2a 73 74 61 74 66 6e |e).Xchar|.*statfn|
|00005030| 61 6d 65 3b 0a 58 7b 0a | 58 09 46 49 4c 45 09 2a |ame;.X{.|X.FILE.*|
|00005040| 69 6e 70 3b 0a 58 0a 58 | 09 69 66 20 28 28 69 6e |inp;.X.X|.if ((in|
|00005050| 70 20 3d 20 66 6f 70 65 | 6e 28 73 74 61 74 66 6e |p = fope|n(statfn|
|00005060| 61 6d 65 2c 20 22 72 22 | 29 29 20 3d 3d 20 4e 55 |ame, "r"|)) == NU|
|00005070| 4c 4c 29 20 7b 0a 58 09 | 09 70 72 69 6e 74 66 28 |LL) {.X.|.printf(|
|00005080| 22 5c 6e 43 61 6e 27 74 | 20 6f 70 65 6e 20 25 73 |"\nCan't| open %s|
|00005090| 20 74 6f 20 72 65 61 64 | 20 6c 65 74 74 65 72 20 | to read| letter |
|000050a0| 73 74 61 74 69 73 74 69 | 63 73 5c 6e 22 2c 20 73 |statisti|cs\n", s|
|000050b0| 74 61 74 66 6e 61 6d 65 | 29 3b 0a 58 09 09 65 78 |tatfname|);.X..ex|
|000050c0| 69 74 28 30 29 3b 0a 58 | 09 09 7d 0a 58 09 6c 6f |it(0);.X|..}.X.lo|
|000050d0| 61 64 5f 31 73 74 61 74 | 73 28 69 6e 70 29 3b 0a |ad_1stat|s(inp);.|
|000050e0| 58 09 66 63 6c 6f 73 65 | 28 69 6e 70 29 3b 0a 58 |X.fclose|(inp);.X|
|000050f0| 7d 0a 58 0a 45 4e 44 5f | 4f 46 5f 73 74 61 74 73 |}.X.END_|OF_stats|
|00005100| 2e 63 0a 69 66 20 74 65 | 73 74 20 31 38 37 39 36 |.c.if te|st 18796|
|00005110| 20 2d 6e 65 20 60 77 63 | 20 2d 63 20 3c 73 74 61 | -ne `wc| -c <sta|
|00005120| 74 73 2e 63 60 3b 20 74 | 68 65 6e 0a 20 20 20 20 |ts.c`; t|hen.    |
|00005130| 65 63 68 6f 20 73 68 61 | 72 3a 20 5c 22 73 74 61 |echo sha|r: \"sta|
|00005140| 74 73 2e 63 5c 22 20 75 | 6e 70 61 63 6b 65 64 20 |ts.c\" u|npacked |
|00005150| 77 69 74 68 20 77 72 6f | 6e 67 20 73 69 7a 65 21 |with wro|ng size!|
|00005160| 0a 66 69 0a 23 20 65 6e | 64 20 6f 66 20 6f 76 65 |.fi.# en|d of ove|
|00005170| 72 77 72 69 74 69 6e 67 | 20 63 68 65 63 6b 0a 66 |rwriting| check.f|
|00005180| 69 0a 69 66 20 74 65 73 | 74 20 2d 66 20 74 65 72 |i.if tes|t -f ter|
|00005190| 6d 69 6e 61 6c 2e 63 20 | 2d 61 20 22 24 7b 31 7d |minal.c |-a "${1}|
|000051a0| 22 20 21 3d 20 22 2d 63 | 22 20 3b 20 74 68 65 6e |" != "-c|" ; then|
|000051b0| 20 0a 20 20 65 63 68 6f | 20 73 68 61 72 3a 20 57 | .  echo| shar: W|
|000051c0| 69 6c 6c 20 6e 6f 74 20 | 6f 76 65 72 2d 77 72 69 |ill not |over-wri|
|000051d0| 74 65 20 65 78 69 73 74 | 69 6e 67 20 66 69 6c 65 |te exist|ing file|
|000051e0| 20 5c 22 74 65 72 6d 69 | 6e 61 6c 2e 63 5c 22 0a | \"termi|nal.c\".|
|000051f0| 65 6c 73 65 0a 65 63 68 | 6f 20 73 68 61 72 3a 20 |else.ech|o shar: |
|00005200| 45 78 74 72 61 63 74 69 | 6e 67 20 5c 22 74 65 72 |Extracti|ng \"ter|
|00005210| 6d 69 6e 61 6c 2e 63 5c | 22 20 5c 28 31 39 32 38 |minal.c\|" \(1928|
|00005220| 35 20 63 68 61 72 61 63 | 74 65 72 73 5c 29 0a 73 |5 charac|ters\).s|
|00005230| 65 64 20 22 73 2f 5e 58 | 2f 2f 22 20 3e 74 65 72 |ed "s/^X|//" >ter|
|00005240| 6d 69 6e 61 6c 2e 63 20 | 3c 3c 27 45 4e 44 5f 4f |minal.c |<<'END_O|
|00005250| 46 5f 74 65 72 6d 69 6e | 61 6c 2e 63 27 0a 58 2f |F_termin|al.c'.X/|
|00005260| 2a 20 54 65 72 6d 69 6e | 61 6c 20 64 65 70 65 6e |* Termin|al depen|
|00005270| 64 65 6e 74 20 72 6f 75 | 74 69 6e 65 73 20 28 77 |dent rou|tines (w|
|00005280| 65 6c 6c 2c 20 6d 6f 73 | 74 20 6f 66 20 74 68 65 |ell, mos|t of the|
|00005290| 6d 20 2e 2e 2e 29 0a 58 | 20 2a 0a 58 20 2a 20 41 |m ...).X| *.X * A|
|000052a0| 75 74 68 6f 72 3a 20 42 | 6f 62 20 42 61 6c 64 77 |uthor: B|ob Baldw|
|000052b0| 69 6e 2c 20 4f 63 74 6f | 62 65 72 20 31 39 38 36 |in, Octo|ber 1986|
|000052c0| 2e 0a 58 20 2a 2f 0a 58 | 0a 58 2f 2a 20 0a 58 20 |..X */.X|.X/* .X |
|000052d0| 2a 20 52 6f 75 74 69 6e | 65 73 20 69 6e 20 74 65 |* Routin|es in te|
|000052e0| 72 6d 69 6e 61 6c 20 61 | 62 73 74 72 61 63 74 69 |rminal a|bstracti|
|000052f0| 6f 6e 3a 0a 58 20 2a 0a | 58 20 2a 09 73 65 74 5f |on:.X *.|X *.set_|
|00005300| 74 65 72 6d 28 29 0a 58 | 20 2a 09 09 49 6e 69 74 |term().X| *..Init|
|00005310| 69 61 6c 69 7a 65 20 74 | 65 72 6d 69 6e 61 6c 2c |ialize t|erminal,|
|00005320| 20 63 6c 65 61 72 20 74 | 68 65 20 73 63 72 65 65 | clear t|he scree|
|00005330| 6e 2e 0a 58 20 2a 0a 58 | 20 2a 09 75 6e 73 65 74 |n..X *.X| *.unset|
|00005340| 5f 74 65 72 6d 28 29 0a | 58 20 2a 09 09 52 65 74 |_term().|X *..Ret|
|00005350| 75 72 6e 20 74 65 72 6d | 69 6e 61 6c 20 74 6f 20 |urn term|inal to |
|00005360| 73 74 61 74 65 20 62 65 | 66 6f 72 65 20 73 65 74 |state be|fore set|
|00005370| 5f 74 65 72 6d 28 29 2e | 0a 58 20 2a 0a 58 20 2a |_term().|.X *.X *|
|00005380| 09 63 68 61 72 32 73 79 | 6d 28 63 68 61 72 29 0a |.char2sy|m(char).|
|00005390| 58 20 2a 09 09 52 65 74 | 75 72 6e 20 74 68 65 20 |X *..Ret|urn the |
|000053a0| 73 79 6d 62 6f 6c 20 75 | 73 65 64 20 74 6f 20 64 |symbol u|sed to d|
|000053b0| 69 73 70 6c 61 79 20 74 | 68 65 20 67 69 76 65 6e |isplay t|he given|
|000053c0| 20 63 68 61 72 20 69 6e | 20 74 68 65 0a 58 20 2a | char in| the.X *|
|000053d0| 09 09 64 65 63 72 79 70 | 74 69 6f 6e 20 77 69 6e |..decryp|tion win|
|000053e0| 64 6f 77 2e 0a 58 20 2a | 0a 58 20 2a 09 70 75 74 |dow..X *|.X *.put|
|000053f0| 73 79 6d 28 73 79 6d 62 | 6f 6c 29 0a 58 20 2a 09 |sym(symb|ol).X *.|
|00005400| 09 44 69 73 70 6c 61 79 | 73 20 74 68 65 20 67 69 |.Display|s the gi|
|00005410| 76 65 6e 20 73 79 6d 62 | 6f 6c 20 6f 6e 20 74 68 |ven symb|ol on th|
|00005420| 65 20 74 65 72 6d 69 6e | 61 6c 2e 20 20 48 61 6e |e termin|al.  Han|
|00005430| 64 6c 65 73 0a 58 20 2a | 20 09 09 65 6e 74 65 72 |dles.X *| ..enter|
|00005440| 69 6e 67 20 61 6e 64 20 | 65 78 69 74 69 6e 67 20 |ing and |exiting |
|00005450| 67 72 61 70 68 69 63 73 | 20 6d 6f 64 65 2e 0a 58 |graphics| mode..X|
|00005460| 20 2a 0a 58 20 2a 09 67 | 65 74 63 6d 64 28 29 0a | *.X *.g|etcmd().|
|00005470| 58 20 2a 09 09 52 65 61 | 64 73 20 73 74 64 69 6e |X *..Rea|ds stdin|
|00005480| 20 66 6f 72 20 61 20 6b | 65 79 73 74 72 6f 6b 65 | for a k|eystroke|
|00005490| 20 61 6e 64 20 72 65 74 | 75 72 6e 73 0a 58 20 2a | and ret|urns.X *|
|000054a0| 09 09 61 20 63 6f 6d 6d | 61 6e 64 20 69 6e 74 65 |..a comm|and inte|
|000054b0| 67 65 72 2e 0a 58 20 2a | 0a 58 20 2a 09 62 65 65 |ger..X *|.X *.bee|
|000054c0| 70 28 29 0a 58 20 2a 09 | 09 43 61 75 73 65 20 74 |p().X *.|.Cause t|
|000054d0| 68 65 20 74 65 72 6d 69 | 6e 61 6c 20 74 6f 20 62 |he termi|nal to b|
|000054e0| 65 65 70 20 6f 72 20 66 | 6c 61 73 68 2e 0a 58 20 |eep or f|lash..X |
|000054f0| 2a 2f 0a 58 0a 58 2f 2a | 20 54 68 65 20 74 72 69 |*/.X.X/*| The tri|
|00005500| 63 6b 20 74 6f 20 64 65 | 76 69 63 65 20 69 6e 64 |ck to de|vice ind|
|00005510| 65 70 65 6e 64 65 6e 63 | 65 20 69 73 20 74 6f 20 |ependenc|e is to |
|00005520| 63 6c 65 61 72 6c 79 20 | 73 65 70 61 72 61 74 65 |clearly |separate|
|00005530| 0a 58 20 2a 20 69 6e 74 | 65 72 6e 61 6c 20 61 6e |.X * int|ernal an|
|00005540| 64 20 65 78 74 65 72 6e | 61 6c 20 72 65 70 72 65 |d extern|al repre|
|00005550| 73 65 6e 74 61 74 69 6f | 6e 73 2e 20 20 4f 6e 20 |sentatio|ns.  On |
|00005560| 74 68 65 20 69 6e 62 6f | 75 6e 64 0a 58 20 2a 20 |the inbo|und.X * |
|00005570| 73 69 64 65 2c 20 77 65 | 20 73 65 70 61 72 61 74 |side, we| separat|
|00005580| 65 20 4b 45 59 53 54 4f | 4b 45 53 20 77 68 69 63 |e KEYSTO|KES whic|
|00005590| 68 20 67 65 6e 65 72 61 | 74 65 20 61 20 73 65 71 |h genera|te a seq|
|000055a0| 75 65 6e 63 65 0a 58 20 | 2a 20 6f 66 20 61 73 63 |uence.X |* of asc|
|000055b0| 69 69 20 63 68 61 72 61 | 63 74 65 72 73 20 28 6f |ii chara|cters (o|
|000055c0| 6e 65 20 70 65 72 20 6b | 65 79 73 74 72 6f 6b 65 |ne per k|eystroke|
|000055d0| 20 69 6e 20 74 68 65 20 | 73 69 6d 70 6c 65 20 63 | in the |simple c|
|000055e0| 61 73 65 29 2c 0a 58 20 | 2a 20 66 72 6f 6d 20 43 |ase),.X |* from C|
|000055f0| 4f 4d 4d 41 4e 44 2d 4b | 45 59 53 20 73 75 63 68 |OMMAND-K|EYS such|
|00005600| 20 61 73 20 6d 6f 76 65 | 2d 63 75 72 73 6f 72 2d | as move|-cursor-|
|00005610| 75 70 2e 20 20 4f 6e 20 | 74 68 65 20 6f 75 74 62 |up.  On |the outb|
|00005620| 6f 75 6e 64 0a 58 20 2a | 20 73 69 64 65 20 77 65 |ound.X *| side we|
|00005630| 20 73 65 70 61 72 61 74 | 65 20 53 59 4d 42 4f 4c | separat|e SYMBOL|
|00005640| 53 20 66 72 6f 6d 20 74 | 68 65 20 73 65 71 75 65 |S from t|he seque|
|00005650| 6e 63 65 20 6f 66 20 61 | 73 63 69 69 20 63 68 61 |nce of a|scii cha|
|00005660| 72 61 63 74 65 72 73 2c | 0a 58 20 2a 20 63 61 6c |racters,|.X * cal|
|00005670| 6c 65 64 20 47 52 41 50 | 48 49 43 53 2c 20 75 73 |led GRAP|HICS, us|
|00005680| 65 64 20 74 6f 20 64 69 | 73 70 6c 61 79 20 74 68 |ed to di|splay th|
|00005690| 65 20 73 79 6d 62 6f 6c | 2e 20 20 45 61 63 68 20 |e symbol|.  Each |
|000056a0| 73 79 6d 62 6f 6c 0a 58 | 20 2a 20 68 61 73 20 61 |symbol.X| * has a|
|000056b0| 20 73 69 6e 67 6c 65 20 | 75 73 65 20 62 79 20 74 | single |use by t|
|000056c0| 68 65 20 63 6f 64 65 2c | 20 74 68 6f 75 67 68 20 |he code,| though |
|000056d0| 74 77 6f 20 73 79 6d 62 | 6f 6c 73 20 6d 69 67 68 |two symb|ols migh|
|000056e0| 74 0a 58 20 2a 20 61 70 | 70 65 61 72 20 74 68 65 |t.X * ap|pear the|
|000056f0| 20 73 61 6d 65 20 6f 6e | 20 74 68 65 20 75 73 65 | same on| the use|
|00005700| 72 27 73 20 74 65 72 6d | 69 6e 61 6c 2e 0a 58 20 |r's term|inal..X |
|00005710| 2a 2f 0a 58 0a 58 2f 2a | 20 53 79 6d 62 6f 6c 73 |*/.X.X/*| Symbols|
|00005720| 20 61 72 65 20 72 65 70 | 72 65 73 65 6e 74 65 64 | are rep|resented|
|00005730| 20 62 79 20 69 6e 74 65 | 67 65 72 73 2e 20 20 49 | by inte|gers.  I|
|00005740| 66 20 74 68 65 20 69 6e | 74 65 67 65 72 20 69 73 |f the in|teger is|
|00005750| 0a 58 20 2a 20 67 72 65 | 61 74 65 72 20 74 68 61 |.X * gre|ater tha|
|00005760| 6e 20 32 35 36 2c 20 74 | 68 65 6e 20 69 74 20 64 |n 256, t|hen it d|
|00005770| 65 6e 6f 74 65 73 20 6f | 6e 65 20 6f 66 20 74 68 |enotes o|ne of th|
|00005780| 65 20 73 79 6d 62 6f 6c | 73 20 64 65 66 69 6e 65 |e symbol|s define|
|00005790| 64 0a 58 20 2a 20 69 6e | 20 74 65 72 6d 69 6e 61 |d.X * in| termina|
|000057a0| 6c 2e 68 2e 0a 58 20 2a | 0a 58 20 2a 20 43 6f 6d |l.h..X *|.X * Com|
|000057b0| 6d 61 6e 64 73 20 61 72 | 65 20 72 65 70 72 65 73 |mands ar|e repres|
|000057c0| 65 6e 74 65 64 20 62 79 | 20 74 77 6f 20 62 79 74 |ented by| two byt|
|000057d0| 65 20 69 6e 74 65 67 65 | 72 73 2e 20 20 54 68 65 |e intege|rs.  The|
|000057e0| 20 68 69 67 68 20 62 79 | 74 65 0a 58 20 2a 20 64 | high by|te.X * d|
|000057f0| 65 73 63 72 69 62 65 73 | 20 74 68 65 20 63 6f 6d |escribes| the com|
|00005800| 6d 61 6e 64 20 28 73 65 | 65 20 74 65 72 6d 69 6e |mand (se|e termin|
|00005810| 61 6c 2e 68 29 2c 20 74 | 68 65 20 6c 6f 77 20 62 |al.h), t|he low b|
|00005820| 79 74 65 20 69 73 20 74 | 68 65 0a 58 20 2a 20 61 |yte is t|he.X * a|
|00005830| 72 67 75 6d 65 6e 74 20 | 74 6f 20 74 68 65 20 63 |rgument |to the c|
|00005840| 6f 6d 6d 61 6e 64 2e 20 | 20 46 6f 72 20 65 78 61 |ommand. | For exa|
|00005850| 6d 70 6c 65 2c 20 74 68 | 65 20 69 6e 73 65 72 74 |mple, th|e insert|
|00005860| 20 63 68 61 72 20 63 6f | 6d 6d 61 6e 64 0a 58 20 | char co|mmand.X |
|00005870| 2a 20 73 65 74 73 20 74 | 68 65 20 6c 6f 77 20 62 |* sets t|he low b|
|00005880| 79 74 65 20 74 6f 20 74 | 68 65 20 63 68 61 72 61 |yte to t|he chara|
|00005890| 63 74 65 72 20 74 6f 20 | 62 65 20 69 6e 73 65 72 |cter to |be inser|
|000058a0| 74 65 64 2e 0a 58 20 2a | 2f 0a 58 0a 58 2f 2a 20 |ted..X *|/.X.X/* |
|000058b0| 49 4e 54 45 52 4e 41 4c | 53 3a 20 73 79 6d 62 6f |INTERNAL|S: symbo|
|000058c0| 6c 73 20 61 6e 64 20 67 | 72 61 70 68 69 63 73 0a |ls and g|raphics.|
|000058d0| 58 20 2a 0a 58 20 2a 20 | 54 68 65 20 74 65 72 6d |X *.X * |The term|
|000058e0| 69 6e 61 6c 20 69 73 20 | 61 73 73 75 6d 65 64 20 |inal is |assumed |
|000058f0| 74 6f 20 62 65 20 69 6e | 20 6f 6e 65 20 6f 66 20 |to be in| one of |
|00005900| 74 68 72 65 65 20 6d 6f | 64 65 73 3a 20 6e 6f 72 |three mo|des: nor|
|00005910| 6d 61 6c 2c 20 67 72 61 | 70 68 69 63 73 2c 0a 58 |mal, gra|phics,.X|
|00005920| 20 2a 20 6f 72 20 73 74 | 61 6e 64 6f 75 74 20 28 | * or st|andout (|
|00005930| 69 6e 76 65 72 73 65 20 | 76 69 64 65 6f 29 2e 20 |inverse |video). |
|00005940| 20 54 68 69 73 20 61 62 | 73 74 72 61 63 74 69 6f | This ab|stractio|
|00005950| 6e 20 74 61 6b 65 73 20 | 63 61 72 65 20 6f 66 20 |n takes |care of |
|00005960| 61 6c 6c 20 74 68 65 0a | 58 20 2a 20 73 77 69 74 |all the.|X * swit|
|00005970| 63 68 69 6e 67 20 62 65 | 74 77 65 65 6e 20 6d 6f |ching be|tween mo|
|00005980| 64 65 73 20 61 6e 64 20 | 74 72 69 65 73 20 74 6f |des and |tries to|
|00005990| 20 61 76 6f 69 64 20 73 | 65 6e 64 69 6e 67 20 72 | avoid s|ending r|
|000059a0| 65 64 75 6e 64 61 6e 74 | 20 65 73 63 61 70 65 0a |edundant| escape.|
|000059b0| 58 20 2a 20 73 65 71 75 | 65 6e 63 65 73 20 74 6f |X * sequ|ences to|
|000059c0| 20 74 68 65 20 74 65 72 | 6d 69 6e 61 6c 2e 0a 58 | the ter|minal..X|
|000059d0| 20 2a 20 0a 58 20 2a 20 | 50 61 72 74 20 6f 66 20 | * .X * |Part of |
|000059e0| 74 68 65 20 74 65 72 6d | 69 6e 61 6c 20 69 6e 69 |the term|inal ini|
|000059f0| 74 69 61 6c 69 7a 61 74 | 69 6f 6e 20 69 73 20 74 |tializat|ion is t|
|00005a00| 6f 20 62 75 69 6c 64 20 | 61 20 74 61 62 6c 65 20 |o build |a table |
|00005a10| 74 68 61 74 20 64 65 73 | 63 72 69 62 65 73 0a 58 |that des|cribes.X|
|00005a20| 20 2a 20 68 6f 77 20 74 | 6f 20 64 69 73 70 6c 61 | * how t|o displa|
|00005a30| 79 20 65 61 63 68 20 73 | 79 6d 62 6f 6c 2e 20 20 |y each s|ymbol.  |
|00005a40| 41 73 63 69 69 20 63 68 | 61 72 61 63 74 65 72 73 |Ascii ch|aracters|
|00005a50| 20 28 69 6e 63 6c 75 64 | 69 6e 67 20 63 6f 6e 74 | (includ|ing cont|
|00005a60| 72 6f 6c 20 63 68 61 72 | 73 29 0a 58 20 2a 20 70 |rol char|s).X * p|
|00005a70| 61 73 73 20 74 68 72 6f | 75 67 68 20 75 6e 63 68 |ass thro|ugh unch|
|00005a80| 61 6e 67 65 64 2c 20 62 | 75 74 20 73 79 6d 62 6f |anged, b|ut symbo|
|00005a90| 6c 73 20 61 72 65 20 6c | 6f 6f 6b 65 64 20 75 70 |ls are l|ooked up|
|00005aa0| 20 69 6e 20 67 72 61 70 | 68 74 61 62 2e 0a 58 20 | in grap|htab..X |
|00005ab0| 2a 20 45 61 63 68 20 65 | 6e 74 72 79 20 69 6e 20 |* Each e|ntry in |
|00005ac0| 74 68 65 20 67 72 61 70 | 68 69 63 73 20 74 61 62 |the grap|hics tab|
|00005ad0| 6c 65 20 73 70 65 63 69 | 66 69 65 73 20 61 20 74 |le speci|fies a t|
|00005ae0| 65 72 6d 69 6e 61 6c 20 | 6d 6f 64 65 20 61 6e 64 |erminal |mode and|
|00005af0| 20 61 0a 58 20 2a 20 73 | 74 72 69 6e 67 20 74 6f | a.X * s|tring to|
|00005b00| 20 70 72 69 6e 74 20 74 | 68 61 74 20 77 69 6c 6c | print t|hat will|
|00005b10| 20 64 69 73 70 6c 61 79 | 20 74 68 61 74 20 73 79 | display| that sy|
|00005b20| 6d 62 6f 6c 20 69 6e 20 | 74 68 65 20 67 69 76 65 |mbol in |the give|
|00005b30| 6e 20 6d 6f 64 65 2e 0a | 58 20 2a 20 49 66 20 6e |n mode..|X * If n|
|00005b40| 65 63 65 73 73 61 72 79 | 20 61 20 6d 6f 64 65 20 |ecessary| a mode |
|00005b50| 73 77 69 74 63 68 20 73 | 65 71 75 65 6e 63 65 20 |switch s|equence |
|00005b60| 77 69 6c 6c 20 62 65 20 | 73 65 6e 74 20 62 65 66 |will be |sent bef|
|00005b70| 6f 72 65 20 70 72 69 6e | 74 69 6e 67 20 74 68 65 |ore prin|ting the|
|00005b80| 0a 58 20 2a 20 73 74 72 | 69 6e 67 2e 0a 58 20 2a |.X * str|ing..X *|
|00005b90| 0a 58 20 2a 20 54 68 65 | 20 67 72 61 70 68 69 63 |.X * The| graphic|
|00005ba0| 73 20 63 6f 72 72 65 73 | 70 6f 6e 64 69 6e 67 20 |s corres|ponding |
|00005bb0| 74 6f 20 73 79 6d 62 6f | 6c 73 20 63 61 6e 20 62 |to symbo|ls can b|
|00005bc0| 65 20 73 65 74 20 62 79 | 20 74 68 65 20 73 68 65 |e set by| the she|
|00005bd0| 6c 6c 20 76 61 72 69 61 | 62 6c 65 0a 58 20 2a 20 |ll varia|ble.X * |
|00005be0| 27 47 52 41 48 50 49 43 | 53 27 20 75 73 69 6e 67 |'GRAHPIC|S' using|
|00005bf0| 20 61 20 66 6f 72 6d 61 | 74 20 73 69 6d 69 6c 61 | a forma|t simila|
|00005c00| 72 20 74 6f 20 74 68 65 | 20 74 65 72 6d 63 61 70 |r to the| termcap|
|00005c10| 20 66 6f 72 6d 61 74 2e | 20 20 54 68 65 20 74 77 | format.|  The tw|
|00005c20| 6f 0a 58 20 2a 20 63 68 | 61 72 61 63 74 65 72 20 |o.X * ch|aracter |
|00005c30| 6e 61 6d 65 73 20 6f 66 | 20 74 68 65 20 73 79 6d |names of| the sym|
|00005c40| 62 6f 6c 73 20 61 72 65 | 20 64 65 66 69 6e 65 64 |bols are| defined|
|00005c50| 20 69 6e 20 74 68 65 20 | 74 61 62 6c 65 20 73 79 | in the |table sy|
|00005c60| 6d 6e 61 6d 65 73 5b 5d | 2c 0a 58 20 2a 20 77 68 |mnames[]|,.X * wh|
|00005c70| 69 63 68 20 63 61 6e 20 | 62 65 20 69 6e 64 65 78 |ich can |be index|
|00005c80| 65 64 20 62 79 20 74 68 | 65 20 73 79 6d 62 6f 6c |ed by th|e symbol|
|00005c90| 20 63 6f 64 65 73 20 28 | 73 65 65 20 74 65 72 6d | codes (|see term|
|00005ca0| 69 6e 61 6c 2e 68 29 2e | 0a 58 20 2a 20 54 68 65 |inal.h).|.X * The|
|00005cb0| 20 6d 61 70 70 69 6e 67 | 20 66 72 6f 6d 20 73 79 | mapping| from sy|
|00005cc0| 6d 62 6f 6c 73 20 74 6f | 20 67 72 61 70 68 69 63 |mbols to| graphic|
|00005cd0| 73 20 69 73 20 66 69 72 | 73 74 20 73 65 74 20 74 |s is fir|st set t|
|00005ce0| 6f 20 61 20 64 65 66 61 | 75 6c 74 20 76 61 6c 75 |o a defa|ult valu|
|00005cf0| 65 2c 0a 58 20 2a 20 61 | 6e 64 20 74 68 65 6e 20 |e,.X * a|nd then |
|00005d00| 75 70 64 61 74 65 64 20 | 62 79 20 74 68 65 20 76 |updated |by the v|
|00005d10| 61 6c 75 65 73 20 66 6f | 75 6e 64 20 69 6e 20 74 |alues fo|und in t|
|00005d20| 68 65 20 73 68 65 6c 6c | 20 76 61 72 69 61 62 6c |he shell| variabl|
|00005d30| 65 20 28 69 66 20 61 6e | 79 29 2e 0a 58 20 2a 0a |e (if an|y)..X *.|
|00005d40| 58 20 2a 20 54 68 65 20 | 47 52 41 50 48 49 43 53 |X * The |GRAPHICS|
|00005d50| 20 73 74 72 69 6e 67 20 | 63 6f 6e 73 69 73 74 73 | string |consists|
|00005d60| 20 6f 66 20 61 20 6e 75 | 6d 62 65 72 20 6f 66 20 | of a nu|mber of |
|00005d70| 65 6e 74 72 69 65 73 20 | 73 65 70 61 72 61 74 65 |entries |separate|
|00005d80| 64 20 62 79 20 63 6f 6c | 6f 6e 0a 58 20 2a 20 63 |d by col|on.X * c|
|00005d90| 68 61 72 61 63 74 65 72 | 73 2e 20 20 45 61 63 68 |haracter|s.  Each|
|00005da0| 20 65 6e 74 72 79 20 68 | 61 76 65 20 74 68 65 20 | entry h|ave the |
|00005db0| 66 6f 72 6d 20 78 78 3d | 5c 4d 63 63 63 2c 20 77 |form xx=|\Mccc, w|
|00005dc0| 68 65 72 65 20 78 78 20 | 69 73 20 74 68 65 20 6e |here xx |is the n|
|00005dd0| 61 6d 65 0a 58 20 2a 20 | 6f 66 20 61 20 73 79 6d |ame.X * |of a sym|
|00005de0| 62 6f 6c 20 28 73 65 65 | 20 73 79 6d 6e 61 6d 65 |bol (see| symname|
|00005df0| 73 5b 5d 20 62 65 6c 6f | 77 29 2c 20 5c 4d 20 69 |s[] belo|w), \M i|
|00005e00| 6e 64 69 63 61 74 65 73 | 20 74 68 65 20 6d 6f 64 |ndicates| the mod|
|00005e10| 65 20 66 6f 72 20 64 69 | 73 70 6c 61 79 69 6e 67 |e for di|splaying|
|00005e20| 0a 58 20 2a 20 74 68 69 | 73 20 73 79 6d 62 6f 6c |.X * thi|s symbol|
|00005e30| 20 28 69 74 20 6d 75 73 | 74 20 62 65 20 6f 6e 65 | (it mus|t be one|
|00005e40| 20 6f 66 20 5c 4e 20 28 | 66 6f 72 20 6e 6f 72 6d | of \N (|for norm|
|00005e50| 61 6c 29 2c 20 5c 47 20 | 28 66 6f 72 20 67 72 61 |al), \G |(for gra|
|00005e60| 70 68 69 63 73 29 2c 20 | 6f 72 20 5c 53 0a 58 20 |phics), |or \S.X |
|00005e70| 2a 20 28 66 6f 72 20 73 | 74 61 6e 64 6f 75 74 20 |* (for s|tandout |
|00005e80| 6f 72 20 69 6e 76 65 72 | 73 65 20 76 69 64 65 6f |or inver|se video|
|00005e90| 29 29 2c 20 61 6e 64 20 | 63 63 63 20 69 73 20 61 |)), and |ccc is a|
|00005ea0| 20 73 65 71 75 65 6e 63 | 65 20 6f 66 20 63 68 61 | sequenc|e of cha|
|00005eb0| 72 61 63 74 65 72 73 0a | 58 20 2a 20 74 6f 20 73 |racters.|X * to s|
|00005ec0| 65 6e 64 2e 20 20 54 6f | 20 69 6e 63 6c 75 64 65 |end.  To| include|
|00005ed0| 20 61 20 63 6f 6c 6f 6e | 20 63 68 61 72 61 63 74 | a colon| charact|
|00005ee0| 65 72 20 69 6e 20 74 68 | 65 20 63 63 63 20 70 6f |er in th|e ccc po|
|00005ef0| 72 74 69 6f 6e 2c 20 70 | 72 65 63 65 65 64 20 69 |rtion, p|receed i|
|00005f00| 74 0a 58 20 2a 20 62 79 | 20 61 20 62 61 63 6b 73 |t.X * by| a backs|
|00005f10| 6c 61 73 68 20 28 69 2e | 65 2e 2c 20 27 5c 27 29 |lash (i.|e., '\')|
|00005f20| 2e 20 20 49 66 20 74 68 | 65 20 5c 4d 20 69 73 20 |.  If th|e \M is |
|00005f30| 6f 6d 6d 69 74 74 65 64 | 2c 20 69 74 20 64 65 66 |ommitted|, it def|
|00005f40| 61 75 6c 74 73 20 74 6f | 20 6e 6f 72 6d 61 6c 2e |aults to| normal.|
|00005f50| 0a 58 20 2a 2f 0a 58 0a | 58 2f 2a 20 49 4e 54 45 |.X */.X.|X/* INTE|
|00005f60| 52 4e 41 4c 53 3a 20 6b | 65 79 73 74 72 6f 6b 65 |RNALS: k|eystroke|
|00005f70| 73 20 61 6e 64 20 63 6f | 6d 6d 61 6e 64 73 0a 58 |s and co|mmands.X|
|00005f80| 20 2a 0a 58 20 2a 20 54 | 68 65 20 74 61 62 6c 65 | *.X * T|he table|
|00005f90| 20 6b 65 79 63 6d 64 74 | 61 62 20 69 73 20 75 73 | keycmdt|ab is us|
|00005fa0| 65 64 20 74 6f 20 63 6f | 6e 76 65 72 74 20 61 20 |ed to co|nvert a |
|00005fb0| 73 65 71 75 65 6e 63 65 | 20 6f 66 20 61 73 63 69 |sequence| of asci|
|00005fc0| 69 20 63 68 61 72 61 63 | 74 65 72 73 0a 58 20 2a |i charac|ters.X *|
|00005fd0| 20 69 6e 74 6f 20 61 20 | 63 6f 6d 6d 61 6e 64 20 | into a |command |
|00005fe0| 69 6e 74 65 67 65 72 2e | 20 20 54 68 65 20 74 61 |integer.|  The ta|
|00005ff0| 62 6c 65 20 73 70 65 63 | 69 66 69 65 73 20 74 68 |ble spec|ifies th|
|00006000| 65 20 65 73 63 61 70 65 | 20 73 65 71 75 65 6e 63 |e escape| sequenc|
|00006010| 65 73 20 74 68 61 74 0a | 58 20 2a 20 6d 69 67 68 |es that.|X * migh|
|00006020| 74 20 62 65 20 67 65 6e | 65 72 61 74 65 64 20 62 |t be gen|erated b|
|00006030| 79 20 74 68 65 20 64 69 | 66 66 65 72 65 6e 74 20 |y the di|fferent |
|00006040| 63 6f 6d 6d 61 6e 64 20 | 6b 65 79 73 2e 20 20 54 |command |keys.  T|
|00006050| 68 65 20 67 65 74 63 6d | 64 20 72 6f 75 74 69 6e |he getcm|d routin|
|00006060| 65 0a 58 20 2a 20 77 69 | 6c 6c 20 72 65 61 64 20 |e.X * wi|ll read |
|00006070| 63 68 61 72 61 63 74 65 | 72 73 20 66 72 6f 6d 20 |characte|rs from |
|00006080| 73 74 64 69 6e 20 75 6e | 74 69 6c 20 69 74 20 68 |stdin un|til it h|
|00006090| 61 73 20 75 6e 69 71 75 | 65 6c 79 20 69 64 65 6e |as uniqu|ely iden|
|000060a0| 74 69 66 69 65 64 20 61 | 0a 58 20 2a 20 63 6f 6d |tified a|.X * com|
|000060b0| 6d 61 6e 64 2c 20 6f 72 | 20 64 65 63 69 64 65 64 |mand, or| decided|
|000060c0| 20 74 68 61 74 20 74 68 | 65 72 65 20 69 73 20 6e | that th|ere is n|
|000060d0| 6f 20 6d 61 74 63 68 2e | 20 20 49 66 20 6e 6f 20 |o match.|  If no |
|000060e0| 6d 61 74 63 68 20 69 73 | 20 66 6f 75 6e 64 2c 0a |match is| found,.|
|000060f0| 58 20 2a 20 74 68 65 20 | 74 65 72 6d 69 6e 61 6c |X * the |terminal|
|00006100| 20 69 73 20 62 65 65 70 | 65 64 20 61 6e 64 20 6d | is beep|ed and m|
|00006110| 61 74 63 68 69 6e 67 20 | 69 73 20 72 65 73 74 61 |atching |is resta|
|00006120| 72 74 65 64 20 77 69 74 | 68 20 74 68 65 20 6e 65 |rted wit|h the ne|
|00006130| 78 74 20 63 68 61 72 61 | 63 74 65 72 2e 0a 58 20 |xt chara|cter..X |
|00006140| 2a 20 42 79 20 64 65 66 | 61 75 6c 74 20 6b 65 79 |* By def|ault key|
|00006150| 73 74 72 6f 6b 65 73 20 | 74 75 72 6e 20 69 6e 74 |strokes |turn int|
|00006160| 6f 20 73 65 6c 66 2d 69 | 6e 73 65 72 74 20 63 6f |o self-i|nsert co|
|00006170| 6d 6d 61 6e 64 73 2e 0a | 58 20 2a 0a 58 20 2a 20 |mmands..|X *.X * |
|00006180| 49 6e 20 67 65 6e 65 72 | 61 6c 20 74 68 65 20 6c |In gener|al the l|
|00006190| 61 73 74 20 63 68 61 72 | 61 63 74 65 72 20 6f 66 |ast char|acter of|
|000061a0| 20 74 68 65 20 73 65 71 | 75 65 6e 63 65 20 69 73 | the seq|uence is|
|000061b0| 20 72 65 74 75 72 6e 65 | 64 20 61 73 20 74 68 65 | returne|d as the|
|000061c0| 20 61 72 67 2e 0a 58 20 | 2a 20 54 68 69 73 20 68 | arg..X |* This h|
|000061d0| 65 6c 70 73 20 77 69 6e | 64 6f 77 73 20 61 73 73 |elps win|dows ass|
|000061e0| 69 67 6e 20 64 69 66 66 | 65 72 65 6e 74 20 6d 65 |ign diff|erent me|
|000061f0| 61 6e 69 6e 67 73 20 74 | 6f 20 6b 65 79 73 74 6f |anings t|o keysto|
|00006200| 6b 65 73 2e 20 20 46 6f | 72 20 65 78 61 6d 70 6c |kes.  Fo|r exampl|
|00006210| 65 0a 58 20 2a 20 54 68 | 65 20 72 65 74 75 72 6e |e.X * Th|e return|
|00006220| 20 6b 65 79 20 63 61 6e | 20 62 65 20 65 69 74 68 | key can| be eith|
|00006230| 65 72 20 61 6e 20 69 6e | 73 65 72 74 2d 72 65 74 |er an in|sert-ret|
|00006240| 75 72 6e 20 63 6f 6d 6d | 61 6e 64 20 6f 72 20 61 |urn comm|and or a|
|00006250| 6e 20 65 78 65 63 75 74 | 65 2d 0a 58 20 2a 20 63 |n execut|e-.X * c|
|00006260| 6f 6d 6d 61 6e 64 2d 6c | 69 6e 65 20 63 6f 6d 6d |ommand-l|ine comm|
|00006270| 61 6e 64 2e 0a 58 20 2a | 2f 0a 58 0a 58 0a 58 23 |and..X *|/.X.X.X#|
|00006280| 69 6e 63 6c 75 64 65 09 | 3c 63 75 72 73 65 73 2e |include.|<curses.|
|00006290| 68 3e 0a 58 23 69 6e 63 | 6c 75 64 65 09 3c 73 67 |h>.X#inc|lude.<sg|
|000062a0| 74 74 79 2e 68 3e 0a 58 | 23 69 6e 63 6c 75 64 65 |tty.h>.X|#include|
|000062b0| 09 3c 73 74 72 69 6e 67 | 73 2e 68 3e 0a 58 23 69 |.<string|s.h>.X#i|
|000062c0| 6e 63 6c 75 64 65 09 22 | 77 69 6e 64 6f 77 2e 68 |nclude."|window.h|
|000062d0| 22 0a 58 23 69 6e 63 6c | 75 64 65 09 22 74 65 72 |".X#incl|ude."ter|
|000062e0| 6d 69 6e 61 6c 2e 68 22 | 0a 58 23 69 6e 63 6c 75 |minal.h"|.X#inclu|
|000062f0| 64 65 09 22 73 70 65 63 | 73 2e 68 22 0a 58 0a 58 |de."spec|s.h".X.X|
|00006300| 0a 58 2f 2a 20 52 6f 75 | 74 69 6e 65 73 20 66 72 |.X/* Rou|tines fr|
|00006310| 6f 6d 20 74 65 72 6d 63 | 61 70 20 6c 69 62 72 61 |om termc|ap libra|
|00006320| 72 79 2e 0a 58 20 2a 2f | 0a 58 65 78 74 65 72 6e |ry..X */|.Xextern|
|00006330| 09 63 68 61 72 20 2a 67 | 65 74 65 6e 76 28 29 2c |.char *g|etenv(),|
|00006340| 20 2a 74 67 65 74 73 74 | 72 28 29 2c 20 2a 74 67 | *tgetst|r(), *tg|
|00006350| 6f 74 6f 28 29 3b 0a 58 | 0a 58 0a 58 2f 2a 20 53 |oto();.X|.X.X/* S|
|00006360| 63 72 65 65 6e 20 63 6f | 6e 74 72 6f 6c 20 73 74 |creen co|ntrol st|
|00006370| 72 69 6e 67 73 20 66 72 | 6f 6d 20 74 65 72 6d 63 |rings fr|om termc|
|00006380| 61 70 20 65 6e 74 72 79 | 2e 0a 58 20 2a 2f 0a 58 |ap entry|..X */.X|
|00006390| 63 68 61 72 09 2a 74 65 | 72 6d 5f 69 73 3b 09 09 |char.*te|rm_is;..|
|000063a0| 2f 2a 20 54 65 72 6d 69 | 6e 61 6c 20 69 6e 69 74 |/* Termi|nal init|
|000063b0| 69 61 6c 69 7a 61 69 6f | 6e 20 73 74 72 69 6e 67 |ializaio|n string|
|000063c0| 2e 20 2a 2f 0a 58 63 68 | 61 72 09 2a 65 72 61 73 |. */.Xch|ar.*eras|
|000063d0| 65 5f 65 6f 6c 3b 09 09 | 2f 2a 20 45 72 61 73 65 |e_eol;..|/* Erase|
|000063e0| 20 74 6f 20 65 6e 64 20 | 6f 66 20 6c 69 6e 65 2e | to end |of line.|
|000063f0| 20 2a 2f 0a 58 63 68 61 | 72 09 2a 65 72 61 73 65 | */.Xcha|r.*erase|
+--------+-------------------------+-------------------------+--------+--------+