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 #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))
35 #define lshift(x,s) ((((x)<<(s))&0xFFFFFFFF) | (((x)>>(32-(s)))&0xFFFFFFFF))
37 #define lshift(x,s) (((x)<<(s)) | ((x)>>(32-(s))))
40 #define ROUND1(a,b,c,d,k,s) a = lshift(a + F(b,c,d) + X[k], s)
41 #define ROUND2(a,b,c,d,k,s) a = lshift(a + G(b,c,d) + X[k] + 0x5A827999,s)
42 #define ROUND3(a,b,c,d,k,s) a = lshift(a + H(b,c,d) + X[k] + 0x6ED9EBA1,s)
44 /* this applies md4 to 64 byte chunks */
45 static void mdfour64(uint32 *M)
48 uint32 AA, BB, CC, DD;
55 A = m->A; B = m->B; C = m->C; D = m->D;
56 AA = A; BB = B; CC = C; DD = D;
58 ROUND1(A,B,C,D, 0, 3); ROUND1(D,A,B,C, 1, 7);
59 ROUND1(C,D,A,B, 2, 11); ROUND1(B,C,D,A, 3, 19);
60 ROUND1(A,B,C,D, 4, 3); ROUND1(D,A,B,C, 5, 7);
61 ROUND1(C,D,A,B, 6, 11); ROUND1(B,C,D,A, 7, 19);
62 ROUND1(A,B,C,D, 8, 3); ROUND1(D,A,B,C, 9, 7);
63 ROUND1(C,D,A,B, 10, 11); ROUND1(B,C,D,A, 11, 19);
64 ROUND1(A,B,C,D, 12, 3); ROUND1(D,A,B,C, 13, 7);
65 ROUND1(C,D,A,B, 14, 11); ROUND1(B,C,D,A, 15, 19);
67 ROUND2(A,B,C,D, 0, 3); ROUND2(D,A,B,C, 4, 5);
68 ROUND2(C,D,A,B, 8, 9); ROUND2(B,C,D,A, 12, 13);
69 ROUND2(A,B,C,D, 1, 3); ROUND2(D,A,B,C, 5, 5);
70 ROUND2(C,D,A,B, 9, 9); ROUND2(B,C,D,A, 13, 13);
71 ROUND2(A,B,C,D, 2, 3); ROUND2(D,A,B,C, 6, 5);
72 ROUND2(C,D,A,B, 10, 9); ROUND2(B,C,D,A, 14, 13);
73 ROUND2(A,B,C,D, 3, 3); ROUND2(D,A,B,C, 7, 5);
74 ROUND2(C,D,A,B, 11, 9); ROUND2(B,C,D,A, 15, 13);
76 ROUND3(A,B,C,D, 0, 3); ROUND3(D,A,B,C, 8, 9);
77 ROUND3(C,D,A,B, 4, 11); ROUND3(B,C,D,A, 12, 15);
78 ROUND3(A,B,C,D, 2, 3); ROUND3(D,A,B,C, 10, 9);
79 ROUND3(C,D,A,B, 6, 11); ROUND3(B,C,D,A, 14, 15);
80 ROUND3(A,B,C,D, 1, 3); ROUND3(D,A,B,C, 9, 9);
81 ROUND3(C,D,A,B, 5, 11); ROUND3(B,C,D,A, 13, 15);
82 ROUND3(A,B,C,D, 3, 3); ROUND3(D,A,B,C, 11, 9);
83 ROUND3(C,D,A,B, 7, 11); ROUND3(B,C,D,A, 15, 15);
85 A += AA; B += BB; C += CC; D += DD;
88 A &= 0xFFFFFFFF; B &= 0xFFFFFFFF;
89 C &= 0xFFFFFFFF; D &= 0xFFFFFFFF;
95 m->A = A; m->B = B; m->C = C; m->D = D;
98 static void copy64(uint32 *M, unsigned char *in)
103 M[i] = (in[i*4+3]<<24) | (in[i*4+2]<<16) |
104 (in[i*4+1]<<8) | (in[i*4+0]<<0);
107 static void copy4(unsigned char *out,uint32 x)
110 out[1] = (x>>8)&0xFF;
111 out[2] = (x>>16)&0xFF;
112 out[3] = (x>>24)&0xFF;
115 void mdfour_begin(struct mdfour *md)
125 static void mdfour_tail(unsigned char *in, int n)
127 unsigned char buf[128];
136 if (n) memcpy(buf, in, n);
152 void mdfour_update(struct mdfour *md, unsigned char *in, int n)
156 if (n == 0) mdfour_tail(in, n);
168 if (n) mdfour_tail(in, n);
172 void mdfour_result(struct mdfour *md, unsigned char *out)
183 void mdfour(unsigned char *out, unsigned char *in, int n)
187 mdfour_update(&md, in, n);
188 mdfour_result(&md, out);
192 static void file_checksum1(char *fname)
196 unsigned char buf[64*1024], sum[16];
198 fd = open(fname,O_RDONLY);
207 int n = read(fd, buf, sizeof(buf));
209 mdfour_update(&md, buf, n);
214 mdfour_result(&md, sum);
217 printf("%02X", sum[i]);
224 static void file_checksum2(char *fname)
228 unsigned char buf[64], sum[16];
230 fd = open(fname,O_RDONLY);
239 int n = read(fd, buf, sizeof(buf));
241 MDupdate(&md, buf, n*8);
245 MDupdate(&md, buf, 0);
250 memcpy(sum, md.buffer, 16);
253 printf("%02X", sum[i]);
258 int main(int argc, char *argv[])
260 file_checksum1(argv[1]);
262 file_checksum2(argv[1]);