test/math/bigint.cpp


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 < ssize(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);
}