結果

問題 No.2896 Monotonic Prime Factors
ユーザー vjudge1vjudge1
提出日時 2024-09-25 14:46:31
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 897 bytes
コンパイル時間 2,151 ms
コンパイル使用メモリ 203,992 KB
実行使用メモリ 26,212 KB
最終ジャッジ日時 2024-09-25 14:46:37
合計ジャッジ時間 6,361 ms
ジャッジサーバーID
(参考情報)
judge1 / judge4
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 25 ms
16,128 KB
testcase_01 AC 21 ms
16,128 KB
testcase_02 AC 21 ms
16,128 KB
testcase_03 AC 22 ms
16,128 KB
testcase_04 AC 216 ms
16,128 KB
testcase_05 RE -
testcase_06 AC 229 ms
24,516 KB
testcase_07 AC 223 ms
16,256 KB
testcase_08 AC 217 ms
16,128 KB
testcase_09 RE -
testcase_10 AC 128 ms
22,340 KB
testcase_11 AC 47 ms
24,392 KB
testcase_12 AC 35 ms
24,516 KB
testcase_13 AC 138 ms
24,388 KB
testcase_14 AC 165 ms
24,516 KB
testcase_15 AC 28 ms
26,212 KB
testcase_16 AC 72 ms
24,260 KB
testcase_17 AC 42 ms
24,384 KB
testcase_18 AC 207 ms
24,384 KB
testcase_19 AC 46 ms
24,388 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<bits/stdc++.h>

using namespace std;

const int N = 2e6 + 5, mod = 998244353;

int f[N], g[N], inv[N], p[N], a, b, cnt, vis[N], q;
vector<int>e;

int C(int n, int m){
  if(m > n)return 0;
  return 1ll * f[n] * g[m] % mod * g[n - m] % mod;
}

int main(){
  f[0] = g[0] = inv[1] = 1;
  for(int i = 1; i <= 1000000; i++){
    f[i] = 1ll * f[i - 1] * i % mod;
    if(i > 1){
      inv[i] = 1ll * inv[mod % i] * (mod - mod / i) % mod;
    }
    g[i] = 1ll * g[i - 1] * inv[i] % mod;
  }
  for(int i = 2; i <= 100000; ++i){
    if(!vis[i]){
      e.push_back(i);
      p[i] = i;
    }
    for(auto v : e){
      if(1ll * v * i > 100000)break;
      vis[v * i] = 1;
      p[v * i] = v;
      if(i % v == 0)break;
    }
  }
  cin >> q;
  while(q--){
    cin >> a >> b;
    for(; a > 1; a /= p[a], cnt++){
    }
    assert(cnt <= 1000000);
    cout << C(cnt - 1, b - 1) << '\n';
  }
  return 0;
}
0