I tried, but I can not spot my mistake. Maybe you guys can find it. I will have to try it some more tomorrow.
Code: Select all
#include <stdlib.h>
#include <stdio.h>
#include <stdint.h>
#include <string.h>
#include <math.h>
#define ROTATE_LEFT(x, n) (((x) << (n)) | ((x) >> (32-(n))))
#define F(x,y,z) (((x) & (y)) | ((~(x)) & (z)))
#define G(x,y,z) (((x) & (z)) | ((y) & (~(z))))
#define H(x,y,z) ((x) ^ (y) ^ (z))
#define I(x,y,z) ((y) ^ ((x) | (~(z))))
void md5_genT(uint32_t *T)
{
uint_fast32_t i;
for(i = 0; i < 64; ++i)
{
T[i] = 4294967296.0 * fabs(sin((float)(i+1)));
//printf("%08x\n", T[i]);
//sleep(1);
}
}
void md5_hash(uint8_t *message, uint_fast32_t mlength, uint32_t *T)
{
uint8_t *nm, nmlength;
uint32_t __s1[] = {7, 12, 17, 22};
uint32_t __s2[] = {5, 9, 14, 20};
uint32_t __s3[] = {4, 11, 16, 23};
uint32_t __s4[] = {6, 10, 15, 21};
uint32_t __k1[] = {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15};
uint32_t __k2[] = {1, 6, 11, 0, 5, 10, 15, 4, 9, 14, 3, 8, 13, 2, 7, 12};
uint32_t __k3[] = {5, 8, 11, 14, 1, 4, 7, 10, 13, 0, 3, 6, 9, 12, 15, 2};
uint32_t __k4[] = {0, 7, 14, 5, 12, 3, 10, 1, 8, 15, 6, 13, 4, 11, 2, 9};
uint32_t X[16];
uint32_t i, j, ii, k;
uint32_t AA, BB, CC, DD;
uint32_t A, B, C, D;
// do padding so message is a multiple of 512 bits (64 bytes)
if((mlength%64) == 0)
{
nmlength = mlength + 64;
}else{
nmlength = mlength + (64 - (mlength%64));
}
nm = (uint8_t*)malloc(nmlength);
memset(nm, 0, nmlength);
memcpy(nm, message, mlength);
printf("nmlength:%u mlength:%u\n", nmlength, mlength);
nm[mlength] = 0x80;
((uint32_t*)&nm[nmlength - 8])[0] = mlength << 3;
((uint32_t*)&nm[nmlength - 8])[1] = 0;
// hash padded message (nmessage)
AA = 0x67452301; BB = 0xefcdab89; CC = 0x98badcfe; DD = 0x10325476;
for(i = 0; i < (nmlength/64); ++i)
{
// create table
for(j = 0; j < 16; ++j)
{
X[j] = nm[i*64+j*4+0] |
(nm[i*64+j*4+1] << 8) |
(nm[i*64+j*4+2] << 16) |
(nm[i*64+j*4+3] << 24);
}
A = AA;
B = BB;
C = CC;
D = DD;
// round one
ii = 0;
for(k = 0; k < 4; ++k)
{
A = B + ROTATE_LEFT((A + F(B, C, D) + X[__k1[k*4+0]] + T[ii+0]), __s1[0]);
D = A + ROTATE_LEFT((D + F(A, B, C) + X[__k1[k*4+1]] + T[ii+1]), __s1[1]);
C = D + ROTATE_LEFT((C + F(D, A, B) + X[__k1[k*4+2]] + T[ii+2]), __s1[2]);
B = C + ROTATE_LEFT((B + F(C, D, A) + X[__k1[k*4+3]] + T[ii+3]), __s1[3]);
printf("X[..]: %u\n", __k1[k*4+0]);
ii += 4;
}
// round two
for(k = 0; k < 4; ++k)
{
A = B + ROTATE_LEFT((A + G(B, C, D) + X[__k2[k*4+0]] + T[ii+0]), __s2[0]);
D = A + ROTATE_LEFT((D + G(A, B, C) + X[__k2[k*4+1]] + T[ii+1]), __s2[1]);
C = D + ROTATE_LEFT((C + G(D, A, B) + X[__k2[k*4+2]] + T[ii+2]), __s2[2]);
B = C + ROTATE_LEFT((B + G(C, D, A) + X[__k2[k*4+3]] + T[ii+3]), __s2[3]);
ii += 4;
}
// round three
for(k = 0; k < 4; ++k)
{
A = B + ROTATE_LEFT((A + H(B, C, D) + X[__k3[k*4+0]] + T[ii+0]), __s3[0]);
D = A + ROTATE_LEFT((D + H(A, B, C) + X[__k3[k*4+1]] + T[ii+1]), __s3[1]);
C = D + ROTATE_LEFT((C + H(D, A, B) + X[__k3[k*4+2]] + T[ii+2]), __s3[2]);
B = C + ROTATE_LEFT((B + H(C, D, A) + X[__k3[k*4+3]] + T[ii+3]), __s3[3]);
ii += 4;
}
// round four
for(k = 0; k < 4; ++k)
{
A = B + ROTATE_LEFT((A + I(B, C, D) + X[__k4[k*4+0]] + T[ii+0]), __s4[0]);
D = A + ROTATE_LEFT((D + I(A, B, C) + X[__k4[k*4+1]] + T[ii+1]), __s4[2]);
C = D + ROTATE_LEFT((C + I(D, A, B) + X[__k4[k*4+2]] + T[ii+2]), __s4[3]);
B = C + ROTATE_LEFT((B + I(C, D, A) + X[__k4[k*4+3]] + T[ii+3]), __s4[4]);
ii += 4;
}
printf("----\n");
printf("ii:%u\n", ii);
AA = AA + A;
BB = BB + B;
CC = CC + C;
DD = DD + D;
}
printf("digest: ");
printf("%02x%02x%02x%02x,", AA>>24, AA>>16&0xFF, AA>>8&0xFF, AA&0xFF);
printf("%02x%02x%02x%02x,", BB>>24, BB>>16&0xFF, BB>>8&0xFF, BB&0xFF);
printf("%02x%02x%02x%02x,", CC>>24, CC>>16&0xFF, CC>>8&0xFF, CC&0xFF);
printf("%02x%02x%02x%02x,", DD>>24, DD>>16&0xFF, DD>>8&0xFF, DD&0xFF);
printf("\n");
return;
}
int main(int argc, char *argv[])
{
uint8_t *message = (uint8_t*)"abc";
uint32_t T[64];
md5_genT(&T[0]);
md5_hash(message, strlen(message), &T[0]);
return 1;
}