結果

問題 No.461 三角形はいくつ？
コンテスト
ユーザー しらっ亭
提出日時 2016-12-12 05:33:10
言語 C++14
(gcc 13.3.0 + boost 1.89.0)
結果
WA  
実行時間 -
コード長 2,706 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,824 ms
コンパイル使用メモリ 175,700 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-11-29 15:08:46
合計ジャッジ時間 5,796 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 28 WA * 13
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;

#define _p(...) (void)printf(__VA_ARGS__)
#define forr(x,arr) for(auto&& x:arr)
#define _overload3(_1,_2,_3,name,...) name
#define _rep2(i,n) _rep3(i,0,n)
#define _rep3(i,a,b) for(int i=int(a);i<int(b);++i)
#define rep(...) _overload3(__VA_ARGS__,_rep3,_rep2,)(__VA_ARGS__)
#define _rrep2(i,n) _rrep3(i,0,n)
#define _rrep3(i,a,b) for(int i=int(b)-1;i>=int(a);i--)
#define rrep(...) _overload3(__VA_ARGS__,_rrep3,_rrep2,)(__VA_ARGS__)
#define all(x) (x).begin(), (x).end()
#define bit(n) (1LL<<(n))
#define sz(x) ((int)(x).size())
#define fst first
#define snd second
using ll=long long;using pii=pair<int,int>;
using vb=vector<bool>;using vs=vector<string>;
using vi=vector<int>;using vvi=vector<vi>;using vvvi=vector<vvi>;
using vl=vector<ll>;using vvl=vector<vl>;using vvvl=vector<vvl>;
using vd=vector<double>;using vvd=vector<vd>;using vvvd=vector<vvd>;
using vpii=vector<pii>;using vvpii=vector<vpii>;using vvvpii=vector<vvpii>;
template<typename T>T read(){T t;cin>>t;return t;}
template<class T,class U>ostream&operator<<(ostream&o,const pair<T,U>&p){o<<'('<<p.fst<<", "<<p.snd<<')';return o;}

// a/(a+b) と足したら1を越えるのが何本あるか
// O(log |S|)
int bs(const long double &ab, const vector<long double> &S) {
  int ok = sz(S);
  int ng = -1;

  while (abs(ok - ng) > 1) {
    int mid = (ok + ng) / 2;
    bool c = ab + S[mid] > 1;
    // cout << "ok: " << ok << ", ng: " << ng << ", ";
    // cout << "check(" << mid << ") -> " << c << endl;
    if (c) ok = mid;
    else ng = mid;
  }
  return sz(S) - ok;
}

void Main() {
  int n = read<int>();

  vector<vector<long double>> R(3);

  rep(i, n) {
    int t = read<int>();
    int a = read<int>();
    int b = read<int>();

    R[t].emplace_back(1.0 * a / (a + b));
  }

  {
    int s0 = sz(R[0]);
    int s1 = sz(R[1]);
    int s2 = sz(R[2]);
    if (s0 > s2) swap(R[0], R[2]);
    if (s1 > s2) swap(R[1], R[2]);
  }

  rep(i, 3) sort(all(R[i]));

  ll ans = 1 + n;

  forr(ab, R[1]) {
    int c0 = bs(ab, R[0]);
    ans += c0;
  }

  forr(ab, R[2]) {
    int c0 = bs(ab, R[0]);
    int c1 = bs(ab, R[1]);
    ans += c0;
    ans += c1;
  }

  // R[0] と R[1]
  forr(ab, R[0]) {
    forr(cd, R[1]) {
      if (ab + cd >= 1) {
        {
          long double z = 1 - (ab < cd ? ab : cd);
          ans += R[2].end() - lower_bound(all(R[2]), z);
        }
        {
          auto e = (1-cd) + (1-ab);
          auto f = cd - (1-ab);
          long double ef(e / (e+f));
          ans -= upper_bound(all(R[2]), ef) - lower_bound(all(R[2]), ef);
        }
      }
    }
  }

  cout << ans << endl;
}
int main() { cin.tie(nullptr); ios::sync_with_stdio(false); Main(); return 0; }
0