Chinaunix首页 | 论坛 | 博客

wangchenxicoolwcx.blog.chinaunix.net

首页　| 　博文目录　| 　关于我

wangchenxicool

  • 博客访问： 4263744
  • 博文数量： 776
  • 博客积分： 13014
  • 博客等级： 上将
  • 技术积分： 10391
  • 用 户 组： 普通用户
  • 注册时间： 2010-02-22 17:00
文章分类

全部博文（776）

文章存档

2015年（55）

2014年（43）

2013年（147）

2012年（20）

2011年（82）

2010年（429）

我的朋友
最近访客
推荐博文
相关博文
md5

分类： C/C++

2010-11-30 16:12:53

 

  1. /* typedef a 32 bit type */
  2. typedef unsigned long int UINT4;
  3. /* Data structure for MD5 (Message Digest) computation */
  4. typedef struct {
  5.         UINT4 i[2]; /* number of _bits_ handled mod 2^64 */
  6.         UINT4 buf[4]; /* scratch buffer */
  7.         unsigned char in[64]; /* input buffer */
  8.         unsigned char digest[16]; /* actual digest after MD5Final call */
  9. } MD5_CTX;

  11. void MD5Init ();
  12. void MD5Update ();
  13. void MD5Final ();
  14. /* forward declaration */
  15. static void Transform ();

  17. static unsigned char PADDING[64] = {
  18.         0x80, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  19.         0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  20.         0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  21.         0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  22.         0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  23.         0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  24.         0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
  25.         0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00
  26. };

  28. /* F, G and H are basic MD5 functions: selection, majority, parity */
  29. #define F(x, y, z) (((x) & (y)) | ((~x) & (z)))
  30. #define G(x, y, z) (((x) & (z)) | ((y) & (~z)))
  31. #define H(x, y, z) ((x) ^ (y) ^ (z))
  32. #define I(x, y, z) ((y) ^ ((x) | (~z)))
  33. /* ROTATE_LEFT rotates x left n bits */
  34. #define ROTATE_LEFT(x, n) (((x) << (n)) | ((x) >> (32-(n))))
  35. /* FF, GG, HH, and II transformations for rounds 1, 2, 3, and 4 */
  36. /* Rotation is separate from addition to prevent recomputation */
  37. #define FF(a, b, c, d, x, s, ac) \
  38.   {(a) += F ((b), (c), (d)) + (x) + (UINT4)(ac); \
  39.    (a) = ROTATE_LEFT ((a), (s)); \
  40.    (a) += (b); \
  41.   }
  42. #define GG(a, b, c, d, x, s, ac) \
  43.   {(a) += G ((b), (c), (d)) + (x) + (UINT4)(ac); \
  44.    (a) = ROTATE_LEFT ((a), (s)); \
  45.    (a) += (b); \
  46.   }
  47. #define HH(a, b, c, d, x, s, ac) \
  48.   {(a) += H ((b), (c), (d)) + (x) + (UINT4)(ac); \
  49.    (a) = ROTATE_LEFT ((a), (s)); \
  50.    (a) += (b); \
  51.   }
  52. #define II(a, b, c, d, x, s, ac) \
  53.   {(a) += I ((b), (c), (d)) + (x) + (UINT4)(ac); \
  54.    (a) = ROTATE_LEFT ((a), (s)); \
  55.    (a) += (b); \
  56.   }

  58. void MD5Init (MD5_CTX *mdContext)
  59. {
  60.         mdContext->i[0] = mdContext->i[1] = (UINT4)0;
  61.         /* Load magic initialization constants.
  62.         */
  63.         mdContext->buf[0] = (UINT4)0x67452301;
  64.         mdContext->buf[1] = (UINT4)0xefcdab89;
  65.         mdContext->buf[2] = (UINT4)0x98badcfe;
  66.         mdContext->buf[3] = (UINT4)0x10325476;
  67. }

  69. void MD5Update (MD5_CTX *mdContext, unsigned char *inBuf, unsigned int inLen)
  70. {
  71.         UINT4 in[16];
  72.         int mdi;
  73.         unsigned int i, ii;
  74.         /* compute number of bytes mod 64 */
  75.         mdi = (int)((mdContext->i[0] >> 3) & 0x3F);
  76.         /* update number of bits */
  77.         if ((mdContext->i[0] + ((UINT4)inLen << 3)) < mdContext->i[0])
  78.         {
  79.                 mdContext->i[1]++;
  80.         }
  81.         mdContext->i[0] += ((UINT4)inLen << 3);
  82.         mdContext->i[1] += ((UINT4)inLen >> 29);
  83.         while (inLen--)
  84.         {
  85.                 /* add new character to buffer, increment mdi */
  86.                 mdContext->in[mdi++] = *inBuf++;
  87.                 /* transform if necessary */
  88.                 if (mdi == 0x40)
  89.                 {
  90.                         for (i = 0, ii = 0; i < 16; i++, ii += 4)
  91.                         {
  92.                                 in[i] = (((UINT4)mdContext->in[ii+3]) << 24) |
  93.                                 (((UINT4)mdContext->in[ii+2]) << 16) |
  94.                                 (((UINT4)mdContext->in[ii+1]) << 8) |
  95.                                 ((UINT4)mdContext->in[ii]);
  96.                         }
  97.                         Transform (mdContext->buf, in);
  98.                         mdi = 0;
  99.                 } // end if

  101.         } // end while

  103. }

  105. void MD5Final (MD5_CTX *mdContext)
  106. {
  107.         UINT4 in[16];
  108.         int mdi;
  109.         unsigned int i, ii;
  110.         unsigned int padLen;
  111.         /* save number of bits */
  112.         in[14] = mdContext->i[0];
  113.         in[15] = mdContext->i[1];
  114.         /* compute number of bytes mod 64 */
  115.         mdi = (int)((mdContext->i[0] >> 3) & 0x3F);
  116.         /* pad out to 56 mod 64 */
  117.         padLen = (mdi < 56) ? (56 - mdi) : (120 - mdi);
  118.         MD5Update (mdContext, PADDING, padLen);
  119.         /* append length in bits and transform */
  120.         for (i = 0, ii = 0; i < 14; i++, ii += 4)
  121.         {
  122.                 in[i] = (((UINT4)mdContext->in[ii+3]) << 24) |
  123.                 (((UINT4)mdContext->in[ii+2]) << 16) |
  124.                 (((UINT4)mdContext->in[ii+1]) << 8) |
  125.                 ((UINT4)mdContext->in[ii]);
  126.         }
  127.         Transform (mdContext->buf, in);
  128.         /* store buffer in digest */
  129.         for (i = 0, ii = 0; i < 4; i++, ii += 4)
  130.         {
  131.                 mdContext->digest[ii] = (unsigned char)(mdContext->buf[i] & 0xFF);
  132.                 mdContext->digest[ii+1] = (unsigned char)((mdContext->buf[i] >> 8) & 0xFF);
  133.                 mdContext->digest[ii+2] = (unsigned char)((mdContext->buf[i] >> 16) & 0xFF);
  134.                 mdContext->digest[ii+3] = (unsigned char)((mdContext->buf[i] >> 24) & 0xFF);
  135.         } // end for

  137. }

  139. /* Basic MD5 step. Transform buf based on in.
  140.  */
  141. static void Transform (UINT4 *buf, UINT4 *in)
  142. {
  143.     UINT4 a = buf[0], b = buf[1], c = buf[2], d = buf[3];
  144.     /* Round 1 */
  145. #define S11 7
  146. #define S12 12
  147. #define S13 17
  148. #define S14 22
  149.     FF ( a, b, c, d, in[ 0], S11, 3614090360); /* 1 */
  150.     FF ( d, a, b, c, in[ 1], S12, 3905402710); /* 2 */
  151.     FF ( c, d, a, b, in[ 2], S13, 606105819); /* 3 */
  152.     FF ( b, c, d, a, in[ 3], S14, 3250441966); /* 4 */
  153.     FF ( a, b, c, d, in[ 4], S11, 4118548399); /* 5 */
  154.     FF ( d, a, b, c, in[ 5], S12, 1200080426); /* 6 */
  155.     FF ( c, d, a, b, in[ 6], S13, 2821735955); /* 7 */
  156.     FF ( b, c, d, a, in[ 7], S14, 4249261313); /* 8 */
  157.     FF ( a, b, c, d, in[ 8], S11, 1770035416); /* 9 */
  158.     FF ( d, a, b, c, in[ 9], S12, 2336552879); /* 10 */
  159.     FF ( c, d, a, b, in[10], S13, 4294925233); /* 11 */
  160.     FF ( b, c, d, a, in[11], S14, 2304563134); /* 12 */
  161.     FF ( a, b, c, d, in[12], S11, 1804603682); /* 13 */
  162.     FF ( d, a, b, c, in[13], S12, 4254626195); /* 14 */
  163.     FF ( c, d, a, b, in[14], S13, 2792965006); /* 15 */
  164.     FF ( b, c, d, a, in[15], S14, 1236535329); /* 16 */
  165.     /* Round 2 */
  166. #define S21 5
  167. #define S22 9
  168. #define S23 14
  169. #define S24 20
  170.     GG ( a, b, c, d, in[ 1], S21, 4129170786); /* 17 */
  171.     GG ( d, a, b, c, in[ 6], S22, 3225465664); /* 18 */
  172.     GG ( c, d, a, b, in[11], S23, 643717713); /* 19 */
  173.     GG ( b, c, d, a, in[ 0], S24, 3921069994); /* 20 */
  174.     GG ( a, b, c, d, in[ 5], S21, 3593408605); /* 21 */
  175.     GG ( d, a, b, c, in[10], S22, 38016083); /* 22 */
  176.     GG ( c, d, a, b, in[15], S23, 3634488961); /* 23 */
  177.     GG ( b, c, d, a, in[ 4], S24, 3889429448); /* 24 */
  178.     GG ( a, b, c, d, in[ 9], S21, 568446438); /* 25 */
  179.     GG ( d, a, b, c, in[14], S22, 3275163606); /* 26 */
  180.     GG ( c, d, a, b, in[ 3], S23, 4107603335); /* 27 */
  181.     GG ( b, c, d, a, in[ 8], S24, 1163531501); /* 28 */
  182.     GG ( a, b, c, d, in[13], S21, 2850285829); /* 29 */
  183.     GG ( d, a, b, c, in[ 2], S22, 4243563512); /* 30 */
  184.     GG ( c, d, a, b, in[ 7], S23, 1735328473); /* 31 */
  185.     GG ( b, c, d, a, in[12], S24, 2368359562); /* 32 */
  186.     /* Round 3 */
  187. #define S31 4
  188. #define S32 11
  189. #define S33 16
  190. #define S34 23
  191.     HH ( a, b, c, d, in[ 5], S31, 4294588738); /* 33 */
  192.     HH ( d, a, b, c, in[ 8], S32, 2272392833); /* 34 */
  193.     HH ( c, d, a, b, in[11], S33, 1839030562); /* 35 */
  194.     HH ( b, c, d, a, in[14], S34, 4259657740); /* 36 */
  195.     HH ( a, b, c, d, in[ 1], S31, 2763975236); /* 37 */
  196.     HH ( d, a, b, c, in[ 4], S32, 1272893353); /* 38 */
  197.     HH ( c, d, a, b, in[ 7], S33, 4139469664); /* 39 */
  198.     HH ( b, c, d, a, in[10], S34, 3200236656); /* 40 */
  199.     HH ( a, b, c, d, in[13], S31, 681279174); /* 41 */
  200.     HH ( d, a, b, c, in[ 0], S32, 3936430074); /* 42 */
  201.     HH ( c, d, a, b, in[ 3], S33, 3572445317); /* 43 */
  202.     HH ( b, c, d, a, in[ 6], S34, 76029189); /* 44 */
  203.     HH ( a, b, c, d, in[ 9], S31, 3654602809); /* 45 */
  204.     HH ( d, a, b, c, in[12], S32, 3873151461); /* 46 */
  205.     HH ( c, d, a, b, in[15], S33, 530742520); /* 47 */
  206.     HH ( b, c, d, a, in[ 2], S34, 3299628645); /* 48 */
  207.     /* Round 4 */
  208. #define S41 6
  209. #define S42 10
  210. #define S43 15
  211. #define S44 21
  212.     II ( a, b, c, d, in[ 0], S41, 4096336452); /* 49 */
  213.     II ( d, a, b, c, in[ 7], S42, 1126891415); /* 50 */
  214.     II ( c, d, a, b, in[14], S43, 2878612391); /* 51 */
  215.     II ( b, c, d, a, in[ 5], S44, 4237533241); /* 52 */
  216.     II ( a, b, c, d, in[12], S41, 1700485571); /* 53 */
  217.     II ( d, a, b, c, in[ 3], S42, 2399980690); /* 54 */
  218.     II ( c, d, a, b, in[10], S43, 4293915773); /* 55 */
  219.     II ( b, c, d, a, in[ 1], S44, 2240044497); /* 56 */
  220.     II ( a, b, c, d, in[ 8], S41, 1873313359); /* 57 */
  221.     II ( d, a, b, c, in[15], S42, 4264355552); /* 58 */
  222.     II ( c, d, a, b, in[ 6], S43, 2734768916); /* 59 */
  223.     II ( b, c, d, a, in[13], S44, 1309151649); /* 60 */
  224.     II ( a, b, c, d, in[ 4], S41, 4149444226); /* 61 */
  225.     II ( d, a, b, c, in[11], S42, 3174756917); /* 62 */
  226.     II ( c, d, a, b, in[ 2], S43, 718787259); /* 63 */
  227.     II ( b, c, d, a, in[ 9], S44, 3951481745); /* 64 */
  228.     buf[0] += a;
  229.     buf[1] += b;
  230.     buf[2] += c;
  231.     buf[3] += d;
  232. }

  234. #include <stdio.h>
  235. #include <sys/types.h>
  236. #include <time.h>
  237. #include <string.h>

  239. static void MDPrint (MD5_CTX *mdContext)
  240. {
  241.         int i;
  242.         for (i = 0; i < 16; i++)
  243.         {
  244.                 printf ("%02x", mdContext->digest[i]);
  245.         }
  246. }
  247. /* size of test block */
  248. #define TEST_BLOCK_SIZE 1000
  249. /* number of blocks to process */
  250. #define TEST_BLOCKS 10000
  251. /* number of test bytes = TEST_BLOCK_SIZE * TEST_BLOCKS */
  252. static long TEST_BYTES = (long)TEST_BLOCK_SIZE * (long)TEST_BLOCKS;
  253. /* A time trial routine, to measure the speed of MD5.
  254.    Measures wall time required to digest TEST_BLOCKS * TEST_BLOCK_SIZE
  255.    characters.
  256.  */
  257. static void MDTimeTrial ()
  258. {
  259.         MD5_CTX mdContext;
  260.         time_t endTime, startTime;
  261.         unsigned char data[TEST_BLOCK_SIZE];
  262.         unsigned int i;
  263.         /* initialize test data */
  264.         for (i = 0; i < TEST_BLOCK_SIZE; i++)
  265.         {
  266.                 data[i] = (unsigned char)(i & 0xFF);
  267.         }
  268.         /* start timer */
  269.         printf ("MD5 time trial. Processing %ld characters...\n", TEST_BYTES);
  270.         time (&startTime);
  271.         /* digest data in TEST_BLOCK_SIZE byte blocks */
  272.         MD5Init (&mdContext);
  273.         for (i = TEST_BLOCKS; i > 0; i--)
  274.         {
  275.                 MD5Update (&mdContext, data, TEST_BLOCK_SIZE);
  276.         }
  277.         MD5Final (&mdContext);
  278.         /* stop timer, get time difference */
  279.         time (&endTime);
  280.         MDPrint (&mdContext);
  281.         printf (" is digest of test input.\n");
  282.         printf("Seconds to process test input: %ld\n", (long)(endTime-startTime));
  283.         printf("Characters processed per second: %ld\n",
  284.         TEST_BYTES / (endTime - startTime));
  285. }
  286. /* Computes the message digest for string inString.
  287.    Prints out message digest, a space, the string (in quotes) and a
  288.    carriage return.
  289.  */
  290. static void MDString(char *inString)
  291. {
  292.         MD5_CTX mdContext;
  293.         unsigned int len = strlen (inString);
  294.         MD5Init (&mdContext);
  295.         MD5Update (&mdContext, inString, len);
  296.         MD5Final (&mdContext);
  297.         MDPrint (&mdContext);
  298.         printf (" \"%s\"\n\n", inString);
  299. }
  300. /* Computes the message digest for a specified file.
  301.    Prints out message digest, a space, the file name, and a carriage
  302.    return.
  303.  */
  304. static void MDFile (char *filename)
  305. {
  306.         FILE *inFile = fopen (filename, "rb");
  307.         MD5_CTX mdContext;
  308.         int bytes;
  309.         unsigned char data[1024];
  310.         if (inFile == NULL)
  311.         {
  312.                 printf ("%s can't be opened.\n", filename);
  313.                 return;
  314.         }
  315.         MD5Init (&mdContext);
  316.         while ((bytes = fread (data, 1, 1024, inFile)) != 0)
  317.         {
  318.                 MD5Update (&mdContext, data, bytes);
  319.         }
  320.         MD5Final (&mdContext);
  321.         MDPrint (&mdContext);
  322.         printf (" %s\n", filename);
  323.         fclose (inFile);
  324. }
  325. /* Writes the message digest of the data from stdin onto stdout,
  326.    followed by a carriage return.
  327.  */
  328. static void MDFilter ()
  329. {
  330.         MD5_CTX mdContext;
  331.         int bytes;
  332.         unsigned char data[16];
  333.         MD5Init (&mdContext);
  334.         while ((bytes = fread (data, 1, 16, stdin)) != 0)
  335.         {
  336.                 MD5Update (&mdContext, data, bytes);
  337.         }
  338.         MD5Final (&mdContext);
  339.         MDPrint (&mdContext);
  340.         printf ("\n");
  341. }
  342. /* Runs a standard suite of test data.
  343.  */
  344. static void MDTestSuite ()
  345. {
  346.         printf ("MD5 test suite results:\n\n");
  347.         MDString ("");
  348.         MDString ("a");
  349.         MDString ("abc");
  350.         MDString ("message digest");
  351.         MDString ("abcdefghijklmnopqrstuvwxyz");
  352.         MDString("ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwxyz0123456789");
  353.         MDString("12345678901234567890123456789012345678901234567890123456789012345678901234567890");
  354.         /* Contents of file foo are "abc" */
  355.         MDFile ("foo");
  356. }
  357. void main (int argc, char *argv[])
  358. {
  359.         int i;
  360.         /* For each command line argument in turn:
  361.         ** filename -- prints message digest and name of file
  362.         ** -sstring -- prints message digest and contents of string 打印消息摘要和字符串内容
  363.         ** -t -- prints time trial statistics for 1M characters
  364.         ** -x -- execute a standard suite of test data
  365.         ** (no args) -- writes messages digest of stdin onto stdout
  366.         */
  367.         if (argc == 1)
  368.         {
  369.                 MDFilter ();
  370.         }
  371.         else
  372.         {
  373.                 for (i = 1; i < argc; i++)
  374.                 {
  375.                         if (argv[i][0] == '-' && argv[i][1] == 's')
  376.                                 MDString (argv[i] + 2);
  377.                         else if (strcmp (argv[i], "-t") == 0)
  378.                                 MDTimeTrial ();
  379.                         else if (strcmp (argv[i], "-x") == 0)
  380.                                 MDTestSuite ();
  381.                         else MDFile (argv[i]);
  382.                 } // end for

  384.         } // end else

  386. }
阅读(1822) | 评论(0) | 转发(0) |
0

上一篇：crc16

下一篇：弹钢琴

给主人留下些什么吧！~~
关于我们 | 关于IT168 | 联系方式 | 广告合作 | 法律声明 | 免费注册

Copyright 2001-2010 ChinaUnix.net All Rights Reserved 北京皓辰网域网络信息技术有限公司. 版权所有

感谢所有关心和支持过ChinaUnix的朋友们

16024965号-6