2 Unix SMB/Netbios implementation.
4 a implementation of MD4 designed for use in the SMB authentication protocol
5 Copyright (C) Andrew Tridgell 1997-1998.
7 This program is free software; you can redistribute it and/or modify
8 it under the terms of the GNU General Public License as published by
9 the Free Software Foundation; either version 2 of the License, or
10 (at your option) any later version.
12 This program is distributed in the hope that it will be useful,
13 but WITHOUT ANY WARRANTY; without even the implied warranty of
14 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 GNU General Public License for more details.
17 You should have received a copy of the GNU General Public License
18 along with this program; if not, write to the Free Software
19 Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
24 /* NOTE: This code makes no attempt to be fast!
26 It assumes that a int is at least 32 bits long
29 static struct mdfour *m;
31 static uint32 F(uint32 X, uint32 Y, uint32 Z)
33 return (X&Y) | ((~X)&Z);
36 static uint32 G(uint32 X, uint32 Y, uint32 Z)
38 return (X&Y) | (X&Z) | (Y&Z);
41 static uint32 H(uint32 X, uint32 Y, uint32 Z)
46 static uint32 lshift(uint32 x, int s)
49 return ((x<<s)&0xFFFFFFFF) | (x>>(32-s));
52 #define ROUND1(a,b,c,d,k,s) a = lshift(a + F(b,c,d) + X[k], s)
53 #define ROUND2(a,b,c,d,k,s) a = lshift(a + G(b,c,d) + X[k] + (uint32)0x5A827999,s)
54 #define ROUND3(a,b,c,d,k,s) a = lshift(a + H(b,c,d) + X[k] + (uint32)0x6ED9EBA1,s)
56 /* this applies md4 to 64 byte chunks */
57 static void mdfour64(uint32 *M)
60 uint32 AA, BB, CC, DD;
67 A = m->A; B = m->B; C = m->C; D = m->D;
68 AA = A; BB = B; CC = C; DD = D;
70 ROUND1(A,B,C,D, 0, 3); ROUND1(D,A,B,C, 1, 7);
71 ROUND1(C,D,A,B, 2, 11); ROUND1(B,C,D,A, 3, 19);
72 ROUND1(A,B,C,D, 4, 3); ROUND1(D,A,B,C, 5, 7);
73 ROUND1(C,D,A,B, 6, 11); ROUND1(B,C,D,A, 7, 19);
74 ROUND1(A,B,C,D, 8, 3); ROUND1(D,A,B,C, 9, 7);
75 ROUND1(C,D,A,B, 10, 11); ROUND1(B,C,D,A, 11, 19);
76 ROUND1(A,B,C,D, 12, 3); ROUND1(D,A,B,C, 13, 7);
77 ROUND1(C,D,A,B, 14, 11); ROUND1(B,C,D,A, 15, 19);
79 ROUND2(A,B,C,D, 0, 3); ROUND2(D,A,B,C, 4, 5);
80 ROUND2(C,D,A,B, 8, 9); ROUND2(B,C,D,A, 12, 13);
81 ROUND2(A,B,C,D, 1, 3); ROUND2(D,A,B,C, 5, 5);
82 ROUND2(C,D,A,B, 9, 9); ROUND2(B,C,D,A, 13, 13);
83 ROUND2(A,B,C,D, 2, 3); ROUND2(D,A,B,C, 6, 5);
84 ROUND2(C,D,A,B, 10, 9); ROUND2(B,C,D,A, 14, 13);
85 ROUND2(A,B,C,D, 3, 3); ROUND2(D,A,B,C, 7, 5);
86 ROUND2(C,D,A,B, 11, 9); ROUND2(B,C,D,A, 15, 13);
88 ROUND3(A,B,C,D, 0, 3); ROUND3(D,A,B,C, 8, 9);
89 ROUND3(C,D,A,B, 4, 11); ROUND3(B,C,D,A, 12, 15);
90 ROUND3(A,B,C,D, 2, 3); ROUND3(D,A,B,C, 10, 9);
91 ROUND3(C,D,A,B, 6, 11); ROUND3(B,C,D,A, 14, 15);
92 ROUND3(A,B,C,D, 1, 3); ROUND3(D,A,B,C, 9, 9);
93 ROUND3(C,D,A,B, 5, 11); ROUND3(B,C,D,A, 13, 15);
94 ROUND3(A,B,C,D, 3, 3); ROUND3(D,A,B,C, 11, 9);
95 ROUND3(C,D,A,B, 7, 11); ROUND3(B,C,D,A, 15, 15);
97 A += AA; B += BB; C += CC; D += DD;
99 A &= 0xFFFFFFFF; B &= 0xFFFFFFFF;
100 C &= 0xFFFFFFFF; D &= 0xFFFFFFFF;
105 m->A = A; m->B = B; m->C = C; m->D = D;
108 static void copy64(uint32 *M, unsigned char *in)
113 M[i] = (in[i*4+3]<<24) | (in[i*4+2]<<16) |
114 (in[i*4+1]<<8) | (in[i*4+0]<<0);
117 static void copy4(unsigned char *out,uint32 x)
120 out[1] = (x>>8)&0xFF;
121 out[2] = (x>>16)&0xFF;
122 out[3] = (x>>24)&0xFF;
125 void mdfour_begin(struct mdfour *md)
135 static void mdfour_tail(unsigned char *in, int n)
137 unsigned char buf[128];
146 if (n) memcpy(buf, in, n);
165 void mdfour_update(struct mdfour *md, unsigned char *in, int n)
169 if (n == 0) mdfour_tail(in, n);
181 if (n) mdfour_tail(in, n);
185 void mdfour_result(struct mdfour *md, unsigned char *out)
196 void mdfour(unsigned char *out, unsigned char *in, int n)
200 mdfour_update(&md, in, n);
201 mdfour_result(&md, out);
205 static void file_checksum1(char *fname)
209 unsigned char buf[64], sum[16];
211 fd = open(fname,O_RDONLY);
220 int n = read(fd, buf, sizeof(buf));
222 mdfour_update(&md, buf, n);
227 mdfour_result(&md, sum);
230 printf("%02X", sum[i]);
236 static void file_checksum2(char *fname)
240 unsigned char buf[64], sum[16];
242 fd = open(fname,O_RDONLY);
251 int n = read(fd, buf, sizeof(buf));
253 MDupdate(&md, buf, n*8);
257 MDupdate(&md, buf, 0);
262 memcpy(sum, md.buffer, 16);
265 printf("%02X", sum[i]);
269 int main(int argc, char *argv[])
271 file_checksum1(argv[1]);
272 file_checksum2(argv[1]);