#line 2 "library/other/template.hpp" #define _CRT_SECURE_NO_WARNINGS #ifdef ONLINE_JUDGE #pragma GCC target("avx512f") #elif defined EVAL #else #pragma GCC target("avx2") #endif #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define rep(i, n) for (int i = 0; i < (n); i++) #define REP(i, n) for (int i = 1; i <= (n); i++) #define all(V) V.begin(), V.end() using i128 = __int128_t; using u128 = __uint128_t; using uint = unsigned int; using lint = long long; using ulint = unsigned long long; using P = std::pair; using LP = std::pair; constexpr int INF = INT_MAX / 2; constexpr lint LINF = LLONG_MAX / 2; constexpr double eps = DBL_EPSILON; constexpr double PI = 3.141592653589793238462643383279; template class prique : public std::priority_queue, std::greater> { }; template inline constexpr decltype(auto) lambda_fix(F&& f) { return [f = std::forward(f)](auto&&... args) { return f(f, std::forward(args)...); }; } template constexpr std::vector make_vec(size_t n) { return std::vector(n); } template constexpr auto make_vec(size_t n, Args&&... args) { return std::vector(args...))>( n, make_vec(std::forward(args)...)); } template constexpr inline bool chmax(T& lhs, const U& rhs) { if (lhs < rhs) { lhs = rhs; return true; } return false; } template constexpr inline bool chmin(T& lhs, const U& rhs) { if (lhs > rhs) { lhs = rhs; return true; } return false; } constexpr inline lint gcd(lint a, lint b) { while (b) { lint c = a; a = b; b = c % b; } return a; } inline lint lcm(lint a, lint b) { return a / gcd(a, b) * b; } constexpr bool isprime(lint n) { if (n == 1) return false; for (int i = 2; i * i <= n; i++) { if (n % i == 0) return false; } return true; } template constexpr T mypow(T a, lint b) { T res(1); while (true) { if (b & 1) res *= a; b >>= 1; if (!b) break; a *= a; } return res; } constexpr lint modpow(lint a, lint b, lint m) { a %= m; lint res(1); while (b) { if (b & 1) { res *= a; res %= m; } a *= a; a %= m; b >>= 1; } return res; } template constexpr void printArray(const std::vector& vec, char split = ' ') { rep(i, vec.size()) { std::cout << vec[i]; std::cout << (i == (int)vec.size() - 1 ? '\n' : split); } } template constexpr void printArray(InputIter l, InputIter r, char split = ' ') { auto rprev = std::prev(r); for (auto i = l; i != r; i++) { std::cout << *i; std::cout << (i == rprev ? '\n' : split); } } LP extGcd(lint a, lint b) { if (b == 0) return {1, 0}; LP s = extGcd(b, a % b); std::swap(s.first, s.second); s.second -= a / b * s.first; return s; } LP ChineseRem(const lint& b1, const lint& m1, const lint& b2, const lint& m2) { lint p = extGcd(m1, m2).first; lint tmp = (b2 - b1) * p % m2; lint r = (b1 + m1 * tmp + m1 * m2) % (m1 * m2); return {r, m1 * m2}; } int LCS(const std::string& a, const std::string& b) { auto dp = make_vec(a.size() + 1, b.size() + 1); rep(i, a.size()) { rep(j, b.size()) { chmax(dp[i + 1][j], dp[i][j]); chmax(dp[i][j + 1], dp[i][j]); if (a[i] == b[j]) chmax(dp[i + 1][j + 1], dp[i][j] + 1); } chmax(dp[i + 1][b.size()], dp[i][b.size()]); } rep(j, b.size()) chmax(dp[a.size()][j + 1], dp[a.size()][j]); return dp[a.size()][b.size()]; } template ::value, std::nullptr_t> = nullptr> void compress(std::vector& vec) { auto tmp = vec; std::sort(all(tmp)); tmp.erase(std::unique(all(tmp)), tmp.end()); for (T& i : vec) i = std::lower_bound(all(tmp), i) - tmp.begin(); } template void compress(T* l, T* r) { std::vector tmp(l, r); std::sort(all(tmp)); tmp.erase(std::unique(all(tmp)), tmp.end()); for (auto i = l; i < r; i++) { *i = std::lower_bound(all(tmp), *i) - tmp.begin(); } } template void compress(InputIter l, InputIter r) { std::vector tmp(l, r); std::sort(all(tmp)); tmp.erase(std::unique(all(tmp)), tmp.end()); for (auto i = l; i < r; i++) { *i = std::lower_bound(all(tmp), *i) - tmp.begin(); } } #line 2 "main.cpp" std::mt19937 mt; int main() { int T, v, x; std::cin >> T; rep(i, T) { std::cin >> v >> x; int p = v * x + 1; std::vector vec; int v_ = v, x_ = x; for (int j = 2; j * j <= v_; j++) { while (v_ % j == 0) { v_ /= j; vec.emplace_back(j); } } if (v_ != 1) vec.emplace_back(v_); for (int j = 2; j * j <= x_; j++) { while (x_ % j == 0) { x_ /= j; vec.emplace_back(j); } } if (x_ != 1) vec.emplace_back(x_); std::sort(all(vec)); vec.erase(std::unique(all(vec)), vec.end()); while (true) { int q = mt() % p; if (!q) continue; bool f = true; for (const auto& j : vec) { if (modpow(q, (p - 1) / j, p) == 1) f = false; } if (f) { lint c = modpow(q, v, p), now = 1; std::vector ans; rep(j, x) { now *= c; now %= p; ans.emplace_back(now); } std::sort(all(ans)); printArray(ans); break; } } } return 0; }