結果

問題 No.1529 Constant Lcm
ユーザー CleyLCleyL
提出日時 2023-11-13 10:20:51
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 16 ms / 3,000 ms
コード長 996 bytes
コンパイル時間 922 ms
コンパイル使用メモリ 76,216 KB
実行使用メモリ 6,676 KB
最終ジャッジ日時 2023-11-13 10:20:55
合計ジャッジ時間 2,369 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,676 KB
testcase_01 AC 2 ms
6,676 KB
testcase_02 AC 2 ms
6,676 KB
testcase_03 AC 2 ms
6,676 KB
testcase_04 AC 2 ms
6,676 KB
testcase_05 AC 2 ms
6,676 KB
testcase_06 AC 2 ms
6,676 KB
testcase_07 AC 2 ms
6,676 KB
testcase_08 AC 2 ms
6,676 KB
testcase_09 AC 2 ms
6,676 KB
testcase_10 AC 14 ms
6,676 KB
testcase_11 AC 7 ms
6,676 KB
testcase_12 AC 3 ms
6,676 KB
testcase_13 AC 13 ms
6,676 KB
testcase_14 AC 9 ms
6,676 KB
testcase_15 AC 10 ms
6,676 KB
testcase_16 AC 7 ms
6,676 KB
testcase_17 AC 5 ms
6,676 KB
testcase_18 AC 5 ms
6,676 KB
testcase_19 AC 9 ms
6,676 KB
testcase_20 AC 15 ms
6,676 KB
testcase_21 AC 16 ms
6,676 KB
testcase_22 AC 15 ms
6,676 KB
testcase_23 AC 15 ms
6,676 KB
testcase_24 AC 15 ms
6,676 KB
testcase_25 AC 15 ms
6,676 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
}
0