結果

問題 No.1396 Giri
ユーザー chocopuu
提出日時 2021-02-21 10:34:00
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 75 ms / 2,000 ms
コード長 1,809 bytes
コンパイル時間 1,921 ms
コンパイル使用メモリ 192,756 KB
最終ジャッジ日時 2025-01-19 02:57:21
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
//#include "atcoder/all"
//using namespace atcoder;
#define int long long
#define REP(i, n) for (int64_t i = 0; i < (int)n; ++i)
#define RREP(i, n) for (int64_t i = (int)n - 1; i >= 0; --i)
#define FOR(i, s, n) for (int64_t i = s; i < (int)n; ++i)
#define RFOR(i, s, n) for (int64_t i = (int)n - 1; i >= s; --i)
#define ALL(a) a.begin(), a.end()
#define IN(a, x, b) (a <= x && x < b)
template<class T>istream&operator >>(istream&is,vector<T>&vec){for(T&x:vec)is>>x;return is;}
template<class T>inline void out(T t){cout << t << "\n";}
template<class T,class... Ts>inline void out(T t,Ts... ts){cout << t << " ";out(ts...);}
template<class T>inline bool CHMIN(T&a,T b){if(a > b){a = b;return true;}return false;}
template<class T>inline bool CHMAX(T&a,T b){if(a < b){a = b;return true;}return false;}
constexpr int64_t INF = 1000000000000000000; // 10^18

const int M = 2750131;
int P[200000];
int invP[M + 1];
int minP[M + 1];
bool isP[M + 1];
 
void Eratosthenes() {
    int p = 0;
    REP(i, M + 1) isP[i] = true, minP[i] = -1;
    isP[0] = isP[1] = false;
    minP[0] = 0; minP[1] = 1;
    FOR(i, 2, M + 1) {
        if (isP[i]) {
            minP[i] = i;
            P[p] = i;
            invP[i] = p++;
            for (int j = i * 2; j < M; j += i) {
                isP[j] = false;
                if (minP[j] == -1) minP[j] = i;
            }
        }
    }
}

signed main() {
    int N;
    cin >> N;
    Eratosthenes();
    int flg = true;
    int ans = 1;
    constexpr int MOD = 998244353;
    RFOR(i, 1, N + 1) {
        if(flg && isP[i]) flg = false;
        else if(isP[i]) {
            int powP = 1;
            while(powP * i <= N) {
                powP *= i;
            }
            ans = (ans * powP) % MOD;
        }
    }
    out(ans);
}
0