結果
問題 | No.2896 Monotonic Prime Factors |
ユーザー | tnakao0123 |
提出日時 | 2024-09-26 11:49:50 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 329 ms / 2,000 ms |
コード長 | 2,745 bytes |
コンパイル時間 | 934 ms |
コンパイル使用メモリ | 63,124 KB |
実行使用メモリ | 17,024 KB |
最終ジャッジ日時 | 2024-09-26 11:50:00 |
合計ジャッジ時間 | 8,129 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 250 ms
16,896 KB |
testcase_01 | AC | 252 ms
16,896 KB |
testcase_02 | AC | 251 ms
16,896 KB |
testcase_03 | AC | 288 ms
16,840 KB |
testcase_04 | AC | 319 ms
17,004 KB |
testcase_05 | AC | 315 ms
16,992 KB |
testcase_06 | AC | 300 ms
17,024 KB |
testcase_07 | AC | 329 ms
16,896 KB |
testcase_08 | AC | 325 ms
16,896 KB |
testcase_09 | AC | 299 ms
16,896 KB |
testcase_10 | AC | 281 ms
16,896 KB |
testcase_11 | AC | 267 ms
16,896 KB |
testcase_12 | AC | 260 ms
16,964 KB |
testcase_13 | AC | 283 ms
16,896 KB |
testcase_14 | AC | 293 ms
16,992 KB |
testcase_15 | AC | 282 ms
16,948 KB |
testcase_16 | AC | 270 ms
16,896 KB |
testcase_17 | AC | 263 ms
16,848 KB |
testcase_18 | AC | 322 ms
17,024 KB |
testcase_19 | AC | 287 ms
16,896 KB |
ソースコード
/* -*- coding: utf-8 -*- * * 2896.cc: No.2896 Monotonic Prime Factors - yukicoder */ #include<cstdio> #include<vector> #include<algorithm> using namespace std; /* constant */ const int MAX_A = 100000; const int MAX_N = MAX_A * 18; const int MAX_P = 320; const int MOD = 998244353; /* typedef */ using vb = vector<bool>; using vi = vector<int>; template<const int MOD> struct MI { int v; MI(): v() {} MI(int _v): v(_v % MOD) { if (v < 0) v += MOD; } MI(long long _v): v(_v % MOD) { if (v < 0) v += MOD; } explicit operator int() const { return v; } MI operator+(const MI m) const { return MI(v + m.v); } MI operator-(const MI m) const { return MI(v + MOD - m.v); } MI operator*(const MI m) const { return MI((long long)v * m.v); } MI &operator+=(const MI m) { return (*this = *this + m); } MI &operator-=(const MI m) { return (*this = *this - m); } MI &operator*=(const MI m) { return (*this = *this * m); } bool operator==(const MI m) const { return v == m.v; } bool operator!=(const MI m) const { return v != m.v; } MI pow(int n) const { // a^n % MOD MI pm = 1, a = *this; while (n > 0) { if (n & 1) pm *= a; a *= a; n >>= 1; } return pm; } MI inv() const { return pow(MOD - 2); } MI operator/(const MI m) const { return *this * m.inv(); } MI &operator/=(const MI m) { return (*this = *this / m); } }; using mi = MI<MOD>; using vmi = vector<mi>; /* global variables */ vb primes; vi pnums; vmi fs, invfs; /* subroutines */ int gen_primes(int maxp) { primes.assign(maxp + 1, true); primes[0] = primes[1] = false; int p; for (p = 2; p * p <= maxp; p++) if (primes[p]) { pnums.push_back(p); for (int q = p * p; q <= maxp; q += p) primes[q] = false; } for (; p <= maxp; p++) if (primes[p]) pnums.push_back(p); return (int)pnums.size(); } int pdnum(int n) { int sum = 0; for (auto pi: pnums) { if (pi * pi > n) { if (n > 1) sum++; return sum; } while (n % pi == 0) n /= pi, sum++; } return -1; } inline mi nck(int n, int k) { // nCk % MOD if (n < k || k < 0) return 0; return fs[n] * invfs[n - k] * invfs[k]; } inline mi nhk(int n, int k) { // nHk % MOD if (n < 0 || k < 0) return 0; return nck(n + k - 1, k); } void prepare_fracs(int n) { fs.resize(n + 1), invfs.resize(n + 1); fs[0] = invfs[0] = 1; for (int i = 1; i <= n; i++) { fs[i] = fs[i - 1] * i; invfs[i] = fs[i].inv(); } } /* main */ int main() { gen_primes(MAX_P); prepare_fracs(MAX_N); int qn; scanf("%d", &qn); int l = 0; while (qn--) { int a, b; scanf("%d%d", &a, &b); l += pdnum(a); mi c = nhk(b, l - b); printf("%d\n", (int)c); } return 0; }