Problem:
Score: 600 points

Problem Statement

M-kun is a brilliant air traffic controller.

On the display of his radar, there are N airplanes numbered 1, 2, ..., N, all flying at the same altitude.

Each of the airplanes flies at a constant speed of 0.1 per second in a constant direction. The current coordinates of the airplane numbered i are (X_i, Y_i), and the direction of the airplane is as follows:

if U_i is U, it flies in the positive y direction;

if U_i is R, it flies in the positive x direction;

if U_i is D, it flies in the negative y direction;

if U_i is L, it flies in the negative x direction.

To help M-kun in his work, determine whether there is a pair of airplanes that will collide with each other if they keep flying as they are now.

If there is such a pair, find the number of seconds after which the first collision will happen.

We assume that the airplanes are negligibly small so that two airplanes only collide when they reach the same coordinates simultaneously.

Constraints

1 \leq N \leq 200000

0 \leq X_i, Y_i \leq 200000

U_i is U, R, D, or L.

The current positions of the N airplanes, (X_i, Y_i), are all distinct.

N, X_i, and Y_i are integers.

Input

Input is given from Standard Input in the following format:

N
X_1 Y_1 U_1
X_2 Y_2 U_2
X_3 Y_3 U_3
:
X_N Y_N U_N

Output

If there is a pair of airplanes that will collide with each other if they keep flying as they are now, print an integer representing the number of seconds after which the first collision will happen.

If there is no such pair, print SAFE.

Sample Input 1

2
11 1 U
11 47 D

Sample Output 1

230

If the airplanes keep flying as they are now, two airplanes numbered 1 and 2 will reach the coordinates (11, 24) simultaneously and collide.

Sample Input 2

4
20 30 U
30 20 R
20 10 D
10 20 L

Sample Output 2

SAFE

No pair of airplanes will collide.

Sample Input 3

8
168 224 U
130 175 R
111 198 D
121 188 L
201 116 U
112 121 R
145 239 D
185 107 L

Sample Output 3

100

Solution in C++ (GCC 9.2.1) (judge verdict: Accepted)

// 20:44 - 21:13
#include <bits/stdc++.h>
using namespace std;

#define fi first
#define se second
#define lch (o << 1)
#define rch (o << 1 | 1)

typedef double db;
typedef long long ll;
typedef unsigned int ui;
typedef pair<int, int> pint;

const int N = 2e5 + 5;
const int INF = 0x3f3f3f3f;
const ll INF_LL = 0x3f3f3f3f3f3f3f3f;

struct Point {
	int x, y, dir;
};
Point p[N];
vector<Point> diag[N * 2][4], sdiag[N * 2][4], xaxis[N * 2][4], yaxis[N * 2][4];

const int DIR[] = {3, 2, 1, 0};

int main() {
	ios::sync_with_stdio(0);

	int n;
	cin >> n;
	for (int i = 1; i <= n; i++) {
		int x, y;
		char c;
		cin >> x >> y >> c;
		int dir = 0;
		if (c == 'U')
			dir = 0;
		else if (c == 'R')
			dir = 1;
		else if (c == 'D')
			dir = 2;
		else
			dir = 3;
		p[i] = (Point){x, y, dir};

		diag[y - x + N][dir].push_back(p[i]);
		sdiag[x + y][dir].push_back(p[i]);
		if (dir % 2 == 1)
			xaxis[y][dir].push_back(p[i]);
		else
			yaxis[x][dir].push_back(p[i]);
	}

	auto CmpByX = [](Point a, Point b) {
		return a.x < b.x;
	};
	auto CmpByY = [](Point a, Point b) {
		return a.y < b.y;
	};

	for (int d = 0; d < 4; d++) {
		for (int i = 0; i < 2 * N; i++) {
			sort(diag[i][d].begin(), diag[i][d].end(), CmpByX);
			sort(sdiag[i][d].begin(), sdiag[i][d].end(), CmpByX);
			sort(xaxis[i][d].begin(), xaxis[i][d].end(), CmpByX);
			sort(yaxis[i][d].begin(), yaxis[i][d].end(), CmpByY);
		}
	}

	ll minT = INF_LL;
	for (int i = 1; i <= n; i++) {
		int x = p[i].x, y = p[i].y;
		int dir = p[i].dir;
		int ndir = (dir + 2) % 4;
		int diagDir = DIR[dir];
		int sdiagDir = (diagDir + 2) % 4;

		auto o = lower_bound(diag[y - x + N][diagDir].begin(), diag[y - x + N][diagDir].end(), p[i], CmpByX);
		auto q = lower_bound(sdiag[x + y][sdiagDir].begin(), sdiag[x + y][sdiagDir].end(), p[i], CmpByX);

		if ((diagDir == 2 || diagDir == 3) && o != diag[y - x + N][diagDir].end())
			minT = min(minT, abs(o->x - x) * 10LL);
		if ((diagDir == 0 || diagDir == 1) && o != diag[y - x + N][diagDir].begin())
			minT = min(minT, abs((o - 1)->x - x) * 10LL);

		if ((sdiagDir == 3 || sdiagDir == 0) && q != sdiag[x + y][sdiagDir].end())
			minT = min(minT, abs(q->x - x) * 10LL);
		if ((sdiagDir == 1 || sdiagDir == 2) && q != sdiag[x + y][sdiagDir].begin())
			minT = min(minT, abs((q - 1)->x - x) * 10LL);

		// parallel to x axis
		if (dir % 2 == 1) {
			auto r = lower_bound(xaxis[y][ndir].begin(), xaxis[y][ndir].end(), p[i], CmpByX);
			if (dir == 1 && r != xaxis[y][ndir].end())
				minT = min(minT, abs(r->x - x) * 5LL);
			if (dir == 3 && r != xaxis[y][ndir].begin())
				minT = min(minT, abs((r - 1)->x - x) * 5LL);
		}

		// parallel to y axis
		if (dir % 2 == 0) {
			auto r = lower_bound(yaxis[x][ndir].begin(), yaxis[x][ndir].end(), p[i], CmpByY);
			if (dir == 0 && r != yaxis[x][ndir].end())
				minT = min(minT, abs(r->y - y) * 5LL);
			if (dir == 2 && r != yaxis[x][ndir].begin())
				minT = min(minT, abs((r - 1)->y - y) * 5LL);
		}
	}

	if (minT == INF_LL)
		cout << "SAFE";
	else
		cout << minT << endl;

	return 0;
}