#include <bits/stdc++.h>

#include <atcoder/all>

template <class T>
long long Inversion(std::vector<T> vec) {
	long long res = 0;
	int n = *max_element(vec.begin(), vec.end()) + 1;
	std::vector<T> d(n + 1, 0);
	for(int i = 0; i < vec.size(); i++) {
		res += i;
		for(int j = vec[i] + 1; j; j -= j & -j) res -= d[j];
		for(int j = vec[i] + 1; j <= n; j += j & -j) d[j] += 1;
	}
	return res;
}

void test(int m) {
	std::vector<int> vec(m);
	std::iota(vec.begin(), vec.end(), 0);
	int g = atcoder::internal::primitive_root_constexpr(m);
	for(auto &v : vec) (v *= g) %= m;
	auto inv = Inversion(vec);
	assert(inv % 2 == 1);
}

int main() {
	for(int i = 3; i <= 40009; i++)
		if(atcoder::internal::is_prime_constexpr(i)) test(i);
	return 0;
}