#955398 (C++14) No.5020 Averaging

提出ソース

結果

問題	No.5020 Averaging
ユーザー	monnu
提出日時	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));
      |                ^~~~~

ソースコード

raw source code

#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';
    }
}

yukicoder

結果

テストケース

コンパイルメッセージ

ソースコード