結果
| 問題 | No.2896 Monotonic Prime Factors |
| コンテスト | |
| ユーザー |
vjudge1
|
| 提出日時 | 2024-09-25 15:31:16 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 1,151 bytes |
| コンパイル時間 | 1,883 ms |
| コンパイル使用メモリ | 193,876 KB |
| 最終ジャッジ日時 | 2025-02-24 12:11:31 |
|
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 2 TLE * 1 -- * 15 |
ソースコード
#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, ct[N];
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;
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;
int sum = 0;
for (int i = 1; i <= cnt; i++) {
while (x % prime[i] == 0) ct[i]++, x /= prime[i];
sum += ct[i];
}
if (sum < y) cout << 0 << '\n';
else cout << C(sum - 1, y - 1) << '\n';
}
return 0;
}
vjudge1