結果
問題 | No.2896 Monotonic Prime Factors |
ユーザー | srjywrdnprkt |
提出日時 | 2024-11-13 15:34:19 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 1,281 bytes |
コンパイル時間 | 2,573 ms |
コンパイル使用メモリ | 214,208 KB |
実行使用メモリ | 13,048 KB |
最終ジャッジ日時 | 2024-11-13 15:34:28 |
合計ジャッジ時間 | 7,793 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 17 ms
11,776 KB |
testcase_01 | AC | 16 ms
11,776 KB |
testcase_02 | AC | 16 ms
11,776 KB |
testcase_03 | AC | 16 ms
11,776 KB |
testcase_04 | AC | 163 ms
11,776 KB |
testcase_05 | RE | - |
testcase_06 | AC | 196 ms
11,776 KB |
testcase_07 | AC | 160 ms
11,864 KB |
testcase_08 | AC | 176 ms
11,776 KB |
testcase_09 | RE | - |
testcase_10 | AC | 111 ms
11,756 KB |
testcase_11 | AC | 37 ms
11,800 KB |
testcase_12 | AC | 27 ms
11,804 KB |
testcase_13 | AC | 112 ms
11,764 KB |
testcase_14 | AC | 141 ms
11,776 KB |
testcase_15 | AC | 21 ms
11,648 KB |
testcase_16 | AC | 56 ms
11,776 KB |
testcase_17 | AC | 32 ms
11,904 KB |
testcase_18 | AC | 162 ms
11,784 KB |
testcase_19 | AC | 34 ms
11,776 KB |
ソースコード
#include <bits/stdc++.h> #include <atcoder/modint> using namespace std; using namespace atcoder; using ll = long long; using mint = modint998244353; vector<ll> spf; map<ll, ll> prime; void osa_k(ll n){ spf.resize(n+1); for (ll i=0; i<=n; i++) spf[i] = i; for (ll i=2; i*i<=n; i++){ if (spf[i] == i){ for (ll j=2; i*j <= n; j++){ spf[i*j] = min(spf[i*j], i); } } } } void prime_factor(ll n){ prime.clear(); while(n != 1){ prime[spf[n]]++; n /= spf[n]; } } const ll MAX=1000000; vector<mint> f, finv; void init(){ f.resize(MAX+1); finv.resize(MAX+1); f[0] = 1; for (int i=1; i<=MAX; i++) f[i] = f[i-1]*i; finv[MAX] = f[MAX].inv(); for (int i=MAX-1; i>=0; i--) finv[i] = finv[i+1] * (i+1); } mint C(ll n, ll k){ if (n < k || k < 0) return 0; return f[n] * finv[k] * finv[n-k] ; } int main(){ cin.tie(nullptr); ios_base::sync_with_stdio(false); init(); osa_k(1e5); ll Q, A, B, cnt=0; cin >> Q; //全素因数を各要素に1つずつ配る。 while(Q--){ cin >> A >> B; prime_factor(A); for (auto [p, e] : prime) cnt += e; cout << C(cnt-1, B-1).val() << endl; } return 0; }