#300546 (C++14) No.5003 物理好きクリッカー

• 問題: No.5003 物理好きクリッカー
• ユーザー: nebukuro09
• 提出日時: 2018-12-01 15:41:53
• 言語: C++14 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 7,962 ms / 10,000 ms
• コード長: 11,315 bytes
• コンパイル時間: 2,236 ms
• 実行使用メモリ: 58,712 KB
• スコア: 149,156,231,963
• 平均クエリ数: 10000.00
• 最終ジャッジ日時: 2021-07-19 07:45:51
• 合計ジャッジ時間: 256,188 ms
• ジャッジサーバーID （参考情報）: judge11 / judge10

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	7,452 ms; 56,428 KB
testcase_01	AC	7,435 ms; 56,312 KB
testcase_02	AC	7,409 ms; 56,468 KB
testcase_03	AC	7,420 ms; 56,096 KB
testcase_04	AC	7,590 ms; 57,500 KB
testcase_05	AC	7,866 ms; 56,436 KB
testcase_06	AC	7,962 ms; 56,536 KB
testcase_07	AC	7,560 ms; 56,260 KB
testcase_08	AC	7,937 ms; 55,952 KB
testcase_09	AC	7,595 ms; 56,208 KB
testcase_10	AC	7,570 ms; 57,004 KB
testcase_11	AC	7,883 ms; 56,712 KB
testcase_12	AC	7,693 ms; 56,464 KB
testcase_13	AC	7,456 ms; 55,972 KB
testcase_14	AC	7,403 ms; 56,356 KB
testcase_15	AC	7,386 ms; 56,452 KB
testcase_16	AC	7,526 ms; 56,208 KB
testcase_17	AC	7,592 ms; 55,748 KB
testcase_18	AC	7,582 ms; 56,712 KB
testcase_19	AC	7,934 ms; 56,288 KB
testcase_20	AC	7,538 ms; 56,940 KB
testcase_21	AC	7,645 ms; 56,888 KB
testcase_22	AC	7,445 ms; 56,940 KB
testcase_23	AC	7,505 ms; 56,932 KB
testcase_24	AC	7,352 ms; 56,376 KB
testcase_25	AC	7,487 ms; 56,536 KB
testcase_26	AC	7,567 ms; 56,352 KB
testcase_27	AC	7,324 ms; 55,536 KB
testcase_28	AC	7,541 ms; 56,568 KB
testcase_29	AC	7,438 ms; 56,260 KB
testcase_30	AC	7,651 ms; 56,632 KB
testcase_31	AC	7,393 ms; 56,068 KB

ソースコード

#include <bits/stdc++.h>
using namespace std;
#define REP(i,n) for (int i=0;i<(n);i++)
#define REP2(i,m,n) for (int i=m;i<(n);i++)
typedef long long ll;
typedef long double ld;

const int N = 10000;
const int F = 5;
const int MAXB = 50;
string STATUS;
string tmp;
int BONUS = N + 10;
int FEVER = N + 10;
int SALE = N + 10;
int FEVER_A[N];
int FEVER_B[N+1];

const string action[F] = {"click", "buy", "sell", "reinforce", "enhclick"};
const string facility[F] = {"hand", "lily", "factory", "casino", "grimoire"};
const uint64_t productivity[F] = {1, 10, 120, 2000, 25000};

uint64_t buy_cost[F][MAXB];
uint64_t rein_cost[F][MAXB];
uint64_t click_cost[MAXB];

const int BEAM_WIDTH = 500;
int prev_id[BEAM_WIDTH * N + 10];
int prev_action[BEAM_WIDTH * N + 10];

const int HASH_BASE = 45;

struct State {
    uint16_t turn;
    uint64_t cookies;
    uint16_t click_level;
    uint16_t facility_count[F];
    uint16_t facility_level[F];
    uint64_t state_score;
    int32_t parent_id;
    int8_t parent_action;

    State() {
        turn = 0;
        cookies = 0;
        click_level = 0;
        memset(facility_count, 0, sizeof(facility_count));
        memset(facility_level, 0, sizeof(facility_level));
        state_score = 0;
        parent_id = -1;
        parent_action = -1;
    }

    uint64_t hash() {
        uint64_t ret = 0;
        ret += click_level;
        REP(i, F) ret *= HASH_BASE, ret += facility_count[i];
        REP(i, F) ret *= HASH_BASE, ret += facility_level[i];
        return ret;
    }

    uint64_t eval() {
        uint64_t ret = 0;
        ret += cookies;
        ret += (N - turn) * (((uint64_t)(1)) << click_level);
        ret += 6 * FEVER_B[turn] * (((uint64_t)(1)) << click_level);
        for (int i = 0; i < F; ++i) {
            ret +=
                (N - turn) *
                productivity[i] *
                (((uint64_t)(1)) << facility_level[i]) *
                facility_count[i];
            ret +=
                6 * FEVER_B[turn] *
                productivity[i] *
                (((uint64_t)(1)) << facility_level[i]) *
                facility_count[i];
        }
        return ret;
    }

    void harvest() {
        for (int i = 0; i < F; ++i) {
            cookies +=
                (FEVER_A[turn] == 1 ? 7 : 1) *
                productivity[i] *
                (((uint64_t)(1)) << facility_level[i]) *
                facility_count[i];
        }
    }

    void harvest_back() {
        for (int i = 0; i < F; ++i) {
            cookies -=
                (FEVER_A[turn] == 1 ? 7 : 1) *
                productivity[i] *
                (((uint64_t)(1)) << facility_level[i]) *
                facility_count[i];
        }
    }

    void click() {
        cookies += (FEVER_A[turn] == 1 ? 7 : 1) * (((uint64_t)(1)) << click_level);
        harvest();
        turn += 1;
    }

    void click_back() {
        turn -= 1;
        harvest_back();
        cookies -= (FEVER_A[turn] == 1 ? 7 : 1) * (((uint64_t)(1)) << click_level);
    }

    bool buy(int f) {
        uint64_t c = buy_cost[f][facility_count[f]];
        if (turn == SALE + 1) c = ceil(c * 0.9);
        if (cookies < c) return false;
        cookies -= c;
        facility_count[f] += 1;
        harvest();
        turn += 1;
        return true;
    }

    void buy_back(int f) {
        turn -= 1;
        harvest_back();
        facility_count[f] -= 1;
        uint64_t c = buy_cost[f][facility_count[f]];
        if (turn == SALE + 1) c = ceil(c * 0.9);
        cookies += c;
    }

    bool sell(int f) {
        if (facility_count[f] == 0) return false;
        cookies += ceil(buy_cost[f][facility_count[f]-1] / 4.0);
        facility_count[f] -= 1;
        harvest();
        turn += 1;
        return true;
    }

    void sell_back(int f) {
        turn -= 1;
        harvest_back();
        facility_count[f] += 1;
        cookies -= ceil(buy_cost[f][facility_count[f]-1] / 4.0);
    }

    bool reinforce(int f) {
        uint64_t c = rein_cost[f][facility_level[f]];
        if (turn == SALE + 1) c = ceil(c * 0.9);
        if (cookies < c) return false;
        cookies -= c;
        facility_level[f] += 1;
        harvest();
        turn += 1;
        return true;
    }

    void reinforce_back(int f) {
        turn -= 1;
        harvest_back();
        facility_level[f] -= 1;
        uint64_t c = rein_cost[f][facility_level[f]];
        if (turn == SALE + 1) c = ceil(c * 0.9);
        cookies += c;
    }

    bool enhclick() {
        uint64_t c = click_cost[click_level];
        if (turn == SALE + 1) c = ceil(c * 0.9);
        if (cookies < c) return false;
        cookies -= c;
        click_level += 1;
        harvest();
        turn += 1;
        return true;
    }

    void enhclick_back() {
        turn -= 1;
        harvest_back();
        click_level -= 1;
        uint64_t c = click_cost[click_level];
        if (turn == SALE + 1) c = ceil(c * 0.9);
        cookies += click_cost[click_level];
    }

    uint64_t eval_click() {
        click();
        uint64_t score = eval();
        click_back();
        return score;
    }

    uint64_t eval_buy(int f) {
        if (!buy(f)) return 0;
        uint64_t score = eval();
        buy_back(f);
        return score;
    }

    uint64_t eval_sell(int f) {
        if (!sell(f)) return 0;
        uint64_t score = eval();
        sell_back(f);
        return score;
    }

    uint64_t eval_reinforce(int f) {
        if (!reinforce(f)) return 0;
        uint64_t score = eval();
        reinforce_back(f);
        return score;
    }

    uint64_t eval_enhclick() {
        if (!enhclick()) return 0;
        uint64_t score = eval();
        enhclick_back();
        return score;
    }

    pair<int, int> best_next_action() { // for greedy search
        uint64_t best_score = 0;
        int x = -1;
        int y = 0;

        uint64_t s;

        s = eval_click();
        if (s > best_score) {
            best_score = s;
            x = 0;
        }

        REP(i, F) {
            s = eval_buy(i);
            if (s > best_score) {
                best_score = s;
                x = 1;
                y = i;
            }
        }

        REP(i, F) {
            s = eval_sell(i);
            if (s > best_score) {
                best_score = s;
                x = 2;
                y = i;
            }
        }

        REP(i, F) {
            s = eval_reinforce(i);
            if (s > best_score) {
                best_score = s;
                x = 3;
                y = i;
            }
        }

        s = eval_enhclick();
        if (s > best_score) {
            best_score = s;
            x = 4;
        }

        return make_pair(x, y);
    }
};

bool operator<(const State& a, const State& b) {
    return a.state_score < b.state_score;
}


void init() {
    buy_cost[0][0] =  150;
    buy_cost[1][0] =  2000;
    buy_cost[2][0] =  30000;
    buy_cost[3][0] =  600000;
    buy_cost[4][0] =  10000000;
    rein_cost[0][0] = buy_cost[0][0] * 10;
    rein_cost[1][0] = buy_cost[1][0] * 10;
    rein_cost[2][0] = buy_cost[2][0] * 10;
    rein_cost[3][0] = buy_cost[3][0] * 10;
    rein_cost[4][0] = buy_cost[4][0] * 10;
    click_cost[0] = 15;

    REP(i, F) REP(j, MAXB-1) {
        buy_cost[i][j+1] = ceil(6.0 * buy_cost[i][j] / 5);
    }

    REP(i, F) REP(j, MAXB-1) {
        rein_cost[i][j+1] = 10 * rein_cost[i][j];
    }

    REP(j, MAXB-1) {
        click_cost[j+1] = click_cost[j] * 10;
    }


    int hoge;
    cin >> hoge >> STATUS;

    REP(i, N) if (STATUS[i] == 'F') {
        for (int j = i + 1; j <= i + 20 && j < N; ++j) {
            FEVER_A[j] = 1;
            FEVER_B[j] = 1;
        }
    }

    for (int i = N - 2; i >= 0; --i) {
        FEVER_B[i] += FEVER_B[i+1];
    }
}

void act(int x, int y) {
    cout << action[x];
    if (x >= 1 && x <= 3) {
        cout << " " << facility[y];
    }
    cout << endl;

    cin >> tmp;
    if (tmp != "ok") {
        exit(EXIT_FAILURE);
    }
}

void greedy_search() {
    State S;

    pair<int, int> p;
    int x, y;

    REP(_, N) {
        p = S.best_next_action();
        x = p.first;
        y = p.second;

        if (x == 0) {
            S.click();
        } else if (x == 1) {
            S.buy(y);
        } else if (x == 2) {
            S.sell(y);
        } else if (x == 3) {
            S.reinforce(y);
        } else {
            S.enhclick();
        }

        act(x, y);
    }
}

void beam_search() {
    State S, T;
    priority_queue<State> pq[2];
    pq[0].push(S);

    int id = -1;
    int cur = 0, tar = 1;

    for (int t = 0; t < N; ++t) {
        pq[tar] = priority_queue<State>();
        set<uint64_t> st;

        for (int k = 0; k < BEAM_WIDTH && !pq[cur].empty(); ++k) {
            S = pq[cur].top();
            pq[cur].pop();

            uint64_t hash = S.hash();
            if (st.find(hash) != st.end()) {
                k -= 1;
                continue;
            }
            st.insert(hash);

            if (BONUS + 1 == t) {
                S.cookies += ceil(1.0 * S.cookies / 100);
            }

            id += 1;
            prev_id[id] = S.parent_id;
            prev_action[id] = S.parent_action;

            T = S;
            T.parent_id = id;

            T.click();
            T.state_score = T.eval();
            T.parent_action = 0;
            pq[tar].push(T);
            T.click_back();

            for (int f = 0; f < F; ++f) {
                if (T.buy(f)) {
                    T.state_score = T.eval();
                    T.parent_action = 10 + f;
                    pq[tar].push(T);
                    T.buy_back(f);
                }
            }

            for (int f = 0; f < F; ++f) {
                if (T.sell(f)) {
                    T.state_score = T.eval();
                    T.parent_action = 20 + f;
                    pq[tar].push(T);
                    T.sell_back(f);
                }
            }

            for (int f = 0; f < F; ++f) {
                if (T.reinforce(f)) {
                    T.state_score = T.eval();
                    T.parent_action = 30 + f;
                    pq[tar].push(T);
                    T.reinforce_back(f);
                }
            }

            if (T.enhclick()) {
                T.state_score = T.eval();
                T.parent_action = 40;
                pq[tar].push(T);
                T.enhclick_back();
            }
        }

        if (STATUS[t] == 'B') {
            BONUS = t;
        } else if (STATUS[t] == 'F') {
            FEVER = t;
        } else if (STATUS[t] == 'S') {
            SALE = t;
        }

        cur ^= 1;
        tar ^= 1;
    }

    State last = pq[cur].top();
    vector<int> actions = {last.parent_action};
    id = last.parent_id;
    //cerr << last.cookies << endl;

    while (prev_action[id] != -1) {
        actions.push_back(prev_action[id]);
        id = prev_id[id];
    }

    reverse(actions.begin(), actions.end());

    int last_buy_or_rein = 0;
    REP(i, N) if (actions[i] / 10 == 1 || actions[i] / 10 == 3) {
        last_buy_or_rein = i;
    }
    sort(actions.begin() + last_buy_or_rein + 1, actions.end());

    for (int i = 0; i < N; ++i) {
        int x = actions[i] / 10;
        int y = actions[i] % 10;
        act(x, y);
    }
}

int main() {
    init();
    beam_search();
}