#include #include #include #include #include using mint = atcoder::modint998244353; std::vector prime_factors(int64_t m) { std::vector pf; for (int64_t p = 2; p * p <= m; ++p) { if (m % p == 0) { do m /= p; while (m % p == 0); pf.push_back(p); } } if (m != 1) { pf.push_back(m); } return pf; } mint solve(int32_t n, int64_t m, std::vector pf, std::vector a, std::vector w) { { std::vector a2; std::vector w2; for (int32_t i = 0; i < n; ++i) if (m % a[i] == 0) { a2.push_back(a[i]); w2.push_back(w[i]); } a.swap(a2); w.swap(w2); n = a.size(); } // A _ i | m for all i const int32_t k = pf.size(); std::vector prd(1 << k, 1); for (int32_t i = 0; i < n; ++i) { int32_t t = 0; for (int32_t j = 0; j < k; ++j) { if ((m / pf[j]) % a[i] == 0) t |= 1 << j; } prd[t] *= 1 + w[i]; } std::vector> prd_vec; for (int32_t s = 0; s < 1 << k; ++s) { if (prd[s] != 1) prd_vec.emplace_back(s, prd[s]); } mint ans = 0; for (int32_t s = 0; s < 1 << k; ++s) { mint add = 1; for (auto [t, pw] : prd_vec) { if ((s & t) == s) { add *= pw; } } if ((k - __builtin_popcount(s)) & 1) { ans -= add; } else { ans += add; } } return ans - (m == 1); } int main() { std::ios::sync_with_stdio(false); std::cin.tie(nullptr); int32_t t; int64_t m; std::cin >> t >> m; auto pf = prime_factors(m); for (int32_t testcase = 0; testcase < t; ++testcase) { int32_t n; std::cin >> n; std::vector a(n); for (auto &e : a) std::cin >> e; std::vector w(n); for (auto &e : w) std::cin >> e; std::cout << solve(n, m, pf, a, w).val() << '\n'; } return 0; }