結果

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

ソースコード

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

#define A first
#define B second

typedef long long Integer;
Integer gcd(Integer a, Integer b) { return a > 0 ? gcd(b % a, a) : b; }
struct rational {
  Integer p, q;
  void normalize() { // keep q positive
    if (q < 0) p *= -1, q *= -1;
    Integer d = gcd(p < 0 ? -p : p, q);
    if (d == 0) p = 0,  q = 1;
    else        p /= d, q /= d;
  }
  rational(Integer p, Integer q = 1) : p(p), q(q) {
    normalize();
  }
  rational &operator += (const rational &a) {
    p = a.q * p + a.p * q; q = a.q * q; normalize();
    return *this;
  }
  rational &operator -= (const rational &a) {
    p = a.q * p - a.p * q; q = a.q * q; normalize();
    return *this;
  }
  rational &operator *= (const rational &a) {
    p *= a.p; q *= a.q; normalize();
    return *this;
  }
  rational &operator /= (const rational &a) {
    p *= a.q; q *= a.p; normalize();
    return *this;
  }
  rational &operator - () {
    p *= -1;
    return *this;
  }
};
rational operator + (const rational &a, const rational &b) {
  return rational(a) += b;
}
rational operator * (const rational &a, const rational &b) {
  return rational(a) *= b;
}
rational operator - (const rational &a, const rational &b) {
  return rational(a) -= b;
}
rational operator / (const rational &a, const rational &b) {
  return rational(a) /= b;
}
bool operator < (const rational &a, const rational &b) { // avoid overflow
  return (long double) a.p * b.q < (long double) a.q * b.p;
}
bool operator <= (const rational &a, const rational &b) {
  return !(b < a);
}
bool operator > (const rational &a, const rational &b) {
  return b < a;
}
bool operator >= (const rational &a, const rational &b) {
  return !(a < b);
}
bool operator == (const rational &a, const rational &b) {
  return !(a < b) && !(b < a);
}
bool operator != (const rational &a, const rational &b) {
  return (a < b) || (b < a);
}

inline int cmp(int a, int b, int c, int d) {
  ll ad = a * d;
  ll cb = c * b;
  ll ret = ad - cb;
  if (ret > 0) return 1;
  if (ret < 0) return -1;
  return 0;
}

inline bool cmp_b(const pii &lhs, decltype(lhs) rhs) {
  return cmp(lhs.A, lhs.B, rhs.A, rhs.B) <= 0;
}

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

  auto check = [&](int s) -> bool {
    ll c = S[s].A;
    ll d = S[s].B;
    return a*(c+d)+c*(a+b) > (a+b)*(c+d);
  };

  while (abs(ok - ng) > 1) {
    int mid = (ok + ng) / 2;
    bool c = check(mid);
    // 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>();

  vpii R[3];

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

    R[t].emplace_back(a, b);
  }

  rep(i, 3) sort(all(R[i]), cmp_b);
  //rep(ii,3) { cout << "R["<<ii<<"]:"; rep(jj, sz(R[ii])) cout << ' ' << R[ii][jj]; cout << endl; }

  ll ans = 1 + n;

  forr(ab, R[1]) {
    int a = ab.A;
    int b = ab.B;
    int c0 = bs(a, b, R[0]);

    //_p("R1: a:%d, b:%d, c0: %d\n", a, b, c0);

    ans += c0;
  }

  forr(ab, R[2]) {
    ll a = ab.A;
    ll b = ab.B;
    int c0 = bs(a, b, R[0]);
    int c1 = bs(a, b, R[1]);

    //_p("R2: a:%lld, b:%lld, c0: %d, c1: %d\n", a, b, c0, c1);

    ans += c0;
    ans += c1;
  }



  // R[0] と R[1] O(N^2)
  vd X;
  vector<rational> Y;
  {
    forr(ab, R[0]) {
      ll a = ab.A;
      ll b = ab.B;
      forr(cd, R[1]) {
        ll c = cd.A;
        ll d = cd.B;
        if (a*(c+d) + c*(a+b) >= (a+b)*(c+d)) {
          if (cmp(a, a+b, c, c+d) <= 0) {
            X.emplace_back(1.0 * a / (a + b));
          }
          else {
            X.emplace_back(1.0 * c / (c + d));
          }
          {
            ll e = (a+b)*d + (c+d)*b;
            ll f = (a+b)*c - (c+d)*b;
            if (f > 0) {
              rational ef(e);
              ef /= f;
              Y.emplace_back(ef);
            }
          }
        }
      }
    }
  }

  // O(N^2 log N^2) TLE?
  sort(all(X));
  sort(all(Y));
  //cout << "X:"; rep(ii,sz(X)) cout << ' ' << X[ii]; cout << endl;

  // 3本
  // O(n * log(N^2))
  forr(ab, R[2]) {
    ll a = ab.A;
    ll b = ab.B;
    double bab = 1.0 * b / (a + b);
    // b/(a+b) 以上の X を数える
    // X は単調増加
    // OK な下限を探す
    int ok = sz(X);
    int ng = -1;

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

    // 交点に重なる場合を除く
    rational a_b(a);
    a_b /= b;
    ans -= upper_bound(all(Y), a_b) - lower_bound(all(Y), a_b);
  }

  cout << ans << endl;

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