結果
問題 | No.461 三角形はいくつ? |
ユーザー | paruki |
提出日時 | 2016-12-13 15:57:27 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 314 ms / 5,000 ms |
コード長 | 1,827 bytes |
コンパイル時間 | 1,841 ms |
コンパイル使用メモリ | 178,536 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-30 01:05:39 |
合計ジャッジ時間 | 6,268 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 41 |
ソースコード
#include "bits/stdc++.h" using namespace std; #define FOR(i,j,k) for(int (i)=(j);(i)<(int)(k);++(i)) #define rep(i,j) FOR(i,0,j) #define each(x,y) for(auto &(x):(y)) #define mp make_pair #define mt make_tuple #define all(x) (x).begin(),(x).end() #define debug(x) cout<<#x<<": "<<(x)<<endl #define smax(x,y) (x)=max((x),(y)) #define smin(x,y) (x)=min((x),(y)) #define MEM(x,y) memset((x),(y),sizeof (x)) #define sz(x) (int)(x).size() typedef long long ll; typedef pair<int, int> pii; typedef vector<int> vi; typedef vector<ll> vll; typedef pair<ll, ll> pll; bool les(pll x, pll y) { return x.first*y.second < y.first*x.second; } bool f(pll x, pll y) { return x.first*y.second + x.second*y.first <= x.second*y.second; } bool g(pll x, pll y) { return x.first*y.second + x.second*y.first < x.second*y.second; } int main(){ ios::sync_with_stdio(false); cin.tie(0); int N; cin >> N; vi p(N), a(N), b(N); rep(i, N)cin >> p[i] >> a[i] >> b[i]; vector<vector<pll>> vv(3); rep(i, N) { vv[p[i]].push_back(mp(b[i], a[i] + b[i])); } rep(i, 3)sort(all(vv[i]), les); ll ans = N + 1; if(sz(vv[0])*sz(vv[1])*sz(vv[2])) each(s, vv[0]) { each(t, vv[1]) if(f(s, t) && f(s,vv[2][0]) && f(t, vv[2][0])){ int lb = 0, ub = sz(vv[2]), m; while(ub - lb > 1) { m = (lb + ub) >> 1; ((f(s, vv[2][m]) && f(t, vv[2][m])) ? lb : ub) = m; } ll d = s.second*t.second; pll ng(d - s.first*t.second - t.first*s.second, d); ans += ub; if(binary_search(vv[2].begin(), vv[2].begin()+ub, ng, les)) --ans; } } rep(i, 3) { int j = (i + 1) % 3; each(s, vv[i])each(t, vv[j])if(g(s, t))ans++; } cout << ans << endl; }