結果
問題 | No.2576 LCM Pattern |
ユーザー |
![]() |
提出日時 | 2023-12-13 16:43:25 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,158 bytes |
コンパイル時間 | 5,783 ms |
コンパイル使用メモリ | 252,468 KB |
最終ジャッジ日時 | 2025-02-18 10:44:58 |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 23 |
ソースコード
#include <bits/stdc++.h> using namespace std; #include <atcoder/all> using namespace atcoder; #define rep(i, n) for(int i=0;i<(n);++i) #define rep1(i, n) for(int i=1;i<=(n);i++) #define ll long long using mint = modint998244353; using P = pair<ll,ll>; using lb = long double; using T = tuple<ll, ll, ll>; #ifdef LOCAL # include <debug_print.hpp> # define dbg(...) debug_print::multi_print(#__VA_ARGS__, __VA_ARGS__) #else # define dbg(...) (static_cast<void>(0)) #endif int main() { int n, m; cin >> n >> m; vector<pair<int, int>> ps; for(int i=2;i*i<=m;i++){ if(m%i) continue; ps.emplace_back(i,0); while(m%i==0){ ps.back().second++; m/=i; } } if(m!=1) ps.emplace_back(m,1); mint ans = 0; int N = ps.size(); rep(i,1<<N){ mint now = 1; rep(j,N) { if(i>>j&1) { now *= ps[j].second; } else{ now *= ps[j].second+1; } } if(__builtin_popcount(i)%2==0) ans += now.pow(n); else ans -= now.pow(n); } cout << ans.val() << endl; return 0; }