結果

問題 No.2497 GCD of LCMs
ユーザー kwm_t
提出日時 2023-10-07 12:04:40
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 21 ms / 2,000 ms
コード長 1,973 bytes
コンパイル時間 4,407 ms
コンパイル使用メモリ 260,628 KB
最終ジャッジ日時 2025-02-17 06:06:33
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 14
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
using mint = modint998244353;
const int mod = 998244353;
//using mint = modint1000000007;
//const int mod = 1000000007;
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 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; }
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];
    vector<vector<int>>to(n);
    rep(i, m) {
        int u, v; cin >> u >> v;
        u--, v--;
        to[u].push_back(v);
        to[v].push_back(u);
    }
    set<int>pst;
    rep(i, n) {
        int x = a[i];
        for (int j = 2; j* j <= x; ++j) {
            if (0 != x % j)continue;
            pst.insert(j);
            while (!(x%j))x /= j;
        }
        if (1 != x)pst.insert(x);
    }
    vector<mint>ans(n, 1);
    for (auto p : pst) {
        vector<int>cnt(n), dp(n, INF);
        rep(i, n) {
            int ca = a[i];
            while (!(ca%p)) {
                cnt[i]++;
                ca /= p;
            }
        }
        std::priority_queue<P, vector<P>, std::greater<P>> pq;
        pq.push({ cnt[0],0 });
        dp[0] = cnt[0];
        while (!pq.empty()) {
            auto top = pq.top();
            pq.pop();
            int val = top.first;
            int pos = top.second;
            if (dp[pos] < val) continue;
            for (auto npos : to[pos]) {
                int nval = max(dp[pos], cnt[npos]);
                if (chmin(dp[npos], nval)) {
                    pq.push({ dp[npos],npos });
                }
            }
        }
        vector<mint>pw(100);
        pw[0] = 1;
        rep2(i, 1, 100)pw[i] = pw[i - 1] * p;
        rep(i, n)ans[i] *= pw[dp[i]];
    }
    rep(i, n)cout << ans[i].val() << endl;
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0