結果
問題 |
No.1529 Constant Lcm
|
ユーザー |
👑 |
提出日時 | 2023-11-13 10:20:51 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 32 ms / 3,000 ms |
コード長 | 996 bytes |
コンパイル時間 | 753 ms |
コンパイル使用メモリ | 75,468 KB |
最終ジャッジ日時 | 2025-02-17 21:52:45 |
ジャッジサーバーID (参考情報) |
judge5 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 24 |
ソースコード
#include <iostream> #include <vector> const long long MOD = 998244353; using namespace std; struct SieveEratos{ vector<bool> t; vector<int> primes; SieveEratos(int n):t(n+1,true){ t[0] = t[1] = false; for(int i = 2; n >= i; i++){ if(t[i]){ primes.emplace_back(i); for(int j = i+i; n >= j; j+=i){ t[j] = false; } } } } bool operator[](int x){return t[x];} }; int main(){ long long n;cin>>n; SieveEratos SE(n+1); long long ans = 1; for(long long i = 2; n >= i; i++){ if(SE[i]){ if(n%i){ long long nw = i; while(nw*i < n)nw *= i; ans = (ans*nw)%MOD; }else{ long long nw = i; while(n%(nw*i) == 0){ nw *= i; } if(nw == n){ nw /= i; } long long rem = n/nw; long long nw2 = 1; while(nw2*i < rem)nw2 *= i; ans = (ans * ((((nw*nw)%MOD)*nw2)%MOD))%MOD; } } } cout << ans << endl; }