結果
問題 | No.2264 Gear Coloring |
ユーザー |
![]() |
提出日時 | 2023-04-07 22:59:19 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 151 ms / 2,000 ms |
コード長 | 2,858 bytes |
コンパイル時間 | 1,486 ms |
コンパイル使用メモリ | 126,128 KB |
最終ジャッジ日時 | 2025-02-12 02:24:07 |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 18 |
ソースコード
#include <iostream> #include <algorithm> #include <iomanip> #include <vector> #include <queue> #include <deque> #include <set> #include <map> #include <tuple> #include <cmath> #include <numeric> #include <functional> #include <cassert> #include <atcoder/modint> #define debug_value(x) cerr << "line" << __LINE__ << ":<" << __func__ << ">:" << #x << "=" << x << endl; #define debug(x) cerr << "line" << __LINE__ << ":<" << __func__ << ">:" << x << endl; template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } using namespace std; typedef long long ll; template<typename T> vector<vector<T>> vec2d(int n, int m, T v){ return vector<vector<T>>(n, vector<T>(m, v)); } template<typename T> vector<vector<vector<T>>> vec3d(int n, int m, int k, T v){ return vector<vector<vector<T>>>(n, vector<vector<T>>(m, vector<T>(k, v))); } template<typename T> void print_vector(vector<T> v, char delimiter=' '){ if(v.empty()) { cout << endl; return; } for(int i = 0; i+1 < v.size(); i++) cout << v[i] << delimiter; cout << v.back() << endl; } using mint = atcoder::modint998244353; ostream& operator<<(ostream& os, const mint& m){ os << m.val(); return os; } class Eratosthenes{ public: int m; vector<bool> is_prime; vector<int> primes; Eratosthenes(int m_){ m = m_; init(); } private: void init(){ is_prime.assign(m+1, true); is_prime[0] = false, is_prime[1] = false; for(int i = 2; i <= m; i++){ if(is_prime[i]){ primes.push_back(i); for(int j = 2; i*j <= m; j++) is_prime[i*j] = false; } } } }; auto et = Eratosthenes(100000); int tot(int n){ int ans = n; int m = n; for(int p: et.primes){ if(n%p == 0){ while(m%p == 0) m /= p; ans /= p; ans *= p-1; } } if(m != 1){ ans /= m; ans *= m-1; } return ans; } vector<int> divs(int n){ vector<int> ans; for(int x = 1; x*x <= n; x++){ if(n%x == 0){ ans.push_back(x); if(x*x != n) ans.push_back(n/x); } } return ans; } int main(){ ios::sync_with_stdio(false); cin.tie(0); cout << setprecision(10) << fixed; int n, m; cin >> n >> m; vector<int> a(n); ll l = 1; for(int i = 0; i < n; i++) { cin >> a[i]; l = lcm(l, (ll)a[i]); } mint ans = 0; auto d = divs(l); for(int g: d){ mint tmp = 1; for(int i = 0; i < n; i++){ int gg = gcd(g, a[i]); tmp *= mint(m).pow(gg); } ans += tmp*tot(l/g); } ans /= l; cout << ans << endl; }