結果

問題 No.1529 Constant Lcm
ユーザー maguroflymagurofly
提出日時 2021-06-04 21:06:05
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 400 ms / 3,000 ms
コード長 1,734 bytes
コンパイル時間 1,984 ms
コンパイル使用メモリ 210,116 KB
実行使用メモリ 15,008 KB
最終ジャッジ日時 2023-08-12 10:46:46
合計ジャッジ時間 6,868 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 2 ms
4,384 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 1 ms
4,376 KB
testcase_05 AC 1 ms
4,384 KB
testcase_06 AC 1 ms
4,376 KB
testcase_07 AC 2 ms
4,380 KB
testcase_08 AC 1 ms
4,380 KB
testcase_09 AC 2 ms
4,380 KB
testcase_10 AC 348 ms
13,708 KB
testcase_11 AC 125 ms
7,368 KB
testcase_12 AC 7 ms
4,380 KB
testcase_13 AC 298 ms
12,488 KB
testcase_14 AC 171 ms
9,088 KB
testcase_15 AC 194 ms
9,428 KB
testcase_16 AC 150 ms
7,844 KB
testcase_17 AC 79 ms
5,840 KB
testcase_18 AC 83 ms
6,272 KB
testcase_19 AC 186 ms
9,316 KB
testcase_20 AC 400 ms
14,884 KB
testcase_21 AC 390 ms
14,920 KB
testcase_22 AC 386 ms
14,768 KB
testcase_23 AC 384 ms
14,724 KB
testcase_24 AC 369 ms
15,008 KB
testcase_25 AC 359 ms
14,732 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define rep(i, l, r) for (auto i = (l); i < (r); i++)
#define chmin(dest, src) if ((dest) > (src)) dest = (src)
#define chmax(dest, src) if ((dest) < (src)) dest = (src)
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const ll MOD = 998244353;

// map<int, int> prime_division(int n) {
//   map<int, int> pd;
//   int i = 2;
//   while (i * i <= n) {
//     while (i % n == 0) {
//       pd[i]++;
//       n /= i;
//     }
//     i++;
//   }
//   if (n > 1) pd[n]++;
//   return pd;
// }

ll pow_mod(ll x, int e) {
  ll r = 1;
  while (e > 0) {
    if (e & 1) r = r * x % MOD;
    x = x * x % MOD;
    e >>= 1;
  }
  return r;
}

int main() {
  int N; cin >> N;

  vector<int> lpf(N + 1);
  vector<int> primes;
  rep (d, 2, N) {
    if (!lpf[d]) {
      lpf[d] = d;
      primes.push_back(d);
    }
    for (int p : primes) {
      if (p * d > N || p > lpf[d]) break;
      lpf[p * d] = p;
    }
  }

  auto prime_division = [lpf](int n, map<int, int> &pd) {
    while (n > 1) {
      int d = lpf[n];
      n /= d;
      pd[d]++;
    }
  };

  // vector<vector<int>> pds(N, map<int, int>());
  // rep (i, 1, N) {
  //   prime_division(i, pds[i]);
  // }

  map<int, int> lcm_pd;
  rep (i, 1, (N + 2) / 2) {
    int j = N - i;
    map<int, int> pd;
    // for (auto &e : pds[i]) pd[e.first] += e.second;
    // for (auto &e : pds[j]) pd[e.first] += e.second;
    // for (int p : pds[i]) pd[p]++;
    // for (int p : pds[j]) pd[p]++;
    prime_division(i, pd);
    prime_division(j, pd);
    for (auto &e : pd) chmax(lcm_pd[e.first], e.second);
  }

  ll ans = 1;
  for (auto &e : lcm_pd) {
    ans *= pow_mod(e.first, e.second);
    ans %= MOD;
  }

  if (ans < 0) ans += MOD;

  cout << ans << endl;

  return 0;
}
0