結果
問題 | No.1262 グラフを作ろう! |
ユーザー | SHIJOU |
提出日時 | 2020-10-17 15:13:22 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 3,543 bytes |
コンパイル時間 | 3,217 ms |
コンパイル使用メモリ | 211,220 KB |
実行使用メモリ | 13,552 KB |
最終ジャッジ日時 | 2023-09-28 08:09:16 |
合計ジャッジ時間 | 67,942 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 765 ms
13,552 KB |
testcase_01 | AC | 526 ms
8,700 KB |
testcase_02 | AC | 838 ms
9,720 KB |
testcase_03 | AC | 667 ms
8,176 KB |
testcase_04 | AC | 232 ms
10,508 KB |
testcase_05 | AC | 450 ms
9,472 KB |
testcase_06 | AC | 243 ms
10,972 KB |
testcase_07 | AC | 104 ms
7,872 KB |
testcase_08 | AC | 627 ms
10,444 KB |
testcase_09 | AC | 431 ms
7,896 KB |
testcase_10 | AC | 387 ms
8,688 KB |
testcase_11 | AC | 677 ms
10,904 KB |
testcase_12 | AC | 720 ms
8,904 KB |
testcase_13 | AC | 699 ms
8,204 KB |
testcase_14 | AC | 365 ms
8,992 KB |
testcase_15 | AC | 345 ms
7,608 KB |
testcase_16 | AC | 671 ms
10,160 KB |
testcase_17 | AC | 589 ms
9,876 KB |
testcase_18 | AC | 902 ms
10,180 KB |
testcase_19 | AC | 234 ms
10,028 KB |
testcase_20 | AC | 450 ms
11,200 KB |
testcase_21 | AC | 908 ms
9,760 KB |
testcase_22 | AC | 172 ms
9,752 KB |
testcase_23 | AC | 418 ms
8,116 KB |
testcase_24 | AC | 653 ms
10,532 KB |
testcase_25 | AC | 150 ms
8,948 KB |
testcase_26 | AC | 227 ms
9,196 KB |
testcase_27 | AC | 460 ms
10,392 KB |
testcase_28 | AC | 383 ms
10,056 KB |
testcase_29 | AC | 316 ms
10,248 KB |
testcase_30 | AC | 315 ms
9,004 KB |
testcase_31 | AC | 830 ms
9,212 KB |
testcase_32 | AC | 508 ms
9,260 KB |
testcase_33 | AC | 200 ms
7,608 KB |
testcase_34 | AC | 564 ms
11,040 KB |
testcase_35 | AC | 642 ms
8,532 KB |
testcase_36 | AC | 223 ms
7,608 KB |
testcase_37 | AC | 932 ms
10,316 KB |
testcase_38 | AC | 539 ms
10,032 KB |
testcase_39 | AC | 398 ms
7,916 KB |
testcase_40 | AC | 392 ms
8,488 KB |
testcase_41 | AC | 377 ms
10,444 KB |
testcase_42 | AC | 805 ms
10,516 KB |
testcase_43 | AC | 425 ms
9,524 KB |
testcase_44 | AC | 462 ms
9,188 KB |
testcase_45 | AC | 751 ms
10,708 KB |
testcase_46 | AC | 355 ms
11,172 KB |
testcase_47 | AC | 459 ms
7,648 KB |
testcase_48 | AC | 694 ms
10,176 KB |
testcase_49 | AC | 910 ms
10,024 KB |
testcase_50 | AC | 790 ms
10,536 KB |
testcase_51 | AC | 528 ms
8,348 KB |
testcase_52 | AC | 382 ms
10,972 KB |
testcase_53 | AC | 170 ms
8,380 KB |
testcase_54 | AC | 795 ms
8,700 KB |
testcase_55 | AC | 573 ms
8,948 KB |
testcase_56 | AC | 374 ms
10,544 KB |
testcase_57 | AC | 1,019 ms
10,732 KB |
testcase_58 | AC | 680 ms
8,564 KB |
testcase_59 | AC | 370 ms
10,512 KB |
testcase_60 | AC | 961 ms
10,480 KB |
testcase_61 | AC | 350 ms
8,688 KB |
testcase_62 | AC | 844 ms
9,792 KB |
testcase_63 | AC | 801 ms
8,376 KB |
testcase_64 | AC | 547 ms
9,148 KB |
testcase_65 | AC | 208 ms
9,196 KB |
testcase_66 | AC | 377 ms
9,144 KB |
testcase_67 | AC | 316 ms
9,460 KB |
testcase_68 | AC | 332 ms
11,192 KB |
testcase_69 | AC | 331 ms
9,172 KB |
testcase_70 | AC | 769 ms
11,056 KB |
testcase_71 | AC | 638 ms
9,920 KB |
testcase_72 | AC | 879 ms
10,840 KB |
testcase_73 | AC | 961 ms
10,008 KB |
testcase_74 | AC | 715 ms
9,216 KB |
testcase_75 | AC | 732 ms
9,916 KB |
testcase_76 | AC | 580 ms
9,920 KB |
testcase_77 | AC | 732 ms
9,632 KB |
testcase_78 | AC | 737 ms
9,984 KB |
testcase_79 | AC | 1,008 ms
11,224 KB |
testcase_80 | AC | 749 ms
11,016 KB |
testcase_81 | AC | 713 ms
9,580 KB |
testcase_82 | AC | 875 ms
9,784 KB |
testcase_83 | AC | 709 ms
10,708 KB |
testcase_84 | AC | 624 ms
9,980 KB |
testcase_85 | AC | 1,083 ms
10,784 KB |
testcase_86 | AC | 781 ms
10,056 KB |
testcase_87 | AC | 970 ms
10,440 KB |
testcase_88 | AC | 714 ms
9,456 KB |
testcase_89 | AC | 653 ms
10,748 KB |
testcase_90 | AC | 439 ms
11,264 KB |
testcase_91 | AC | 438 ms
11,448 KB |
testcase_92 | AC | 441 ms
11,308 KB |
testcase_93 | AC | 1,286 ms
11,392 KB |
testcase_94 | TLE | - |
testcase_95 | -- | - |
ソースコード
//#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; } } rep(i, m) { int a; cin >> a; ans -= a; auto div = pp.div(a); for(int d:div) { ans += ll(d)*phi[a/d]; } } cout << ans << endl; }