Submission #554758

# Submission time Handle Problem Language Result Execution time Memory
554758 2022-04-29T11:16:14 Z alextodoran The Forest of Fangorn (CEOI14_fangorn) C++17
80 / 100
1252 ms 65536 KB
/**
 ____ ____ ____ ____ ____
||a |||t |||o |||d |||o ||
||__|||__|||__|||__|||__||
|/__\|/__\|/__\|/__\|/__\|

**/

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef long double ld;

const int T_MAX = 2000;
const int C_MAX = 10000;

struct Point {
    ll x, y;
};

ll W, H;
int T, C;
Point me;
Point trees[T_MAX + 2];
Point camps[C_MAX + 2];

ld getlen (const Point &a, const Point &b) {
    ld x, y;
    if (a.y == b.y) {
        x = -1;
    } else if (a.y > b.y) {
        x = a.x + (ld) (b.x - a.x) / (b.y - a.y) * (0 - a.y);
    } else {
        x = a.x + (ld) (b.x - a.x) / (b.y - a.y) * (H - a.y);
    }
    if (a.x == b.x) {
        y = -1;
    } else if (a.x > b.x) {
        y = a.y + (ld) (b.y - a.y) / (b.x - a.x) * (0 - a.x);
    } else {
        y = a.y + (ld) (b.y - a.y) / (b.x - a.x) * (W - a.x);
    }
    if (!(0 <= x && x <= W)) {
        return (a.x < b.x ? W + y : W + H + W + (H - y));
    } else {
        return (a.y < b.y ? W + H + (W - x) : x);
    }
}

ll clen[C_MAX + 2];
int cperm[C_MAX + 2];

int id[T_MAX + 2][C_MAX + 2];
int idme[T_MAX + 2];

int tperm[T_MAX + 2];
void splitBy (int t) {
    for (int t1 = 1; t1 <= T; t1++) {
        trees[t1].x = trees[t].x * 2 - trees[t1].x;
        trees[t1].y = trees[t].y * 2 - trees[t1].y;
    }
    iota(tperm + 1, tperm + T + 1, 1);
    swap(tperm[t], tperm[T]);
    sort(tperm + 1, tperm + (T - 1) + 1, [&] (const int &t1, const int &t2) {
        return getlen(trees[t], trees[t1]) < getlen(trees[t], trees[t2]);
    });
    for (int ci = 1, ti = 0; ci <= C; ci++) {
        while (ti + 1 < T && getlen(trees[t], trees[tperm[ti + 1]]) < clen[cperm[ci]]) {
            ti++;
        }
        id[t][cperm[ci]] = (ti > 0 ? ti : T - 1);
    }
    for (int ti = 0; ti + 1 < T; ti++) {
        if (getlen(trees[t], trees[tperm[ti + 1]]) < getlen(trees[t], me)) {
            idme[t] = ti + 1;
        } else {
            break;
        }
    }
    if (idme[t] == 0) {
        idme[t] = T - 1;
    }
    for (int t1 = 1; t1 <= T; t1++) {
        trees[t1].x = trees[t].x * 2 - trees[t1].x;
        trees[t1].y = trees[t].y * 2 - trees[t1].y;
    }
}

int main () {
    ios_base::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);

    cin >> W >> H;
    cin >> me.x >> me.y;
    cin >> C;
    for (int c = 1; c <= C; c++) {
        cin >> camps[c].x >> camps[c].y;
    }
    cin >> T;
    for (int t = 1; t <= T; t++) {
        cin >> trees[t].x >> trees[t].y;
    }

    for (int c = 1; c <= C; c++) {
        if (camps[c].y == 0) {
            clen[c] = camps[c].x;
        } else if (camps[c].x == W) {
            clen[c] = W + camps[c].y;
        } else if (camps[c].y == H) {
            clen[c] = W + H + (W - camps[c].x);
        } else {
            clen[c] = W + H + W + (H - camps[c].y);
        }
    }
    iota(cperm + 1, cperm + C + 1, 1);
    sort(cperm + 1, cperm + C + 1, [&] (const int &c1, const int &c2) {
        return clen[c1] < clen[c2];
    });

    for (int t = 1; t <= T; t++) {
        splitBy(t);
    }

    vector <int> reach;
    for (int c = 1; c <= C; c++) {
        bool ok = true;
        for (int t = 1; t <= T; t++) {
            if (id[t][c] != idme[t]) {
                ok = false;
                break;
            }
        }
        if (ok == true) {
            reach.push_back(c);
        }
    }
    cout << (int) reach.size() << "\n";
    for (int c : reach) {
        cout << c << " ";
    }
    cout << "\n";

    return 0;
}
# Verdict Execution time Memory Grader output
1 Correct 1 ms 340 KB Output is correct
2 Correct 1 ms 340 KB Output is correct
3 Correct 1 ms 340 KB Output is correct
4 Correct 1 ms 340 KB Output is correct
5 Correct 1 ms 340 KB Output is correct
6 Correct 1 ms 324 KB Output is correct
7 Correct 2 ms 596 KB Output is correct
8 Correct 2 ms 596 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 460 KB Output is correct
2 Correct 1 ms 468 KB Output is correct
3 Correct 2 ms 596 KB Output is correct
4 Correct 2 ms 596 KB Output is correct
5 Correct 3 ms 596 KB Output is correct
6 Correct 5 ms 852 KB Output is correct
7 Correct 1 ms 280 KB Output is correct
8 Correct 1 ms 340 KB Output is correct
9 Correct 2 ms 584 KB Output is correct
10 Correct 8 ms 1084 KB Output is correct
11 Correct 13 ms 1100 KB Output is correct
12 Correct 10 ms 1096 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 340 KB Output is correct
2 Correct 1 ms 444 KB Output is correct
3 Correct 1 ms 468 KB Output is correct
4 Correct 284 ms 4348 KB Output is correct
5 Correct 58 ms 2240 KB Output is correct
6 Correct 1216 ms 8504 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1252 ms 17844 KB Output is correct
2 Runtime error 1153 ms 65536 KB Execution killed with signal 9
3 Halted 0 ms 0 KB -