結果
問題 | No.2497 GCD of LCMs |
ユーザー |
|
提出日時 | 2023-09-03 12:07:59 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 225 ms / 2,000 ms |
コード長 | 1,562 bytes |
コンパイル時間 | 4,556 ms |
コンパイル使用メモリ | 261,248 KB |
最終ジャッジ日時 | 2025-02-16 18:22:29 |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 14 |
ソースコード
#include <bits/stdc++.h>#include <atcoder/all>using namespace std;using namespace atcoder;struct Fast {Fast() {std::cin.tie(nullptr);ios::sync_with_stdio(false);cout << setprecision(10);}} fast;#define rep(i, a, b) for (int(i) = (a); (i) < (int)(b); (i)++)using ll = long long;using mint = modint998244353;int main() {int n, m;cin >> n >> m;vector<ll> a(n);rep(i, 0, n) cin >> a[i];vector<vector<int>> g(n, vector<int>());rep(i, 0, m) {int u, v;cin >> u >> v;u--;v--;g[u].push_back(v);g[v].push_back(u);}set<ll> primes;rep(i, 0, n) {ll x = a[i];for (ll d = 2; d * d <= x; d += (d & 1) + 1) {if (x % d) continue;primes.insert(d);while (!(x % d)) x /= d;}if (x > 1) primes.insert(x);}ll amax = 0;rep(i, 0, n) amax = max(amax, a[i]);vector<mint> ans(n, 1);for (auto p : primes) {vector<bool> seen(n, false);for (ll pow = 1; pow <= amax; pow *= p) {if (a[0] % (pow * p) == 0) continue;vector<bool> reachable(n, false);stack<int> st;st.push(0);while (!st.empty()) {auto x = st.top();st.pop();reachable[x] = true;for (auto y : g[x]) {if (a[y] % (pow * p) && !reachable[y]) {reachable[y] = true;st.push(y);}}}rep(i, 0, n) if (reachable[i] && !seen[i]) {ans[i] *= pow;seen[i] = true;}}}rep(i, 0, n) cout << ans[i].val() << "\n";}