md5c.c 10 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338
  1. #include "md5c.h"
  2. #include <string.h>
  3. #include <stdio.h>
  4. #define S11 7
  5. #define S12 12
  6. #define S13 17
  7. #define S14 22
  8. #define S21 5
  9. #define S22 9
  10. #define S23 14
  11. #define S24 20
  12. #define S31 4
  13. #define S32 11
  14. #define S33 16
  15. #define S34 23
  16. #define S41 6
  17. #define S42 10
  18. #define S43 15
  19. #define S44 21
  20. static void MD5_memcpy (POINTER output, POINTER input, unsigned int len);
  21. static void MD5Transform (UINT4 state[4], unsigned char block[64]);
  22. static void Encode (unsigned char *output, UINT4 *input, unsigned int len);
  23. static void MD5_memset (POINTER output, int value, unsigned int len);
  24. static void Decode (UINT4 *output, unsigned char *input, unsigned int len);
  25. static unsigned char PADDING[64] = {
  26. 0x80, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
  27. 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
  28. 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
  29. };
  30. /* F, G, H and I are basic MD5 functions.
  31. */
  32. #define F(x, y, z) (((x) & (y)) | ((~x) & (z)))
  33. #define G(x, y, z) (((x) & (z)) | ((y) & (~z)))
  34. #define H(x, y, z) ((x) ^ (y) ^ (z))
  35. #define I(x, y, z) ((y) ^ ((x) | (~z)))
  36. /* ROTATE_LEFT rotates x left n bits.
  37. */
  38. #define ROTATE_LEFT(x, n) (((x) << (n)) | ((x) >> (32-(n))))
  39. /* FF, GG, HH, and II transformations for rounds 1, 2, 3, and 4.
  40. Rotation is separate from addition to prevent recomputation.
  41. */
  42. #define FF(a, b, c, d, x, s, ac) { \
  43. (a) += F ((b), (c), (d)) + (x) + (UINT4)(ac); \
  44. (a) = ROTATE_LEFT ((a), (s)); \
  45. (a) += (b); \
  46. }
  47. #define GG(a, b, c, d, x, s, ac) { \
  48. (a) += G ((b), (c), (d)) + (x) + (UINT4)(ac); \
  49. (a) = ROTATE_LEFT ((a), (s)); \
  50. (a) += (b); \
  51. }
  52. #define HH(a, b, c, d, x, s, ac) { \
  53. (a) += H ((b), (c), (d)) + (x) + (UINT4)(ac); \
  54. (a) = ROTATE_LEFT ((a), (s)); \
  55. (a) += (b); \
  56. }
  57. #define II(a, b, c, d, x, s, ac) { \
  58. (a) += I ((b), (c), (d)) + (x) + (UINT4)(ac); \
  59. (a) = ROTATE_LEFT ((a), (s)); \
  60. (a) += (b); \
  61. }
  62. /* MD5 initialization. Begins an MD5 operation, writing a new context.
  63. */
  64. void MD5Init (MD5_CTX *context) /* context */
  65. {
  66. context->count[0] = context->count[1] = 0;
  67. /* Load magic initialization constants.
  68. */
  69. context->state[0] = 0x67452301;
  70. context->state[1] = 0xefcdab89;
  71. context->state[2] = 0x98badcfe;
  72. context->state[3] = 0x10325476;
  73. }
  74. /* MD5 block update operation. Continues an MD5 message-digest
  75. operation, processing another message block, and updating the
  76. context.
  77. */
  78. void MD5Update (MD5_CTX *context, unsigned char *input, unsigned int inputLen)
  79. {
  80. unsigned int i, index, partLen;
  81. /* Compute number of bytes mod 64 */
  82. index = (unsigned int)((context->count[0] >> 3) & 0x3F);
  83. /* Update number of bits */
  84. if ((context->count[0] += ((UINT4)inputLen << 3))
  85. < ((UINT4)inputLen << 3))
  86. context->count[1]++;
  87. context->count[1] += ((UINT4)inputLen >> 29);
  88. partLen = 64 - index;
  89. /* Transform as many times as possible.
  90. */
  91. if (inputLen >= partLen) {
  92. MD5_memcpy((POINTER)&context->buffer[index], (POINTER)input, partLen);
  93. MD5Transform (context->state, context->buffer);
  94. for (i = partLen; i + 63 < inputLen; i += 64)
  95. MD5Transform (context->state, &input[i]);
  96. index = 0;
  97. }
  98. else
  99. i = 0;
  100. /* Buffer remaining input */
  101. MD5_memcpy((POINTER)&context->buffer[index], (POINTER)&input[i],inputLen-i);
  102. }
  103. /* MD5 finalization. Ends an MD5 message-digest operation, writing the
  104. the message digest and zeroizing the context.
  105. */
  106. void MD5Final (unsigned char digest[16], MD5_CTX *context)
  107. {
  108. unsigned char bits[8];
  109. unsigned int index, padLen;
  110. /* Save number of bits */
  111. Encode (bits, context->count, 8);
  112. /* Pad out to 56 mod 64.
  113. */
  114. index = (unsigned int)((context->count[0] >> 3) & 0x3f);
  115. padLen = (index < 56) ? (56 - index) : (120 - index);
  116. MD5Update (context, PADDING, padLen);
  117. /* Append length (before padding) */
  118. MD5Update (context, bits, 8);
  119. /* Store state in digest */
  120. Encode (digest, context->state, 16);
  121. /* Zeroize sensitive information.
  122. */
  123. MD5_memset ((POINTER)context, 0, sizeof (*context));
  124. }
  125. /* MD5 basic transformation. Transforms state based on block.
  126. */
  127. static void MD5Transform (UINT4 state[4], unsigned char block[64])
  128. {
  129. UINT4 a = state[0], b = state[1], c = state[2], d = state[3], x[16];
  130. Decode (x, block, 64);
  131. /* Round 1 */
  132. FF (a, b, c, d, x[ 0], S11, 0xd76aa478); /* 1 */
  133. FF (d, a, b, c, x[ 1], S12, 0xe8c7b756); /* 2 */
  134. FF (c, d, a, b, x[ 2], S13, 0x242070db); /* 3 */
  135. FF (b, c, d, a, x[ 3], S14, 0xc1bdceee); /* 4 */
  136. FF (a, b, c, d, x[ 4], S11, 0xf57c0faf); /* 5 */
  137. FF (d, a, b, c, x[ 5], S12, 0x4787c62a); /* 6 */
  138. FF (c, d, a, b, x[ 6], S13, 0xa8304613); /* 7 */
  139. FF (b, c, d, a, x[ 7], S14, 0xfd469501); /* 8 */
  140. FF (a, b, c, d, x[ 8], S11, 0x698098d8); /* 9 */
  141. FF (d, a, b, c, x[ 9], S12, 0x8b44f7af); /* 10 */
  142. FF (c, d, a, b, x[10], S13, 0xffff5bb1); /* 11 */
  143. FF (b, c, d, a, x[11], S14, 0x895cd7be); /* 12 */
  144. FF (a, b, c, d, x[12], S11, 0x6b901122); /* 13 */
  145. FF (d, a, b, c, x[13], S12, 0xfd987193); /* 14 */
  146. FF (c, d, a, b, x[14], S13, 0xa679438e); /* 15 */
  147. FF (b, c, d, a, x[15], S14, 0x49b40821); /* 16 */
  148. /* Round 2 */
  149. GG (a, b, c, d, x[ 1], S21, 0xf61e2562); /* 17 */
  150. GG (d, a, b, c, x[ 6], S22, 0xc040b340); /* 18 */
  151. GG (c, d, a, b, x[11], S23, 0x265e5a51); /* 19 */
  152. GG (b, c, d, a, x[ 0], S24, 0xe9b6c7aa); /* 20 */
  153. GG (a, b, c, d, x[ 5], S21, 0xd62f105d); /* 21 */
  154. GG (d, a, b, c, x[10], S22, 0x2441453); /* 22 */
  155. GG (c, d, a, b, x[15], S23, 0xd8a1e681); /* 23 */
  156. GG (b, c, d, a, x[ 4], S24, 0xe7d3fbc8); /* 24 */
  157. GG (a, b, c, d, x[ 9], S21, 0x21e1cde6); /* 25 */
  158. GG (d, a, b, c, x[14], S22, 0xc33707d6); /* 26 */
  159. GG (c, d, a, b, x[ 3], S23, 0xf4d50d87); /* 27 */
  160. GG (b, c, d, a, x[ 8], S24, 0x455a14ed); /* 28 */
  161. GG (a, b, c, d, x[13], S21, 0xa9e3e905); /* 29 */
  162. GG (d, a, b, c, x[ 2], S22, 0xfcefa3f8); /* 30 */
  163. GG (c, d, a, b, x[ 7], S23, 0x676f02d9); /* 31 */
  164. GG (b, c, d, a, x[12], S24, 0x8d2a4c8a); /* 32 */
  165. /* Round 3 */
  166. HH (a, b, c, d, x[ 5], S31, 0xfffa3942); /* 33 */
  167. HH (d, a, b, c, x[ 8], S32, 0x8771f681); /* 34 */
  168. HH (c, d, a, b, x[11], S33, 0x6d9d6122); /* 35 */
  169. HH (b, c, d, a, x[14], S34, 0xfde5380c); /* 36 */
  170. HH (a, b, c, d, x[ 1], S31, 0xa4beea44); /* 37 */
  171. HH (d, a, b, c, x[ 4], S32, 0x4bdecfa9); /* 38 */
  172. HH (c, d, a, b, x[ 7], S33, 0xf6bb4b60); /* 39 */
  173. HH (b, c, d, a, x[10], S34, 0xbebfbc70); /* 40 */
  174. HH (a, b, c, d, x[13], S31, 0x289b7ec6); /* 41 */
  175. HH (d, a, b, c, x[ 0], S32, 0xeaa127fa); /* 42 */
  176. HH (c, d, a, b, x[ 3], S33, 0xd4ef3085); /* 43 */
  177. HH (b, c, d, a, x[ 6], S34, 0x4881d05); /* 44 */
  178. HH (a, b, c, d, x[ 9], S31, 0xd9d4d039); /* 45 */
  179. HH (d, a, b, c, x[12], S32, 0xe6db99e5); /* 46 */
  180. HH (c, d, a, b, x[15], S33, 0x1fa27cf8); /* 47 */
  181. HH (b, c, d, a, x[ 2], S34, 0xc4ac5665); /* 48 */
  182. /* Round 4 */
  183. II (a, b, c, d, x[ 0], S41, 0xf4292244); /* 49 */
  184. II (d, a, b, c, x[ 7], S42, 0x432aff97); /* 50 */
  185. II (c, d, a, b, x[14], S43, 0xab9423a7); /* 51 */
  186. II (b, c, d, a, x[ 5], S44, 0xfc93a039); /* 52 */
  187. II (a, b, c, d, x[12], S41, 0x655b59c3); /* 53 */
  188. II (d, a, b, c, x[ 3], S42, 0x8f0ccc92); /* 54 */
  189. II (c, d, a, b, x[10], S43, 0xffeff47d); /* 55 */
  190. II (b, c, d, a, x[ 1], S44, 0x85845dd1); /* 56 */
  191. II (a, b, c, d, x[ 8], S41, 0x6fa87e4f); /* 57 */
  192. II (d, a, b, c, x[15], S42, 0xfe2ce6e0); /* 58 */
  193. II (c, d, a, b, x[ 6], S43, 0xa3014314); /* 59 */
  194. II (b, c, d, a, x[13], S44, 0x4e0811a1); /* 60 */
  195. II (a, b, c, d, x[ 4], S41, 0xf7537e82); /* 61 */
  196. II (d, a, b, c, x[11], S42, 0xbd3af235); /* 62 */
  197. II (c, d, a, b, x[ 2], S43, 0x2ad7d2bb); /* 63 */
  198. II (b, c, d, a, x[ 9], S44, 0xeb86d391); /* 64 */
  199. state[0] += a;
  200. state[1] += b;
  201. state[2] += c;
  202. state[3] += d;
  203. /* Zeroize sensitive information.
  204. */
  205. MD5_memset ((POINTER)x, 0, sizeof (x));
  206. }
  207. /* Encodes input (UINT4) into output (unsigned char). Assumes len is
  208. a multiple of 4.
  209. */
  210. static void Encode (unsigned char *output, UINT4 *input, unsigned int len)
  211. {
  212. unsigned int i, j;
  213. for (i = 0, j = 0; j < len; i++, j += 4) {
  214. output[j] = (unsigned char)(input[i] & 0xff);
  215. output[j+1] = (unsigned char)((input[i] >> 8) & 0xff);
  216. output[j+2] = (unsigned char)((input[i] >> 16) & 0xff);
  217. output[j+3] = (unsigned char)((input[i] >> 24) & 0xff);
  218. }
  219. }
  220. /* Decodes input (unsigned char) into output (UINT4). Assumes len is
  221. a multiple of 4.
  222. */
  223. static void Decode (UINT4 *output, unsigned char *input, unsigned int len)
  224. {
  225. unsigned int i, j;
  226. for (i = 0, j = 0; j < len; i++, j += 4)
  227. output[i] = ((UINT4)input[j]) | (((UINT4)input[j+1]) << 8) |
  228. (((UINT4)input[j+2]) << 16) | (((UINT4)input[j+3]) << 24);
  229. }
  230. /* Note: Replace "for loop" with standard memcpy if possible.
  231. */
  232. static void MD5_memcpy (POINTER output, POINTER input, unsigned int len)
  233. {
  234. unsigned int i;
  235. for (i = 0; i < len; i++)
  236. output[i] = input[i];
  237. }
  238. /* Note: Replace "for loop" with standard memset if possible.
  239. */
  240. static void MD5_memset (POINTER output, int value, unsigned int len)
  241. {
  242. unsigned int i;
  243. for (i = 0; i < len; i++)
  244. ((char *)output)[i] = (char)value;
  245. }
  246. /* Digests a string and prints the result.
  247. */
  248. void MDString (char *string,unsigned char digest[16])
  249. {
  250. MD5_CTX context;
  251. unsigned int len = strlen (string);
  252. MD5Init (&context);
  253. MD5Update (&context, (unsigned char *)string, len);
  254. MD5Final (digest, &context);
  255. }
  256. /* Digests a file and prints the result.
  257. */
  258. int MD5File (char *filename,unsigned char digest[16])
  259. {
  260. FILE *file;
  261. MD5_CTX context;
  262. int len;
  263. unsigned char buffer[1024];
  264. if ((file = fopen (filename, "rb")) == NULL)
  265. return -1;
  266. else {
  267. MD5Init (&context);
  268. while (len == fread (buffer, 1, 1024, file))
  269. MD5Update (&context, buffer, len);
  270. MD5Final (digest, &context);
  271. fclose (file);
  272. }
  273. return 0;
  274. }
  275. void MD5UpdaterString(MD5_CTX *context,const char *string)
  276. {
  277. unsigned int len = strlen (string);
  278. MD5Update (context, (unsigned char *)string, len);
  279. }
  280. int MD5FileUpdateFile (MD5_CTX *context,char *filename)
  281. {
  282. FILE *file;
  283. int len;
  284. unsigned char buffer[1024];
  285. if ((file = fopen (filename, "rb")) == NULL)
  286. return -1;
  287. else {
  288. while (len == fread (buffer, 1, 1024, file))
  289. MD5Update (context, buffer, len);
  290. fclose (file);
  291. }
  292. return 0;
  293. }