結果

問題 No.461 三角形はいくつ？
コンテスト
ユーザー FF256grhy
提出日時 2016-12-12 02:57:59
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 88 ms / 5,000 ms
コード長 1,281 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,097 ms
コンパイル使用メモリ 45,540 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-11-29 13:38:35
合計ジャッジ時間 3,292 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 41
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:27:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   27 |         scanf("%d", &n);
      |         ~~~~~^~~~~~~~~~
main.cpp:28:26: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   28 |         inc(i, n) { scanf("%d%d%d", &p[i], &a[i], &b[i]); }
      |                     ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

ソースコード

diff #
raw source code 

#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;
	}
};

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] ].push_back(Frac(a[i], a[i] + b[i])); }
	inc(i, 3) { vec[   i  ].push_back(Frac(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(e1.q - e1.p, e1.q);
		Frac y(e2.q - e2.p, e2.q);
		Frac z((e1.q - e1.p) * e2.q + (e2.q - e2.p) * e1.q, e1.q * e2.q);
		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