sha256

0

Witam,

mam pewien problem z kodem dla funckji skrótu sha256.
Poniżej znajduje się mój spreparowany na potrzeby testowania kod programu.
Wszystko oparte jest na publikacji: http://csrc.nist.gov/publications/fips/fips180-2/fips180-2.pdf

Problem w działaniu programu jest iż dla t=2 w pętli

  for(int t=0; t<64; t++)
  {

    T1=h+S1(e)+Ch(e,f,g)+K[t]+W[t];
    T2=S0(a)+Maj(a,b,c);

    h=g;
    g=f;
    f=e;
    e=(d+T1);
    d=c;
    c=b;
    b=a;
    a=(T1+T2);

   }

przypisuje dla wartości a i e złe wyniki. Nie potrafię zlokalizować w czym jest błąd i dlatego zwracam się z prośbą do was. Jeśli będą potrzebne inne informacje w celu rozwiazania problemu proszę piszcie.

Cały kod:


#include <iostream.h>
#include <math.h>

using namespace std;

#define Ch(x, y, z)     ((x & y)^(~x & z))
#define Maj(x, y, z)    ((x & y) ^ (x & z) ^ (z & y))
#define SHR(x, n)       (x >> n)
#define ROTR(x, n)      ((x >> n) | (x << (32 - n)))
#define S0(x)           (ROTR(x, 2) ^ ROTR(x, 13) ^ ROTR(x, 22))
#define S1(x)           (ROTR(x, 6) ^ ROTR(x, 11) ^ ROTR(x, 25))
#define s0(x)           (ROTR(x, 7) ^ ROTR(x, 18) ^ SHR(x, 3))
#define s1(x)           (ROTR(x, 17) ^ ROTR(x, 19) ^ SHR(x, 10))

int main()
{

int K[64] = {
                0x428a2f98, 0x71374491, 0xb5c0fbcf, 0xe9b5dba5, 0x3956c25b, 0x59f111f1, 0x923f82a4, 0xab1c5ed5,
                0xd807aa98, 0x12835b01, 0x243185be, 0x550c7dc3, 0x72be5d74, 0x80deb1fe, 0x9bdc06a7, 0xc19bf174,
                0xe49b69c1, 0xefbe4786, 0x0fc19dc6, 0x240ca1cc, 0x2de92c6f, 0x4a7484aa, 0x5cb0a9dc, 0x76f988da,
                0x983e5152, 0xa831c66d, 0xb00327c8, 0xbf597fc7, 0xc6e00bf3, 0xd5a79147, 0x06ca6351, 0x14292967,
                0x27b70a85, 0x2e1b2138, 0x4d2c6dfc, 0x53380d13, 0x650a7354, 0x766a0abb, 0x81c2c92e, 0x92722c85,
                0xa2bfe8a1, 0xa81a664b, 0xc24b8b70, 0xc76c51a3, 0xd192e819, 0xd6990624, 0xf40e3585, 0x106aa070,
                0x19a4c116, 0x1e376c08, 0x2748774c, 0x34b0bcb5, 0x391c0cb3, 0x4ed8aa4a, 0x5b9cca4f, 0x682e6ff3,
                0x748f82ee, 0x78a5636f, 0x84c87814, 0x8cc70208, 0x90befffa, 0xa4506ceb, 0xbef9a3f7, 0xc67178f
            };

int H[8] =  {
                0x6a09e667, 0xbb67ae85, 0x3c6ef372, 0xa54ff53a,
                0x510e527f, 0x9b05688c, 0x1f83d9ab, 0x5be0cd19
            };

int longe = 24;
int bit_text[24];

cout.setf(ios::hex,   ios::basefield);


//abc
bit_text[0]=0;
bit_text[1]=1;
bit_text[2]=1;
bit_text[3]=0;
bit_text[4]=0;
bit_text[5]=0;
bit_text[6]=0;
bit_text[7]=1;
bit_text[8]=0;
bit_text[9]=1;
bit_text[10]=1;
bit_text[11]=0;
bit_text[12]=0;
bit_text[13]=0;
bit_text[14]=1;
bit_text[15]=0;
bit_text[16]=0;
bit_text[17]=1;
bit_text[18]=1;
bit_text[19]=0;
bit_text[20]=0;
bit_text[21]=0;
bit_text[22]=1;
bit_text[23]=1;

int pow_dopelnienia         =       int((longe+65)/512)+1;
int wielomian_dopelnienia   =       512*pow_dopelnienia;
int ile_zer                 =       wielomian_dopelnienia-(longe+65);

int* bit_sha    =   new int[wielomian_dopelnienia];

for(int i=0; i<longe; i++)
{
  bit_sha[i]=bit_text[i];
}

bit_sha[longe]=1;

for(int i=longe+1; i<longe+1+ile_zer; i++)
{
  bit_sha[i]=0;
}

int tmp     =   longe;
int p       =   wielomian_dopelnienia-1 ;

    while(p>=longe+1+ile_zer)
    {
        bit_sha[p]   =   tmp%2;
        tmp          =   tmp/2;
        p--;
    };


int M[1][16];

int tmp_32  =   0;
int many    =   0;

for(int i=0; i<pow_dopelnienia; i++)
{
  for(int j=0; j<16; j++)
  {
      tmp_32=0;
    for(int k=31; k>=0; k--)
    {
      tmp_32+=int(bit_sha[many]*pow(2,k));
      many++;
    }
    M[i][j]=tmp_32;
    cout << M[i][j] << " ";
  }
}

int W[64];

int a,b,c,d,e,f,g,h;
int T1=0;
int T2=0;

for(int i=0; i<pow_dopelnienia;i++)
{
  for(int t=0; t<64; t++)
  {
      W[t]=M[i][t];
    if(t>15)
      W[t]=s1(W[t-2])+W[t-7]+s0(W[t-15])+W[t-16];
  }

a=H[0];
b=H[1];
c=H[2];
d=H[3];
e=H[4];
f=H[5];
g=H[6];
h=H[7];

      cout << "\n" << a << " " << b << " " << c << " " << d << " " << e << " " << f << " " << g << " " << h << " \n";

  for(int t=0; t<64; t++)
  {

    T1=h+S1(e)+Ch(e,f,g)+K[t]+W[t];
    T2=S0(a)+Maj(a,b,c);

    h=g;
    g=f;
    f=e;
    e=(d+T1);
    d=c;
    c=b;
    b=a;
    a=(T1+T2);

       cout << "\n" << a << " " << b << " " << c << " " << d << " " << e << " " << f << " " << g << " " << h << " ";

  }

H[0]+=a;
H[1]+=b;
H[2]+=c;
H[3]+=d;
H[4]+=e;
H[5]+=f;
H[6]+=g;
H[7]+=h;

}

cout << "\n\n" << H[0] << H[1] << H[2] << H[3] << H[4] << H[5] << H[6] << H[7] <<"\n";

return 0;
}

Przepraszam za drugi post ale bład jest już dla t=1, dla t=0 wszystko działa dobrze.

0

Znalazłem taki wyciek pamięci: jest zadeklarowana tablica M[1][16], a w tej pętli jest odczyt indeksów 0..63 (tablica ma tylko 16)

for(int t=0; t<64; t++)
{
    W[t]=M[i][t];
// ...
0

Ale to za wiele nie zmienia. Po zabezpieczeniu przed wyciekiem pamięci w tym miejscu do postaci:

  for(int t=0; t<64; t++)
  {
    if(t<16)
      W[t]=M[i][t];
    if(t>15)
      W[t]=s1(W[t-2])+W[t-7]+s0(W[t-15])+W[t-16];
  } 

Nic w wyniku się nie zmienia.

A założenie w teorii dla tej pętli jest:

W[t] = M[i][t]                                           0<=t<=15
W[t] = s1(W[t-2])+W[t-7]+s0(W[t-15])+W[t-16]            16<=t<=63
0

Temat możemy zamknąć, rozwiązaniem okazała się praca na niewłaściwym typie zmiennej.

Zamiast "int" należało użyć "uint32_t", dlatego że liczby w zapisie hex'a były zbyt duże do przechowywania w zmiennej "int".

Dzięki wszystkim, którzy podjeli próbę pomocy.

1 użytkowników online, w tym zalogowanych: 0, gości: 1