結果

問題 No.3250 最小公倍数
コンテスト
ユーザー miya145592
提出日時 2025-08-29 23:21:35
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
MLE  
実行時間 -
コード長 2,646 bytes
コンパイル時間 3,938 ms
コンパイル使用メモリ 300,196 KB
実行使用メモリ 814,516 KB
最終ジャッジ日時 2025-10-16 16:40:10
合計ジャッジ時間 11,038 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample -- * 1
other MLE * 1 -- * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const int MOD = 998244353;

int N;
vector<int> A;
vector<vector<int>> G;
vector<int> prime_list;
unordered_map<int,int> rev_list;
int L;

struct NodeData {
    vector<int> exps;        // 各素数の指数
    unordered_set<int> big;  // 1000より大きな素数
};

vector<NodeData> ans;

// エラトステネスの篩
vector<bool> sieve(int n) {
    vector<bool> is_prime(n+1, true);
    is_prime[0] = is_prime[1] = false;
    for (int i = 2; i * i <= n; i++) {
        if (is_prime[i]) {
            for (int j = 2*i; j <= n; j += i) {
                is_prime[j] = false;
            }
        }
    }
    return is_prime;
}

// 素因数分解
void factorization(int n, NodeData &arr) {
    int temp = n;
    for (int p : prime_list) {
        if (temp % p == 0) {
            int cnt = 0;
            while (temp % p == 0) {
                cnt++;
                temp /= p;
            }
            int idx = rev_list[p];
            arr.exps[idx] = cnt;
        }
    }
    if (temp != 1) {
        arr.big.insert(temp);
    }
}

void dfs(int v, int p) {
    factorization(A[v], ans[v]);
    for (int nv : G[v]) {
        if (nv == p) continue;
        dfs(nv, v);
        for (int i = 0; i < L; i++) {
            ans[v].exps[i] = max(ans[v].exps[i], ans[nv].exps[i]);
        }
        for (int x : ans[nv].big) {
            ans[v].big.insert(x);
        }
    }
}

ll modpow(ll a, ll e, ll mod) {
    ll res = 1;
    while (e > 0) {
        if (e & 1) res = (res * a) % mod;
        a = (a * a) % mod;
        e >>= 1;
    }
    return res;
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);

    cin >> N;
    A.resize(N);
    for (int i = 0; i < N; i++) cin >> A[i];

    G.assign(N, {});
    for (int i = 0; i < N-1; i++) {
        int u, v;
        cin >> u >> v;
        u--; v--;
        G[u].push_back(v);
        G[v].push_back(u);
    }

    // 素数リスト作成
    auto P = sieve(1000);
    for (int i = 0; i < (int)P.size(); i++) {
        if (P[i]) {
            rev_list[i] = (int)prime_list.size();
            prime_list.push_back(i);
        }
    }
    L = (int)prime_list.size();

    ans.assign(N, NodeData{vector<int>(L,0), {}});

    dfs(0, -1);

    for (int v = 0; v < N; v++) {
        ll tmp = 1;
        for (int i = 0; i < L; i++) {
            int e = ans[v].exps[i];
            if (e > 0) {
                tmp = (tmp * modpow(prime_list[i], e, MOD)) % MOD;
            }
        }
        for (int p : ans[v].big) {
            tmp = (tmp * (p % MOD)) % MOD;
        }
        cout << tmp << "\n";
    }
}
0