結果

問題 No.5016 Worst Mayor
ユーザー assy1028assy1028
提出日時 2023-04-29 15:43:06
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 546 ms / 2,000 ms
コード長 6,579 bytes
コンパイル時間 1,932 ms
コンパイル使用メモリ 130,788 KB
実行使用メモリ 30,768 KB
スコア 10,910,284,428
平均クエリ数 400.00
最終ジャッジ日時 2023-04-29 15:43:40
合計ジャッジ時間 26,772 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge13
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 489 ms
30,060 KB
testcase_01 AC 364 ms
30,008 KB
testcase_02 AC 431 ms
30,156 KB
testcase_03 AC 411 ms
29,904 KB
testcase_04 AC 420 ms
30,432 KB
testcase_05 AC 475 ms
30,092 KB
testcase_06 AC 447 ms
30,272 KB
testcase_07 AC 442 ms
30,560 KB
testcase_08 AC 417 ms
30,304 KB
testcase_09 AC 449 ms
30,196 KB
testcase_10 AC 494 ms
30,260 KB
testcase_11 AC 395 ms
30,632 KB
testcase_12 AC 415 ms
30,568 KB
testcase_13 AC 337 ms
30,708 KB
testcase_14 AC 417 ms
30,432 KB
testcase_15 AC 423 ms
30,044 KB
testcase_16 AC 405 ms
30,584 KB
testcase_17 AC 407 ms
30,640 KB
testcase_18 AC 432 ms
29,972 KB
testcase_19 AC 546 ms
30,532 KB
testcase_20 AC 401 ms
30,072 KB
testcase_21 AC 411 ms
30,360 KB
testcase_22 AC 450 ms
30,128 KB
testcase_23 AC 399 ms
30,392 KB
testcase_24 AC 456 ms
30,696 KB
testcase_25 AC 449 ms
30,184 KB
testcase_26 AC 443 ms
30,512 KB
testcase_27 AC 405 ms
30,160 KB
testcase_28 AC 469 ms
30,028 KB
testcase_29 AC 457 ms
30,060 KB
testcase_30 AC 458 ms
30,696 KB
testcase_31 AC 474 ms
30,168 KB
testcase_32 AC 427 ms
30,720 KB
testcase_33 AC 435 ms
30,296 KB
testcase_34 AC 437 ms
30,768 KB
testcase_35 AC 414 ms
30,676 KB
testcase_36 AC 373 ms
30,244 KB
testcase_37 AC 476 ms
30,024 KB
testcase_38 AC 378 ms
30,108 KB
testcase_39 AC 438 ms
30,696 KB
testcase_40 AC 427 ms
30,620 KB
testcase_41 AC 411 ms
30,216 KB
testcase_42 AC 403 ms
30,112 KB
testcase_43 AC 433 ms
30,116 KB
testcase_44 AC 377 ms
30,644 KB
testcase_45 AC 431 ms
30,676 KB
testcase_46 AC 450 ms
30,692 KB
testcase_47 AC 458 ms
30,116 KB
testcase_48 AC 466 ms
30,452 KB
testcase_49 AC 411 ms
30,512 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