結果
問題 | 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 longusing 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))#endifint 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;}