結果

問題 No.5016 Worst Mayor
ユーザー tempura_pptempura_pp
提出日時 2023-04-29 16:56:52
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,842 ms / 2,000 ms
コード長 5,198 bytes
コンパイル時間 2,412 ms
コンパイル使用メモリ 157,340 KB
実行使用メモリ 41,076 KB
スコア 20,950,504,453
平均クエリ数 400.00
最終ジャッジ日時 2023-04-29 17:01:14
合計ジャッジ時間 97,628 ms
ジャッジサーバーID
（参考情報） 		judge16 / judge14
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1,822 ms
40,224 KB
testcase_01 AC 1,800 ms
40,456 KB
testcase_02 AC 1,787 ms
40,584 KB
testcase_03 AC 1,820 ms
40,320 KB
testcase_04 AC 1,760 ms
40,392 KB
testcase_05 AC 1,788 ms
40,468 KB
testcase_06 AC 1,782 ms
40,064 KB
testcase_07 AC 1,795 ms
40,800 KB
testcase_08 AC 1,781 ms
40,680 KB
testcase_09 AC 1,793 ms
40,060 KB
testcase_10 AC 1,781 ms
40,036 KB
testcase_11 AC 1,791 ms
40,784 KB
testcase_12 AC 1,747 ms
40,032 KB
testcase_13 AC 1,770 ms
40,200 KB
testcase_14 AC 1,792 ms
40,352 KB
testcase_15 AC 1,817 ms
40,540 KB
testcase_16 AC 1,763 ms
40,744 KB
testcase_17 AC 1,791 ms
40,600 KB
testcase_18 AC 1,818 ms
40,104 KB
testcase_19 AC 1,773 ms
40,556 KB
testcase_20 AC 1,784 ms
40,508 KB
testcase_21 AC 1,784 ms
40,468 KB
testcase_22 AC 1,791 ms
40,224 KB
testcase_23 AC 1,806 ms
40,500 KB
testcase_24 AC 1,776 ms
40,164 KB
testcase_25 AC 1,781 ms
40,104 KB
testcase_26 AC 1,801 ms
40,756 KB
testcase_27 AC 1,796 ms
40,492 KB
testcase_28 AC 1,814 ms
39,804 KB
testcase_29 AC 1,786 ms
40,532 KB
testcase_30 AC 1,842 ms
39,848 KB
testcase_31 AC 1,815 ms
40,680 KB
testcase_32 AC 1,810 ms
40,384 KB
testcase_33 AC 1,827 ms
40,280 KB
testcase_34 AC 1,819 ms
40,136 KB
testcase_35 AC 1,798 ms
40,008 KB
testcase_36 AC 1,782 ms
40,496 KB
testcase_37 AC 1,781 ms
41,076 KB
testcase_38 AC 1,794 ms
40,488 KB
testcase_39 AC 1,778 ms
39,896 KB
testcase_40 AC 1,782 ms
40,512 KB
testcase_41 AC 1,780 ms
40,640 KB
testcase_42 AC 1,781 ms
39,944 KB
testcase_43 AC 1,791 ms
40,460 KB
testcase_44 AC 1,818 ms
40,112 KB
testcase_45 AC 1,815 ms
40,620 KB
testcase_46 AC 1,807 ms
40,728 KB
testcase_47 AC 1,779 ms
39,924 KB
testcase_48 AC 1,818 ms
40,376 KB
testcase_49 AC 1,814 ms
40,460 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <algorithm>
#include <assert.h>
#include <bitset>
#include <complex>
#include <deque>
#include <functional>
#include <iomanip>
#include <iostream>
#include <map>
#include <math.h>
#include <numeric>
#include <queue>
#include <set>
#include <stack>
#include <string>
#include <vector>
using namespace std;
#define REP(i, m, n) for(int i = (int)(m); i < (int)(n); ++i)
#define rep(i, n) REP(i, 0, n)
using ll = long long;
constexpr int inf = 1e9 + 7;
constexpr ll longinf = 1LL << 60;
constexpr ll mod = 1e9 + 7;

struct Xor128 {
    unsigned x, y, z, w;

    Xor128(unsigned w_ = 88675123) : x(123456789), y(362436069), z(521288629), w(w_){};

    inline unsigned xor128() {
        unsigned t;
        t = x ^ (x << 11);
        x = y;
        y = z;
        z = w;
        return w = (w ^ (w >> 19)) ^ (t ^ (t >> 8));
    }

    int nextInt(int x, int y) {
        if(x > y)
            swap(x, y);
        return xor128() % (y - x) + x;
    }

    double nextDouble(double a, double b) {
        return (double)(xor128() & 0xffff) / 0xffff * (b - a) + a;
    }
};

int dist[200][200];
int N = 3000, T = 400, M = 196;
int cost[3030];
int s[3030], t[3030];
int NC = 100000, NH = 22301;

void input() {
    int _;
    cin >> _ >> _;
    rep(i, N) {
        int a, b, c, d;
        cin >> a >> b >> c >> d;
        --a;
        --b;
        --c;
        --d;
        s[i] = a * 14 + b;
        t[i] = c * 14 + d;
    }
}

pair<int, int> etov(int e) {
    int x, y;
    if(e < 13 * 14) {
        x = (e / 13) * 14 + e % 13;
        y = x + 1;
    } else {
        e -= 13 * 14;
        x = e;
        y = e + 14;
    }
    return {x, y};
}

struct State {
    int score;
    int dist[196][196];
    vector<int> result;
    vector<int> up;
    State() : result(0), score(0), up(0) {
        rep(i, 196) rep(j, 196) {
            dist[i][j] = NC * 50;
        }
        rep(i, 14) rep(j, 14) {
            int x = i * 14 + j;
            if(j > 0) {
                dist[x][x - 1] = NC;
            }
            if(j < 13) {
                dist[x][x + 1] = NC;
            }
            if(i > 0) {
                dist[x][x - 14] = NC;
            }
            if(i < 13) {
                dist[x][x + 14] = NC;
            }
            dist[x][x] = 0;
        }
        rep(i, M) rep(j, M) rep(k, M) {
            dist[i][j] = min(dist[i][j], dist[i][k] + dist[k][j]);
        }
    }
    static State getNewState(State &s, int e, int score) {
        s.result.push_back(e);
        s.up.push_back(score - s.score);
        auto [v1, v2] = etov(e);
        s.dist[v1][v2] = NH;
        s.dist[v2][v1] = NH;
        rep(i, M) rep(j, M) {
            s.dist[i][j] = min(s.dist[i][j], min(s.dist[i][v1] + s.dist[v2][j] + NH, s.dist[i][v2] + s.dist[v1][j] + NH));
        }
        s.score = score;
        return s;
    }
    int getNextScore(int e) {
        int res = 0;
        auto [v1, v2] = etov(e);
        rep(i, N) {
            int nc = min(dist[s[i]][t[i]], min(dist[s[i]][v1] + dist[v2][t[i]] + NH, dist[s[i]][v2] + dist[v1][t[i]] + NH));
            res += nc % 100;
        }
        return res + score;
    }
};

auto rng = Xor128();
State search() {
    State s;
    vector<State> bests;
    bests.emplace_back(s);
    rep(i, 90) {
        int t = i == 0 ? 400 : 40;
        priority_queue<pair<int, pair<int, int>>> pq;
        rep(c, bests.size()) {
            rep(_, t) {
                int e = rng.nextInt(0, 2 * 14 * 13);
                int ns = bests[c].getNextScore(e);
                pq.push({ns, {-c, e}});
            }
        }
        vector<State> nb;
        rep(_, 27) {
            auto p = pq.top();
            pq.pop();
            int s = p.first, id = -p.second.first, e = p.second.second;
            State state = bests[id];
            auto ns = State::getNewState(state, e, s);
            nb.emplace_back(ns);
        }
        nb.swap(bests);
    }
    return bests[0];
}

int num = 0;
void output(int e) {
    auto [x, y] = etov(e);
    cout << 1 << " ";
    cout << x / 14 + 1 << " " << x % 14 + 1 << " ";
    cout << y / 14 + 1 << " " << y % 14 + 1 << endl;
}

int cur = 0;
int U = 1000000, V = 1;
pair<int, int> read() {
    // cerr << cur << " " << U << " " << V << endl;
    // return {U, V};
    int u, v;
    cin >> u >> v;
    return {u, v};
}

void turn(int t, State &state) {
    auto p = read();
    int u = p.first, v = p.second;
    int ma = num < 90 ? state.up[num] : cur;
    int ben = 60 * (ma - cur) * (T - t);
    int pay = (int)(10000000 / sqrt(v));
    cerr << ben << " " << pay << endl;
    if((ben > pay || num < 75) && u > pay) {
        output(state.result[num]);
        ++num;
        U = u - pay;
        cur = ma;
    } else {
        int npay = (int)(10000000 / sqrt(v + 1));
        int dpay = max(0, min(90 - num, T - t - 1)) * (pay - npay);
        if(dpay > 50000) {
            cout << 2 << endl;
            V = v + 1;
        } else {
            cout << 3 << endl;
            U = u + 50000;
        }
    }
    U += cur * 60;
}

int main() {
    input();
    State state = search();
    rep(i, T) {
        turn(i, state);
    }
    cerr << U << " " << V << endl;
    cerr << num << endl;
    return 0;
}
0