結果
問題 | No.2249 GCDistance |
ユーザー |
|
提出日時 | 2023-03-17 22:08:00 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,118 ms / 5,000 ms |
コード長 | 994 bytes |
コンパイル時間 | 773 ms |
コンパイル使用メモリ | 77,596 KB |
最終ジャッジ日時 | 2025-02-11 13:15:49 |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 10 |
ソースコード
#include<iostream>#include<vector>#include<tuple>#include<algorithm>using namespace std;using ll=long long;#define rep(i,n) for(int i=0;i<n;i++)#define rrep(i,n) for(int i=(n)-1;i>=0;i--)#define all(v) v.begin(),v.end()#define rall(v) v.rbegin(),v.rend()template<class T> bool chmax(T &a, T b){if (a < b){a = b;return true;} else return false;}template<class T> bool chmin(T &a, T b){if (a > b){a = b;return true;} else return false;}int main(){ios::sync_with_stdio(false);cin.tie(nullptr);int t;cin>>t;vector<ll>V(1e7+1),S(1e7+1);vector<bool>flag(1e7+1);for(int x=2;x<=1e7;x++){V[x]=x;}for(int x=2;x<=1e7;x++){if(flag[x])continue;flag[x]=true;V[x]=x-1;for(int y=2*x;y<=1e7;y+=x){flag[y]=true;V[y]/=x;V[y]*=x-1;}}rep(i,1e7)S[i+1]=V[i+1]+S[i];while(t--){ll n;cin>>n;cout<<n*(n-1)-S[n]<<"\n";}}