結果
問題 | No.2249 GCDistance |
ユーザー | roaris |
提出日時 | 2023-03-17 22:56:13 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 2,846 ms / 5,000 ms |
コード長 | 1,541 bytes |
コンパイル時間 | 3,347 ms |
コンパイル使用メモリ | 208,800 KB |
実行使用メモリ | 120,392 KB |
最終ジャッジ日時 | 2023-10-18 15:56:40 |
合計ジャッジ時間 | 36,938 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge15 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2,514 ms
120,392 KB |
testcase_01 | AC | 2,821 ms
120,392 KB |
testcase_02 | AC | 2,810 ms
120,392 KB |
testcase_03 | AC | 2,816 ms
120,392 KB |
testcase_04 | AC | 2,636 ms
120,392 KB |
testcase_05 | AC | 2,824 ms
120,392 KB |
testcase_06 | AC | 2,846 ms
120,392 KB |
testcase_07 | AC | 2,827 ms
120,392 KB |
testcase_08 | AC | 2,529 ms
120,392 KB |
testcase_09 | AC | 2,810 ms
120,392 KB |
testcase_10 | AC | 2,781 ms
120,392 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; int64_t euler_phi(int64_t n) { int64_t ret = n; for(int64_t i = 2; i * i <= n; i++) { if(n % i == 0) { ret -= ret / i; while(n % i == 0) n /= i; } } if(n > 1) ret -= ret / n; return ret; } 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); cout << ans-phi[N]+1 << endl; } }