結果

問題 No.202 1円玉投げ
コンテスト
ユーザー Yang33
提出日時 2018-08-02 11:29:02
言語 C++14
(gcc 13.3.0 + boost 1.89.0)
結果
WA  
実行時間 -
コード長 4,811 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,015 ms
コンパイル使用メモリ 174,564 KB
実行使用メモリ 14,976 KB
最終ジャッジ日時 2024-12-22 11:30:31
合計ジャッジ時間 4,206 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 24 WA * 14
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;

using VS = vector<string>;    using LL = long long;
using VI = vector<int>;       using VVI = vector<VI>;
using PII = pair<int, int>;   using PLL = pair<LL, LL>;
using VL = vector<LL>;        using VVL = vector<VL>;

#define ALL(a)  begin((a)),end((a))
#define RALL(a) (a).rbegin(), (a).rend()
#define SZ(a) int((a).size())
#define SORT(c) sort(ALL((c)))
#define RSORT(c) sort(RALL((c)))
#define UNIQ(c) (c).erase(unique(ALL((c))), end((c)))
#define FOR(i, s, e) for (int(i) = (s); (i) < (e); (i)++)
#define FORR(i, s, e) for (int(i) = (s); (i) > (e); (i)--)
#define debug(x) cerr << #x << ": " << x << endl
const int INF = 1e9;                          const LL LINF = 1e16;
const LL MOD = 1000000007;                    const double PI = acos(-1.0);
int DX[9] = { 0,0, 0, 1, -1, 1, 1, -1, -1 };    int DY[9] = { 0,1, -1, 0, 0, 1, -1, 1, -1 };

/* -----  2018/08/02  Problem: yukicoder 202  / Link: http://yukicoder.me/problems/no/202  ----- */
/* ------問題------

ユキさんは、手元にある大量の1円玉を1枚ずつ部屋の床に投げて遊んでいます。

投げた1円玉が床の上で止まったとき、それが既に投げられた床の上の1円玉と
少しでも重なって止まった場合には、次の1円玉を投げる前に、その投げた1円玉を
取り除きます。
そうではない場合、すなわち、投げた1円玉が床の上の1円玉と全く重ならずに
止まった場合には、その投げた1円玉をそのまま床の上に残します。
（この場合、その投げた1円玉はそれ以降ずっと移動しないものとします）。

ユキさんは、すでに1円玉をN回投げ終えました。さて、床の上には何枚の
1円玉があるでしょうか？

なお、
投げた1円玉は、床の上（または床の上の1円玉の上）に、必ず倒れた状態で
止まるものとします。
1円玉は円形とし、半径は10とします。厚みは無視します。
2つの1円玉の円がちょうど互いに外接する状態では、それら2つの1円玉は
重なっていないものとします。
N回投げたそれぞれの1円玉が止まった時の円の中心位置が、
床面上の座標として記録されています。

-----問題ここまで----- */
/* -----解説等-----



----解説ここまで---- */
using Double = long double;
const double EPS = 1e-9;

typedef struct Point {
	Double x, y;
	Point(Double x_, Double y_) : x(x_), y(y_) {}

	bool operator < (const Point& a) const {
		return fabs(x - a.x) < EPS ? y + EPS < a.y : x + EPS < a.x;
	}
	bool operator > (const Point& a) const {
		return a < *this;
	}
	bool operator == (const Point& a) const {
		return !(a < *this) && !(a > *this);
	}
	Point operator + (const Point& a) const {
		return Point(x + a.x, y + a.y);
	}
	Point operator - (const Point& a) const {
		return Point(x - a.x, y - a.y);
	}
	Point operator * (const Point& a) const {
		return Point(x * a.x, y * a.y);
	}
	Point operator / (const Point& a) const {
		return Point(x / a.x, y / a.y);
	}
	Point operator * (const Double& d) const {
		return Point(x * d, y * d);
	}
	Point operator / (const Double& d) const {
		return Point(x / d, y / d);
	}
} Vector;

// ノルム
Double norm(const Point& a) {
	return a.x * a.x + a.y * a.y;
}
// サイズ
Double abs(const Point& a) {
	return sqrt(norm(a));
}

// 円構造体
struct Circle {
	Point c; Double r;
	Circle(Point c_, Double r_) : c(c_), r(r_) {}

	Double area() const { return PI * r * r; }
};

// 円の衝突
int circle_intersect(const Circle& a, const Circle& b) {
	Double d = abs(a.c - b.c);
	if (a.r + b.r + EPS < d) return 0; // 0 : outside
	if (b.r + d + EPS < a.r) return -3; // 0 : b in a
	if (a.r + d + EPS < b.r) return -4; // 0 : a in b
	if (fabs(a.r + b.r - d) < EPS) return 1; // 1 : outside
	if (fabs(b.r + d - a.r) < EPS) return 3; // 1 : b in a
	if (fabs(a.r + d - b.r) < EPS) return 4; // 1 : a in b
	return 2; // 2
}

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

	int N; cin >> N;
	// 50*50
	vector<vector<vector<Circle>>>block(410, vector<vector<Circle>>(410));
	auto check = [&block](const Circle& c) {
		int X = c.c.x / 50 + 2, Y = c.c.y / 50 + 2;

		FOR(k, 0, 9) {
			int nx = X + DX[k], ny = Y + DY[k];
			for (const auto& cc : block[nx][ny]) {
				int ci = circle_intersect(cc, c);
				if (ci==2) {
					return false;
				}
			}
		}

		return true;
	};
	auto PUSH = [&block](const Circle& c) {
		int X = c.c.x / 50 + 2, Y = c.c.y / 50 + 2;

		block[X][Y].push_back(c);
	};
	FOR(i, 0, N) {
		int x, y; cin >> x >> y;
		Circle c(Point(x, y), 10);
		if (check(c)) {
			PUSH(c);
		}
	}
	int ans = 0;
	FOR(i, 0, 410) {
		FOR(j, 0, 410) {
			ans += SZ(block[i][j]);
		}
	}
	cout << ans << "\n";

	return 0;
}
0