1*6ae1554aSColin Percival /*
2*6ae1554aSColin Percival * Copyright (c) 1989, 1993
3*6ae1554aSColin Percival * The Regents of the University of California. All rights reserved.
4*6ae1554aSColin Percival *
5*6ae1554aSColin Percival * This code is derived from software contributed to Berkeley by
6*6ae1554aSColin Percival * Rick Adams.
7*6ae1554aSColin Percival *
8*6ae1554aSColin Percival * Authors:
9*6ae1554aSColin Percival * Stan King, John Eldridge, based on algorithm suggested by
10*6ae1554aSColin Percival * Bob Morris
11*6ae1554aSColin Percival * 29-Sep-82
12*6ae1554aSColin Percival *
13*6ae1554aSColin Percival * Redistribution and use in source and binary forms, with or without
14*6ae1554aSColin Percival * modification, are permitted provided that the following conditions
15*6ae1554aSColin Percival * are met:
16*6ae1554aSColin Percival * 1. Redistributions of source code must retain the above copyright
17*6ae1554aSColin Percival * notice, this list of conditions and the following disclaimer.
18*6ae1554aSColin Percival * 2. Redistributions in binary form must reproduce the above copyright
19*6ae1554aSColin Percival * notice, this list of conditions and the following disclaimer in the
20*6ae1554aSColin Percival * documentation and/or other materials provided with the distribution.
21*6ae1554aSColin Percival * 3. Neither the name of the University nor the names of its contributors
22*6ae1554aSColin Percival * may be used to endorse or promote products derived from this software
23*6ae1554aSColin Percival * without specific prior written permission.
24*6ae1554aSColin Percival *
25*6ae1554aSColin Percival * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
26*6ae1554aSColin Percival * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
27*6ae1554aSColin Percival * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
28*6ae1554aSColin Percival * ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
29*6ae1554aSColin Percival * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
30*6ae1554aSColin Percival * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
31*6ae1554aSColin Percival * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
32*6ae1554aSColin Percival * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
33*6ae1554aSColin Percival * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
34*6ae1554aSColin Percival * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
35*6ae1554aSColin Percival * SUCH DAMAGE.
36*6ae1554aSColin Percival */
37*6ae1554aSColin Percival
38*6ae1554aSColin Percival #include <errno.h>
39*6ae1554aSColin Percival #include <math.h>
40*6ae1554aSColin Percival #include <stdio.h>
41*6ae1554aSColin Percival #include <stdlib.h>
42*6ae1554aSColin Percival #include <string.h>
43*6ae1554aSColin Percival #include <ctype.h>
44*6ae1554aSColin Percival #include <unistd.h>
45*6ae1554aSColin Percival
46*6ae1554aSColin Percival #define LINELENGTH 2048
47*6ae1554aSColin Percival #define ROTATE(ch, perm) \
48*6ae1554aSColin Percival isascii(ch) ? ( \
49*6ae1554aSColin Percival isupper(ch) ? ('A' + (ch - 'A' + perm) % 26) : \
50*6ae1554aSColin Percival islower(ch) ? ('a' + (ch - 'a' + perm) % 26) : ch) : ch
51*6ae1554aSColin Percival
52*6ae1554aSColin Percival /*
53*6ae1554aSColin Percival * letter frequencies (taken from some unix(tm) documentation)
54*6ae1554aSColin Percival * (unix is a trademark of Bell Laboratories)
55*6ae1554aSColin Percival */
56*6ae1554aSColin Percival static double stdf[26] = {
57*6ae1554aSColin Percival 7.97, 1.35, 3.61, 4.78, 12.37, 2.01, 1.46, 4.49, 6.39, 0.04,
58*6ae1554aSColin Percival 0.42, 3.81, 2.69, 5.92, 6.96, 2.91, 0.08, 6.63, 8.77, 9.68,
59*6ae1554aSColin Percival 2.62, 0.81, 1.88, 0.23, 2.07, 0.06,
60*6ae1554aSColin Percival };
61*6ae1554aSColin Percival
62*6ae1554aSColin Percival static void printit(char *);
63*6ae1554aSColin Percival
64*6ae1554aSColin Percival int
main(int argc,char ** argv)65*6ae1554aSColin Percival main(int argc, char **argv)
66*6ae1554aSColin Percival {
67*6ae1554aSColin Percival int ch, dot, i, nread, winnerdot = 0;
68*6ae1554aSColin Percival char *inbuf;
69*6ae1554aSColin Percival int obs[26], try, winner;
70*6ae1554aSColin Percival
71*6ae1554aSColin Percival if (argc > 1)
72*6ae1554aSColin Percival printit(argv[1]);
73*6ae1554aSColin Percival
74*6ae1554aSColin Percival if (!(inbuf = malloc((size_t)LINELENGTH))) {
75*6ae1554aSColin Percival (void)fprintf(stderr, "caesar: out of memory.\n");
76*6ae1554aSColin Percival exit(1);
77*6ae1554aSColin Percival }
78*6ae1554aSColin Percival
79*6ae1554aSColin Percival /* adjust frequency table to weight low probs REAL low */
80*6ae1554aSColin Percival for (i = 0; i < 26; ++i)
81*6ae1554aSColin Percival stdf[i] = log(stdf[i]) + log(26.0 / 100.0);
82*6ae1554aSColin Percival
83*6ae1554aSColin Percival /* zero out observation table */
84*6ae1554aSColin Percival bzero(obs, 26 * sizeof(int));
85*6ae1554aSColin Percival
86*6ae1554aSColin Percival if ((nread = read(STDIN_FILENO, inbuf, (size_t)LINELENGTH)) < 0) {
87*6ae1554aSColin Percival (void)fprintf(stderr, "caesar: %s\n", strerror(errno));
88*6ae1554aSColin Percival exit(1);
89*6ae1554aSColin Percival }
90*6ae1554aSColin Percival for (i = nread; i--;) {
91*6ae1554aSColin Percival ch = (unsigned char) inbuf[i];
92*6ae1554aSColin Percival if (isascii(ch)) {
93*6ae1554aSColin Percival if (islower(ch))
94*6ae1554aSColin Percival ++obs[ch - 'a'];
95*6ae1554aSColin Percival else if (isupper(ch))
96*6ae1554aSColin Percival ++obs[ch - 'A'];
97*6ae1554aSColin Percival }
98*6ae1554aSColin Percival }
99*6ae1554aSColin Percival
100*6ae1554aSColin Percival /*
101*6ae1554aSColin Percival * now "dot" the freqs with the observed letter freqs
102*6ae1554aSColin Percival * and keep track of best fit
103*6ae1554aSColin Percival */
104*6ae1554aSColin Percival for (try = winner = 0; try < 26; ++try) { /* += 13) { */
105*6ae1554aSColin Percival dot = 0;
106*6ae1554aSColin Percival for (i = 0; i < 26; i++)
107*6ae1554aSColin Percival dot += obs[i] * stdf[(i + try) % 26];
108*6ae1554aSColin Percival /* initialize winning score */
109*6ae1554aSColin Percival if (try == 0)
110*6ae1554aSColin Percival winnerdot = dot;
111*6ae1554aSColin Percival if (dot > winnerdot) {
112*6ae1554aSColin Percival /* got a new winner! */
113*6ae1554aSColin Percival winner = try;
114*6ae1554aSColin Percival winnerdot = dot;
115*6ae1554aSColin Percival }
116*6ae1554aSColin Percival }
117*6ae1554aSColin Percival
118*6ae1554aSColin Percival for (;;) {
119*6ae1554aSColin Percival for (i = 0; i < nread; ++i) {
120*6ae1554aSColin Percival ch = (unsigned char) inbuf[i];
121*6ae1554aSColin Percival putchar(ROTATE(ch, winner));
122*6ae1554aSColin Percival }
123*6ae1554aSColin Percival if (nread < LINELENGTH)
124*6ae1554aSColin Percival break;
125*6ae1554aSColin Percival if ((nread = read(STDIN_FILENO, inbuf, (size_t)LINELENGTH)) < 0) {
126*6ae1554aSColin Percival (void)fprintf(stderr, "caesar: %s\n", strerror(errno));
127*6ae1554aSColin Percival exit(1);
128*6ae1554aSColin Percival }
129*6ae1554aSColin Percival }
130*6ae1554aSColin Percival exit(0);
131*6ae1554aSColin Percival }
132*6ae1554aSColin Percival
133*6ae1554aSColin Percival static void
printit(char * arg)134*6ae1554aSColin Percival printit(char *arg)
135*6ae1554aSColin Percival {
136*6ae1554aSColin Percival int ch, rot;
137*6ae1554aSColin Percival
138*6ae1554aSColin Percival if ((rot = atoi(arg)) < 0) {
139*6ae1554aSColin Percival (void)fprintf(stderr, "caesar: bad rotation value.\n");
140*6ae1554aSColin Percival exit(1);
141*6ae1554aSColin Percival }
142*6ae1554aSColin Percival while ((ch = getchar()) != EOF)
143*6ae1554aSColin Percival putchar(ROTATE(ch, rot));
144*6ae1554aSColin Percival exit(0);
145*6ae1554aSColin Percival }
146