結果
問題 | No.2249 GCDistance |
ユーザー |
|
提出日時 | 2023-03-21 19:09:26 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 642 ms / 5,000 ms |
コード長 | 952 bytes |
コンパイル時間 | 2,819 ms |
コンパイル使用メモリ | 248,384 KB |
実行使用メモリ | 162,296 KB |
最終ジャッジ日時 | 2024-09-18 14:27:34 |
合計ジャッジ時間 | 10,237 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 10 |
ソースコード
#include <bits/stdc++.h>#define rep(i, n) for (int i = 0; i < (n); ++i)using namespace std;using ll = long long;// totient tablevector<int> ps, pf, phi;void sieve(int mx) {pf.resize(mx+1);phi.resize(mx+1);phi[1] = 1;rep(i, mx+1) pf[i] = i;for (int i = 2; i <= mx; ++i) {if (pf[i] == i) {ps.push_back(i);phi[i] = i-1;}rep(j, ps.size()) {int x = ps[j]*i;if (x > mx) break;pf[x] = ps[j];if (i%ps[j] == 0) {phi[x] = phi[i]*ps[j];break;}phi[x] = phi[i]*(ps[j]-1);}}}int main() {sieve(1e7);vector<ll> s(1e7+1);rep(i, 1e7) s[i+1] = s[i]+phi[i+1];int t;cin >> t;while (t--) {int n;cin >> n;ll ans = (ll)n*(n-1) - s[n] + 1;cout << ans << '\n';}return 0;}