Home | History | Annotate | Line # | Download | only in caesar
caesar.c revision 1.20
      1 /*	$NetBSD: caesar.c,v 1.20 2005/11/19 14:22:21 rillig Exp $	*/
      2 
      3 /*
      4  * Copyright (c) 1989, 1993
      5  *	The Regents of the University of California.  All rights reserved.
      6  *
      7  * This code is derived from software contributed to Berkeley by
      8  * Rick Adams.
      9  *
     10  * Authors:
     11  *	Stan King, John Eldridge, based on algorithm suggested by
     12  *	Bob Morris
     13  * 29-Sep-82
     14  *      Roland Illig, 2005
     15  *
     16  * Redistribution and use in source and binary forms, with or without
     17  * modification, are permitted provided that the following conditions
     18  * are met:
     19  * 1. Redistributions of source code must retain the above copyright
     20  *    notice, this list of conditions and the following disclaimer.
     21  * 2. Redistributions in binary form must reproduce the above copyright
     22  *    notice, this list of conditions and the following disclaimer in the
     23  *    documentation and/or other materials provided with the distribution.
     24  * 3. Neither the name of the University nor the names of its contributors
     25  *    may be used to endorse or promote products derived from this software
     26  *    without specific prior written permission.
     27  *
     28  * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
     29  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
     30  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
     31  * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
     32  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
     33  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
     34  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
     35  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
     36  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
     37  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
     38  * SUCH DAMAGE.
     39  */
     40 
     41 #include <sys/cdefs.h>
     42 #ifndef lint
     43 __COPYRIGHT("@(#) Copyright (c) 1989, 1993\n\
     44 	The Regents of the University of California.  All rights reserved.\n");
     45 #endif /* not lint */
     46 
     47 #ifndef lint
     48 #if 0
     49 static char sccsid[] = "@(#)caesar.c	8.1 (Berkeley) 5/31/93";
     50 #else
     51 __RCSID("$NetBSD: caesar.c,v 1.20 2005/11/19 14:22:21 rillig Exp $");
     52 #endif
     53 #endif /* not lint */
     54 
     55 #include <err.h>
     56 #include <errno.h>
     57 #include <limits.h>
     58 #include <math.h>
     59 #include <stdio.h>
     60 #include <string.h>
     61 #include <stdlib.h>
     62 
     63 #define NCHARS			(1 << CHAR_BIT)
     64 #define LETTERS			(26)
     65 
     66 /*
     67  * letter frequencies (taken from some unix(tm) documentation)
     68  * (unix is a trademark of Bell Laboratories)
     69  */
     70 static const unsigned char upper[LETTERS] = "ABCDEFGHIJKLMNOPQRSTUVWXYZ";
     71 static const unsigned char lower[LETTERS] = "abcdefghijklmnopqrstuvwxyz";
     72 static double stdf[LETTERS] = {
     73 	7.97, 1.35, 3.61, 4.78, 12.37, 2.01, 1.46, 4.49, 6.39, 0.04,
     74 	0.42, 3.81, 2.69, 5.92,  6.96, 2.91, 0.08, 6.63, 8.77, 9.68,
     75 	2.62, 0.81, 1.88, 0.23,  2.07, 0.06
     76 };
     77 
     78 static unsigned char rottbl[NCHARS];
     79 
     80 
     81 static void
     82 init_rottbl(int rot)
     83 {
     84 	size_t i;
     85 
     86 	rot %= LETTERS;		/* prevent integer overflow */
     87 
     88 	for (i = 0; i < NCHARS; i++)
     89 		rottbl[i] = (unsigned char)i;
     90 
     91 	for (i = 0; i < LETTERS; i++)
     92 		rottbl[upper[i]] = upper[(i + rot) % LETTERS];
     93 
     94 	for (i = 0; i < LETTERS; i++)
     95 		rottbl[lower[i]] = lower[(i + rot) % LETTERS];
     96 }
     97 
     98 static void
     99 print_file(void)
    100 {
    101 	int ch;
    102 
    103 	while ((ch = getchar()) != EOF) {
    104 		if (putchar(rottbl[ch]) == EOF) {
    105 			err(EXIT_FAILURE, "<stdout>");
    106 			/* NOTREACHED */
    107 		}
    108 	}
    109 }
    110 
    111 static void
    112 print_array(const unsigned char *a, size_t len)
    113 {
    114 	size_t i;
    115 
    116 	for (i = 0; i < len; i++) {
    117 		if (putchar(rottbl[a[i]]) == EOF) {
    118 			err(EXIT_FAILURE, "<stdout>");
    119 			/* NOTREACHED */
    120 		}
    121 	}
    122 }
    123 
    124 static int
    125 get_rotation(const char *arg)
    126 {
    127 	long rot;
    128 	char *endp;
    129 
    130 	errno = 0;
    131 	rot = strtol(arg, &endp, 10);
    132 	if (!(errno == 0 && *endp == '\0' && 0 <= rot && rot <= INT_MAX)) {
    133 		errx(EXIT_FAILURE, "bad rotation value: %s", arg);
    134 		/* NOTREACHED */
    135 	}
    136 	return (int) rot;
    137 }
    138 
    139 static void
    140 guess_and_rotate(void)
    141 {
    142 	unsigned char inbuf[2048];
    143 	unsigned int obs[NCHARS];
    144 	size_t i, nread;
    145 	double dot, winnerdot;
    146 	int try, winner;
    147 	int ch;
    148 
    149 	/* adjust frequency table to weight low probs REAL low */
    150 	for (i = 0; i < LETTERS; i++)
    151 		stdf[i] = log(stdf[i]) + log(LETTERS / 100.0);
    152 
    153 	/* zero out observation table */
    154 	(void)memset(obs, 0, sizeof(obs));
    155 
    156 	for (nread = 0; nread < sizeof(inbuf); nread++) {
    157 		if ((ch = getchar()) == EOF)
    158 			break;
    159 		inbuf[nread] = (unsigned char) ch;
    160 	}
    161 
    162 	for (i = 0; i < nread; i++)
    163 		obs[inbuf[i]]++;
    164 
    165 	/*
    166 	 * now "dot" the freqs with the observed letter freqs
    167 	 * and keep track of best fit
    168 	 */
    169 	winner = 0;
    170 	winnerdot = 0.0;
    171 	for (try = 0; try < LETTERS; try++) {
    172 		dot = 0.0;
    173 		for (i = 0; i < LETTERS; i++)
    174 			dot += (obs[upper[i]] + obs[lower[i]])
    175 			     * stdf[(i + try) % LETTERS];
    176 
    177 		if (try == 0 || dot > winnerdot) {
    178 			/* got a new winner! */
    179 			winner = try;
    180 			winnerdot = dot;
    181 		}
    182 	}
    183 
    184 	init_rottbl(winner);
    185 	print_array(inbuf, nread);
    186 	print_file();
    187 }
    188 
    189 int
    190 main(int argc, char **argv)
    191 {
    192 
    193 	if (argc == 1) {
    194 		guess_and_rotate();
    195 	} else if (argc == 2) {
    196 		init_rottbl(get_rotation(argv[1]));
    197 		print_file();
    198 	} else {
    199 		(void)fprintf(stderr, "usage: caesar [rotation]\n");
    200 		exit(EXIT_FAILURE);
    201 		/* NOTREACHED */
    202 	}
    203 
    204 	if (ferror(stdin)) {
    205 		errx(EXIT_FAILURE, "<stdin>");
    206 		/* NOTREACHED */
    207 	}
    208 
    209 	(void)fflush(stdout);
    210 	if (ferror(stdout)) {
    211 		errx(EXIT_FAILURE, "<stdout>");
    212 		/* NOTREACHED */
    213 	}
    214 
    215 	return 0;
    216 }
    217