結果

問題 No.3250 最小公倍数
コンテスト
ユーザー jiangxinyang
提出日時 2025-08-04 01:03:53
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 3,085 bytes
コンパイル時間 2,299 ms
コンパイル使用メモリ 214,304 KB
実行使用メモリ 22,996 KB
最終ジャッジ日時 2025-10-16 16:05:55
合計ジャッジ時間 9,704 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample -- * 1
other TLE * 1 -- * 21
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:52:16: warning: ignoring return value of ‘FILE* freopen(const char*, const char*, FILE*)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   52 |         freopen("lcm.in", "r", stdin);
      |         ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~
main.cpp:53:16: warning: ignoring return value of ‘FILE* freopen(const char*, const char*, FILE*)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   53 |         freopen("lcm.out", "w", stdout);
      |         ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~

ソースコード

diff # 

#include <bits/stdc++.h>
#define F first
#define S second

using namespace std;

inline int read();
void write(int);
void writeln(int);

const int P = 998244353;
const int A = 1e6;

int n;
vector<int> prime;

void getprime() {
    prime.resize(A + 1);
    iota(prime.begin(), prime.end(), 0);
    for(int i = 2; i * i <= A; i++) 
        if(prime[i] == i) 
			for(int j = i * i; j <= A; j += i) 
				if(prime[j] == j) 
					prime[j] = i;
	return ;
}

void fact(int x, unordered_map<int, int>& f) {
    f.clear();
    if(x == 1) return;
    while(x > 1) {
        int p = prime[x];
        int cnt = 0;
        while(x % p == 0) x /= p, cnt++;
        f[p] = cnt;
    }
}

long long qpow(long long p, int e) {
    long long ret = 1;
    p %= P;
    while(e) {
        if(e & 1) ret = (ret * p) % P;
        p = (p * p) % P;
        e >>= 1;
    }
    return ret;
}

int main() {
	
	freopen("lcm.in", "r", stdin);
	freopen("lcm.out", "w", stdout);
    
    getprime();

    n = read();
    vector<int> a(n + 1);
    for(int i = 1; i <= n; i++) a[i] = read();

    vector<vector<int> > G(n + 1);
    for(int i = 0; i < n - 1; i++) {
        int u = read(), v = read();
        G[u].push_back(v);
        G[v].push_back(u);
    }

    vector<int> fa(n + 1, -1);
    queue<int> que;
    
    que.push(1);
    fa[1] = 0;
    while(!que.empty()) {
        int u = que.front();
        que.pop();
        for(int v : G[u]) if(fa[v] == -1) fa[v] = u, que.push(v);
    }
    vector<vector<int> > cd(n + 1);
    for(int v = 2; v <= n; ++v) cd[fa[v]].push_back(v);
    vector<unordered_map<int, int> > f(n + 1);
    for(int i = 1; i <= n; i++) fact(a[i], f[i]);
    vector<int> cct(n + 1, 0);
    for(int u = 1; u <= n; u++) cct[u] = cd[u].size();
    queue<int> q;
    for(int u = 1; u <= n; u++) if(cct[u] == 0) q.push(u);
    vector<int> ans(n + 1);
    while(!q.empty()) {
        int u = q.front();
        q.pop();
        for(int v : cd[u]) {
            if(f[u].size() < f[v].size()) swap(f[u], f[v]);
            for(const auto& kv : f[v]) {
                int p = kv.F, e = kv.S;
                auto it = f[u].find(p);
                if(it == f[u].end()) f[u][p] = e;
                else it->S = max(it->S, e);
            }
            f[v].clear();
        }
        long long res = 1;
        for(const auto& kv : f[u]) res = res * qpow(kv.F, kv.S) % P;
        ans[u] = res;
        int p = fa[u];
        if(p != 0) { if(--cct[p] == 0) q.push(p); }
    }

    for(int i = 1; i <= n; i++) printf("%d\n", ans[i]);

//	printf("\nThe time used: ");
//	printf("%.2lfs",(double)clock()/CLOCKS_PER_SEC);

	return 0;

}

inline int read() {
	int res = 0, f = 1;
	char ch = getchar();
	while( !(ch >= '0' && ch <= '9') ) {
		if(ch == '-') f = -1;
		ch = getchar();
	}
	while(ch >= '0' && ch <= '9') {
		res = (res << 1) + (res << 3) + (ch ^ 48);
		ch = getchar();
	}
	return res * f;
}

void write(int x) {
	static int sta[35];
	int top = 0;
	do {
		sta[top++] = x % 10;
		x /= 10;
	} while(x);
	while(top) putchar(sta[--top] ^ 48);
}

void writeln(int x) {
	write(x);
	putchar('\n');
}
0