結果

問題 No.5016 Worst Mayor
ユーザー tempura_pptempura_pp
提出日時 2023-04-29 16:30:20
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 875 ms / 2,000 ms
コード長 5,141 bytes
コンパイル時間 2,374 ms
コンパイル使用メモリ 156,716 KB
実行使用メモリ 38,388 KB
スコア 20,247,431,335
平均クエリ数 400.00
最終ジャッジ日時 2023-04-29 16:31:12
合計ジャッジ時間 48,225 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 831 ms
36,220 KB
testcase_01 AC 800 ms
37,076 KB
testcase_02 AC 815 ms
36,776 KB
testcase_03 AC 846 ms
35,624 KB
testcase_04 AC 836 ms
35,268 KB
testcase_05 AC 827 ms
35,928 KB
testcase_06 AC 855 ms
35,208 KB
testcase_07 AC 856 ms
36,204 KB
testcase_08 AC 835 ms
36,916 KB
testcase_09 AC 851 ms
36,540 KB
testcase_10 AC 867 ms
35,048 KB
testcase_11 AC 828 ms
35,712 KB
testcase_12 AC 839 ms
35,908 KB
testcase_13 AC 827 ms
36,792 KB
testcase_14 AC 838 ms
35,692 KB
testcase_15 AC 851 ms
35,312 KB
testcase_16 AC 825 ms
36,184 KB
testcase_17 AC 823 ms
35,688 KB
testcase_18 AC 852 ms
35,024 KB
testcase_19 AC 845 ms
36,016 KB
testcase_20 AC 827 ms
35,304 KB
testcase_21 AC 848 ms
35,632 KB
testcase_22 AC 825 ms
36,892 KB
testcase_23 AC 848 ms
36,480 KB
testcase_24 AC 859 ms
34,996 KB
testcase_25 AC 835 ms
34,776 KB
testcase_26 AC 848 ms
35,996 KB
testcase_27 AC 850 ms
36,392 KB
testcase_28 AC 854 ms
37,164 KB
testcase_29 AC 838 ms
36,756 KB
testcase_30 AC 855 ms
35,848 KB
testcase_31 AC 826 ms
36,972 KB
testcase_32 AC 860 ms
34,952 KB
testcase_33 AC 860 ms
36,688 KB
testcase_34 AC 839 ms
35,180 KB
testcase_35 AC 858 ms
35,124 KB
testcase_36 AC 861 ms
35,324 KB
testcase_37 AC 823 ms
36,372 KB
testcase_38 AC 842 ms
36,140 KB
testcase_39 AC 857 ms
34,900 KB
testcase_40 AC 808 ms
36,568 KB
testcase_41 AC 847 ms
35,420 KB
testcase_42 AC 875 ms
38,388 KB
testcase_43 AC 831 ms
35,812 KB
testcase_44 AC 859 ms
34,972 KB
testcase_45 AC 856 ms
37,260 KB
testcase_46 AC 826 ms
35,072 KB
testcase_47 AC 864 ms
36,772 KB
testcase_48 AC 871 ms
36,184 KB
testcase_49 AC 840 ms
35,524 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, 85) {
        int t = i == 0 ? 400 : 20;
        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(_, 20) {
            auto p = pq.top();
            pq.pop();
            int s = p.first, id = p.second.first, e = p.second.second;
            auto ns = State::getNewState(bests[id], 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 < 85 ? state.up[num] : cur;
    int ben = 60 * (ma - cur) * (T - t);
    int pay = (int)(10000000 / sqrt(v));
    cerr << ben << " " << pay << endl;
    if(num < 83 && 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(85 - 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