結果
問題 | No.2249 GCDistance |
ユーザー | roaris |
提出日時 | 2023-03-17 22:58:05 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 2,599 ms / 5,000 ms |
コード長 | 1,323 bytes |
コンパイル時間 | 2,308 ms |
コンパイル使用メモリ | 208,572 KB |
実行使用メモリ | 120,408 KB |
最終ジャッジ日時 | 2023-10-18 15:59:29 |
合計ジャッジ時間 | 35,015 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2,524 ms
120,408 KB |
testcase_01 | AC | 2,594 ms
120,408 KB |
testcase_02 | AC | 2,599 ms
120,408 KB |
testcase_03 | AC | 2,595 ms
120,408 KB |
testcase_04 | AC | 2,552 ms
120,408 KB |
testcase_05 | AC | 2,594 ms
120,408 KB |
testcase_06 | AC | 2,583 ms
120,408 KB |
testcase_07 | AC | 2,579 ms
120,408 KB |
testcase_08 | AC | 2,518 ms
120,408 KB |
testcase_09 | AC | 2,594 ms
120,408 KB |
testcase_10 | AC | 2,580 ms
120,408 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; #define rep(i, n) for (int i=0; i<n; i++) #define pb push_back typedef long long ll; struct Fast_factorize { vector<int> p; Fast_factorize(int n) { p.resize(n+1); fill(p.begin(), p.end(), -1); for (int i=2; i<=n; i++) if (p[i]==-1) { for (int j=i; j<=n; j+=i) { if (p[j]==-1) p[j] = i; } } } vector<int> factorize(int n) { vector<int> res; while (n>1) { res.pb(p[n]); n /= p[n]; } return res; } }; int main() { cin.tie(0); ios::sync_with_stdio(false); int MAX = 10000010; Fast_factorize ff(MAX); ll phi[MAX]; phi[0] = 0ll; for (int i=1; i<MAX; i++) { phi[i] = (ll)(i); vector<int> ps = ff.factorize(i); sort(ps.begin(), ps.end()); ps.erase(unique(ps.begin(), ps.end()), ps.end()); int x = i; for (int p : ps) { phi[i] -= phi[i]/p; while (x%p==0) x /= p; } if (x>1) phi[i] -= phi[i]/x; } for (int i=1; i<MAX; i++) phi[i] += phi[i-1]; int T; cin >> T; while (T--) { int N; cin >> N; ll ans = (ll)(N)*(N-1)-phi[N]+1; printf("%lld\n", ans); } }