blob: 53e18dd1d8461186a70470f0233803c9377b4c1d (
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
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
|
#include "../util.h"
#include <math/bigint.cpp>
template<ll MOD>
struct modInt {
ll value = 0;
modInt() {}
modInt(const bigint& x) {
stringstream a;
a << x;
string b = a.str();
for (ll i = b[0] == '-' ? 1 : 0; i < sz(b); i++) {
value *= 10;
value += b[i] - '0';
value %= MOD;
}
if (b[0] == '-') value = (MOD - value) % MOD;
}
modInt(ll x) : value(((x % MOD) + MOD) % MOD) {}
modInt operator+(modInt o) const {return value + o.value;}
modInt operator-(modInt o) const {return value - o.value;}
modInt operator*(modInt o) const {return value * o.value;}
modInt& operator+=(modInt o) {return *this = *this + o;}
modInt& operator-=(modInt o) {return *this = *this - o;}
modInt& operator*=(modInt o) {return *this = *this * o;}
ll& operator*() {return value;}
bool operator==(const modInt& o) const {return value == o.value;}
bool operator!=(const modInt& o) const {return value != o.value;}
};
constexpr ll MOD = 1'394'633'899;
constexpr ll POOL = 8;
void stress_test(int LIM) {
int queries = 0;
for (int tries = 0; tries < LIM; tries++) {
vector<modInt<MOD>> expectedPool(POOL);
vector<bigint> gotPool(POOL);
for (int i = 0; i < POOL; i++) {
ll x = Random::integer<ll>(-1'000'000'000'000'000'000ll, 1'000'000'000'000'000'000ll);
expectedPool[i] = x;
gotPool[i] = x;
if (expectedPool[i] != modInt<MOD>(gotPool[i])) cerr << "error: 0" << FAIL;
}
for (int i = 0; i < 200; i++) {
int a = Random::integer<int>(0, POOL);
int b = Random::integer<int>(0, POOL);
int o = Random::integer<int>(0, 3);
if (Random::integer<int>(0, 2) == 0) {//x=
auto tmpExpected = expectedPool[a];
auto tmpGot = gotPool[a];
if (o == 0) {
tmpExpected += expectedPool[b];
tmpGot += gotPool[b];
}
if (o == 1) {
tmpExpected -= expectedPool[b];
tmpGot -= gotPool[b];
}
if (o == 2) {
tmpExpected -= expectedPool[b];
tmpGot -= gotPool[b];
}
if (tmpExpected != modInt<MOD>(tmpGot)) {
cerr << gotPool[a];
if (o == 0) cerr << "+";
if (o == 1) cerr << "-";
if (o == 2) cerr << "*";
cerr << gotPool[b] << "=" << tmpGot << endl;
cerr << "error: 1" << FAIL;
}
expectedPool[b] = tmpExpected;
gotPool[b] = tmpGot;
} else {//x
int c = Random::integer<int>(0, POOL);
modInt<MOD> tmpExpected;
bigint tmpGot;
if (o == 0) {
tmpExpected = expectedPool[a] + expectedPool[b];
tmpGot = gotPool[a] + gotPool[b];
}
if (o == 1) {
tmpExpected = expectedPool[a] - expectedPool[b];
tmpGot = gotPool[a] - gotPool[b];
}
if (o == 2) {
tmpExpected = expectedPool[a] * expectedPool[b];
tmpGot = gotPool[a] * gotPool[b];
}
if (tmpExpected != modInt<MOD>(tmpGot)) {
cerr << gotPool[a];
if (o == 0) cerr << "+";
if (o == 1) cerr << "-";
if (o == 2) cerr << "*";
cerr << gotPool[b] << "=" << tmpGot << endl;
cerr << "error: 2" << FAIL;
}
expectedPool[c] = tmpExpected;
gotPool[c] = tmpGot;
}
queries++;
}
}
cerr << "tested random queries: " << queries << endl;
}
int main() {
stress_test(100);
if (!sanitize) stress_test(1000);
}
|
