結果

問題 No.461 三角形はいくつ？
コンテスト
ユーザー anta
提出日時 2016-12-12 00:39:06
言語 C++14
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 46 ms / 5,000 ms
コード長 2,374 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,875 ms
コンパイル使用メモリ 174,416 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-11-29 04:16:14
合計ジャッジ時間 3,685 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include "bits/stdc++.h"
using namespace std;
#define rep(i,n) for(int (i)=0;(i)<(int)(n);++(i))
#define rer(i,l,u) for(int (i)=(int)(l);(i)<=(int)(u);++(i))
#define reu(i,l,u) for(int (i)=(int)(l);(i)<(int)(u);++(i))
static const int INF = 0x3f3f3f3f; static const long long INFL = 0x3f3f3f3f3f3f3f3fLL;
typedef vector<int> vi; typedef pair<int, int> pii; typedef vector<pair<int, int> > vpii; typedef long long ll;
template<typename T, typename U> static void amin(T &x, U y) { if(y < x) x = y; }
template<typename T, typename U> static void amax(T &x, U y) { if(x < y) x = y; }

struct Ratio {
	long long num, den;
	Ratio() : num(0), den(1) {}
	Ratio(long long x, long long y) : num(x), den(y) { if(den < 0 || (den == 0 && num < 0)) num = -num, den = -den; }
	bool operator<(const Ratio &that) const { return (ll)num * that.den < (ll)that.num * den; }
	bool operator<=(const Ratio &that) const { return (ll)num * that.den <= (ll)that.num * den; }
	bool operator==(const Ratio &that) const { return (ll)num * that.den == (ll)that.num * den; }
	Ratio operator+(const Ratio &that) const {
		return Ratio(num * that.den + that.num * den, den * that.den);
	}
	Ratio operator-(const Ratio &that) const {
		return Ratio(num * that.den - that.num * den, den * that.den);
	}
};

int main() {
	int N;
	while(~scanf("%d", &N)) {
		vector<Ratio> segments[3];
		rep(p, 3)
			segments[p].emplace_back(1, 1);
		rep(i, N) {
			int p; int a; int b;
			scanf("%d%d%d", &p, &a, &b);
			segments[p].emplace_back(a, a + b);
		}
		rep(p, 3)
			sort(segments[p].begin(), segments[p].end());
		auto ge1 = [](Ratio a, Ratio b) {
			return (ll)a.num * b.den + (ll)b.num * a.den >= (ll)a.den * b.den;
		};
		ll ans = 0;
		for(auto s : segments[0]) {
			int n = (int)segments[2].size();
			int jR = n;
			int jL = n;
			int j1 = n, j2 = n;
			Ratio val1 = Ratio(1, 1) - s;
			for(; j1 > 0 && val1 <= segments[2][j1 - 1]; -- j1);
			for(auto t : segments[1]) {
				Ratio val = Ratio(2, 1) - (s + t);
				for(; jR > 0 && val < segments[2][jR - 1]; -- jR);
				ans += (int)segments[2].size() - jR;

				Ratio val2 = Ratio(1, 1) - t;
				for(; jL > 0 && val <= segments[2][jL - 1]; -- jL);
				for(; j2 > 0 && val2 <= segments[2][j2 - 1]; -- j2);
				if(Ratio(1, 1) <= s + t) {
					int L = max(j1, j2), R = jL;
					if(L < R)
						ans += R - L;
				}
			}
		}
		printf("%lld\n", ans);
	}
	return 0;
}
0