結果

問題 No.3717 GCD LCM GCD
コンテスト
ユーザー KEYBO
提出日時 2026-09-18 22:53:47
言語 C++23(gnu拡張gcc16)
(gcc 16.1.0 + boost 1.92.0 + ACL)
コンパイル:
g++-16 -O2 -lm -std=gnu++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
WA  
実行時間 -
コード長 1,803 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 4,804 ms
コンパイル使用メモリ 391,928 KB
実行使用メモリ 92,160 KB
最終ジャッジ日時 2026-09-18 22:53:54
合計ジャッジ時間 6,709 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 7 WA * 1
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
#include <atcoder/all>
using namespace atcoder;

using ll = int64_t;
using ul = uint64_t;
using ld = long double;
using vi = vector<int>;
using vc = vector<char>;
using vs = vector<string>;
using vb = vector<bool>;
using vl = vector<ll>;
using vvi = vector<vi>;
using vvc = vector<vc>;
using vvb = vector<vb>;
using vvl = vector<vl>;
using mint = modint998244353;
using vm = vector<mint>;

void sieve(int N, vb &del, vi &div) {
  for (int p = 2; p <= N; p++) {
    if (del[p]) continue;
    div[p] = p;
    for (int a = p*2; a <= N; a += p) {
      del[a] = true;
      div[a] = min(div[a], p);
    }
  }
  return;
}

// pow_mod
ll pwmd(ll b, ll e, ll m) {
  if (e == 0) return 1;
  ll bef;
  if (e%2 == 0) bef = pwmd(b, e/2, m)%m;
  else bef = pwmd(b, e - 1, m)%m;
  if (e%2 == 0) return bef*bef%m;
  else return bef*b%m;
}

int main() {
  int lim = 1000000;
  vb del(lim + 1, false);
  vi div(lim + 1, 1e9);
  sieve(lim, del, div);
  ll N,K;
  cin >> N >> K;
  vl A(N);
  for (int i = 0; i < N; i++) {
    cin >> A[i];
  }
  vvi exps(lim + 1);
  for (int i = 0; i < N; i++) {
    if (A[i] == 1) continue;
    pair<int,int> cur = {-1, 0};
    while(A[i] > 1) {
      int pri = div[A[i]];
      if (cur.first != -1 && cur.first != pri) exps[cur.first].push_back(cur.second);
      if (cur.first != pri) cur.second = 0;
      cur.first = pri;
      cur.second++;
      A[i] /= pri;
    }
    exps[cur.first].push_back(cur.second);
  }
  ll mod = 998244353;
  mint ans = 1;
  for (ll i = 2; i <= lim; i++) {
    if (exps[i].size() == 0) continue;
    while(exps[i].size() < N) exps[i].push_back(0);
    sort(exps[i].begin(), exps[i].end());
    ll exp = exps[i][(N + K - 1)/K - 1];
    ans *= pwmd(i, exp, mod);
  }
  cout << ans.val() << endl;
  return 0;
}
0