結果
問題 | No.1396 Giri |
ユーザー | Example0911 |
提出日時 | 2021-02-14 21:42:24 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 34 ms / 2,000 ms |
コード長 | 1,540 bytes |
コンパイル時間 | 2,624 ms |
コンパイル使用メモリ | 199,656 KB |
最終ジャッジ日時 | 2025-01-18 20:49:39 |
ジャッジサーバーID (参考情報) |
judge1 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 23 |
ソースコード
#include "bits/stdc++.h" using namespace std; using ll = long long; using P = pair<ll, ll>; const ll INF = (1LL << 61); ll mod = 998244353; bool prime(ll n) { if (n == 1)return true; for (ll i = 2; i * i <= n; i++) { if (n % i == 0)return false; } return true; } struct Sieve { int n; vector<int> f, primes; // n以下の素数を列挙する Sieve(int n = 1) :n(n), f(n + 1) { f[0] = f[1] = -1; for (ll i = 2; i <= n; ++i) { if (f[i]) continue; primes.push_back(i); f[i] = i; for (ll j = i * i; j <= n; j += i) { if (!f[j]) f[j] = i; } } } // xが素数かどうか判定する bool isPrime(int x) { return f[x] == x; } // 素因数を全列挙 vector<int> factorList(int x) { vector<int> res; while (x != 1) { res.push_back(f[x]); x /= f[x]; } return res; } // ランレングス圧縮 vector<P> factor(int x) { vector<int> fl = factorList(x); if (fl.size() == 0) return {}; vector<P> res(1, P(fl[0], 0)); for (int p : fl) { if (res.back().first == p) { res.back().second++; } else { res.emplace_back(p, 1); } } return res; } }; signed main() { ios::sync_with_stdio(false); cin.tie(0); ll N; cin >> N; map<ll, ll>mp; ll tmp = 0; for (int i = N; ; i--) { if (prime(i)) { tmp = i; break; } } Sieve s(1000010); vector<int> now = s.primes; ll tmp2 = N; ll ans = 1; for (auto x : now) { if (x == tmp)continue; N = tmp2; N /= x; while (N != 0) { ans *= x; ans %= mod; N /= x; } } cout << ans << endl; return 0; }