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 inline uint32 F(uint32 X, uint32 Y, uint32 Z)
33 return (X&Y) | ((~X)&Z);
36 static inline uint32 G(uint32 X, uint32 Y, uint32 Z)
38 return (X&Y) | (X&Z) | (Y&Z);
41 static inline uint32 H(uint32 X, uint32 Y, uint32 Z)
46 static inline uint32 lshift(uint32 x, int s)
50 return ((x<<s)&0xFFFFFFFF) | (x>>(32-s));
52 return ((x<<s) | (x>>(32-s)));
56 #define ROUND1(a,b,c,d,k,s) a = lshift(a + F(b,c,d) + X[k], s)
57 #define ROUND2(a,b,c,d,k,s) a = lshift(a + G(b,c,d) + X[k] + 0x5A827999,s)
58 #define ROUND3(a,b,c,d,k,s) a = lshift(a + H(b,c,d) + X[k] + 0x6ED9EBA1,s)
60 /* this applies md4 to 64 byte chunks */
61 static void mdfour64(uint32 *M)
64 uint32 AA, BB, CC, DD;
71 A = m->A; B = m->B; C = m->C; D = m->D;
72 AA = A; BB = B; CC = C; DD = D;
74 ROUND1(A,B,C,D, 0, 3); ROUND1(D,A,B,C, 1, 7);
75 ROUND1(C,D,A,B, 2, 11); ROUND1(B,C,D,A, 3, 19);
76 ROUND1(A,B,C,D, 4, 3); ROUND1(D,A,B,C, 5, 7);
77 ROUND1(C,D,A,B, 6, 11); ROUND1(B,C,D,A, 7, 19);
78 ROUND1(A,B,C,D, 8, 3); ROUND1(D,A,B,C, 9, 7);
79 ROUND1(C,D,A,B, 10, 11); ROUND1(B,C,D,A, 11, 19);
80 ROUND1(A,B,C,D, 12, 3); ROUND1(D,A,B,C, 13, 7);
81 ROUND1(C,D,A,B, 14, 11); ROUND1(B,C,D,A, 15, 19);
83 ROUND2(A,B,C,D, 0, 3); ROUND2(D,A,B,C, 4, 5);
84 ROUND2(C,D,A,B, 8, 9); ROUND2(B,C,D,A, 12, 13);
85 ROUND2(A,B,C,D, 1, 3); ROUND2(D,A,B,C, 5, 5);
86 ROUND2(C,D,A,B, 9, 9); ROUND2(B,C,D,A, 13, 13);
87 ROUND2(A,B,C,D, 2, 3); ROUND2(D,A,B,C, 6, 5);
88 ROUND2(C,D,A,B, 10, 9); ROUND2(B,C,D,A, 14, 13);
89 ROUND2(A,B,C,D, 3, 3); ROUND2(D,A,B,C, 7, 5);
90 ROUND2(C,D,A,B, 11, 9); ROUND2(B,C,D,A, 15, 13);
92 ROUND3(A,B,C,D, 0, 3); ROUND3(D,A,B,C, 8, 9);
93 ROUND3(C,D,A,B, 4, 11); ROUND3(B,C,D,A, 12, 15);
94 ROUND3(A,B,C,D, 2, 3); ROUND3(D,A,B,C, 10, 9);
95 ROUND3(C,D,A,B, 6, 11); ROUND3(B,C,D,A, 14, 15);
96 ROUND3(A,B,C,D, 1, 3); ROUND3(D,A,B,C, 9, 9);
97 ROUND3(C,D,A,B, 5, 11); ROUND3(B,C,D,A, 13, 15);
98 ROUND3(A,B,C,D, 3, 3); ROUND3(D,A,B,C, 11, 9);
99 ROUND3(C,D,A,B, 7, 11); ROUND3(B,C,D,A, 15, 15);
101 A += AA; B += BB; C += CC; D += DD;
104 A &= 0xFFFFFFFF; B &= 0xFFFFFFFF;
105 C &= 0xFFFFFFFF; D &= 0xFFFFFFFF;
111 m->A = A; m->B = B; m->C = C; m->D = D;
114 static void copy64(uint32 *M, unsigned char *in)
119 M[i] = (in[i*4+3]<<24) | (in[i*4+2]<<16) |
120 (in[i*4+1]<<8) | (in[i*4+0]<<0);
123 static void copy4(unsigned char *out,uint32 x)
126 out[1] = (x>>8)&0xFF;
127 out[2] = (x>>16)&0xFF;
128 out[3] = (x>>24)&0xFF;
131 void mdfour_begin(struct mdfour *md)
141 static void mdfour_tail(unsigned char *in, int n)
143 unsigned char buf[128];
152 if (n) memcpy(buf, in, n);
168 void mdfour_update(struct mdfour *md, unsigned char *in, int n)
172 if (n == 0) mdfour_tail(in, n);
184 if (n) mdfour_tail(in, n);
188 void mdfour_result(struct mdfour *md, unsigned char *out)
199 void mdfour(unsigned char *out, unsigned char *in, int n)
203 mdfour_update(&md, in, n);
204 mdfour_result(&md, out);
208 static void file_checksum1(char *fname)
212 unsigned char buf[64], sum[16];
214 fd = open(fname,O_RDONLY);
223 int n = read(fd, buf, sizeof(buf));
225 mdfour_update(&md, buf, n);
230 mdfour_result(&md, sum);
233 printf("%02X", sum[i]);
239 static void file_checksum2(char *fname)
243 unsigned char buf[64], sum[16];
245 fd = open(fname,O_RDONLY);
254 int n = read(fd, buf, sizeof(buf));
256 MDupdate(&md, buf, n*8);
260 MDupdate(&md, buf, 0);
265 memcpy(sum, md.buffer, 16);
268 printf("%02X", sum[i]);
272 int main(int argc, char *argv[])
274 file_checksum1(argv[1]);
275 file_checksum2(argv[1]);