- //mpz_sub_ui(t_exp, p, 1);
- mpz_divexact_ui(t_exp, p_minus1, 2);
-
-#if 0
- //printf("prime check a=2\n");
- mpz_powm(t, t_p1, t_exp, proth);
- mpz_add_ui(t, t, 1);
- //printf("cmp1\n");
- if(mpz_cmp(t, proth) == 0)
- {
- printf("prime with a=2!\n");
- break;
- }
-#endif