#829261 (C++17) No.5013 セクスタプル (open)

提出ソース

結果

問題	No.5013 セクスタプル (open)
ユーザー	NAMIDAIRO
提出日時	2022-12-29 22:16:29
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	4 ms / 2,000 ms
コード長	3,832 bytes
コンパイル時間	1,454 ms
実行使用メモリ	4,084 KB
スコア	13,479
最終ジャッジ日時	2022-12-29 22:16:36
合計ジャッジ時間	6,796 ms
ジャッジサーバーID （参考情報）	judge13 / judge12
	純コード判定しない問題か言語

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	4 ms 3,936 KB
testcase_01	AC	4 ms 3,844 KB
testcase_02	AC	4 ms 3,828 KB
testcase_03	AC	3 ms 3,776 KB
testcase_04	AC	3 ms 3,824 KB
testcase_05	AC	3 ms 3,880 KB
testcase_06	AC	3 ms 3,776 KB
testcase_07	AC	3 ms 3,976 KB
testcase_08	AC	3 ms 3,984 KB
testcase_09	AC	4 ms 3,776 KB
testcase_10	AC	4 ms 4,080 KB
testcase_11	AC	3 ms 3,880 KB
testcase_12	AC	3 ms 4,084 KB
testcase_13	AC	4 ms 3,844 KB
testcase_14	AC	4 ms 3,976 KB
testcase_15	AC	4 ms 4,084 KB
testcase_16	AC	4 ms 3,932 KB
testcase_17	AC	3 ms 3,924 KB
testcase_18	AC	3 ms 4,032 KB
testcase_19	AC	4 ms 3,924 KB
testcase_20	AC	4 ms 3,884 KB
testcase_21	AC	4 ms 3,880 KB
testcase_22	AC	4 ms 3,880 KB
testcase_23	AC	4 ms 3,920 KB
testcase_24	AC	3 ms 4,080 KB
testcase_25	AC	3 ms 3,776 KB
testcase_26	AC	3 ms 3,936 KB
testcase_27	AC	4 ms 3,896 KB
testcase_28	AC	3 ms 3,900 KB
testcase_29	AC	4 ms 3,780 KB
testcase_30	AC	4 ms 3,896 KB
testcase_31	AC	4 ms 3,884 KB
testcase_32	AC	4 ms 3,808 KB
testcase_33	AC	3 ms 4,080 KB
testcase_34	AC	3 ms 3,840 KB
testcase_35	AC	4 ms 3,896 KB
testcase_36	AC	4 ms 3,876 KB
testcase_37	AC	3 ms 4,028 KB
testcase_38	AC	4 ms 3,776 KB
testcase_39	AC	4 ms 3,828 KB
testcase_40	AC	3 ms 3,900 KB
testcase_41	AC	4 ms 3,980 KB
testcase_42	AC	4 ms 4,080 KB
testcase_43	AC	4 ms 3,828 KB
testcase_44	AC	4 ms 3,936 KB
testcase_45	AC	3 ms 3,844 KB
testcase_46	AC	4 ms 3,776 KB
testcase_47	AC	3 ms 3,840 KB
testcase_48	AC	4 ms 3,928 KB
testcase_49	AC	4 ms 3,904 KB
testcase_50	AC	4 ms 3,880 KB
testcase_51	AC	3 ms 3,772 KB
testcase_52	AC	3 ms 3,924 KB
testcase_53	AC	4 ms 3,924 KB
testcase_54	AC	3 ms 3,924 KB
testcase_55	AC	4 ms 3,908 KB
testcase_56	AC	4 ms 4,028 KB
testcase_57	AC	3 ms 3,896 KB
testcase_58	AC	3 ms 3,928 KB
testcase_59	AC	4 ms 3,880 KB
testcase_60	AC	3 ms 3,936 KB
testcase_61	AC	4 ms 3,844 KB
testcase_62	AC	3 ms 3,908 KB
testcase_63	AC	4 ms 3,932 KB
testcase_64	AC	4 ms 4,032 KB
testcase_65	AC	4 ms 3,844 KB
testcase_66	AC	4 ms 3,936 KB
testcase_67	AC	3 ms 3,884 KB
testcase_68	AC	4 ms 3,932 KB
testcase_69	AC	4 ms 3,824 KB
testcase_70	AC	3 ms 3,824 KB
testcase_71	AC	3 ms 3,844 KB
testcase_72	AC	3 ms 3,840 KB
testcase_73	AC	3 ms 3,844 KB
testcase_74	AC	3 ms 3,772 KB
testcase_75	AC	3 ms 3,880 KB
testcase_76	AC	4 ms 3,840 KB
testcase_77	AC	4 ms 3,824 KB
testcase_78	AC	4 ms 3,876 KB
testcase_79	AC	3 ms 3,976 KB
testcase_80	AC	3 ms 3,928 KB
testcase_81	AC	4 ms 3,884 KB
testcase_82	AC	4 ms 4,028 KB
testcase_83	AC	4 ms 3,936 KB
testcase_84	AC	4 ms 3,932 KB
testcase_85	AC	4 ms 3,884 KB
testcase_86	AC	4 ms 4,080 KB
testcase_87	AC	4 ms 3,936 KB
testcase_88	AC	3 ms 3,908 KB
testcase_89	AC	4 ms 3,880 KB
testcase_90	AC	4 ms 3,904 KB
testcase_91	AC	4 ms 3,904 KB
testcase_92	AC	3 ms 3,884 KB
testcase_93	AC	3 ms 3,844 KB
testcase_94	AC	3 ms 3,904 KB
testcase_95	AC	3 ms 3,776 KB
testcase_96	AC	3 ms 3,924 KB
testcase_97	AC	4 ms 3,824 KB
testcase_98	AC	3 ms 4,084 KB
testcase_99	AC	4 ms 4,032 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <iostream>
#include <vector>
#include <algorithm>
#include <map>
#include <queue>
#include <set>
#include <ctime>
#include <random>
#include <cmath>
using namespace std;
typedef long long ll;
#define rep(i, n) for(int i = 0; i < (n); i++)

template<class T>
using vi = vector<T>;

template<class T>
using vii = vector<vi<T>>;

template<class T>
using viii = vector<vii<T>>;

using pii = pair<int, int>;

int n = 36;
const int len = 6;
const int TL = 1 * 1000;
const int TL2 = 2 * 1000;

struct dice {
    vi<int> vec;
    int idx = 0;
    dice() : vec(len), idx(0){}
};

vi<dice> d(n);
vii<int> board(len, vi<int>(len));

int score() {
    const int inf = 1e9;
    int res = 0;
    rep(i, len){
        vi<int> cnt(len, 0);
        rep(j, len) {
            rep(k, len) { //どの値
                int pos = board[i][j];
                int num = d[pos].vec[k];
                if (num) cnt[k] += num;
                else cnt[k] = -inf;
            }
        }

        rep(k, len) {
            if (cnt[k] > 0) res += 3 + cnt[k] - 6;
        }
    }

    rep(j, len){
        vi<int> cnt(len, 0);
        rep(i, len) {
            rep(k, len) { //どの値
                int pos = board[i][j];
                int num = d[pos].vec[k];
                if (num) cnt[k] += num;
                else cnt[k] = -inf;
            }
        }

        rep(k, len) {
            if (cnt[k] > 0) res += 3 + cnt[k] - 6;
        }
    }

    return res * 1000;
}

void inv() {
    rep(i, len) rep(j, i) {
        swap(board[i][j], board[j][i]);
    }
    return;
}

int main()
{
    random_device rnd;
    clock_t start = clock();

    rep(i, n) {
        rep(j, len) {
            int t;
            cin >> t;
            t--;
            d[i].vec[t]++;
        }
        d[i].idx = i;
    }

    int scr = 0;
    {
        vi<bool> used(n);        
        rep(i, len) {
            scr = 0;
            rep(j, len) {
                vi<pii> cnt;
                rep(k, n) {
                    if (used[k]) continue;
                    cnt.push_back({ d[k].vec[j], k });
                }
                sort(cnt.begin(), cnt.end(), greater<pii>{});
                int ts = 0;
                rep(k2, len) ts += cnt[k2].first;
                if (ts > scr) {
                    rep(k2, len) board[i][k2] = cnt[k2].second;
                    scr = ts;
                }
            }

            rep(j, len) used[board[i][j]] = true;
        }
    }

    scr = score();
    double t0 = 1000, t1 = 900, tmp = t0;
    int cn = 0;
    while (clock() - start < TL2) {
        cn++;
        if (cn % 100 == 0) {
            double lt = (double)(clock() - start) / TL2;
            tmp = pow(t0, 1 - lt) * pow(t1, lt);
        }

        if (clock() - start > TL && (rnd() & 1)) inv();

        if (rnd() & 1) {
            int idx0 = rnd() % len;
            int idx1 = (idx0 + rnd() % (len - 1) + 1) % len;
            int row = rnd() % len;     

            swap(board[row][idx0], board[row][idx1]);
            int ts = score();
            bernoulli_distribution bd(min<double>(1, exp((ts - scr) / tmp)));
            if (ts >= scr || bd(rnd)) scr = ts;
            //if (ts >= scr) scr = ts;
            else swap(board[row][idx0], board[row][idx1]);
        }
        else {
            int row = rnd() % len;
            vi<int> td = board[row];
            shuffle(board[row].begin(), board[row].end(), rnd);
            int ts = score();
            bernoulli_distribution bd(min<double>(1, exp((ts - scr) / tmp)));
            if (ts >= scr || bd(rnd)) scr = ts;
            else swap(board[row], td);
        }
    }    

    vi<pii> ans(n);
    rep(i, len) rep(j, len) {
        ans[board[i][j]] = { i + 1, j + 1 };
    }
    for (pii elem : ans) cout << elem.first << " " << elem.second << endl;
    return 0;
}

yukicoder

結果

テストケース

ソースコード