結果
問題 | No.1262 グラフを作ろう! |
ユーザー | SHIJOU |
提出日時 | 2020-10-17 15:17:53 |
言語 | C++17 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 922 ms / 3,000 ms |
コード長 | 3,657 bytes |
コンパイル時間 | 2,299 ms |
コンパイル使用メモリ | 213,472 KB |
実行使用メモリ | 19,256 KB |
最終ジャッジ日時 | 2023-09-28 08:11:29 |
合計ジャッジ時間 | 49,863 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 538 ms
12,700 KB |
testcase_01 | AC | 399 ms
11,696 KB |
testcase_02 | AC | 686 ms
15,100 KB |
testcase_03 | AC | 311 ms
9,852 KB |
testcase_04 | AC | 241 ms
16,756 KB |
testcase_05 | AC | 424 ms
13,940 KB |
testcase_06 | AC | 253 ms
18,560 KB |
testcase_07 | AC | 102 ms
9,296 KB |
testcase_08 | AC | 610 ms
17,192 KB |
testcase_09 | AC | 177 ms
9,080 KB |
testcase_10 | AC | 314 ms
11,044 KB |
testcase_11 | AC | 645 ms
18,240 KB |
testcase_12 | AC | 473 ms
11,780 KB |
testcase_13 | AC | 341 ms
10,320 KB |
testcase_14 | AC | 314 ms
12,060 KB |
testcase_15 | AC | 119 ms
8,216 KB |
testcase_16 | AC | 601 ms
16,352 KB |
testcase_17 | AC | 522 ms
14,408 KB |
testcase_18 | AC | 736 ms
16,184 KB |
testcase_19 | AC | 247 ms
15,392 KB |
testcase_20 | AC | 471 ms
19,028 KB |
testcase_21 | AC | 705 ms
14,712 KB |
testcase_22 | AC | 180 ms
14,884 KB |
testcase_23 | AC | 259 ms
9,736 KB |
testcase_24 | AC | 621 ms
16,848 KB |
testcase_25 | AC | 155 ms
12,328 KB |
testcase_26 | AC | 236 ms
12,888 KB |
testcase_27 | AC | 457 ms
16,504 KB |
testcase_28 | AC | 395 ms
15,064 KB |
testcase_29 | AC | 339 ms
16,288 KB |
testcase_30 | AC | 301 ms
12,548 KB |
testcase_31 | AC | 581 ms
12,820 KB |
testcase_32 | AC | 449 ms
12,992 KB |
testcase_33 | AC | 105 ms
8,400 KB |
testcase_34 | AC | 561 ms
18,420 KB |
testcase_35 | AC | 382 ms
10,780 KB |
testcase_36 | AC | 92 ms
8,432 KB |
testcase_37 | AC | 772 ms
16,316 KB |
testcase_38 | AC | 488 ms
15,208 KB |
testcase_39 | AC | 218 ms
9,200 KB |
testcase_40 | AC | 304 ms
10,920 KB |
testcase_41 | AC | 386 ms
16,996 KB |
testcase_42 | AC | 699 ms
16,776 KB |
testcase_43 | AC | 405 ms
14,136 KB |
testcase_44 | AC | 426 ms
13,484 KB |
testcase_45 | AC | 715 ms
17,800 KB |
testcase_46 | AC | 368 ms
18,920 KB |
testcase_47 | AC | 138 ms
8,524 KB |
testcase_48 | AC | 623 ms
16,024 KB |
testcase_49 | AC | 729 ms
15,592 KB |
testcase_50 | AC | 693 ms
17,040 KB |
testcase_51 | AC | 194 ms
10,012 KB |
testcase_52 | AC | 391 ms
18,196 KB |
testcase_53 | AC | 158 ms
10,576 KB |
testcase_54 | AC | 426 ms
11,584 KB |
testcase_55 | AC | 414 ms
12,440 KB |
testcase_56 | AC | 377 ms
16,508 KB |
testcase_57 | AC | 857 ms
17,776 KB |
testcase_58 | AC | 345 ms
11,208 KB |
testcase_59 | AC | 379 ms
17,268 KB |
testcase_60 | AC | 797 ms
17,032 KB |
testcase_61 | AC | 278 ms
11,456 KB |
testcase_62 | AC | 638 ms
14,620 KB |
testcase_63 | AC | 370 ms
11,056 KB |
testcase_64 | AC | 408 ms
12,220 KB |
testcase_65 | AC | 215 ms
12,976 KB |
testcase_66 | AC | 325 ms
11,896 KB |
testcase_67 | AC | 310 ms
13,704 KB |
testcase_68 | AC | 345 ms
18,848 KB |
testcase_69 | AC | 311 ms
13,252 KB |
testcase_70 | AC | 620 ms
18,172 KB |
testcase_71 | AC | 374 ms
15,128 KB |
testcase_72 | AC | 643 ms
17,832 KB |
testcase_73 | AC | 434 ms
15,476 KB |
testcase_74 | AC | 210 ms
13,356 KB |
testcase_75 | AC | 405 ms
15,200 KB |
testcase_76 | AC | 337 ms
14,880 KB |
testcase_77 | AC | 302 ms
14,144 KB |
testcase_78 | AC | 442 ms
15,724 KB |
testcase_79 | AC | 675 ms
17,972 KB |
testcase_80 | AC | 620 ms
18,592 KB |
testcase_81 | AC | 266 ms
14,024 KB |
testcase_82 | AC | 345 ms
14,536 KB |
testcase_83 | AC | 571 ms
18,068 KB |
testcase_84 | AC | 357 ms
15,096 KB |
testcase_85 | AC | 699 ms
17,852 KB |
testcase_86 | AC | 390 ms
15,256 KB |
testcase_87 | AC | 589 ms
16,708 KB |
testcase_88 | AC | 254 ms
13,960 KB |
testcase_89 | AC | 539 ms
17,520 KB |
testcase_90 | AC | 322 ms
19,152 KB |
testcase_91 | AC | 323 ms
19,044 KB |
testcase_92 | AC | 321 ms
19,256 KB |
testcase_93 | AC | 332 ms
19,028 KB |
testcase_94 | AC | 324 ms
18,948 KB |
testcase_95 | AC | 922 ms
19,044 KB |
ソースコード
//#define _GLIBCXX_DEBUG #include <bits/stdc++.h> #define rep(i, n) for(int i=0; i<n; ++i) #define all(v) v.begin(), v.end() #define rall(v) v.rbegin(), v.rend() using namespace std; using ll = int64_t; using ld = long double; using P = pair<int, int>; using vs = vector<string>; using vi = vector<int>; using vvi = vector<vi>; template<class T> using PQ = priority_queue<T>; template<class T> using PQG = priority_queue<T, vector<T>, greater<T> >; const int INF = 0xccccccc; const ll LINF = 0xcccccccccccccccLL; template<typename T1, typename T2> inline bool chmax(T1 &a, T2 b) {return a < b && (a = b, true);} template<typename T1, typename T2> inline bool chmin(T1 &a, T2 b) {return a > b && (a = b, true);} template<typename T1, typename T2> istream &operator>>(istream &is, pair<T1, T2> &p) { return is >> p.first >> p.second;} template<typename T1, typename T2> ostream &operator<<(ostream &os, const pair<T1, T2> &p) { return os << p.first << ' ' << p.second;} struct Sieve { int n; vector<int> f, primes; Sieve(int n=1):n(n), f(n+1) { f[0] = f[1] = -1; for(int i=2; i<=n; ++i) { if(!f[i]) { primes.emplace_back(i); f[i] = i; } for(int p:primes) { int u = p*i; if(u > n or p > f[i]) break; f[u] = p; } } } //素数判定 bool isPrime(int x) {return f[x] == x;} //素数判定(1つ用)(宣言はsqrt(x)のサイズで) bool isPrime1(int64_t x) { bool m = true; for(int z:primes) { if(!(x%z)) { m = false; break; } } return m; } //素数列挙 vector<int > factorList(int x) { assert(x <= n); vector<int> res; while(x != 1) { res.emplace_back(f[x]); x /= f[x]; } return res; } //素因数分解 pair(素数, 素因数の数) vector<pair<int, int> > factor(int x) { vector<int> fl = factorList(x); if(fl.size() == 0) return {}; vector<pair<int, int> > res(1, pair<int, int>(fl[0], 0)); for(int p : fl) { if(res.back().first == p) { ++res.back().second; } else { res.emplace_back(p, 1); } } return res; } //素数列挙(1つ用)(宣言はsqrt(x)のサイズで) vector<int64_t> factorList1(int64_t x) { vector<int64_t> vec; for(int z:primes) { while(!(x%z)) { vec.emplace_back(z); x /= z; } } if(x != 1) vec.emplace_back(x); return vec; } //素数列挙(1つ用)(宣言はsqrt(x)のサイズで) pair(素数, 素因数の数) vector<pair<int64_t, int> > factor1(int64_t x) { vector<int64_t> vec = factorList1(x); vector<pair<int64_t, int> > vecc; for(vector<int64_t>::iterator itr = vec.begin(); itr != vec.end(); ++itr) { if(itr != vec.begin() && *itr == *(itr-1)) { (vecc.end()-1)->second++; } else { vecc.emplace_back(pair<int64_t, int>(*itr, 1)); } } return vecc; } //約数列挙 vector<int> div(int x) { vector<int> res(1, 1); vector<pair<int, int> > fac = factor(x); for(int i = 0; i < int(fac.size()); ++i) { int si = res.size(); for(int j = 0; j < fac[i].second; ++j) { for(int k = 0; k < si; ++k) { res.emplace_back(res[k+j*si]*fac[i].first); } } } sort(res.begin(), res.end()); return res; } } pp(1000000); //head int main() { ios::sync_with_stdio(false); cin.tie(0); int n, m; cin >> n >> m; ll ans = 0; vi phi(n+1); for(int i = 1; i <= n; i++) { int &u = phi[i] = i; auto f = pp.factor(i); for(P x:f) { u /= x.first; u *= x.first-1; } } vector<ll> sup(n+1, -1); rep(i, m) { int a; cin >> a; ans -= a; if(sup[a] != -1) { ans += sup[a]; continue; } auto div = pp.div(a); ll &z = sup[a] = 0; for(int d:div) { z += ll(d)*phi[a/d]; } ans += z; } cout << ans << endl; }