結果
問題 | No.2249 GCDistance |
ユーザー | pointN |
提出日時 | 2023-03-17 22:24:45 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,421 ms / 5,000 ms |
コード長 | 1,698 bytes |
コンパイル時間 | 1,292 ms |
コンパイル使用メモリ | 132,248 KB |
実行使用メモリ | 120,396 KB |
最終ジャッジ日時 | 2023-10-18 15:21:00 |
合計ジャッジ時間 | 19,765 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1,318 ms
120,396 KB |
testcase_01 | AC | 1,392 ms
120,396 KB |
testcase_02 | AC | 1,421 ms
120,396 KB |
testcase_03 | AC | 1,389 ms
120,396 KB |
testcase_04 | AC | 1,367 ms
120,396 KB |
testcase_05 | AC | 1,407 ms
120,396 KB |
testcase_06 | AC | 1,400 ms
120,396 KB |
testcase_07 | AC | 1,392 ms
120,396 KB |
testcase_08 | AC | 1,353 ms
120,396 KB |
testcase_09 | AC | 1,387 ms
120,396 KB |
testcase_10 | AC | 1,381 ms
120,396 KB |
ソースコード
#include <iostream> #include <string> #include <vector> #include <algorithm> #include <functional> #include <cmath> #include <iomanip> #include <stack> #include <queue> #include <numeric> #include <map> #include <unordered_map> #include <set> #include <fstream> #include <chrono> #include <random> #include <bitset> //#include <atcoder/all> #define rep(i,n) for(int i=0;i<(n);i++) #define all(x) x.begin(), x.end() #define rall(x) x.rbegin(), x.rend() #define sz(x) ((int)(x).size()) #define pb push_back using ll = long long; using namespace std; template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; } template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; } ll gcd(ll a, ll b) {return b?gcd(b,a%b):a;} ll lcm(ll a, ll b) {return a/gcd(a,b)*b;} ll euler_phi(ll n) { ll ret = n; for(ll 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; } int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); int e7 = 10000000; vector<int> prime(e7+1, 0); prime[0] = prime[1] = 1; for(int i=2; i <= e7; i++){ if(prime[i]==0){ for(int j=i*2; j <= e7; j+=i){ if(prime[j]==0) prime[j]=i; } prime[i]=i; } } vector<ll> ps(e7+2,0); for(int i=2;i<=e7;i++){ int now = i; ll ret = now; while(now!=1){ int d = prime[now]; ret -= ret / d; while(now%d==0) now/=d; } if(now>1) ret -= ret/now; ps[i+1] = ps[i]+ret; } int T; cin >> T; while(T--){ ll N; cin >> N; ll ans = N*(N-1); ans -= ps[N+1]; cout << ans << '\n'; } return 0; }