HASH算法，sha系列的基本写完

sha1 ,见我之前的一篇博客 http://andydhu.blog.51cto.com/3337368/824735

sha256

 
  
  1. #include<stdio.h> 
  2. #include<stdlib.h> 
  3. #include<time.h> 
  4. //#define DEBUG 1 
  5. #define RoundNum 64 
  6. #define BlockWidthNum 16 
  7. #define ConstNum 8 
  8. unsigned int H[ConstNum] = {0x6a09e667, 0xbb67ae85, 0x3c6ef372, 0xa54ff53a, 0x510e527f, 0x9b05688c, 0x1f83d9ab, 0x5be0cd19}; 
  9. unsigned int K[RoundNum] = { 
  10. 0x428a2f98, 0x71374491, 0xb5c0fbcf, 0xe9b5dba5, 0x3956c25b, 0x59f111f1, 0x923f82a4, 0xab1c5ed5, 
  11. 0xd807aa98, 0x12835b01, 0x243185be, 0x550c7dc3, 0x72be5d74, 0x80deb1fe, 0x9bdc06a7, 0xc19bf174, 
  12. 0xe49b69c1, 0xefbe4786, 0x0fc19dc6, 0x240ca1cc, 0x2de92c6f, 0x4a7484aa, 0x5cb0a9dc, 0x76f988da, 
  13. 0x983e5152, 0xa831c66d, 0xb00327c8, 0xbf597fc7, 0xc6e00bf3, 0xd5a79147, 0x06ca6351, 0x14292967, 
  14. 0x27b70a85, 0x2e1b2138, 0x4d2c6dfc, 0x53380d13, 0x650a7354, 0x766a0abb, 0x81c2c92e, 0x92722c85, 
  15. 0xa2bfe8a1, 0xa81a664b, 0xc24b8b70, 0xc76c51a3, 0xd192e819, 0xd6990624, 0xf40e3585, 0x106aa070, 
  16. 0x19a4c116, 0x1e376c08, 0x2748774c, 0x34b0bcb5, 0x391c0cb3, 0x4ed8aa4a, 0x5b9cca4f, 0x682e6ff3, 
  17. 0x748f82ee, 0x78a5636f, 0x84c87814, 0x8cc70208, 0x90befffa, 0xa4506ceb, 0xbef9a3f7, 0xc67178f2}; 
  18.  
  19. unsigned int M[BlockWidthNum]; 
  20. unsigned int W[RoundNum]; 
  21. unsigned char LastString[BlockWidthNum*4*2]; 
  22. FILE * fp; 
  23. fpos_t pos; //文件长度 
  24. int flag ;  
  25. int LastRoundCount; 
  26.  
  27. int GetMessage(unsigned long long rLen); 
  28. void SetW(); 
  29. void Round(); 
  30. unsigned int ShiftNum3(unsigned int input, int OpOne, int OpTwo, int Opthree); 
  31.  
  32. int main(int argc, char *argv[]) 
  33. { 
  34. //  fp = fopen(argv[1],"rb"); 
  35.     if(!(fp = fopen(argv[1],"rb"))) 
  36.     { 
  37.         exit(1); 
  38.     } 
  39.     clock_t start, finish;   
  40.     start = clock();   
  41.      
  42.     fseek(fp, 0, SEEK_END); 
  43.     fgetpos(fp, &pos); 
  44.  
  45.     unsigned long long realLen = pos; 
  46.  
  47.     memset(LastString,0x00,strlen(LastString)*sizeof(unsigned char)); 
  48.      
  49.     int LastOut = realLen%(BlockWidthNum*4); 
  50.     fseek(fp,-1*LastOut,SEEK_END); 
  51.  
  52.     int i; 
  53.     fread(LastString,LastOut,1,fp); 
  54.     LastString[LastOut] = 0x80;  
  55.      
  56.     LastRoundCount = 1; 
  57.     flag = 0; 
  58.  
  59.     if( LastOut <= BlockWidthNum*4-1-8) 
  60.     { 
  61.         flag = 1; 
  62.     //  memcpy(LastString+512/8-8,&copyLen,8); 
  63.         for( i = 8; i >= 1; i--) 
  64.             LastString[512/8 - i] = (char)((realLen<<3)>>((i-1)*8));  
  65.     } 
  66.     else 
  67.     { 
  68.         flag = 2; 
  69.     //  memcpy(LastString+512/8*2-8,&copyLen,8); 
  70.         for( i = 1; i <= 8; i++) 
  71.             LastString[512/8*2 - i] = (char)((realLen<<3)>>((i-1)*8));  
  72.     } 
  73.  
  74.     realLen -= LastOut; 
  75.  
  76.     fseek(fp,0,SEEK_SET); 
  77.      
  78.     while(GetMessage(realLen)) 
  79.     { 
  80.         SetW(); 
  81.         Round(); 
  82.     } 
  83.      
  84.     printf("SHA256:  "); 
  85.     for(i = 0 ; i < ConstNum; i++) 
  86.         printf("%08x",H[i]); 
  87.     printf("\n");  
  88.      
  89.     finish = clock();   
  90.     double duration = (double)(finish - start) / CLOCKS_PER_SEC;  
  91.     printf("Cal time: %f seconds\n",duration); 
  92.     fclose(fp); 
  93.     return 0; 
  94. } 
  95.  
  96. int GetMessage(unsigned long long rLen) 
  97. { 
  98.     fgetpos(fp, &pos); 
  99.  
  100.     unsigned char temp[4]; 
  101.     if(pos + 16*4 <= rLen ) 
  102.     { 
  103.          
  104.         fread(M,4,16,fp); 
  105.         return 1; 
  106.     } 
  107.     if(LastRoundCount <= flag) 
  108.     { 
  109.         memcpy(M,(LastString+512/8*(LastRoundCount - 1)),64);//最后的一段或是两段 
  110.         LastRoundCount++; 
  111.         return 1; 
  112.     } 
  113.     else 
  114.         return 0; //false;   
  115. } 
  116. void SetW() 
  117. { 
  118.     int index = 0; 
  119. //  memcpy(W,M,16*4); 
  120.  
  121.     for(index = 0; index <= 15;index++) 
  122.     { 
  123.         W[index] = (M[index]>>24)|(M[index]<<24)|((M[index]&0x00ff0000)>>8)|((M[index]&0x0000ff00)<<8); 
  124.     }  
  125.     unsigned int s0 = 0, s1 = 0; 
  126.     for(index = 16; index < RoundNum; index++) 
  127.     { 
  128.         s0 = 0; 
  129.         s1 = 0; 
  130.         s0 = ShiftNum3(W[index-2],17,19,-10); 
  131.         s1 = ShiftNum3(W[index-15],7,18,-3); 
  132.         W[index] = (s0 + W[index-7] + s1 + W[index-16]); 
  133. //      printf(" %08x ", W[index]); 
  134.     } 
  135. } 
  136.  
  137. unsigned int ShiftNum3(unsigned int input, int OpOne, int OpTwo, int Opthree) 
  138. { 
  139.     unsigned int Ans = 0 ;  
  140.     Ans ^= ((input >> OpOne)|(input<<(32-OpOne))); 
  141.     Ans ^= ((input >> OpTwo)|(input<<(32-OpTwo))); 
  142.     if(Opthree > 0) 
  143.         Ans ^= ((input >> Opthree)|(input<<(32-Opthree))); 
  144.     else 
  145.     { 
  146.     //  printf("%d",Opthree); 
  147.         Opthree = 0 - Opthree; 
  148.     //  printf("%d",Opthree); 
  149.         Ans ^= (input >> Opthree); 
  150.     } 
  151.     return Ans; 
  152. } 
  153.  
  154. void Round() 
  155. { 
  156.     unsigned int A,B,C,D,E,F,G,I; 
  157.     A = H[0]; 
  158.     B = H[1]; 
  159.     C = H[2]; 
  160.     D = H[3]; 
  161.     E = H[4]; 
  162.     F = H[5]; 
  163.     G = H[6]; 
  164.     I = H[7]; 
  165.      
  166.     int index = 0; 
  167.     unsigned int t1 = 0, t2= 0, s0 = 0, s1 = 0; 
  168.     for(index = 0;index < RoundNum ;index++) 
  169.     { 
  170.         s0 = 0; 
  171.         s1 = 0; 
  172.         s0 = ShiftNum3(A, 2, 13, 22); 
  173.         t2 = s0 + ((A&B)^(A&C)^(B&C)); 
  174.         s1 = ShiftNum3(E, 6, 11, 25); 
  175.         t1 = I + s1 + ((E&F)^((~E)&G)) + K[index] + W[index]; 
  176.          
  177.         I = G; 
  178.         G = F; 
  179.         F = E; 
  180.         E = D + t1; 
  181.         D = C; 
  182.         C = B; 
  183.         B = A; 
  184.         A = t1 + t2; 
  185. #ifdef DEBUG 
  186.     printf("ABCDE  %x %x %x %x %x \n",A,B,C,D,E); 
  187.     system("pause"); 
  188. #endif 
  189.     } 
  190.  
  191.     H[0]+=A; 
  192.     H[1]+=B; 
  193.     H[2]+=C; 
  194.     H[3]+=D; 
  195.     H[4]+=E; 
  196.     H[5]+=F; 
  197.     H[6]+=G; 
  198.     H[7]+=I; 
  199. } 
 
  
sha384
 
  
 
  
  1. #include<stdio.h> 
  2. #include<stdlib.h> 
  3. #include<time.h> 
  4. //#define DEBUG 1 
  5. #define RoundNum 80 
  6. #define BlockWidthNum 16 
  7. #define ConstNum 8 
  8. unsigned long long H[ConstNum] = { 
  9. 0xcbbb9d5dc1059ed8, 0x629a292a367cd507, 0x9159015a3070dd17, 0x152fecd8f70e5939, 
  10. 0x67332667ffc00b31, 0x8eb44a8768581511, 0xdb0c2e0d64f98fa7, 0x47b5481dbefa4fa4}; 
  11. unsigned long long K[RoundNum] = { 
  12. 0x428a2f98d728ae22, 0x7137449123ef65cd, 0xb5c0fbcfec4d3b2f, 0xe9b5dba58189dbbc,  
  13. 0x3956c25bf348b538, 0x59f111f1b605d019, 0x923f82a4af194f9b, 0xab1c5ed5da6d8118,  
  14. 0xd807aa98a3030242, 0x12835b0145706fbe, 0x243185be4ee4b28c, 0x550c7dc3d5ffb4e2,  
  15. 0x72be5d74f27b896f, 0x80deb1fe3b1696b1, 0x9bdc06a725c71235, 0xc19bf174cf692694,  
  16. 0xe49b69c19ef14ad2, 0xefbe4786384f25e3, 0x0fc19dc68b8cd5b5, 0x240ca1cc77ac9c65,  
  17. 0x2de92c6f592b0275, 0x4a7484aa6ea6e483, 0x5cb0a9dcbd41fbd4, 0x76f988da831153b5,  
  18. 0x983e5152ee66dfab, 0xa831c66d2db43210, 0xb00327c898fb213f, 0xbf597fc7beef0ee4,  
  19. 0xc6e00bf33da88fc2, 0xd5a79147930aa725, 0x06ca6351e003826f, 0x142929670a0e6e70,  
  20. 0x27b70a8546d22ffc, 0x2e1b21385c26c926, 0x4d2c6dfc5ac42aed, 0x53380d139d95b3df,  
  21. 0x650a73548baf63de, 0x766a0abb3c77b2a8, 0x81c2c92e47edaee6, 0x92722c851482353b,  
  22. 0xa2bfe8a14cf10364, 0xa81a664bbc423001, 0xc24b8b70d0f89791, 0xc76c51a30654be30,  
  23. 0xd192e819d6ef5218, 0xd69906245565a910, 0xf40e35855771202a, 0x106aa07032bbd1b8,  
  24. 0x19a4c116b8d2d0c8, 0x1e376c085141ab53, 0x2748774cdf8eeb99, 0x34b0bcb5e19b48a8,  
  25. 0x391c0cb3c5c95a63, 0x4ed8aa4ae3418acb, 0x5b9cca4f7763e373, 0x682e6ff3d6b2b8a3,  
  26. 0x748f82ee5defb2fc, 0x78a5636f43172f60, 0x84c87814a1f0ab72, 0x8cc702081a6439ec,  
  27. 0x90befffa23631e28, 0xa4506cebde82bde9, 0xbef9a3f7b2c67915, 0xc67178f2e372532b,  
  28. 0xca273eceea26619c, 0xd186b8c721c0c207, 0xeada7dd6cde0eb1e, 0xf57d4f7fee6ed178,  
  29. 0x06f067aa72176fba, 0x0a637dc5a2c898a6, 0x113f9804bef90dae, 0x1b710b35131c471b,  
  30. 0x28db77f523047d84, 0x32caab7b40c72493, 0x3c9ebe0a15c9bebc, 0x431d67c49c100d4c,  
  31. 0x4cc5d4becb3e42b6, 0x597f299cfc657e2a, 0x5fcb6fab3ad6faec, 0x6c44198c4a475817}; 
  32.  
  33. unsigned long long M[BlockWidthNum]; 
  34. unsigned long long W[RoundNum]; 
  35. unsigned char LastString[BlockWidthNum*8*2]; 
  36. FILE * fp; 
  37. fpos_t pos; //文件长度 
  38. int flag ;  
  39. int LastRoundCount; 
  40.  
  41. int GetMessage(unsigned long long rLen); 
  42. void SetW(); 
  43. void Round(); 
  44. unsigned long long ShiftNum3(unsigned long long input, int OpOne, int OpTwo, int Opthree); 
  45.  
  46. int main(int argc, char *argv[]) 
  47. { 
  48. //  fp = fopen(argv[1],"rb"); 
  49.     if(!(fp = fopen(argv[1],"rb"))) 
  50.     { 
  51.         exit(1); 
  52.     } 
  53.     clock_t start, finish;   
  54.     start = clock();   
  55.      
  56.     fseek(fp, 0, SEEK_END); 
  57.     fgetpos(fp, &pos); 
  58.  
  59.     unsigned long long realLen = pos;//暂时不支持hash超大文件 
  60.  
  61.     memset(LastString,0x00,strlen(LastString)*sizeof(unsigned char)); 
  62.      
  63.     int LastOut = realLen%(BlockWidthNum*8); 
  64.     fseek(fp,-1*LastOut,SEEK_END); 
  65.  
  66.     int i; 
  67.     fread(LastString,LastOut,1,fp); 
  68.     LastString[LastOut] = 0x80;  
  69.      
  70.     LastRoundCount = 1; 
  71.     flag = 0; 
  72. //暂时不支持hash大小超过2^64的文件， 
  73. //sha512更长了，这样更加安全了，但是也不一定要hash这种超大的文件 
  74.     if( LastOut <= BlockWidthNum*8-1-16) 
  75.     { 
  76.         flag = 1; 
  77.     //  memcpy(LastString+512/8-8,&copyLen,8); 
  78.         for( i = 8; i >= 1; i--) 
  79.             LastString[BlockWidthNum*8 - i] = (char)((realLen<<3)>>((i-1)*8));  
  80.     } 
  81.     else 
  82.     { 
  83.         flag = 2; 
  84.     //  memcpy(LastString+512/8*2-8,&copyLen,8); 
  85.         for( i = 1; i <= 8; i++) 
  86.             LastString[BlockWidthNum*8*2 - i] = (char)((realLen<<3)>>((i-1)*8));  
  87.     } 
  88.  
  89.     realLen -= LastOut; 
  90.  
  91.     fseek(fp,0,SEEK_SET); 
  92.      
  93.     while(GetMessage(realLen)) 
  94.     { 
  95.         SetW(); 
  96.         Round(); 
  97.     } 
  98.      
  99.     printf("SHA384:  "); 
  100.     for(i = 0 ; i < 6; i++) 
  101.         printf("%016llx",H[i]); 
  102.     printf("\n");  
  103.      
  104.     finish = clock();   
  105.     double duration = (double)(finish - start) / CLOCKS_PER_SEC;  
  106.     printf("Cal time: %f seconds\n",duration); 
  107.     fclose(fp); 
  108.     return 0; 
  109. } 
  110.  
  111. int GetMessage(unsigned long long rLen) 
  112. { 
  113.     fgetpos(fp, &pos); 
  114.  
  115.     if(pos + 16 * 8 <= rLen ) 
  116.     { 
  117.         fread(M,8,16,fp); 
  118.         return 1; 
  119.     } 
  120.     if(LastRoundCount <= flag) 
  121.     { 
  122.         memcpy(M,(LastString+BlockWidthNum*8*(LastRoundCount - 1)),BlockWidthNum*8);//最后的一段或是两段 
  123.         LastRoundCount++; 
  124.         return 1; 
  125.     } 
  126.     else 
  127.         return 0; //false;   
  128. } 
  129. void SetW() 
  130. { 
  131.     int index = 0; 
  132. //  memcpy(W,M,16*4); 
  133.  
  134.     for(index = 0; index <= 15;index++) 
  135.     { 
  136.         W[index] = (M[index]>>56)|(M[index]<<56)|((M[index]&0x00ff000000000000)>>40)|((M[index]&0x000000000000ff00)<<40) 
  137.                     |((M[index]&0x0000ff0000000000)>>24)|((M[index]&0x0000000000ff0000)<<24)|((M[index]&0x000000ff00000000)>>8)|((M[index]&0x00000000ff000000)<<8); 
  138. //      printf("%016llx ",W[index]); 
  139.     }  
  140.     unsigned long long s0 = 0, s1 = 0; 
  141.     for(index = 16; index < RoundNum; index++) 
  142.     { 
  143.         s0 = 0; 
  144.         s1 = 0; 
  145.         s0 = ShiftNum3(W[index-2],19,61,-6); 
  146.         s1 = ShiftNum3(W[index-15],1,8,-7); 
  147.         W[index] = (s0 + W[index-7] + s1 + W[index-16]); 
  148. //      printf("%016llx ",W[index]); 
  149.     } 
  150. } 
  151.  
  152. unsigned long long ShiftNum3(unsigned long long input, int OpOne, int OpTwo, int Opthree) 
  153. { 
  154.     unsigned long long Ans = 0 ;  
  155.     Ans ^= ((input >> OpOne)|(input<<(64-OpOne))); 
  156.     Ans ^= ((input >> OpTwo)|(input<<(64-OpTwo))); 
  157.     if(Opthree > 0) 
  158.         Ans ^= ((input >> Opthree)|(input<<(64-Opthree))); 
  159.     else 
  160.     { 
  161.     //  printf("%d",Opthree); 
  162.         Opthree = 0 - Opthree; 
  163.     //  printf("%d",Opthree); 
  164.         Ans ^= (input >> Opthree); 
  165.     } 
  166.     return Ans; 
  167. } 
  168.  
  169. void Round() 
  170. { 
  171.     unsigned long long A,B,C,D,E,F,G,I; 
  172.     A = H[0]; 
  173.     B = H[1]; 
  174.     C = H[2]; 
  175.     D = H[3]; 
  176.     E = H[4]; 
  177.     F = H[5]; 
  178.     G = H[6]; 
  179.     I = H[7]; 
  180.      
  181.     int index = 0; 
  182.     unsigned long long t1 = 0, t2= 0, s0 = 0, s1 = 0; 
  183.     for(index = 0;index < RoundNum ;index++) 
  184.     { 
  185.         s0 = 0; 
  186.         s1 = 0; 
  187.         s0 = ShiftNum3(A, 28, 34, 39); 
  188.         t2 = s0 + ((A&B)^(A&C)^(B&C)); 
  189.         s1 = ShiftNum3(E, 14, 18, 41); 
  190.         t1 = I + s1 + ((E&F)^((~E)&G)) + K[index] + W[index]; 
  191.          
  192.         I = G; 
  193.         G = F; 
  194.         F = E; 
  195.         E = D + t1; 
  196.         D = C; 
  197.         C = B; 
  198.         B = A; 
  199.         A = t1 + t2; 
  200. #ifdef DEBUG 
  201.     printf("ABCDE  %x %x %x %x %x \n",A,B,C,D,E); 
  202.     system("pause"); 
  203. #endif 
  204.     } 
  205.  
  206.     H[0]+=A; 
  207.     H[1]+=B; 
  208.     H[2]+=C; 
  209.     H[3]+=D; 
  210.     H[4]+=E; 
  211.     H[5]+=F; 
  212.     H[6]+=G; 
  213.     H[7]+=I; 
  214. } 
 
  
sha512
 
  
 
  
  1. #include<stdio.h> 
  2. #include<stdlib.h> 
  3. #include<time.h> 
  4. //#define DEBUG 1 
  5. #define RoundNum 80 
  6. #define BlockWidthNum 16 
  7. #define ConstNum 8 
  8. unsigned long long H[ConstNum] = { 
  9. 0x6a09e667f3bcc908, 0xbb67ae8584caa73b, 0x3c6ef372fe94f82b, 0xa54ff53a5f1d36f1,  
  10. 0x510e527fade682d1, 0x9b05688c2b3e6c1f, 0x1f83d9abfb41bd6b, 0x5be0cd19137e2179}; 
  11. unsigned long long K[RoundNum] = { 
  12. 0x428a2f98d728ae22, 0x7137449123ef65cd, 0xb5c0fbcfec4d3b2f, 0xe9b5dba58189dbbc,  
  13. 0x3956c25bf348b538, 0x59f111f1b605d019, 0x923f82a4af194f9b, 0xab1c5ed5da6d8118,  
  14. 0xd807aa98a3030242, 0x12835b0145706fbe, 0x243185be4ee4b28c, 0x550c7dc3d5ffb4e2,  
  15. 0x72be5d74f27b896f, 0x80deb1fe3b1696b1, 0x9bdc06a725c71235, 0xc19bf174cf692694,  
  16. 0xe49b69c19ef14ad2, 0xefbe4786384f25e3, 0x0fc19dc68b8cd5b5, 0x240ca1cc77ac9c65,  
  17. 0x2de92c6f592b0275, 0x4a7484aa6ea6e483, 0x5cb0a9dcbd41fbd4, 0x76f988da831153b5,  
  18. 0x983e5152ee66dfab, 0xa831c66d2db43210, 0xb00327c898fb213f, 0xbf597fc7beef0ee4,  
  19. 0xc6e00bf33da88fc2, 0xd5a79147930aa725, 0x06ca6351e003826f, 0x142929670a0e6e70,  
  20. 0x27b70a8546d22ffc, 0x2e1b21385c26c926, 0x4d2c6dfc5ac42aed, 0x53380d139d95b3df,  
  21. 0x650a73548baf63de, 0x766a0abb3c77b2a8, 0x81c2c92e47edaee6, 0x92722c851482353b,  
  22. 0xa2bfe8a14cf10364, 0xa81a664bbc423001, 0xc24b8b70d0f89791, 0xc76c51a30654be30,  
  23. 0xd192e819d6ef5218, 0xd69906245565a910, 0xf40e35855771202a, 0x106aa07032bbd1b8,  
  24. 0x19a4c116b8d2d0c8, 0x1e376c085141ab53, 0x2748774cdf8eeb99, 0x34b0bcb5e19b48a8,  
  25. 0x391c0cb3c5c95a63, 0x4ed8aa4ae3418acb, 0x5b9cca4f7763e373, 0x682e6ff3d6b2b8a3,  
  26. 0x748f82ee5defb2fc, 0x78a5636f43172f60, 0x84c87814a1f0ab72, 0x8cc702081a6439ec,  
  27. 0x90befffa23631e28, 0xa4506cebde82bde9, 0xbef9a3f7b2c67915, 0xc67178f2e372532b,  
  28. 0xca273eceea26619c, 0xd186b8c721c0c207, 0xeada7dd6cde0eb1e, 0xf57d4f7fee6ed178,  
  29. 0x06f067aa72176fba, 0x0a637dc5a2c898a6, 0x113f9804bef90dae, 0x1b710b35131c471b,  
  30. 0x28db77f523047d84, 0x32caab7b40c72493, 0x3c9ebe0a15c9bebc, 0x431d67c49c100d4c,  
  31. 0x4cc5d4becb3e42b6, 0x597f299cfc657e2a, 0x5fcb6fab3ad6faec, 0x6c44198c4a475817}; 
  32.  
  33. unsigned long long M[BlockWidthNum]; 
  34. unsigned long long W[RoundNum]; 
  35. unsigned char LastString[BlockWidthNum*8*2]; 
  36. FILE * fp; 
  37. fpos_t pos; //文件长度 
  38. int flag ;  
  39. int LastRoundCount; 
  40.  
  41. int GetMessage(unsigned long long rLen); 
  42. void SetW(); 
  43. void Round(); 
  44. unsigned long long ShiftNum3(unsigned long long input, int OpOne, int OpTwo, int Opthree); 
  45.  
  46. int main(int argc, char *argv[]) 
  47. { 
  48. //  fp = fopen(argv[1],"rb"); 
  49.     if(!(fp = fopen(argv[1],"rb"))) 
  50.     { 
  51.         exit(1); 
  52.     } 
  53.     clock_t start, finish;   
  54.     start = clock();   
  55.      
  56.     fseek(fp, 0, SEEK_END); 
  57.     fgetpos(fp, &pos); 
  58.  
  59.     unsigned long long realLen = pos;//暂时不支持hash超大文件 
  60.  
  61.     memset(LastString,0x00,strlen(LastString)*sizeof(unsigned char)); 
  62.      
  63.     int LastOut = realLen%(BlockWidthNum*8); 
  64.     fseek(fp,-1*LastOut,SEEK_END); 
  65.  
  66.     int i; 
  67.     fread(LastString,LastOut,1,fp); 
  68.     LastString[LastOut] = 0x80;  
  69.      
  70.     LastRoundCount = 1; 
  71.     flag = 0; 
  72. //暂时不支持hash大小超过2^64的文件， 
  73. //sha512更长了，这样更加安全了，但是也不一定要hash这种超大的文件 
  74.     if( LastOut <= BlockWidthNum*8-1-16) 
  75.     { 
  76.         flag = 1; 
  77.     //  memcpy(LastString+512/8-8,&copyLen,8); 
  78.         for( i = 8; i >= 1; i--) 
  79.             LastString[BlockWidthNum*8 - i] = (char)((realLen<<3)>>((i-1)*8));  
  80.     } 
  81.     else 
  82.     { 
  83.         flag = 2; 
  84.     //  memcpy(LastString+512/8*2-8,&copyLen,8); 
  85.         for( i = 1; i <= 8; i++) 
  86.             LastString[BlockWidthNum*8*2 - i] = (char)((realLen<<3)>>((i-1)*8));  
  87.     } 
  88.  
  89.     realLen -= LastOut; 
  90.  
  91.     fseek(fp,0,SEEK_SET); 
  92.      
  93.     while(GetMessage(realLen)) 
  94.     { 
  95.         SetW(); 
  96.         Round(); 
  97.     } 
  98.      
  99.     printf("SHA512:  "); 
  100.     for(i = 0 ; i < ConstNum; i++) 
  101.         printf("%016llx",H[i]); 
  102.     printf("\n");  
  103.      
  104.     finish = clock();   
  105.     double duration = (double)(finish - start) / CLOCKS_PER_SEC;  
  106.     printf("Cal time: %f seconds\n",duration); 
  107.     fclose(fp); 
  108.     return 0; 
  109. } 
  110.  
  111. int GetMessage(unsigned long long rLen) 
  112. { 
  113.     fgetpos(fp, &pos); 
  114.  
  115.     if(pos + 16 * 8 <= rLen ) 
  116.     { 
  117.         fread(M,8,16,fp); 
  118.         return 1; 
  119.     } 
  120.     if(LastRoundCount <= flag) 
  121.     { 
  122.         memcpy(M,(LastString+BlockWidthNum*8*(LastRoundCount - 1)),BlockWidthNum*8);//最后的一段或是两段 
  123.         LastRoundCount++; 
  124.         return 1; 
  125.     } 
  126.     else 
  127.         return 0; //false;   
  128. } 
  129. void SetW() 
  130. { 
  131.     int index = 0; 
  132. //  memcpy(W,M,16*4); 
  133.  
  134.     for(index = 0; index <= 15;index++) 
  135.     { 
  136.         W[index] = (M[index]>>56)|(M[index]<<56)|((M[index]&0x00ff000000000000)>>40)|((M[index]&0x000000000000ff00)<<40) 
  137.                     |((M[index]&0x0000ff0000000000)>>24)|((M[index]&0x0000000000ff0000)<<24)|((M[index]&0x000000ff00000000)>>8)|((M[index]&0x00000000ff000000)<<8); 
  138. //      printf("%016llx ",W[index]); 
  139.     }  
  140.     unsigned long long s0 = 0, s1 = 0; 
  141.     for(index = 16; index < RoundNum; index++) 
  142.     { 
  143.         s0 = 0; 
  144.         s1 = 0; 
  145.         s0 = ShiftNum3(W[index-2],19,61,-6); 
  146.         s1 = ShiftNum3(W[index-15],1,8,-7); 
  147.         W[index] = (s0 + W[index-7] + s1 + W[index-16]); 
  148. //      printf("%016llx ",W[index]); 
  149.     } 
  150. } 
  151.  
  152. unsigned long long ShiftNum3(unsigned long long input, int OpOne, int OpTwo, int Opthree) 
  153. { 
  154.     unsigned long long Ans = 0 ;  
  155.     Ans ^= ((input >> OpOne)|(input<<(64-OpOne))); 
  156.     Ans ^= ((input >> OpTwo)|(input<<(64-OpTwo))); 
  157.     if(Opthree > 0) 
  158.         Ans ^= ((input >> Opthree)|(input<<(64-Opthree))); 
  159.     else 
  160.     { 
  161.     //  printf("%d",Opthree); 
  162.         Opthree = 0 - Opthree; 
  163.     //  printf("%d",Opthree); 
  164.         Ans ^= (input >> Opthree); 
  165.     } 
  166.     return Ans; 
  167. } 
  168.  
  169. void Round() 
  170. { 
  171.     unsigned long long A,B,C,D,E,F,G,I; 
  172.     A = H[0]; 
  173.     B = H[1]; 
  174.     C = H[2]; 
  175.     D = H[3]; 
  176.     E = H[4]; 
  177.     F = H[5]; 
  178.     G = H[6]; 
  179.     I = H[7]; 
  180.      
  181.     int index = 0; 
  182.     unsigned long long t1 = 0, t2= 0, s0 = 0, s1 = 0; 
  183.     for(index = 0;index < RoundNum ;index++) 
  184.     { 
  185.         s0 = 0; 
  186.         s1 = 0; 
  187.         s0 = ShiftNum3(A, 28, 34, 39); 
  188.         t2 = s0 + ((A&B)^(A&C)^(B&C)); 
  189.         s1 = ShiftNum3(E, 14, 18, 41); 
  190.         t1 = I + s1 + ((E&F)^((~E)&G)) + K[index] + W[index]; 
  191.          
  192.         I = G; 
  193.         G = F; 
  194.         F = E; 
  195.         E = D + t1; 
  196.         D = C; 
  197.         C = B; 
  198.         B = A; 
  199.         A = t1 + t2; 
  200. #ifdef DEBUG 
  201.     printf("ABCDE  %x %x %x %x %x \n",A,B,C,D,E); 
  202.     system("pause"); 
  203. #endif 
  204.     } 
  205.  
  206.     H[0]+=A; 
  207.     H[1]+=B; 
  208.     H[2]+=C; 
  209.     H[3]+=D; 
  210.     H[4]+=E; 
  211.     H[5]+=F; 
  212.     H[6]+=G; 
  213.     H[7]+=I; 
  214. }