結果

問題 No.1262 グラフを作ろう！
ユーザー kyoprounokyoprouno
提出日時 2021-01-04 01:56:47
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 552 ms / 3,000 ms
コード長 7,137 bytes
コンパイル時間 2,225 ms
コンパイル使用メモリ 182,416 KB
実行使用メモリ 15,620 KB
最終ジャッジ日時 2024-04-21 21:29:14
合計ジャッジ時間 28,678 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 277 ms
14,708 KB
testcase_01 AC 263 ms
12,504 KB
testcase_02 AC 334 ms
14,464 KB
testcase_03 AC 405 ms
15,104 KB
testcase_04 AC 28 ms
8,064 KB
testcase_05 AC 141 ms
11,008 KB
testcase_06 AC 20 ms
7,856 KB
testcase_07 AC 37 ms
8,448 KB
testcase_08 AC 190 ms
11,888 KB
testcase_09 AC 217 ms
12,800 KB
testcase_10 AC 143 ms
11,068 KB
testcase_11 AC 197 ms
12,032 KB
testcase_12 AC 343 ms
14,336 KB
testcase_13 AC 308 ms
15,036 KB
testcase_14 AC 143 ms
10,496 KB
testcase_15 AC 165 ms
12,160 KB
testcase_16 AC 198 ms
12,252 KB
testcase_17 AC 226 ms
12,120 KB
testcase_18 AC 317 ms
14,464 KB
testcase_19 AC 41 ms
8,320 KB
testcase_20 AC 97 ms
9,776 KB
testcase_21 AC 358 ms
15,232 KB
testcase_22 AC 23 ms
7,936 KB
testcase_23 AC 197 ms
12,288 KB
testcase_24 AC 200 ms
12,392 KB
testcase_25 AC 34 ms
8,192 KB
testcase_26 AC 59 ms
8,960 KB
testcase_27 AC 119 ms
10,412 KB
testcase_28 AC 103 ms
9,984 KB
testcase_29 AC 71 ms
9,216 KB
testcase_30 AC 98 ms
9,932 KB
testcase_31 AC 306 ms
15,224 KB
testcase_32 AC 168 ms
11,904 KB
testcase_33 AC 91 ms
9,984 KB
testcase_34 AC 140 ms
10,836 KB
testcase_35 AC 260 ms
14,336 KB
testcase_36 AC 107 ms
10,448 KB
testcase_37 AC 437 ms
15,104 KB
testcase_38 AC 167 ms
11,464 KB
testcase_39 AC 179 ms
12,200 KB
testcase_40 AC 154 ms
11,392 KB
testcase_41 AC 81 ms
9,528 KB
testcase_42 AC 255 ms
13,648 KB
testcase_43 AC 129 ms
10,684 KB
testcase_44 AC 151 ms
11,264 KB
testcase_45 AC 221 ms
12,928 KB
testcase_46 AC 58 ms
8,960 KB
testcase_47 AC 223 ms
13,568 KB
testcase_48 AC 210 ms
12,720 KB
testcase_49 AC 302 ms
14,868 KB
testcase_50 AC 345 ms
13,312 KB
testcase_51 AC 230 ms
13,568 KB
testcase_52 AC 74 ms
9,196 KB
testcase_53 AC 54 ms
8,812 KB
testcase_54 AC 367 ms
15,232 KB
testcase_55 AC 274 ms
12,672 KB
testcase_56 AC 95 ms
9,472 KB
testcase_57 AC 324 ms
15,232 KB
testcase_58 AC 271 ms
14,444 KB
testcase_59 AC 77 ms
9,472 KB
testcase_60 AC 304 ms
14,848 KB
testcase_61 AC 129 ms
10,624 KB
testcase_62 AC 280 ms
14,588 KB
testcase_63 AC 303 ms
15,488 KB
testcase_64 AC 186 ms
12,248 KB
testcase_65 AC 47 ms
8,704 KB
testcase_66 AC 130 ms
10,748 KB
testcase_67 AC 85 ms
9,472 KB
testcase_68 AC 49 ms
8,632 KB
testcase_69 AC 94 ms
9,856 KB
testcase_70 AC 210 ms
12,416 KB
testcase_71 AC 199 ms
12,472 KB
testcase_72 AC 258 ms
13,560 KB
testcase_73 AC 331 ms
15,584 KB
testcase_74 AC 435 ms
14,792 KB
testcase_75 AC 283 ms
13,248 KB
testcase_76 AC 189 ms
12,160 KB
testcase_77 AC 259 ms
14,068 KB
testcase_78 AC 227 ms
12,984 KB
testcase_79 AC 297 ms
14,468 KB
testcase_80 AC 197 ms
12,032 KB
testcase_81 AC 252 ms
14,080 KB
testcase_82 AC 318 ms
15,488 KB
testcase_83 AC 275 ms
12,032 KB
testcase_84 AC 225 ms
12,416 KB
testcase_85 AC 394 ms
15,104 KB
testcase_86 AC 258 ms
13,824 KB
testcase_87 AC 308 ms
14,716 KB
testcase_88 AC 260 ms
14,192 KB
testcase_89 AC 179 ms
11,796 KB
testcase_90 AC 155 ms
15,556 KB
testcase_91 AC 154 ms
15,488 KB
testcase_92 AC 154 ms
15,528 KB
testcase_93 AC 419 ms
15,488 KB
testcase_94 AC 552 ms
15,444 KB
testcase_95 AC 306 ms
15,620 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC optimize("O3")

#include <bits/stdc++.h>
using namespace std;
#ifdef LOCAL
#define _GLIBCXX_DEBUG
#endif

#pragma region Macros

#define FOR(i, l, r) for(int i = (l) ;i < (r) ;i++)
#define REP(i, n) FOR(i, 0, n)
#define REPS(i, n) FOR(i, 1, n+1)
#define RFOR(i, l, r) for(int i = (l) ; i >= (r) ; i--)
#define RREP(i, n) RFOR(i, n - 1, 0)
#define RREPS(i, n) RFOR(i, n , 1)
#define SZ(x) ((int)(x).size())
#define all(x) (x).begin(), (x).end()
#define rall(x) (x).rbegin(), (x).rend()
template <class T = int>
using V = vector<T>;
template <class T = int>
using VV = V<V<T>>;
#define ll long long
using ld = long double;
#define pll pair<ll,ll>
#define pii pair<int,int>
#define pq priority_queue
#define pb push_back
#define eb emplace_back
#define fi first
#define se second
#define endl '\n'
#define ios ios_base::sync_with_stdio(0),cin.tie(0),cout.tie(0);
#define lb(c,x) distance(c.begin(),lower_bound(all(c),x))
#define ub(c,x) distance(c.begin(),upper_bound(all(c),x))

#define VEC(type, name, size) \
	V<type> name(size);       \
	INPUT(name)
#define VVEC(type, name, h, w)    \
	VV<type> name(h, V<type>(w)); \
	INPUT(name)
#define INT(...)     \
	int __VA_ARGS__; \
	INPUT(__VA_ARGS__)
#define LL(...)     \
	ll __VA_ARGS__; \
	INPUT(__VA_ARGS__)
#define STR(...)        \
	string __VA_ARGS__; \
	INPUT(__VA_ARGS__)
#define CHAR(...)     \
	char __VA_ARGS__; \
	INPUT(__VA_ARGS__)
#define DOUBLE(...)     \
	DOUBLE __VA_ARGS__; \
	INPUT(__VA_ARGS__)
#define LD(...)     \
	LD __VA_ARGS__; \
	INPUT(__VA_ARGS__)

void scan(int &a) { cin >> a;}
void scan(long long &a) { cin >> a;}
void scan(char &a) { cin >> a;}
void scan(double &a) { cin >> a;}
void scan(long double &a) { cin >> a;}
void scan(char a[]) { scanf("%s", a);}
void scan(string &a) { cin >> a;}
template <class T>
void scan(V<T> &);
template <class T, class L>
void scan(pair<T, L> &);
template <class T>
void scan(V<T> &a){
	for(auto &i : a) scan(i);
}
template <class T, class L>
void scan(pair<T, L> &p){
	scan(p.first);
	scan(p.second);
}
template <class T>
void scan(T &a) { cin >> a;}
void INPUT() {}
template <class Head, class... Tail>
void INPUT(Head &head, Tail &... tail){
	scan(head);
	INPUT(tail...);
}
template <class T>
inline void print(T x) { cout << x << '\n'; }
template <typename T1, typename T2>
istream &operator>>(istream &is, pair<T1, T2> &p) {
	is >> p.first >> p.second;
	return is;
}
template <typename T1, typename T2>
ostream &operator<<(ostream &os, const pair<T1, T2> &p) {
	os << p.first << " " << p.second;
	return os;
}
template <typename T>
istream &operator>>(istream &is, vector<T> &v) {
	for(T &in : v) is >> in;
	return is;
}
template <class T>
ostream &operator<<(ostream &os, const V<T> &v) {
	REP(i, SZ(v)) {
		if(i) os << " ";
		os << v[i];
	}
	return os;
}
//debug
template <typename T>
void view(const V<T> &v) {
	cerr << "{ ";
	for(const auto &e : v) {
		cerr << e << ", ";
	}
	cerr << "\b\b }";
}
template <typename T>
void view(const VV<T> &vv) {
	cerr << "{\n";
	for(const auto &v : vv) {
		cerr << "\t";
		view(v);
		cerr << "\n";
	}
	cerr << "}";
}
template <typename T, typename U>
void view(const V<pair<T, U>> &v) {
	cerr << "{\n";
	for(const auto &c : v) cerr << "\t(" << c.first << ", " << c.second << ")\n";
	cerr << "}";
}
template <typename T, typename U>
void view(const map<T, U> &m) {
	cerr << "{\n";
	for(auto &t : m) cerr << "\t[" << t.first << "] : " << t.second << "\n";
	cerr << "}";
}
template <typename T, typename U>
void view(const pair<T, U> &p) { cerr << "(" << p.first << ", " << p.second << ")"; }
template <typename T>
void view(const set<T> &s) {
	cerr << "{ ";
	for(auto &t : s) {
		view(t);
		cerr << ", ";
	}
	cerr << "\b\b }";
}
template <typename T>
void view(T e) { cerr << e; }
#ifdef LOCAL
void debug_out() {}
template <typename Head, typename... Tail>
void debug_out(Head H, Tail... T) {
	view(H);
	cerr << ", ";
	debug_out(T...);
}
#define debug(...)                                           \
	do {                                                     \
		cerr << __LINE__ << " [" << #__VA_ARGS__ << "] : ["; \
		debug_out(__VA_ARGS__);                              \
		cerr << "\b\b]\n";                                   \
	} while(0)
#else
#define debug(...) (void(0))
#endif
template <class T>
V<T> press(V<T> &x) {
	V<T> res = x;
	sort(all(res));
	res.erase(unique(all(res)), res.end());
	REP(i, SZ(x)) {
		x[i] = lower_bound(all(res), x[i]) - res.begin();
	}
	return res;
}
template <class T>
inline bool chmin(T &a, T b) {
	if(a > b) {
		a = b;
		return true;
	}
	return false;
}
template <class T>
inline bool chmax(T &a, T b) {
	if(a < b) {
		a = b;
		return true;
	}
	return false;
}
inline void Yes(bool b = true) { cout << (b ? "Yes" : "No") << '\n'; }
inline void YES(bool b = true) { cout << (b ? "YES" : "NO") << '\n'; }
inline void err(bool b = true) {
	if(b) {
		cout << -1 << '\n';
		exit(0);
	}
}
template <class T>
inline void fin(bool b = true, T e = 0) {
	if(b) {
		cout << e << '\n';
		exit(0);
	}
}
template <class T>
T divup(T x, T y) { return (x + (y - 1)) / y; }
template <typename T>
T pow(T a, long long n, T e = 1) {
	T ret = e;
	while(n) {
		if(n & 1) ret *= a;
		a *= a;
		n >>= 1;
	}
	return ret;
}
struct iofast {
	iofast() {
		ios_base::sync_with_stdio(false);
		cin.tie(nullptr);
		cout << fixed << setprecision(15);
	}
} iofast_;
const int inf = 1e9;
const ll INF = 1e18;
#pragma endregion


inline int topbit(unsigned long long x){
	return x?63-__builtin_clzll(x):-1;
}

inline int popcount(unsigned long long x){
	return __builtin_popcountll(x);
}

inline int parity(unsigned long long x){//popcount%2
	return __builtin_parity(x);
}

struct Sieve{
    int n;
    vector<int> f, primes;
    //コンストラクタを定義
    Sieve(int n=1):n(n), f(n+1){//メンバ変数の初期化
        f[0]=f[1]=-1;
        for(ll i=2;i<=n;i++){
            if(f[i]) continue;
            primes.push_back(i);
            f[i]=i;
            for(ll j=i*i;j<=n;j+=i){
                if(!f[j]) f[j]=i;
            }
        }
    }
    bool isPrime(int x) { return f[x]==x;} //xが素数かどうかを判定
    vector<int> factorList(int x){//xの素因数を重複ありで小さい順に列挙
        vector<int> res;
        while(x!=1){
            res.push_back(f[x]);
            x/=f[x];
        }
        return res;
    }
    vector<pair<int,int>> factor(int x){//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;
    }
};


int main(){
    LL(N,M);
    VEC(ll, a, M);
    ll ans=0;
    
    Sieve sieve=Sieve(1000005);
    
    
    REP(i,M){
        V<pair<int,int>> p=sieve.factor(a[i]);
        ll sum=a[i];
        ll s=SZ(p);
        REP(i,s){
            sum/=p[i].fi;
            sum*=(p[i].fi-1)*(p[i].se+1)+1;
        }
        ans+=sum-a[i];
    }
    cout << ans << endl;
    return 0;
}
0