結果

問題 No.1917 LCMST
ユーザー Georgi PetkovGeorgi Petkov
提出日時 2023-02-11 00:00:40
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
MLE  
実行時間 -
コード長 2,382 bytes
コンパイル時間 2,170 ms
コンパイル使用メモリ 172,252 KB
実行使用メモリ 814,236 KB
最終ジャッジ日時 2023-09-22 00:50:35
合計ジャッジ時間 7,742 ms
ジャッジサーバーID
(参考情報)
judge12 / judge11
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
5,700 KB
testcase_01 AC 3 ms
5,868 KB
testcase_02 AC 3 ms
5,692 KB
testcase_03 MLE -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
testcase_41 -- -
testcase_42 -- -
testcase_43 -- -
testcase_44 -- -
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
#define endl '\n'

using namespace std;

const int maxn = 1e5 + 3;
const long long inf = 1e18;

int n;
bool is[maxn];
long long ans;
vector <int> mins[maxn];

vector <int> get_divisors(int x) {
    vector <int> ans;
    for (int i = 1; i * i <= x; i++)
        if (x % i == 0) {
            ans.push_back(i);
            if (i * i != x)
                ans.push_back(x / i);
        }

    return ans;
}

int get_minimum(int idx) { 
    while (!mins[idx].empty() && !is[mins[idx].back()])
        mins[idx].pop_back();

    return mins[idx].empty() ? -1 : mins[idx].back();
}

void add_num(int x) {
    vector <int> d = get_divisors(x); 
    for (auto i: d)
        mins[i].push_back(x);
}

void read() {
    cin >> n;
    for (int i = 1; i <= n; i++) {
        int x;
        cin >> x;

        if (is[x]) {
            ans += x;
            continue;
        }  

        is[x] = true;
    } 
}

struct edge {
    int from, to;
    long long w;

    edge(int _from, int _to, long long _w) {from = _from, to = _to, w = _w;}
    bool operator < (const edge &other) const {
        return w > other.w;
    }
};

void add_next(int x, int d, priority_queue <edge> &pq) {
    auto temp = get_minimum(d);
    if (temp == -1)
        return;

    pq.push(edge(x, temp, (long long)x * temp / d));
}

void add_edges(int x, priority_queue <edge> &pq) {
    auto d = get_divisors(x);
    for (auto i: d)
        add_next(x, i, pq);
}

edge get_best_edge(priority_queue <edge> &pq) {
    for (;;) {
        auto e = pq.top();
        pq.pop();

        if (!is[e.to]) {
            add_edges(e.from, pq);
            continue;
        }

        return e;
    } 
}

void solve() {
    int left = 0;
    for (int i = maxn-1; i >= 1; i--)
        if (is[i])
            add_num(i), left++;

    for (int i = 1; i < maxn; i++)
        if (is[i]) {
            is[i] = false;
            priority_queue <edge> pq;
            add_edges(i, pq);
            left--;
            while (left) {
                auto e = get_best_edge(pq);
                ans += e.w;
                is[e.to] = false;
                left--;

                add_edges(e.from, pq);
                add_edges(e.to, pq);
            }

            break;
        }

    cout << ans << endl;
}

int main() {
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);

    read();
    solve();
}
0