結果

問題 No.1917 LCMST
ユーザー Georgi PetkovGeorgi Petkov
提出日時 2023-02-11 04:55:10
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 3,850 ms / 4,000 ms
コード長 3,005 bytes
コンパイル時間 1,788 ms
コンパイル使用メモリ 178,168 KB
実行使用メモリ 147,324 KB
最終ジャッジ日時 2024-07-07 21:07:39
合計ジャッジ時間 50,330 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
5,760 KB
testcase_01 AC 4 ms
5,760 KB
testcase_02 AC 4 ms
5,632 KB
testcase_03 AC 10 ms
5,888 KB
testcase_04 AC 11 ms
6,144 KB
testcase_05 AC 10 ms
5,888 KB
testcase_06 AC 10 ms
5,888 KB
testcase_07 AC 9 ms
6,016 KB
testcase_08 AC 4 ms
5,888 KB
testcase_09 AC 2,060 ms
79,456 KB
testcase_10 AC 2,114 ms
79,576 KB
testcase_11 AC 2,091 ms
79,476 KB
testcase_12 AC 1,536 ms
46,004 KB
testcase_13 AC 433 ms
13,816 KB
testcase_14 AC 1,763 ms
46,500 KB
testcase_15 AC 245 ms
9,444 KB
testcase_16 AC 457 ms
14,008 KB
testcase_17 AC 725 ms
19,336 KB
testcase_18 AC 938 ms
28,368 KB
testcase_19 AC 245 ms
8,840 KB
testcase_20 AC 165 ms
7,548 KB
testcase_21 AC 255 ms
9,888 KB
testcase_22 AC 169 ms
7,552 KB
testcase_23 AC 171 ms
7,424 KB
testcase_24 AC 242 ms
9,820 KB
testcase_25 AC 175 ms
7,552 KB
testcase_26 AC 191 ms
8,472 KB
testcase_27 AC 203 ms
8,476 KB
testcase_28 AC 218 ms
8,580 KB
testcase_29 AC 1,902 ms
145,320 KB
testcase_30 AC 1,883 ms
145,280 KB
testcase_31 AC 1,902 ms
145,288 KB
testcase_32 AC 1,880 ms
145,340 KB
testcase_33 AC 1,892 ms
145,448 KB
testcase_34 AC 90 ms
5,632 KB
testcase_35 AC 86 ms
5,760 KB
testcase_36 AC 87 ms
5,632 KB
testcase_37 AC 3,530 ms
147,324 KB
testcase_38 AC 3,850 ms
146,216 KB
testcase_39 AC 3,460 ms
147,268 KB
testcase_40 AC 3,453 ms
145,960 KB
testcase_41 AC 3,426 ms
147,284 KB
testcase_42 AC 3,161 ms
145,220 KB
testcase_43 AC 94 ms
6,016 KB
testcase_44 AC 93 ms
6,016 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#define endl '\n'

using namespace std;

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

int n;
bool is[maxn];
long long ans;
vector <int> not_taken[maxn];
int taken[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;
}

void add_num(int x) {
    for (auto i: get_divisors(x))
        not_taken[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;
    }
};

ostream &operator << (ostream &out, const edge &e) {
    out << "edge from " << e.from << " to " << e.to << " with weight " << e.w << endl;
    return out;
}

void add_edges(int x, priority_queue <edge> &pq) {
    for (auto i: get_divisors(x)) {
        if (i > bucket) {
            for (auto j: not_taken[i])
                if (is[j])
                    pq.push(edge(x, j, (long long)x * j / i));         
        }

        else {
            if (!taken[i] || x < taken[i])
                taken[i] = x;
        }
    }
}

edge get_best_edge(priority_queue <edge> &pq) { 
    while (!pq.empty()) {
        auto e = pq.top();
        if (!is[e.to]) {
            pq.pop();
            continue;
        }

        return e;
    } 

    return edge(-1, -1, inf);
}

edge get_single(int idx) {
    while (!not_taken[idx].empty() && !is[not_taken[idx].back()])
        not_taken[idx].pop_back();

    if (not_taken[idx].empty() || !taken[idx])
        return edge(-1, -1, inf);

    int a = taken[idx], b = not_taken[idx].back();
    return edge(a, b, (long long)a * b / idx);
}

edge get_multiple() {
    edge ans = {-1, -1, inf};
    for (int i = 1; i <= bucket; i++)
        ans = max(ans, get_single(i));

    return ans;
}

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);
                e = max(e, get_multiple());
                ans += e.w;
                is[e.to] = is[e.from] = false;
                add_edges(e.to, pq);
                left--;
            }

            break;
        }

    cout << ans << endl;
}

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

    read();
    solve();
}
0