結果
問題 | No.461 三角形はいくつ? |
ユーザー |
|
提出日時 | 2016-12-11 18:31:23 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,799 bytes |
コンパイル時間 | 845 ms |
コンパイル使用メモリ | 79,776 KB |
実行使用メモリ | 10,280 KB |
最終ジャッジ日時 | 2024-11-29 03:49:35 |
合計ジャッジ時間 | 147,626 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 31 TLE * 10 |
ソースコード
#include <iostream>#include <vector>#include <algorithm>#include <array>#include <cassert>#define repeat(i,n) for (int i = 0; (i) < int(n); ++(i))#define whole(f,x,...) ([&](decltype((x)) whole) { return (f)(begin(whole), end(whole), ## __VA_ARGS__); })(x)typedef long long ll;using namespace std;struct rational { ll p, q; };rational make_rational(ll p) { return { p, 1 }; }rational make_rational(ll p, ll q) { if (p == 0) return { 0, 1 }; ll gcd = __gcd(p, q); return { p/gcd, q/gcd }; }rational operator + (rational a, rational b) { return make_rational(a.p * b.q + b.p * a.q, a.q * b.q); }rational operator - (ll a, rational b) { return make_rational(a * b.q - b.p, b.q); }bool operator < (rational a, rational b) { return a.p * b.q < a.q * b.p; }bool operator == (rational a, rational b) { return a.p * b.q == a.q * b.p; }bool operator <= (rational a, ll b) { return a.p <= a.q * b; }int main() {int n; cin >> n;array<vector<rational>, 3> qs = {};repeat (i,n) {int p, a, b; cin >> p >> a >> b;qs[p].push_back(make_rational(a, a + b));}repeat (i,3) {qs[i].emplace_back(make_rational(1));whole(sort, qs[i]);}ll cnt = 0;repeat (i,3) {for (rational q1 : qs[(i + 1) % 3]) {for (rational q2 : qs[(i + 2) % 3]) {rational qm = max(1 - q1, 1 - q2);rational qc = (1 - q1) + (1 - q2);if (qc <= 1) {int lm = whole(lower_bound, qs[i], qm) - qs[i].begin();int lc = whole(lower_bound, qs[i], qc) - qs[i].begin();cnt += qs[i].size() - lm - (qs[i][lc] == qc);}}}}assert (cnt % 3 == 0);cout << cnt / 3 << endl;return 0;}