結果
問題 | No.2896 Monotonic Prime Factors |
ユーザー | momoyuu |
提出日時 | 2024-09-20 21:40:27 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 224 ms / 2,000 ms |
コード長 | 980 bytes |
コンパイル時間 | 992 ms |
コンパイル使用メモリ | 107,744 KB |
実行使用メモリ | 19,328 KB |
最終ジャッジ日時 | 2024-09-20 21:40:40 |
合計ジャッジ時間 | 4,274 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 18 |
ソースコード
#include<iostream> #include<vector> #include<algorithm> using namespace std; using ll = long long; #include<atcoder/modint> using mint = atcoder::modint998244353; const int B = 2e6; mint fac[B+1],ifac[B+1]; #include<map> int main(){ cin.tie(nullptr); ios::sync_with_stdio(false); int q; cin>>q; fac[0] = 1; for(int i = 1;i<=B;i++) fac[i] = fac[i-1] * i; ifac[B] = fac[B].inv(); for(int i = B-1;i>=0;i--) ifac[i] = ifac[i+1] * (i+1); const int want = 1e5 + 10; vector<int> p(want+1,1); for(int i = 2;i<=want;i++){ if(p[i]!=1) continue; for(int j = i;j<=want;j+=i) p[j] = i; } ll all = 0; for(int i = 0;i<q;i++){ ll a,b; cin>>a>>b; while(a!=1){ all++; a /= p[a]; } mint ans = 0; if(all>=b){ ll use = all - b; ans = fac[use+b-1] * ifac[use] * ifac[b-1]; } cout<<ans.val()<<endl; } }