結果

問題 No.1529 Constant Lcm
ユーザー maguroflymagurofly
提出日時 2021-06-04 21:06:05
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 455 ms / 3,000 ms
コード長 1,734 bytes
コンパイル時間 2,342 ms
コンパイル使用メモリ 211,428 KB
実行使用メモリ 14,996 KB
最終ジャッジ日時 2024-04-30 00:59:14
合計ジャッジ時間 7,867 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 384 ms
14,012 KB
testcase_11 AC 140 ms
7,680 KB
testcase_12 AC 8 ms
5,376 KB
testcase_13 AC 341 ms
12,544 KB
testcase_14 AC 194 ms
9,344 KB
testcase_15 AC 222 ms
9,472 KB
testcase_16 AC 162 ms
8,192 KB
testcase_17 AC 82 ms
6,016 KB
testcase_18 AC 90 ms
6,400 KB
testcase_19 AC 202 ms
9,472 KB
testcase_20 AC 445 ms
14,868 KB
testcase_21 AC 433 ms
14,988 KB
testcase_22 AC 438 ms
14,992 KB
testcase_23 AC 455 ms
14,996 KB
testcase_24 AC 446 ms
14,992 KB
testcase_25 AC 430 ms
14,868 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