結果
問題 |
No.2380 Sylow P-subgroup
|
ユーザー |
![]() |
提出日時 | 2023-07-14 21:41:27 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,860 bytes |
コンパイル時間 | 1,808 ms |
コンパイル使用メモリ | 192,880 KB |
最終ジャッジ日時 | 2025-02-15 13:48:58 |
ジャッジサーバーID (参考情報) |
judge2 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 14 |
ソースコード
#include<bits/stdc++.h> using namespace std; using ll = long long; #include <cstdint> template <std::uint_fast64_t Modulus> class modint { using u64 = std::uint_fast64_t; public: u64 a; constexpr modint(const u64 x = 0) noexcept : a(x % Modulus) {} constexpr u64 &value() noexcept { return a; } constexpr const u64 &value() const noexcept { return a; } constexpr modint operator+(const modint rhs) const noexcept { return modint(*this) += rhs; } constexpr modint operator-(const modint rhs) const noexcept { return modint(*this) -= rhs; } constexpr modint operator*(const modint rhs) const noexcept { return modint(*this) *= rhs; } constexpr modint operator/(const modint rhs) const noexcept { return modint(*this) /= rhs; } constexpr modint &operator+=(const modint rhs) noexcept { a += rhs.a; if (a >= Modulus) { a -= Modulus; } return *this; } constexpr modint &operator-=(const modint rhs) noexcept { if (a < rhs.a) { a += Modulus; } a -= rhs.a; return *this; } constexpr modint &operator*=(const modint rhs) noexcept { a = a * rhs.a % Modulus; return *this; } constexpr modint &operator/=(modint rhs) noexcept { u64 exp = Modulus - 2; while (exp) { if (exp % 2) { *this *= rhs; } rhs *= rhs; exp /= 2; } return *this; } }; const ll mod = 998244353LL; long long power(ll x, ll n){ long long ret = 1; while(n){ if(n&1){ ret *= x; ret %= mod; } x *= x; x %= mod; n >>= 1; } return ret; } int main(){ ll n, p; cin >> n >> p; modint <mod> ans; ans.a = 0; modint <mod> tmp = 1; while(n > 0){ n /= p; ans.a += n; } ans.a = power(p, ans.a); cout << ans.a << endl; }