結果

問題 No.461 三角形はいくつ？
ユーザー FF256grhy
提出日時 2016-12-14 09:59:14
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 97 ms / 5,000 ms
コード長 1,405 bytes
コンパイル時間 685 ms
コンパイル使用メモリ 46,720 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-11-30 07:04:59
合計ジャッジ時間 3,737 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 41
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:33:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   33 |         scanf("%d", &n);
      |         ~~~~~^~~~~~~~~~
main.cpp:34:26: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   34 |         inc(i, n) { scanf("%d%d%d", &p[i], &a[i], &b[i]); }
      |                     ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

ソースコード

diff # 

#include <cstdio>

#include <vector>
#include <utility>
#include <algorithm>

typedef long long   signed int LL;
#define incID(i, l, r) for(int i = (l)    ; i <  (r); i++)
#define inc( i, n) incID(i, 0, n)

struct Frac {
	LL p, q;
	Frac(LL pp, LL qq) { p = pp; q = qq; }
	bool operator< (const Frac & obj) const {
		return p * obj.q <  obj.p * q;
	}
	bool operator==(const Frac & obj) const {
		return p * obj.q == obj.p * q;
	}
	Frac operator+ (const Frac & obj) const {
		return Frac(p * obj.q + obj.p * q, q * obj.q);
	}
	Frac operator- (const Frac & obj) const {
		return Frac(p * obj.q - obj.p * q, q * obj.q);
	}
};

int n, p[4000], a[4000], b[4000];
LL ans;
std::vector<Frac> vec[3];

int main() {
	scanf("%d", &n);
	inc(i, n) { scanf("%d%d%d", &p[i], &a[i], &b[i]); }
	
	inc(i, n) { vec[ p[i] ].emplace_back(a[i], a[i] + b[i]); }
	inc(i, 3) { vec[   i  ].emplace_back(1, 1); }
	
	inc(i, 3) { std::sort(vec[i].begin(), vec[i].end()); }
	
	for(auto && e1: vec[1]) {
	for(auto && e2: vec[2]) {
		Frac x = Frac(1, 1) - e1;
		Frac y = Frac(1, 1) - e2;
		Frac z = x + y;
		if(Frac(1, 1) < z) { continue; }
		if(y < x) { std::swap(x, y); }
		
		auto lby = std::lower_bound(vec[0].begin(), vec[0].end(), y);
		auto lbz = std::lower_bound(vec[0].begin(), vec[0].end(), z);
		int v = vec[0].end() - lby, w = (lbz != vec[0].end() && *lbz == z);
		ans += v - w;
	}
	}
	
	printf("%lld\n", ans);
	
	return 0;
}
0