結果
| 問題 |
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;
}