#include using namespace std; const long long MOD = 998244353; long long modpow(long long a, long long b){ long long ans = 1; while (b > 0){ if (b % 2 == 1){ ans *= a; ans %= MOD; } a *= a; a %= MOD; b /= 2; } return ans; } long long modinv(long long a){ return modpow(a, MOD - 2); } int main(){ int N; long long M; cin >> N >> M; vector A(N); for (int i = 0; i < N; i++){ cin >> A[i]; } int L = 1; for (int i = 0; i < N; i++){ L = lcm(L, A[i]); } vector F; for (int i = 1; i * i <= L; i++){ if (L % i == 0){ F.push_back(i); if (i * i < L){ F.push_back(L / i); } } } sort(F.begin(), F.end()); int cnt = F.size(); vector C1(cnt, 0); for (int i = cnt - 1; i >= 0; i--){ C1[i] = L / F[i]; for (int j = i + 1; j < cnt; j++){ if (F[j] % F[i] == 0){ C1[i] -= C1[j]; } } } vector C2(cnt, 1); for (int i = 0; i < cnt; i++){ for (int j = 0; j < N; j++){ C2[i] *= modpow(M, gcd(A[j], F[i])); C2[i] %= MOD; } } long long ans = 0; for (int i = 0; i < cnt; i++){ ans += C1[i] * C2[i]; ans %= MOD; } ans *= modinv(L); ans %= MOD; cout << ans << endl; }