結果
問題 |
No.2896 Monotonic Prime Factors
|
ユーザー |
![]() |
提出日時 | 2024-09-25 15:33:18 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 487 ms / 2,000 ms |
コード長 | 1,134 bytes |
コンパイル時間 | 1,860 ms |
コンパイル使用メモリ | 194,064 KB |
最終ジャッジ日時 | 2025-02-24 12:11:46 |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 18 |
ソースコード
#include <bits/stdc++.h> using namespace std; using ll = long long; const int N = 1e5 + 10, M = 1e7 + 10, mod = 998244353; int q, prime[N], cnt, id[M], sum; ll inv[M], fact[M], factinv[M]; bool f[N]; ll C(int n, int m) { return fact[n] * factinv[m] % mod * factinv[n - m] % mod; } int main() { ios::sync_with_stdio(0), cin.tie(0); f[1] = 1; for (int i = 2; i < N; i++) { if (!f[i]) prime[++cnt] = i, id[i] = cnt; for (int j = 1; j <= cnt && 1ll * prime[j] * i < N; j++) { f[i * prime[j]] = 1; } } inv[1] = 1; for (int i = 2; i < M; i++) { inv[i] = (mod - mod / i) * inv[mod % i] % mod; } fact[0] = factinv[0] = 1; for (int i = 1; i < M; i++) { fact[i] = fact[i - 1] * i % mod; factinv[i] = factinv[i - 1] * inv[i] % mod; } cin >> q; while (q--) { int x, y; cin >> x >> y; for (int i = 2; i * i <= x; i++) { while (x % i == 0) x /= i, sum++; } if (x > 1) sum++; if (sum < y) cout << 0 << '\n'; else cout << C(sum - 1, y - 1) << '\n'; } return 0; }