結果

問題 No.5013 セクスタプル (open)
ユーザー NAMIDAIRONAMIDAIRO
提出日時 2022-12-29 23:02:12
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 3,532 bytes
コンパイル時間 1,343 ms
実行使用メモリ 6,952 KB
スコア 17,041
最終ジャッジ日時 2022-12-29 23:05:40
合計ジャッジ時間 207,519 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge15
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 99 TLE * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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 = 1000 * 1000;
const int TL2 = 1920 * 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();

         {
            int idx0 = rnd() % len;
            int idx1 = (idx0 + rnd() % (len - 1) + 1) % len;
            int row0 = rnd() % len;   
            int row1 = (row0 + rnd() % (len - 1) + 1) % len;

            swap(board[row0][idx0], board[row1][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[row0][idx0], board[row1][idx1]);
        }
    }    

    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;
}
0