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