Задачи из книги Молдовяна

Проблема с задачей из книги Практикум по криптосистемам с открытым ключом, параграф 5.1, самая первая задача.
Написал на C++ библиотека GMP
#include <iostream> #include <gmp.h> #include <gmpxx.h> #include <time.h> #include <QDebug> #define FALSE false #define TRUE true #define BASE 10 using namespace std; void generate_rand_autoinit(mpz_t r_out, unsigned long base) { srand(time(0)+5*rand()-7); mpz_t r; mpz_t r_max; mpz_init_set_ui(r,rand()); mpz_init(r_out); mpz_init_set_ui(r_max,2); mpz_pow_ui(r_max,r_max,base); //int length_of_range = strlen(mpz_get_str(NULL,10,r_max)); while(base) { mpz_mul_ui(r,r,rand()); --base; } gmp_randstate_t alg; gmp_randinit_mt(alg); gmp_randseed(alg,r); do{ mpz_urandomm(r_out,alg,r_max); }while( mpz_sgn(r_out)==0 ); gmp_randclear(alg); mpz_clear(r); mpz_clear(r_max); } void generate_prime_autoinit_mill_rabin(mpz_t r_out, unsigned long base) { mpz_t r; mpz_init(r); mpz_init(r_out); //Генерируем нечётное число do{ mpz_clear(r_out); generate_rand_autoinit(r_out,base); mpz_mod_ui(r,r_out,2); }while(mpz_cmp_ui(r,1)!=0); //Проводим тест Миллера-Рабина, пока //тест не выдаст положительный резуельтат int i = mpz_probab_prime_p(r_out,50); while(i==0) { mpz_add_ui(r_out,r_out,2); // mpz_nextprime(r_out,r_out); i = mpz_probab_prime_p(r_out,50); } mpz_clear(r); } void _y(mpz_t y, mpz_t a, mpz_t x, mpz_t p) { mpz_init(y); mpz_powm(y,a,x,p); } void _Z(mpz_t Z, mpz_t a, mpz_t U, mpz_t p) { mpz_init(Z); mpz_powm(Z,a,U,p); } void _k(mpz_t k, mpz_t U, mpz_t x, mpz_t H, mpz_t Z, mpz_t gamma) { mpz_init(k); mpz_mul(k,x,H); mpz_mod(k,k,gamma); mpz_mul(k,k,Z); mpz_mod(k,k,gamma); mpz_sub(k,U,k); mpz_mod(k,k,gamma); mpz_div_ui(k,k,2); mpz_mod(k,k,gamma); } void _g(mpz_t g, mpz_t U, mpz_t x, mpz_t H, mpz_t Z, mpz_t gamma) { mpz_init(g); mpz_mul(g,x,H); mpz_mod(g,g,gamma); mpz_mul(g,g,Z); mpz_mod(g,g,gamma); mpz_add(g,U,g); mpz_mod(g,g,gamma); mpz_div_ui(g,g,2); mpz_mod(g,g,gamma); } bool _cmp(mpz_t a, mpz_t k, mpz_t y, mpz_t H, mpz_t S, mpz_t p) { mpz_t tmp,tmp2; mpz_init(tmp); mpz_init(tmp2); mpz_powm(tmp,a,k,p); mpz_mul(tmp,S,tmp); mpz_mod(tmp,tmp,p); mpz_mul(tmp,H,tmp); mpz_mod(tmp,tmp,p); mpz_powm(tmp,y,tmp,p); mpz_powm(tmp2,a,k,p); mpz_mul(tmp,tmp2,tmp); mpz_mod(tmp,tmp,p); if(mpz_cmp(tmp,S)!=0) return FALSE; return TRUE; } void _S_via_a(mpz_t S, mpz_t a, mpz_t g, mpz_t p) { mpz_init(S); mpz_powm(S,a,g,p); } void _S_via_y(mpz_t S, mpz_t y, mpz_t g, mpz_t p) { mpz_init(S); mpz_powm(S,y,g,p); } void generate_a(mpz_t a, mpz_t gamma, mpz_t p) { mpz_t _gamma; mpz_init(_gamma); mpz_sub_ui(_gamma,p,1); gmp_printf("p-1=%Zd ",_gamma); mpz_div(_gamma,_gamma,gamma); gmp_printf("(p-1)/gamma=%Zd\n",_gamma); //mpz_add_ui(_gamma,_gamma,1); mpz_t b; while(true) { //Генерируем 1<b<p generate_rand_autoinit(b,BASE); while(mpz_cmp(b,p) >= 0) { mpz_clear(b); generate_rand_autoinit(b,BASE); } gmp_printf("b=%Zd ",b); //z:=b^gamma' mod p mpz_powm(b,b,_gamma,p); gmp_printf("z=%Zd\n",b); //z=1 ? if(mpz_cmp_ui(b,1)!=0) { mpz_init_set(a,b); mpz_clear(_gamma); mpz_clear(b); break; } else { mpz_clear(b); continue; } }//while(true) } int main() { mpz_t x,y,a,p,Z,U,H,k,g,gamma,S; generate_prime_autoinit_mill_rabin(p,BASE); generate_rand_autoinit(x,BASE); while(mpz_cmp(x,p)>=0) { mpz_clear(x); generate_rand_autoinit(x,BASE); } generate_rand_autoinit(H,BASE); while(mpz_cmp(H,p)>=0) { mpz_clear(H); generate_rand_autoinit(H,BASE); } generate_rand_autoinit(U,BASE); while(mpz_cmp(U,p)>=0) { mpz_clear(U); generate_rand_autoinit(U,BASE); } //mpz_init(gamma); //mpz_sub_ui(gamma,p,1); generate_prime_autoinit_mill_rabin(gamma,BASE); while(mpz_cmp(gamma,p)>=0) { mpz_clear(gamma); generate_prime_autoinit_mill_rabin(gamma,BASE); } generate_a(a,gamma,p); _y(y,a,x,p); _Z(Z,a,U,p); _k(k,U,x,H,Z,gamma); _g(g,U,x,H,Z,gamma); _S_via_a(S,a,g,p); gmp_printf("p=%Zd ",p); gmp_printf("gamma=%Zd\n",gamma); gmp_printf("x=%Zd ",x); gmp_printf("y=%Zd\n",y); gmp_printf("U=%Zd ",U); gmp_printf("Z=%Zd\n",Z); gmp_printf("a=%Zd\n",a); gmp_printf("H=%Zd\n",H); gmp_printf("k=%Zd ",k); gmp_printf("g=%Zd\n",g); gmp_printf("S=%Zd\n",S); cout << _cmp(a,k,y,H,S,p) << "\n"; return EXIT_SUCCESS; }

Ну в итоге части проверочного уравнения не равны...где ошибку допустил?