結果

問題 No.461 三角形はいくつ？
ユーザー kimiyuki
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

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