結果

問題 No.5020 Averaging
ユーザー monnumonnu
提出日時 2024-02-25 14:59:02
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 934 ms / 1,000 ms
コード長 5,708 bytes
コンパイル時間 1,897 ms
コンパイル使用メモリ 179,808 KB
実行使用メモリ 6,676 KB
スコア 14,265,049
最終ジャッジ日時 2024-02-25 14:59:53
合計ジャッジ時間 51,081 ms
ジャッジサーバーID
(参考情報)
judge11 / judge15
純コード判定しない問題か言語
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 907 ms
6,676 KB
testcase_01 AC 914 ms
6,676 KB
testcase_02 AC 931 ms
6,676 KB
testcase_03 AC 908 ms
6,676 KB
testcase_04 AC 908 ms
6,676 KB
testcase_05 AC 931 ms
6,676 KB
testcase_06 AC 913 ms
6,676 KB
testcase_07 AC 906 ms
6,676 KB
testcase_08 AC 931 ms
6,676 KB
testcase_09 AC 930 ms
6,676 KB
testcase_10 AC 906 ms
6,676 KB
testcase_11 AC 906 ms
6,676 KB
testcase_12 AC 931 ms
6,676 KB
testcase_13 AC 930 ms
6,676 KB
testcase_14 AC 909 ms
6,676 KB
testcase_15 AC 906 ms
6,676 KB
testcase_16 AC 930 ms
6,676 KB
testcase_17 AC 921 ms
6,676 KB
testcase_18 AC 906 ms
6,676 KB
testcase_19 AC 930 ms
6,676 KB
testcase_20 AC 931 ms
6,676 KB
testcase_21 AC 907 ms
6,676 KB
testcase_22 AC 906 ms
6,676 KB
testcase_23 AC 930 ms
6,676 KB
testcase_24 AC 930 ms
6,676 KB
testcase_25 AC 906 ms
6,676 KB
testcase_26 AC 906 ms
6,676 KB
testcase_27 AC 930 ms
6,676 KB
testcase_28 AC 934 ms
6,676 KB
testcase_29 AC 906 ms
6,676 KB
testcase_30 AC 907 ms
6,676 KB
testcase_31 AC 930 ms
6,676 KB
testcase_32 AC 930 ms
6,676 KB
testcase_33 AC 906 ms
6,676 KB
testcase_34 AC 905 ms
6,676 KB
testcase_35 AC 931 ms
6,676 KB
testcase_36 AC 906 ms
6,676 KB
testcase_37 AC 931 ms
6,676 KB
testcase_38 AC 930 ms
6,676 KB
testcase_39 AC 906 ms
6,676 KB
testcase_40 AC 908 ms
6,676 KB
testcase_41 AC 930 ms
6,676 KB
testcase_42 AC 931 ms
6,676 KB
testcase_43 AC 906 ms
6,676 KB
testcase_44 AC 906 ms
6,676 KB
testcase_45 AC 929 ms
6,676 KB
testcase_46 AC 930 ms
6,676 KB
testcase_47 AC 907 ms
6,676 KB
testcase_48 AC 906 ms
6,676 KB
testcase_49 AC 930 ms
6,676 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'int main()':
main.cpp:191:15: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
  191 |     for (auto [u, v] : ans)
      |               ^
main.cpp:57:37: warning: 'score' may be used uninitialized [-Wmaybe-uninitialized]
   57 |         double score = log10(max<ll>(1, score));
      |                              ~~~~~~~^~~~~~~~~~
main.cpp:57:16: note: 'score' was declared here
   57 |         double score = log10(max<ll>(1, score));
      |                ^~~~~

ソースコード

diff #

#include <bits/stdc++.h>
#include <random>
using namespace std;
using ll = long long;
#define TARGET (500000000000000000)

random_device seed_gen;
mt19937 engine(seed_gen());
uniform_real_distribution<> uniform_distribution(0.0, 1.0);
double getRandomValue()
{
    return uniform_distribution(engine);
}

int main()
{
    int N;
    cin >> N;
    vector<ll> A(N), B(N);
    for (int i = 0; i < N; i++)
    {
        cin >> A[i] >> B[i];
    }

    vector<int> contribution(N, -1);
    contribution[0] = 0;
    for (int i = 1; i < N; i++)
    {
        int j = getRandomValue() * i;
        contribution[j]++;
        contribution[i] = contribution[j];
    }

    ll sumA = 0;
    ll sumB = 0;
    for (int i = 0; i < N; i++)
    {
        sumA += A[i] / (1LL << (contribution[i]));
        sumB += B[i] / (1LL << (contribution[i]));
    }

    clock_t start = clock();
    double T_init = 1.0e2;
    double T_end = 1.0e-5;
    double limit = 0.9;
    while (true)
    {
        clock_t end = clock();
        double elapsed = (double)(end - start) / CLOCKS_PER_SEC;
        if (elapsed > limit)
        {
            break;
        }
        double T = T_init * pow(T_end / T_init, elapsed / limit);

        ll tmp = max(abs(sumA - TARGET), abs(sumB - TARGET));
        double score = log10(max<ll>(1, score));
        if (false && getRandomValue() < 0.3)
        {
            // add
            vector<int> in, out;
            for (int i = 0; i < N; i++)
            {
                if (contribution[i] == -1)
                {
                    out.push_back(i);
                }
                else
                {
                    in.push_back(i);
                }
            }
            int i = getRandomValue() * out.size();
            i = out[i];
            int j = getRandomValue() * in.size();
            j = in[j];
            sumA -= A[j] / (1LL << contribution[j]);
            sumB -= B[j] / (1LL << contribution[j]);
            contribution[j]++;
            contribution[i] = contribution[j];
            sumA += A[j] / (1LL << contribution[j]);
            sumB += B[j] / (1LL << contribution[j]);
            sumA += A[i] / (1LL << contribution[i]);
            sumB += B[i] / (1LL << contribution[i]);

            ll tmp = max(abs(sumA - TARGET), abs(sumB - TARGET));
            double score_tmp = log10(max<ll>(1, score));
            if (score_tmp > score)
            {
                sumA -= A[j] / (1LL << contribution[j]);
                sumB -= B[j] / (1LL << contribution[j]);
                sumA -= A[i] / (1LL << contribution[i]);
                sumB -= B[i] / (1LL << contribution[i]);
                contribution[j]--;
                contribution[i] = 0;
                sumA += A[j] / (1LL << contribution[j]);
                sumB += B[j] / (1LL << contribution[j]);
            }
        }
        else if (false && getRandomValue() < 0.5)
        {
            // delete
            vector<int> in;
            for (int i = 1; i < N; i++)
            {
                if (contribution[i] != -1)
                {
                    in.push_back(i);
                }
            }
            if (in.size() == 0)
            {
                continue;
            }
            int i = getRandomValue() * in.size();
            i = in[i];
            vector<int> same;
            for (int j = 0; j < N; j++)
            {
                if (i == j)
                {
                    continue;
                }
                if (contribution[i] == contribution[j])
                {
                    same.push_back(j);
                }
            }
            if (same.size() == 0)
            {
                continue;
            }
        }
        else
        {
            int i = getRandomValue() * N;
            int j = getRandomValue() * N;
            if (contribution[i] == contribution[j])
            {
                continue;
            }
            sumA -= A[i] / (1LL << contribution[i]);
            sumB -= B[i] / (1LL << contribution[i]);
            sumA -= A[j] / (1LL << contribution[j]);
            sumB -= B[j] / (1LL << contribution[j]);
            swap(contribution[i], contribution[j]);
            sumA += A[i] / (1LL << contribution[i]);
            sumB += B[i] / (1LL << contribution[i]);
            sumA += A[j] / (1LL << contribution[j]);
            sumB += B[j] / (1LL << contribution[j]);

            ll score_tmp = max(abs(sumA - TARGET), abs(sumB - TARGET));
            double diff = score - score_tmp;
            if (exp(diff / T) >= getRandomValue())
            {
                sumA -= A[i] / (1LL << contribution[i]);
                sumB -= B[i] / (1LL << contribution[i]);
                sumA -= A[j] / (1LL << contribution[j]);
                sumB -= B[j] / (1LL << contribution[j]);
                swap(contribution[i], contribution[j]);
                sumA += A[i] / (1LL << contribution[i]);
                sumB += B[i] / (1LL << contribution[i]);
                sumA += A[j] / (1LL << contribution[j]);
                sumB += B[j] / (1LL << contribution[j]);
            }
        }
    }

    vector<pair<int, int>> ans;
    priority_queue<pair<int, int>> pq;
    for (int i = 0; i < N; i++)
    {
        pq.push({contribution[i], i});
    }
    while (pq.size() > 1)
    {
        int u = pq.top().second;
        pq.pop();
        int v = pq.top().second;
        int d = pq.top().first;
        pq.pop();
        if (u > v)
        {
            swap(u, v);
        }
        ans.push_back({u, v});
        pq.push({d, u});
    }

    cout << ans.size() << '\n';
    for (auto [u, v] : ans)
    {
        cout << u + 1 << ' ' << v + 1 << '\n';
    }
}
0