結果

問題 No.1529 Constant Lcm
ユーザー maguroflymagurofly
提出日時 2021-06-04 21:06:05
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 439 ms / 3,000 ms
コード長 1,734 bytes
コンパイル時間 2,273 ms
コンパイル使用メモリ 213,068 KB
実行使用メモリ 15,164 KB
最終ジャッジ日時 2024-11-19 11:16:58
合計ジャッジ時間 7,437 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,816 KB
testcase_02 AC 2 ms
6,820 KB
testcase_03 AC 2 ms
6,820 KB
testcase_04 AC 2 ms
6,816 KB
testcase_05 AC 2 ms
6,816 KB
testcase_06 AC 2 ms
6,816 KB
testcase_07 AC 2 ms
6,816 KB
testcase_08 AC 2 ms
6,816 KB
testcase_09 AC 2 ms
6,816 KB
testcase_10 AC 380 ms
14,000 KB
testcase_11 AC 136 ms
7,680 KB
testcase_12 AC 7 ms
6,820 KB
testcase_13 AC 325 ms
12,672 KB
testcase_14 AC 184 ms
9,344 KB
testcase_15 AC 214 ms
9,728 KB
testcase_16 AC 162 ms
8,320 KB
testcase_17 AC 82 ms
6,816 KB
testcase_18 AC 88 ms
6,820 KB
testcase_19 AC 203 ms
9,472 KB
testcase_20 AC 438 ms
14,964 KB
testcase_21 AC 429 ms
14,960 KB
testcase_22 AC 422 ms
15,164 KB
testcase_23 AC 439 ms
15,040 KB
testcase_24 AC 432 ms
15,092 KB
testcase_25 AC 414 ms
15,100 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