2 a implementation of MD4 designed for use in the SMB authentication protocol
3 Copyright (C) Andrew Tridgell 1997-1998.
5 This program is free software; you can redistribute it and/or modify
6 it under the terms of the GNU General Public License as published by
7 the Free Software Foundation; either version 2 of the License, or
8 (at your option) any later version.
10 This program is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 GNU General Public License for more details.
15 You should have received a copy of the GNU General Public License
16 along with this program; if not, write to the Free Software
17 Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
22 /* NOTE: This code makes no attempt to be fast!
24 It assumes that a int is at least 32 bits long
27 static struct mdfour *m;
29 #define MASK32 (0xffffffff)
31 #define F(X,Y,Z) ((((X)&(Y)) | ((~(X))&(Z))))
32 #define G(X,Y,Z) ((((X)&(Y)) | ((X)&(Z)) | ((Y)&(Z))))
33 #define H(X,Y,Z) (((X)^(Y)^(Z)))
34 #define lshift(x,s) (((((x)<<(s))&MASK32) | (((x)>>(32-(s)))&MASK32)))
36 #define ROUND1(a,b,c,d,k,s) a = lshift((a + F(b,c,d) + M[k])&MASK32, s)
37 #define ROUND2(a,b,c,d,k,s) a = lshift((a + G(b,c,d) + M[k] + 0x5A827999)&MASK32,s)
38 #define ROUND3(a,b,c,d,k,s) a = lshift((a + H(b,c,d) + M[k] + 0x6ED9EBA1)&MASK32,s)
40 /* this applies md4 to 64 byte chunks */
41 static void mdfour64(uint32 *M)
43 uint32 AA, BB, CC, DD;
46 A = m->A; B = m->B; C = m->C; D = m->D;
47 AA = A; BB = B; CC = C; DD = D;
49 ROUND1(A,B,C,D, 0, 3); ROUND1(D,A,B,C, 1, 7);
50 ROUND1(C,D,A,B, 2, 11); ROUND1(B,C,D,A, 3, 19);
51 ROUND1(A,B,C,D, 4, 3); ROUND1(D,A,B,C, 5, 7);
52 ROUND1(C,D,A,B, 6, 11); ROUND1(B,C,D,A, 7, 19);
53 ROUND1(A,B,C,D, 8, 3); ROUND1(D,A,B,C, 9, 7);
54 ROUND1(C,D,A,B, 10, 11); ROUND1(B,C,D,A, 11, 19);
55 ROUND1(A,B,C,D, 12, 3); ROUND1(D,A,B,C, 13, 7);
56 ROUND1(C,D,A,B, 14, 11); ROUND1(B,C,D,A, 15, 19);
59 ROUND2(A,B,C,D, 0, 3); ROUND2(D,A,B,C, 4, 5);
60 ROUND2(C,D,A,B, 8, 9); ROUND2(B,C,D,A, 12, 13);
61 ROUND2(A,B,C,D, 1, 3); ROUND2(D,A,B,C, 5, 5);
62 ROUND2(C,D,A,B, 9, 9); ROUND2(B,C,D,A, 13, 13);
63 ROUND2(A,B,C,D, 2, 3); ROUND2(D,A,B,C, 6, 5);
64 ROUND2(C,D,A,B, 10, 9); ROUND2(B,C,D,A, 14, 13);
65 ROUND2(A,B,C,D, 3, 3); ROUND2(D,A,B,C, 7, 5);
66 ROUND2(C,D,A,B, 11, 9); ROUND2(B,C,D,A, 15, 13);
68 ROUND3(A,B,C,D, 0, 3); ROUND3(D,A,B,C, 8, 9);
69 ROUND3(C,D,A,B, 4, 11); ROUND3(B,C,D,A, 12, 15);
70 ROUND3(A,B,C,D, 2, 3); ROUND3(D,A,B,C, 10, 9);
71 ROUND3(C,D,A,B, 6, 11); ROUND3(B,C,D,A, 14, 15);
72 ROUND3(A,B,C,D, 1, 3); ROUND3(D,A,B,C, 9, 9);
73 ROUND3(C,D,A,B, 5, 11); ROUND3(B,C,D,A, 13, 15);
74 ROUND3(A,B,C,D, 3, 3); ROUND3(D,A,B,C, 11, 9);
75 ROUND3(C,D,A,B, 7, 11); ROUND3(B,C,D,A, 15, 15);
80 A &= MASK32; B &= MASK32;
81 C &= MASK32; D &= MASK32;
83 m->A = A; m->B = B; m->C = C; m->D = D;
86 static void copy64(uint32 *M, const unsigned char *in)
91 M[i] = (in[i*4+3]<<24) | (in[i*4+2]<<16) |
92 (in[i*4+1]<<8) | (in[i*4+0]<<0);
95 static void copy4(unsigned char *out,uint32 x)
99 out[2] = (x>>16)&0xFF;
100 out[3] = (x>>24)&0xFF;
103 void mdfour_begin(struct mdfour *md)
114 static void mdfour_tail(const unsigned char *in, int n)
116 unsigned char buf[128];
125 if (n) memcpy(buf, in, n);
141 void mdfour_update(struct mdfour *md, const unsigned char *in, int n)
148 mdfour_tail(md->tail, md->tail_len);
153 int len = 64 - md->tail_len;
154 if (len > n) len = n;
155 memcpy(md->tail+md->tail_len, in, len);
159 if (md->tail_len == 64) {
176 memcpy(md->tail, in, n);
182 void mdfour_result(struct mdfour *md, unsigned char *out)
193 void mdfour(unsigned char *out, const unsigned char *in, int n)
197 mdfour_update(&md, in, n);
198 mdfour_update(&md, NULL, 0);
199 mdfour_result(&md, out);
203 static void file_checksum1(char *fname)
207 unsigned char buf[1024], sum[16];
210 fd = open(fname,O_RDONLY|O_BINARY);
216 chunk = 1 + random() % (sizeof(buf) - 1);
221 int n = read(fd, buf, chunk);
223 mdfour_update(&md, buf, n);
225 if (n < chunk) break;
230 mdfour_update(&md, NULL, 0);
232 mdfour_result(&md, sum);
235 printf("%02x", sum[i]);
242 static void file_checksum2(char *fname)
246 unsigned char buf[64], sum[16];
248 fd = open(fname,O_RDONLY|O_BINARY);
257 int n = read(fd, buf, sizeof(buf));
259 MDupdate(&md, buf, n*8);
263 MDupdate(&md, buf, 0);
268 memcpy(sum, md.buffer, 16);
271 printf("%02x", sum[i]);
276 int main(int argc, char *argv[])
278 file_checksum1(argv[1]);
280 file_checksum2(argv[1]);