結果

問題 No.2249 GCDistance
ユーザー Meet BrahmbhattMeet Brahmbhatt
提出日時 2023-03-18 05:52:32
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,397 ms / 5,000 ms
コード長 1,010 bytes
コンパイル時間 2,252 ms
コンパイル使用メモリ 203,788 KB
実行使用メモリ 160,768 KB
最終ジャッジ日時 2024-09-18 13:04:47
合計ジャッジ時間 20,440 ms
ジャッジサーバーID
(参考情報)
judge2 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1,334 ms
160,768 KB
testcase_01 AC 1,389 ms
160,640 KB
testcase_02 AC 1,397 ms
160,640 KB
testcase_03 AC 1,374 ms
160,640 KB
testcase_04 AC 1,338 ms
160,640 KB
testcase_05 AC 1,392 ms
160,768 KB
testcase_06 AC 1,377 ms
160,768 KB
testcase_07 AC 1,394 ms
160,640 KB
testcase_08 AC 1,316 ms
160,768 KB
testcase_09 AC 1,385 ms
160,640 KB
testcase_10 AC 1,378 ms
160,640 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

/**
 *  - Meet Brahmbhatt
 *  - Hard work always pays off
**/
#include"bits/stdc++.h"
using namespace std;

#ifdef MeetBrahmbhatt
#include "debug.h"
#else
#define dbg(...) 72
#endif

#define endl "\n"
#define int long long
#define sz(x) (int)(x.size())

const long long INF = 4e18;
const int32_t M = 1e9 + 7;
const int32_t MM = 998244353;

const int32_t N = 1e7 + 5;
vector<int> phi(N);
vector<bool> P(N, 1);
vector<int> dp(N);

void pre() {
  iota(phi.begin(), phi.end(), 0);
  for (int i = 2; i < N; i++) {
    if (P[i]) {
      phi[i]--;
      for (int j = 2 * i; j < N; j += i) {
        P[j] = 0;
        phi[j] /= i;
        phi[j] *= (i - 1);
      }
    }
  }
  dp[2] = 1;
  for (int i = 3; i < N; i++) {
    dp[i] = dp[i - 1] + (i - 1 - phi[i]) * 2 + phi[i];
  }
}

void solve() {
  int n;
  cin >> n;
  cout << dp[n] << endl;
}

signed main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);
  cout << fixed << setprecision(9);
  int tt = 1;
  pre();
  cin >> tt;
  while (tt--) solve();
  return 0;
}
0