結果
問題 | No.1262 グラフを作ろう! |
ユーザー |
![]() |
提出日時 | 2020-10-16 23:27:55 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,037 ms / 3,000 ms |
コード長 | 1,233 bytes |
コンパイル時間 | 1,734 ms |
コンパイル使用メモリ | 174,328 KB |
実行使用メモリ | 88,960 KB |
最終ジャッジ日時 | 2024-07-21 00:06:25 |
合計ジャッジ時間 | 64,349 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 96 |
ソースコード
#include <bits/stdc++.h> //#include <chrono> //#pragma GCC optimize("O3") using namespace std; #define reps(i,s,n) for(int i = s; i < n; i++) #define rep(i,n) reps(i,0,n) #define Rreps(i,n,e) for(int i = n - 1; i >= e; --i) #define Rrep(i,n) Rreps(i,n,0) #define ALL(a) a.begin(), a.end() #define fi first #define se second using ll = long long; using vec = vector<ll>; using mat = vector<vec>; ll N,M,H,W,Q,K,A,B; string S; typedef pair<ll, ll> P; const ll INF = (1LL<<48); int main() { cin>>N>>M; vec euler(N + 1, 1), res(N + 1, 0); mat prime(N + 1, vec(0)); reps(i, 2, N+1){ if(prime[i].empty()){ for(int j = i; j <= N; j += i){ prime[j].push_back(i); } } } reps(i, 1, N + 1){ int temp(i); for(ll v : prime[i]){ int k(0), tt(temp); while(temp%v == 0){ temp /= v; ++k; } tt /= temp; euler[i] *= tt - tt / v; } } reps(i, 1, N + 1){ for(int j = 1; j * i <= N; ++j) res[j * i] += euler[i] * j; res[i] -= i; } ll Res(0); rep(i, M){ cin>>A; Res += res[A]; } cout<<Res<<endl; }