md5.cpp 9.9 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268
  1. /*
  2. **********************************************************************
  3. ** md5.c **
  4. ** RSA Data Security, Inc. MD5 Message Digest Algorithm **
  5. ** Created: 2/17/90 RLR **
  6. ** Revised: 1/91 SRD,AJ,BSK,JT Reference C Version **
  7. **********************************************************************
  8. */
  9. /*
  10. **********************************************************************
  11. ** Copyright (C) 1990, RSA Data Security, Inc. All rights reserved. **
  12. ** **
  13. ** License to copy and use this software is granted provided that **
  14. ** it is identified as the "RSA Data Security, Inc. MD5 Message **
  15. ** Digest Algorithm" in all material mentioning or referencing this **
  16. ** software or this function. **
  17. ** **
  18. ** License is also granted to make and use derivative works **
  19. ** provided that such works are identified as "derived from the RSA **
  20. ** Data Security, Inc. MD5 Message Digest Algorithm" in all **
  21. ** material mentioning or referencing the derived work. **
  22. ** **
  23. ** RSA Data Security, Inc. makes no representations concerning **
  24. ** either the merchantability of this software or the suitability **
  25. ** of this software for any particular purpose. It is provided "as **
  26. ** is" without express or implied warranty of any kind. **
  27. ** **
  28. ** These notices must be retained in any copies of any part of this **
  29. ** documentation and/or software. **
  30. **********************************************************************
  31. */
  32. /* -- include the following line if the md5.h header file is separate -- */
  33. #include "md5.h"
  34. /* forward declaration */
  35. static void Transform (uint32_t *buf, uint32_t *in);
  36. static unsigned char PADDING[64] = {
  37. 0x80, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  38. 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  39. 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  40. 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  41. 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  42. 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  43. 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  44. 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00
  45. };
  46. /* F, G and H are basic MD5 functions: selection, majority, parity */
  47. #define F(x, y, z) (((x) & (y)) | ((~x) & (z)))
  48. #define G(x, y, z) (((x) & (z)) | ((y) & (~z)))
  49. #define H(x, y, z) ((x) ^ (y) ^ (z))
  50. #define I(x, y, z) ((y) ^ ((x) | (~z)))
  51. /* ROTATE_LEFT rotates x left n bits */
  52. #define ROTATE_LEFT(x, n) (((x) << (n)) | ((x) >> (32-(n))))
  53. /* FF, GG, HH, and II transformations for rounds 1, 2, 3, and 4 */
  54. /* Rotation is separate from addition to prevent recomputation */
  55. #define FF(a, b, c, d, x, s, ac) \
  56. {(a) += F ((b), (c), (d)) + (x) + (uint32_t)(ac); \
  57. (a) = ROTATE_LEFT ((a), (s)); \
  58. (a) += (b); \
  59. }
  60. #define GG(a, b, c, d, x, s, ac) \
  61. {(a) += G ((b), (c), (d)) + (x) + (uint32_t)(ac); \
  62. (a) = ROTATE_LEFT ((a), (s)); \
  63. (a) += (b); \
  64. }
  65. #define HH(a, b, c, d, x, s, ac) \
  66. {(a) += H ((b), (c), (d)) + (x) + (uint32_t)(ac); \
  67. (a) = ROTATE_LEFT ((a), (s)); \
  68. (a) += (b); \
  69. }
  70. #define II(a, b, c, d, x, s, ac) \
  71. {(a) += I ((b), (c), (d)) + (x) + (uint32_t)(ac); \
  72. (a) = ROTATE_LEFT ((a), (s)); \
  73. (a) += (b); \
  74. }
  75. void MD5Init (MD5_CTX *mdContext)
  76. {
  77. mdContext->i[0] = mdContext->i[1] = (uint32_t)0;
  78. /* Load magic initialization constants.
  79. */
  80. mdContext->buf[0] = (uint32_t)0x67452301;
  81. mdContext->buf[1] = (uint32_t)0xefcdab89;
  82. mdContext->buf[2] = (uint32_t)0x98badcfe;
  83. mdContext->buf[3] = (uint32_t)0x10325476;
  84. }
  85. void MD5Update (MD5_CTX *mdContext,unsigned char *inBuf,unsigned int inLen) {
  86. uint32_t in[16];
  87. int mdi;
  88. unsigned int i, ii;
  89. /* compute number of bytes mod 64 */
  90. mdi = (int)((mdContext->i[0] >> 3) & 0x3F);
  91. /* update number of bits */
  92. if ((mdContext->i[0] + ((uint32_t)inLen << 3)) < mdContext->i[0])
  93. mdContext->i[1]++;
  94. mdContext->i[0] += ((uint32_t)inLen << 3);
  95. mdContext->i[1] += ((uint32_t)inLen >> 29);
  96. while (inLen--) {
  97. /* add new character to buffer, increment mdi */
  98. mdContext->in[mdi++] = *inBuf++;
  99. /* transform if necessary */
  100. if (mdi == 0x40) {
  101. for (i = 0, ii = 0; i < 16; i++, ii += 4)
  102. in[i] = (((uint32_t)mdContext->in[ii+3]) << 24) |
  103. (((uint32_t)mdContext->in[ii+2]) << 16) |
  104. (((uint32_t)mdContext->in[ii+1]) << 8) |
  105. ((uint32_t)mdContext->in[ii]);
  106. Transform (mdContext->buf, in);
  107. mdi = 0;
  108. }
  109. }
  110. }
  111. void MD5Final (MD5_CTX *mdContext) {
  112. uint32_t in[16];
  113. int mdi;
  114. unsigned int i, ii;
  115. unsigned int padLen;
  116. /* save number of bits */
  117. in[14] = mdContext->i[0];
  118. in[15] = mdContext->i[1];
  119. /* compute number of bytes mod 64 */
  120. mdi = (int)((mdContext->i[0] >> 3) & 0x3F);
  121. /* pad out to 56 mod 64 */
  122. padLen = (mdi < 56) ? (56 - mdi) : (120 - mdi);
  123. MD5Update (mdContext, PADDING, padLen);
  124. /* append length in bits and transform */
  125. for (i = 0, ii = 0; i < 14; i++, ii += 4)
  126. in[i] = (((uint32_t)mdContext->in[ii+3]) << 24) |
  127. (((uint32_t)mdContext->in[ii+2]) << 16) |
  128. (((uint32_t)mdContext->in[ii+1]) << 8) |
  129. ((uint32_t)mdContext->in[ii]);
  130. Transform (mdContext->buf, in);
  131. /* store buffer in digest */
  132. for (i = 0, ii = 0; i < 4; i++, ii += 4) {
  133. mdContext->digest[ii] = (unsigned char)(mdContext->buf[i] & 0xFF);
  134. mdContext->digest[ii+1] =
  135. (unsigned char)((mdContext->buf[i] >> 8) & 0xFF);
  136. mdContext->digest[ii+2] =
  137. (unsigned char)((mdContext->buf[i] >> 16) & 0xFF);
  138. mdContext->digest[ii+3] =
  139. (unsigned char)((mdContext->buf[i] >> 24) & 0xFF);
  140. }
  141. }
  142. /* Basic MD5 step. Transform buf based on in.
  143. */
  144. static void Transform (uint32_t *buf, uint32_t *in) {
  145. uint32_t a = buf[0], b = buf[1], c = buf[2], d = buf[3];
  146. /* Round 1 */
  147. #define S11 7
  148. #define S12 12
  149. #define S13 17
  150. #define S14 22
  151. FF ( a, b, c, d, in[ 0], S11, 3614090360); /* 1 */
  152. FF ( d, a, b, c, in[ 1], S12, 3905402710); /* 2 */
  153. FF ( c, d, a, b, in[ 2], S13, 606105819); /* 3 */
  154. FF ( b, c, d, a, in[ 3], S14, 3250441966); /* 4 */
  155. FF ( a, b, c, d, in[ 4], S11, 4118548399); /* 5 */
  156. FF ( d, a, b, c, in[ 5], S12, 1200080426); /* 6 */
  157. FF ( c, d, a, b, in[ 6], S13, 2821735955); /* 7 */
  158. FF ( b, c, d, a, in[ 7], S14, 4249261313); /* 8 */
  159. FF ( a, b, c, d, in[ 8], S11, 1770035416); /* 9 */
  160. FF ( d, a, b, c, in[ 9], S12, 2336552879); /* 10 */
  161. FF ( c, d, a, b, in[10], S13, 4294925233); /* 11 */
  162. FF ( b, c, d, a, in[11], S14, 2304563134); /* 12 */
  163. FF ( a, b, c, d, in[12], S11, 1804603682); /* 13 */
  164. FF ( d, a, b, c, in[13], S12, 4254626195); /* 14 */
  165. FF ( c, d, a, b, in[14], S13, 2792965006); /* 15 */
  166. FF ( b, c, d, a, in[15], S14, 1236535329); /* 16 */
  167. /* Round 2 */
  168. #define S21 5
  169. #define S22 9
  170. #define S23 14
  171. #define S24 20
  172. GG ( a, b, c, d, in[ 1], S21, 4129170786); /* 17 */
  173. GG ( d, a, b, c, in[ 6], S22, 3225465664); /* 18 */
  174. GG ( c, d, a, b, in[11], S23, 643717713); /* 19 */
  175. GG ( b, c, d, a, in[ 0], S24, 3921069994); /* 20 */
  176. GG ( a, b, c, d, in[ 5], S21, 3593408605); /* 21 */
  177. GG ( d, a, b, c, in[10], S22, 38016083); /* 22 */
  178. GG ( c, d, a, b, in[15], S23, 3634488961); /* 23 */
  179. GG ( b, c, d, a, in[ 4], S24, 3889429448); /* 24 */
  180. GG ( a, b, c, d, in[ 9], S21, 568446438); /* 25 */
  181. GG ( d, a, b, c, in[14], S22, 3275163606); /* 26 */
  182. GG ( c, d, a, b, in[ 3], S23, 4107603335); /* 27 */
  183. GG ( b, c, d, a, in[ 8], S24, 1163531501); /* 28 */
  184. GG ( a, b, c, d, in[13], S21, 2850285829); /* 29 */
  185. GG ( d, a, b, c, in[ 2], S22, 4243563512); /* 30 */
  186. GG ( c, d, a, b, in[ 7], S23, 1735328473); /* 31 */
  187. GG ( b, c, d, a, in[12], S24, 2368359562); /* 32 */
  188. /* Round 3 */
  189. #define S31 4
  190. #define S32 11
  191. #define S33 16
  192. #define S34 23
  193. HH ( a, b, c, d, in[ 5], S31, 4294588738); /* 33 */
  194. HH ( d, a, b, c, in[ 8], S32, 2272392833); /* 34 */
  195. HH ( c, d, a, b, in[11], S33, 1839030562); /* 35 */
  196. HH ( b, c, d, a, in[14], S34, 4259657740); /* 36 */
  197. HH ( a, b, c, d, in[ 1], S31, 2763975236); /* 37 */
  198. HH ( d, a, b, c, in[ 4], S32, 1272893353); /* 38 */
  199. HH ( c, d, a, b, in[ 7], S33, 4139469664); /* 39 */
  200. HH ( b, c, d, a, in[10], S34, 3200236656); /* 40 */
  201. HH ( a, b, c, d, in[13], S31, 681279174); /* 41 */
  202. HH ( d, a, b, c, in[ 0], S32, 3936430074); /* 42 */
  203. HH ( c, d, a, b, in[ 3], S33, 3572445317); /* 43 */
  204. HH ( b, c, d, a, in[ 6], S34, 76029189); /* 44 */
  205. HH ( a, b, c, d, in[ 9], S31, 3654602809); /* 45 */
  206. HH ( d, a, b, c, in[12], S32, 3873151461); /* 46 */
  207. HH ( c, d, a, b, in[15], S33, 530742520); /* 47 */
  208. HH ( b, c, d, a, in[ 2], S34, 3299628645); /* 48 */
  209. /* Round 4 */
  210. #define S41 6
  211. #define S42 10
  212. #define S43 15
  213. #define S44 21
  214. II ( a, b, c, d, in[ 0], S41, 4096336452); /* 49 */
  215. II ( d, a, b, c, in[ 7], S42, 1126891415); /* 50 */
  216. II ( c, d, a, b, in[14], S43, 2878612391); /* 51 */
  217. II ( b, c, d, a, in[ 5], S44, 4237533241); /* 52 */
  218. II ( a, b, c, d, in[12], S41, 1700485571); /* 53 */
  219. II ( d, a, b, c, in[ 3], S42, 2399980690); /* 54 */
  220. II ( c, d, a, b, in[10], S43, 4293915773); /* 55 */
  221. II ( b, c, d, a, in[ 1], S44, 2240044497); /* 56 */
  222. II ( a, b, c, d, in[ 8], S41, 1873313359); /* 57 */
  223. II ( d, a, b, c, in[15], S42, 4264355552); /* 58 */
  224. II ( c, d, a, b, in[ 6], S43, 2734768916); /* 59 */
  225. II ( b, c, d, a, in[13], S44, 1309151649); /* 60 */
  226. II ( a, b, c, d, in[ 4], S41, 4149444226); /* 61 */
  227. II ( d, a, b, c, in[11], S42, 3174756917); /* 62 */
  228. II ( c, d, a, b, in[ 2], S43, 718787259); /* 63 */
  229. II ( b, c, d, a, in[ 9], S44, 3951481745); /* 64 */
  230. buf[0] += a;
  231. buf[1] += b;
  232. buf[2] += c;
  233. buf[3] += d;
  234. }
  235. /*
  236. **********************************************************************
  237. ** End of md5.c **
  238. ******************************* (cut) ********************************
  239. */