結果
問題 | No.1529 Constant Lcm |
ユーザー | outline |
提出日時 | 2021-06-04 20:16:40 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 813 ms / 3,000 ms |
コード長 | 4,567 bytes |
コンパイル時間 | 1,484 ms |
コンパイル使用メモリ | 145,648 KB |
実行使用メモリ | 10,880 KB |
最終ジャッジ日時 | 2024-04-29 23:25:20 |
合計ジャッジ時間 | 10,422 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 3 ms
5,376 KB |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | AC | 2 ms
5,376 KB |
testcase_05 | AC | 2 ms
5,376 KB |
testcase_06 | AC | 2 ms
5,376 KB |
testcase_07 | AC | 2 ms
5,376 KB |
testcase_08 | AC | 2 ms
5,376 KB |
testcase_09 | AC | 2 ms
5,376 KB |
testcase_10 | AC | 705 ms
10,240 KB |
testcase_11 | AC | 254 ms
6,144 KB |
testcase_12 | AC | 13 ms
5,376 KB |
testcase_13 | AC | 594 ms
9,112 KB |
testcase_14 | AC | 344 ms
7,040 KB |
testcase_15 | AC | 397 ms
7,424 KB |
testcase_16 | AC | 310 ms
6,656 KB |
testcase_17 | AC | 155 ms
5,376 KB |
testcase_18 | AC | 170 ms
5,376 KB |
testcase_19 | AC | 372 ms
7,296 KB |
testcase_20 | AC | 804 ms
10,716 KB |
testcase_21 | AC | 808 ms
10,880 KB |
testcase_22 | AC | 788 ms
10,812 KB |
testcase_23 | AC | 813 ms
10,764 KB |
testcase_24 | AC | 787 ms
10,752 KB |
testcase_25 | AC | 749 ms
10,828 KB |
ソースコード
#include <iostream> #include <vector> #include <algorithm> #include <cmath> #include <queue> #include <string> #include <map> #include <set> #include <stack> #include <tuple> #include <deque> #include <array> #include <numeric> #include <bitset> #include <iomanip> #include <cassert> #include <chrono> #include <random> #include <limits> #include <iterator> #include <functional> #include <sstream> #include <fstream> #include <complex> #include <cstring> #include <unordered_map> #include <unordered_set> using namespace std; // #pragma GCC optimize("O3") // #pragma GCC optimize("unroll-loops") // #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,avx2,tune=native") // #pragma GCC target("avx512f,avx512dq,avx512cd,avx512bw,avx512vl") using ll = long long; constexpr int INF = 1001001001; // constexpr int mod = 1000000007; constexpr int mod = 998244353; template<class T> inline bool chmax(T& x, T y){ if(x < y){ x = y; return true; } return false; } template<class T> inline bool chmin(T& x, T y){ if(x > y){ x = y; return true; } return false; } struct mint { int x; mint() : x(0) {} mint(int64_t y) : x(y >= 0 ? y % mod : (mod - (-y) % mod) % mod) {} mint& operator+=(const mint& p){ if((x += p.x) >= mod) x -= mod; return *this; } mint& operator-=(const mint& p){ if((x -= p.x) < 0) x += mod; return *this; } mint& operator*=(const mint& p){ x = (int)(1LL * x * p.x % mod); return *this; } mint& operator/=(const mint& p){ *this *= p.inverse(); return *this; } mint operator-() const { return mint(-x); } mint operator+(const mint& p) const { return mint(*this) += p; } mint operator-(const mint& p) const { return mint(*this) -= p; } mint operator*(const mint& p) const { return mint(*this) *= p; } mint operator/(const mint& p) const { return mint(*this) /= p; } bool operator==(const mint& p) const { return x == p.x; } bool operator!=(const mint& p) const { return x != p.x; } mint pow(int64_t n) const { mint res = 1, mul = x; while(n > 0){ if(n & 1) res *= mul; mul *= mul; n >>= 1; } return res; } // x^(a^b) // warning : x と mod は互いに素 // x % mod == 0 かつ a > 0 はこれを呼び出さずに 0 を返すように処理 mint pow2(int64_t a, int64_t b) const { if(b == 0) return *this; if((a %= mod - 1) == 0) return mint(1); int64_t n = 1; while(b > 0){ if(b & 1) (n *= a) %= mod - 1; (a *= a) %= mod - 1; b >>= 1; } return pow(n); } mint inverse() const { return pow(mod - 2); } friend ostream& operator<<(ostream& os, const mint& p){ return os << p.x; } friend istream& operator>>(istream& is, mint& p){ int64_t val; is >> val; p = mint(val); return is; } }; struct ExtendedSieve{ vector<int> table; const int n; ExtendedSieve(int n) : n(n) { table.assign(n + 1, 0); table[0] = table[1] = 1; for(int i = 2; i <= n; ++i){ if(table[i] > 0) continue; table[i] = i; for(int j = i + i; j <= n; j += i){ if(table[j] == 0) table[j] = i; } } } int operator[](const int& k){ return table[k]; } bool is_prime(int x){ if(x < 2) return false; return x == table[x]; } map<int, int> prime_factorization(int x){ map<int, int> res; while(x > 1){ ++res[table[x]]; x /= table[x]; } return res; } vector<int> prime_factorization2(int x){ vector<int> res; while(x > 1){ if(!res.size() || res.back() != table[x]){ res.emplace_back(table[x]); } x /= table[x]; } return res; } }; int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); int N; cin >> N; ExtendedSieve sieve(N); map<int, int> primes; for(int i = 1; i < N; ++i){ auto p1 = sieve.prime_factorization(i); auto p2 = sieve.prime_factorization(N - i); for(auto& [p, c] : p2){ chmax(primes[p], (p1[p] += c)); } } mint ans = 1; for(auto [p, c] : primes){ ans *= mint(p).pow(c); } cout << ans << endl; return 0; }