結果
| 問題 |
No.2264 Gear Coloring
|
| コンテスト | |
| ユーザー |
kwm_t
|
| 提出日時 | 2023-04-08 18:11:35 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 163 ms / 2,000 ms |
| コード長 | 1,553 bytes |
| コンパイル時間 | 4,235 ms |
| コンパイル使用メモリ | 254,652 KB |
| 最終ジャッジ日時 | 2025-02-12 03:57:35 |
|
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 18 |
ソースコード
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
//using mint = modint1000000007;
//const int mod = 1000000007;
using mint = modint998244353;
const int mod = 998244353;
//const int INF = 1e9;
//const long long LINF = 1e18;
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define rep2(i,l,r)for(int i=(l);i<(r);++i)
#define rrep(i, n) for (int i = (n-1); i >= 0; --i)
#define rrep2(i,l,r)for(int i=(r-1);i>=(l);--i)
#define all(x) (x).begin(),(x).end()
#define allR(x) (x).rbegin(),(x).rend()
#define endl "\n"
#define P pair<int,int>
template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; }
template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; }
//約数列挙
vector<long long> DivisorEnumeration(long long n) {
vector<long long> x;
for (long long i = 1; i*i <= n; i++) {
if (0 == n % i) {
x.push_back(i);
if (n != i * i) x.push_back(n / i);
}
}
sort(all(x));
return x;
}
int main() {
ios::sync_with_stdio(false);
cin.tie(nullptr);
int n, m; cin >> n >> m;
vector<int>a(n);
rep(i, n)cin >> a[i];
int l = 1;
rep(i, n)l = lcm(l, a[i]);
auto d = DivisorEnumeration(l);
vector<int>c(d.size());
rep(i, c.size())c[i] = l / d[i];
mint ans = 0;
rrep(i, d.size()) {
mint tmp = c[i];
rep(j, n) tmp *= pow_mod(m, gcd(d[i],a[j]), mod);
ans += tmp;
rep(j, i) if (0 == d[i] % d[j]) c[j] -= c[i];
}
ans /= l;
cout << ans.val() << endl;
return 0;
}
kwm_t