結果
問題 | 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;}