結果
問題 | No.2249 GCDistance |
ユーザー |
![]() |
提出日時 | 2023-03-06 17:30:37 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 482 ms / 5,000 ms |
コード長 | 1,925 bytes |
コンパイル時間 | 2,023 ms |
コンパイル使用メモリ | 198,844 KB |
最終ジャッジ日時 | 2025-02-11 05:51:50 |
ジャッジサーバーID (参考情報) |
judge1 / judge6 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 10 |
ソースコード
#include<bits/stdc++.h>#pragma GCC optimize("Ofast")//#pragma GCC target("avx2")#pragma GCC optimize("unroll-loops")using namespace std;//#include<boost/multiprecision/cpp_int.hpp>//#include<boost/multiprecision/cpp_dec_float.hpp>//namespace mp=boost::multiprecision;//#define mulint mp::cpp_int//#define mulfloat mp::cpp_dec_float_100struct __INIT{__INIT(){cin.tie(0);ios::sync_with_stdio(false);cout<<fixed<<setprecision(15);}} __init;//#define INF (1<<30)#define LINF (lint)(1LL<<56)#define MINF (lint)(2e18)#define endl "\n"#define rep(i,n) for(lint (i)=0;(i)<(n);(i)++)#define reprev(i,n) for(lint (i)=(n-1);(i)>=0;(i)--)#define flc(x) __builtin_popcountll(x)#define pint pair<int,int>#define pdouble pair<double,double>#define plint pair<lint,lint>#define fi first#define se second#define all(x) x.begin(),x.end()//#define vec vector<lint>#define nep(x) next_permutation(all(x))typedef long long lint;int dx[8]={1,1,0,-1,-1,-1,0,1};int dy[8]={0,1,1,1,0,-1,-1,-1};const int MAX_N=3e5+5;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;}//vector<int> bucket[MAX_N/1000];//constexpr int MOD=1000000007;constexpr int MOD=998244353;//typedef __int128_t llint;//#include<atcoder/all>//using namespace atcoder;vector< int > euler_phi_table(int n) {vector< int > euler(n + 1);for(int i = 0; i <= n; i++) {euler[i] = i;}for(int i = 2; i <= n; i++) {if(euler[i] == i) {for(int j = i; j <= n; j += i) {euler[j] = euler[j] / i * (i - 1);}}}return euler;}int main(void){int Q;cin >> Q;lint N=10000010;vector<int> res=euler_phi_table(N+1);res[1]--;lint ans[N+2]={};rep(i,N+1) ans[i+1]=ans[i]+2*(i)-res[i+1];while(Q--){int n;cin >> n;cout << ans[n] << endl;}}