結果

問題 No.1396 Giri
ユーザー bayashiko
提出日時 2021-02-14 21:30:55
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 2,460 bytes
コンパイル時間 2,961 ms
コンパイル使用メモリ 235,880 KB
最終ジャッジ日時 2025-01-18 20:41:10
ジャッジサーバーID
（参考情報） 		judge2 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2 TLE * 1
other AC * 15 TLE * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC optimize("Ofast")
#include<bits/stdc++.h>
using namespace std;
//#include<boost/multiprecision/cpp_int.hpp>
//#include<boost/multiprecision/cpp_dec_float.hpp>
//namespace mp=boost::multiprecision;
//#define mulint mp::cpp_int
//#define mulfloat mp::cpp_dec_float_100
//#include<atcoder/all>
//using namespace atcoder;
struct __INIT{__INIT(){cin.tie(0);ios::sync_with_stdio(false);cout<<fixed<<setprecision(15);}} __init;
#define max3(a,b,c) max(a,max(b,c))
#define min3(a,b,c) min(a,min(b,c))
#define INF (1<<30)
#define LINF (lint)(1LL<<56)
#define endl "\n"
#define rep(i,n) for(lint (i)=0;(i)<(n);(i)++)
#define reprev(i,n) for(lint (i)=(n-1);(i)>=0;(i)--)
#define flc(x) __builtin_popcountll(x)
#define pint pair<int,int>
#define pdouble pair<double,double>
#define plint pair<lint,lint>
#define fi first
#define se second
#define all(x) x.begin(),x.end()
#define vec vector<lint>
#define nep(x) next_permutation(all(x))
typedef long long lint;
typedef __int128_t llint;
int dx[8]={1,1,0,-1,-1,-1,0,1};
int dy[8]={0,1,1,1,0,-1,-1,-1};
const int MAX_N=3e5+5;
//vector<int> bucket[MAX_N/1000];
//constexpr int MOD=1000000007;
constexpr int MOD=998244353;

int cnt[1000001];
 
multiset<lint> decomp(lint N){
    multiset<lint> res;
    lint num=N;
    for(int i=2;i*i<=num;i++){
        while(N%i==0) N/=i,res.insert(i);
        if(N==1) break;
    }
    if(N!=1) res.insert(N);
    return res; 
}
 
lint powmod(lint a,lint b){
  return b?powmod(a*a%MOD,b/2)*(b%2?a:1)%MOD:1;
}

vector<bool> prime;
void primeinit(int N){
    if(N+1>prime.size()){
        prime.resize(N+1,true);
    } 
    prime[0] = false;
    prime[1] = false;
    for(int i=2; i*i<=N;i++) if(prime[i]) for(int j=2; i*j<=N;j++) prime[i*j]=false;
}

map<int,int> res;
set<int> apr;


int main(void){
    lint N;
    cin >> N;
    primeinit(N+1);
    lint cs=N;
    while(!prime[cs]) cs--;
    vector<lint> A;
    for(int i=1;i<=N;i++) if(i!=cs) A.push_back(i);
    N--;
    rep(i,N){
        res.clear();
        multiset<lint> cmp=decomp(A[i]);
        while(!cmp.empty()){
            auto itr=cmp.begin();
            res[*itr]=cmp.count(*itr);
            apr.insert(*itr);
            cmp.erase(*itr);
        }
        while(!apr.empty()){
            auto itr=apr.begin();
            cnt[*itr]=max(cnt[*itr],res[*itr]);
            apr.erase(itr);
        }
    }
    lint lcm=1;
    for(int i=2;i<=1000000;i++) lcm=(lcm*powmod(i,cnt[i]))%MOD;
    cout << lcm << endl;
}
0