結果

問題 No.461 三角形はいくつ？
ユーザー koba-e964koba-e964
提出日時 2016-12-12 12:10:30
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 2,317 ms / 5,000 ms
コード長 1,932 bytes
コンパイル時間 736 ms
コンパイル使用メモリ 70,228 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-29 15:41:57
合計ジャッジ時間 35,100 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 1 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 6 ms
5,248 KB
testcase_05 AC 1,438 ms
5,248 KB
testcase_06 AC 835 ms
5,248 KB
testcase_07 AC 971 ms
5,248 KB
testcase_08 AC 659 ms
5,248 KB
testcase_09 AC 4 ms
5,248 KB
testcase_10 AC 205 ms
5,248 KB
testcase_11 AC 1,098 ms
5,248 KB
testcase_12 AC 487 ms
5,248 KB
testcase_13 AC 1,654 ms
5,248 KB
testcase_14 AC 1,589 ms
5,248 KB
testcase_15 AC 1,661 ms
5,248 KB
testcase_16 AC 648 ms
5,248 KB
testcase_17 AC 638 ms
5,248 KB
testcase_18 AC 2,311 ms
5,248 KB
testcase_19 AC 2,317 ms
5,248 KB
testcase_20 AC 1,659 ms
5,248 KB
testcase_21 AC 1,736 ms
5,248 KB
testcase_22 AC 1,700 ms
5,248 KB
testcase_23 AC 554 ms
5,248 KB
testcase_24 AC 533 ms
5,248 KB
testcase_25 AC 7 ms
5,248 KB
testcase_26 AC 7 ms
5,248 KB
testcase_27 AC 7 ms
5,248 KB
testcase_28 AC 7 ms
5,248 KB
testcase_29 AC 491 ms
5,248 KB
testcase_30 AC 228 ms
5,248 KB
testcase_31 AC 941 ms
5,248 KB
testcase_32 AC 1,250 ms
5,248 KB
testcase_33 AC 6 ms
5,248 KB
testcase_34 AC 5 ms
5,248 KB
testcase_35 AC 6 ms
5,248 KB
testcase_36 AC 666 ms
5,248 KB
testcase_37 AC 660 ms
5,248 KB
testcase_38 AC 673 ms
5,248 KB
testcase_39 AC 621 ms
5,248 KB
testcase_40 AC 665 ms
5,248 KB
testcase_41 AC 683 ms
5,248 KB
testcase_42 AC 1,175 ms
5,248 KB
testcase_43 AC 1,168 ms
5,248 KB
testcase_44 AC 1,203 ms
5,248 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <algorithm>
#include <cassert>
#include <iostream>
#include <vector>

#define REP(i,s,n) for(int i=(int)(s);i<(int)(n);i++)

using namespace std;
typedef long long int ll;
typedef vector<ll> VL;
typedef pair<ll, ll> PL;

ll gcd(ll x, ll y) {
  return y == 0 ? x : gcd(y, x % y);
}

PL add(PL x, PL y) {
  PL w = PL(x.first * y.second + x.second * y.first,
		x.second * y.second);
  ll g = gcd(w.first, w.second);
  w.first /= g;
  w.second /= g;
  return w;
}

PL neg(PL x) { return PL(- x.first, x.second); }

bool fr_le(PL x, PL y) {
  return x.first * y.second <= y.first * x.second;
}
bool fr_lt(PL x, PL y) {
  return x.first * y.second < y.first * x.second;
}
int main(void){
  int n;
  cin >> n;
  vector<PL> pool[3];
  REP(i, 0, 3) {
    pool[i].push_back(PL(0, 1));
  }
  REP(i, 0, n) {
    int p, a, b;
    cin >> p >> a >> b;
    int g = __gcd(a, b);
    a /= g;
    b /= g;
    pool[p].push_back(PL(b, a + b));
  }
  if (pool[0].size() > pool[2].size()) {
    swap(pool[0], pool[2]);
  }
  if (pool[1].size() > pool[2].size()) {
    swap(pool[1], pool[2]);
  }
  ll tot = 0;
  sort(pool[2].begin(), pool[2].end(), fr_lt);
  REP(i, 0, pool[0].size()) {
    PL u = pool[0][i];
    REP(j, 0, pool[1].size()) {
      PL v = pool[1][j];
      if (not fr_le(add(u, v), PL(1, 1))) {
	continue;
      }
      PL ma = u;
      if (fr_le(u, v)) { ma = v; }
      // find #{x | x + u <= 1 && x + v <= 1 && x + u + v != 1}
      PL x_ma = add(PL(1, 1), neg(ma));
      int idx = upper_bound(pool[2].begin(), pool[2].end(), x_ma, fr_lt)
	- pool[2].begin();
      tot += idx;
      PL excl = add(PL(1, 1), neg(add(u, v)));
      {
	int idx1 = upper_bound(pool[2].begin(), pool[2].end(), excl, fr_lt)
	- pool[2].begin();
	int idx2 = lower_bound(pool[2].begin(), pool[2].end(), excl, fr_lt)
	- pool[2].begin();
	assert (idx1 - idx2 <= 1);
	assert (idx1 >= idx2);
	tot -= idx1 - idx2;
      }
    }
  }
  cout << tot << endl;
}
0