blob: 9fedc8115e846f3c4d2c99370b3893d99a00c81b (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
|
#include <botan/numthry.h>
using namespace Botan;
#include <iostream>
void test_ressol(const BigInt& p, RandomNumberGenerator& rng)
{
std::cout << p << std::endl;
// const BigInt p_16 = p / 16;
int noroot = 0, false_result = 0;
for(int j = 0; j != 1000; ++j)
{
BigInt x = random_integer(rng, 0, p);
//if(x % p_16 == 0)
//std::cout << "p = " << p << " x = " << x << "\n";
BigInt sqrt_x = ressol(x, p);
if(sqrt_x < 0)
{
++noroot;
continue;
}
BigInt check = square(sqrt_x) % p;
if(check != x % p)
{
std::cout << "FAIL "
<< "x = " << x << "; "
<< "p = " << p << "; "
<< "s = " << sqrt_x << "; "
<< "s^2%p = " << check << "\n";
++false_result;
}
}
/*
std::cout << "nomatch=" << nomatch << " "
<< "noroot=" << noroot << " "
<< "false=" << false_result << "\n";
*/
}
int main()
{
RandomNumberGenerator* rng = RandomNumberGenerator::make_rng();
#if 0
std::cout << ressol(8, 17) << "\n";
std::cout << ressol_orig(8, 17) << "\n";
#endif
#if 1
for(int j = 16; j != 1024; ++j)
{
std::cout << "Round " << j << "\n";
BigInt p = random_prime(*rng, j);
test_ressol(p, *rng);
//printf("%d\n", j);
}
#endif
/*
for(u32bit j = 9; j != PRIME_TABLE_SIZE; ++j)
{
std::cout << "PRIME[" << j << "] == " << PRIMES[j] << std::endl;
//printf("%d - ", PRIMES[j]);
test_ressol(PRIMES[j], *rng);
//printf("\n");
}
*/
}
|