結果

問題 No.5016 Worst Mayor
ユーザー assy1028assy1028
提出日時 2023-04-29 15:35:50
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,181 ms / 2,000 ms
コード長 6,582 bytes
コンパイル時間 1,371 ms
コンパイル使用メモリ 131,340 KB
実行使用メモリ 30,384 KB
スコア 8,591,891,545
平均クエリ数 400.00
最終ジャッジ日時 2023-04-29 15:36:48
合計ジャッジ時間 47,674 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge13
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 853 ms
30,156 KB
testcase_01 AC 772 ms
29,900 KB
testcase_02 AC 796 ms
30,296 KB
testcase_03 AC 748 ms
30,284 KB
testcase_04 AC 906 ms
29,972 KB
testcase_05 AC 908 ms
29,704 KB
testcase_06 AC 787 ms
29,684 KB
testcase_07 AC 952 ms
30,196 KB
testcase_08 AC 876 ms
30,120 KB
testcase_09 AC 820 ms
30,080 KB
testcase_10 AC 988 ms
29,652 KB
testcase_11 AC 860 ms
29,956 KB
testcase_12 AC 855 ms
30,236 KB
testcase_13 AC 642 ms
30,020 KB
testcase_14 AC 777 ms
30,244 KB
testcase_15 AC 936 ms
30,316 KB
testcase_16 AC 807 ms
29,756 KB
testcase_17 AC 779 ms
29,844 KB
testcase_18 AC 859 ms
30,032 KB
testcase_19 AC 1,181 ms
30,280 KB
testcase_20 AC 801 ms
29,904 KB
testcase_21 AC 829 ms
29,640 KB
testcase_22 AC 834 ms
29,752 KB
testcase_23 AC 831 ms
29,848 KB
testcase_24 AC 935 ms
30,328 KB
testcase_25 AC 837 ms
30,332 KB
testcase_26 AC 905 ms
30,224 KB
testcase_27 AC 918 ms
29,948 KB
testcase_28 AC 960 ms
30,028 KB
testcase_29 AC 855 ms
29,892 KB
testcase_30 AC 840 ms
29,832 KB
testcase_31 AC 929 ms
30,312 KB
testcase_32 AC 807 ms
29,648 KB
testcase_33 AC 863 ms
29,556 KB
testcase_34 AC 904 ms
29,856 KB
testcase_35 AC 862 ms
29,740 KB
testcase_36 AC 666 ms
29,628 KB
testcase_37 AC 897 ms
29,636 KB
testcase_38 AC 725 ms
29,976 KB
testcase_39 AC 803 ms
29,968 KB
testcase_40 AC 869 ms
29,712 KB
testcase_41 AC 751 ms
30,228 KB
testcase_42 AC 763 ms
30,284 KB
testcase_43 AC 914 ms
29,932 KB
testcase_44 AC 749 ms
29,988 KB
testcase_45 AC 820 ms
30,384 KB
testcase_46 AC 917 ms
29,864 KB
testcase_47 AC 831 ms
29,804 KB
testcase_48 AC 966 ms
30,116 KB
testcase_49 AC 845 ms
30,208 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <algorithm>
#include <iostream>
#include <cstdio>
#include <map>
#include <numeric>
#include <cmath>
#include <set>
#include <sstream>
#include <string>
#include <vector>
#include <queue>
#include <stack>
#include <complex>
#include <string.h>
#include <unordered_set>
#include <unordered_map>
#include <bitset>
#include <iomanip>
#include <sys/time.h>
#include <tuple>
#include <random>
using namespace std;

#define endl '\n'
#define ALL(v) (v).begin(), (v).end()
#define RALL(v) (v).rbegin(), (v).rend()
#define UNIQ(v) (v).erase(unique((v).begin(), (v).end()), (v).end())

typedef long long ll;
typedef long double ld;
typedef pair<int, int> P;
typedef complex<double> comp;
typedef vector< vector<ld> > matrix;
struct pairhash {
public:
    template<typename T, typename U>
    size_t operator()(const pair<T, U> &x) const {
    size_t seed = hash<T>()(x.first);
    return hash<U>()(x.second) + 0x9e3779b9 + (seed<<6) + (seed>>2);
    }
};
const int INF = 1e9 + 9;
//const ll INF = 1e18 + 9;
const ll MOD = 1e9 + 7;
//const ll MOD = 998244353;
const double EPS = 1e-8;
const double PI = acos(-1);

unsigned long xor128(void) {
    static unsigned long x=123456789,y=362436069,z=521288629,w=88675123;
    unsigned long t;
    t=(x^(x<<11));x=y;y=z;z=w; return( w=(w^(w>>19))^(t^(t>>8)) );
}

// [0, k)
int rand(int k) {
    return xor128() % k;
}

// [a, b)
int rand(int a, int b) {
    return a + xor128() % (b - a);
}

int n, t;
int a[3010], b[3010], c[3010], d[3010];
const int MAX_Y = 14;
const int MAX_X = 14;

double route_weight[15][15][2];
double weight_each[3010][15][15][2];
bool built[15][15][2];

struct road {
    double count;
    int x, y, t;
};
bool operator< (const road& r1, const road& r2) {
    return r1.count < r2.count;
}

vector<road> high_way;
bool prev_build = false;

void build(const road& r) {
    high_way.push_back(r);
    built[r.x][r.y][r.t] = true;
    prev_build = true;
    cout << "1 " << r.x+1 << " " << r.y+1 << " " << r.x+2-r.t << " " << r.y+r.t+1 << endl; cout.flush();
}

void get_coop() {
    prev_build = false;
    cout << "2" << endl; cout.flush();
}

void get_money() {
    prev_build = false;
    cout << "3" << endl; cout.flush();
}

bool in_area(const int idx, const road& r) {
    const int minY = min(a[idx], c[idx]);
    const int maxY = max(a[idx], c[idx]);
    const int minX = min(b[idx], d[idx]);
    const int maxX = max(b[idx], d[idx]);
    return minY <= r.x && r.x <= maxY && minY <= r.x+1-r.t && r.x+1-r.t <= maxY &&
           minX <= r.y && r.y <= maxX && minX <= r.y+r.t && r.y+r.t <= maxX;
}

void add_weight(const int idx) {
    const int minY = min(a[idx], c[idx]);
    const int maxY = max(a[idx], c[idx]);
    const int minX = min(b[idx], d[idx]);
    const int maxX = max(b[idx], d[idx]);
    for (int y = minY; y <= maxY; y++) {
        for (int x = minX; x <= maxX; x++) {
            route_weight[y][x][0] += weight_each[idx][y][x][0];
            route_weight[y][x][1] += weight_each[idx][y][x][1];
        }
    }
}

void subtract_weight(const int idx) {
    const int minY = min(a[idx], c[idx]);
    const int maxY = max(a[idx], c[idx]);
    const int minX = min(b[idx], d[idx]);
    const int maxX = max(b[idx], d[idx]);
    for (int y = minY; y <= maxY; y++) {
        for (int x = minX; x <= maxX; x++) {
            route_weight[y][x][0] -= weight_each[idx][y][x][0];
            route_weight[y][x][1] -= weight_each[idx][y][x][1];
        }
    }
}

void add_weight_each(const int idx, const int minY, const int maxY, const int minX, const int maxX, const double weight) {
    for (int y = minY; y <= maxY; y++) {
        for (int x = minX; x <= maxX; x++) {
            if (y < maxY) weight_each[idx][y][x][0] += weight;
            if (x < maxX) weight_each[idx][y][x][1] += weight;
        }
    }
}

void calc_weight_each(const int idx) {
    memset(weight_each[idx], 0, sizeof(weight_each[idx]));
    int cnt = 0;
    vector<road> h;
    for (const auto& r : high_way) {
        if (in_area(idx, r)) {
            cnt++;
            h.push_back(r);
        }
    }
    const double weight = cnt == 0 ? 1.0 : 1.0 / cnt;
    const int minY = min(a[idx], c[idx]);
    const int maxY = max(a[idx], c[idx]);
    const int minX = min(b[idx], d[idx]);
    const int maxX = max(b[idx], d[idx]);
    if (cnt == 0) {
        add_weight_each(idx, minY, maxY, minX, maxX, weight);
    } else if ((a[idx] - c[idx]) * (b[idx] - d[idx]) > 0){
        for (const auto& r : h) {
            add_weight_each(idx, minY, r.x, minX, r.y, weight);
            add_weight_each(idx, r.x+1-r.t, maxY, r.y+r.t, maxX, weight);
        }
    } else {
        for (const auto& r : h) {
            add_weight_each(idx, r.x+1-r.t, maxY, minX, r.y, weight);
            add_weight_each(idx, minY, r.x, r.y+r.t, maxX, weight);
        }
    }
}

void init() {
    for (int i = 0; i < n; i++) {
        calc_weight_each(i);
        add_weight(i);
    }
}

road get_candidate() {
    if (prev_build) {
        const auto& r = high_way.back();
        for (int i = 0; i < n; i++) {
            if (in_area(i, r)) {
                subtract_weight(i);
                calc_weight_each(i);
                add_weight(i);
            }
        }
    }
    
    double max_weight = 0;
    road ret;
    for (int y = 0; y < MAX_Y; y++) {
        for (int x = 0; x < MAX_X; x++) {
            if (y < MAX_Y-1 && !built[y][x][0] && route_weight[y][x][0] > max_weight) {
                ret = road{route_weight[y][x][0], y, x, 0};
                max_weight = route_weight[y][x][0];
            }
            if (x < MAX_X-1 && !built[y][x][1] && route_weight[y][x][1] > max_weight) {
                ret = road{route_weight[y][x][1], y, x, 1};
                max_weight = route_weight[y][x][1];
            }
        }
    }
    return ret;
}

void solve() {
    init();
    int u, v;
    for (int iter = 0; iter < t; iter++) {
        cin >> u >> v;
        const int cost = 1e7 / sqrt(v);
        const auto& r = get_candidate();
        const int gain = r.count * (t-iter) * 60;// * pow(0.995, iter);
        //cerr << "??? " << u << " " << v << " " << cost << " " << gain << endl;
        if (u < cost || gain < cost) {
            get_coop();
        } else {
            build(r);
        }
    }
}

void input() {
    cin >> n >> t;
    for (int i = 0; i < n; i++) {
        cin >> a[i] >> b[i] >> c[i] >> d[i];
        a[i]--;
        b[i]--;
        c[i]--;
        d[i]--;
    }
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout << fixed << setprecision(15);

    input();
    solve();
}
0