結果

問題 No.461 三角形はいくつ？
ユーザー fumofumofuni
提出日時 2022-08-30 00:04:33
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,362 ms / 5,000 ms
コード長 2,334 bytes
コンパイル時間 3,278 ms
コンパイル使用メモリ 209,084 KB
最終ジャッジ日時 2025-02-06 23:45:08
ジャッジサーバーID
（参考情報） 		judge5 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include<bits/stdc++.h>
using namespace std;
#pragma GCC optimize("O3")
#define rep(i,n) for(ll i=0;i<n;i++)
#define repl(i,l,r) for(ll i=(l);i<(r);i++)
#define per(i,n) for(ll i=(n)-1;i>=0;i--)
#define perl(i,r,l) for(ll i=r-1;i>=l;i--)
#define fi first
#define se second
#define pb push_back
#define ins insert
#define pqueue(x) priority_queue<x,vector<x>,greater<x>>
#define all(x) (x).begin(),(x).end()
#define CST(x) cout<<fixed<<setprecision(x)
#define vtpl(x,y,z) vector<tuple<x,y,z>>
#define rev(x) reverse(x);
using ll=long long;
using vl=vector<ll>;
using vvl=vector<vector<ll>>;
using pl=pair<ll,ll>;
using vpl=vector<pl>;
using vvpl=vector<vpl>;
const ll MOD=1000000007;
const ll MOD9=998244353;
const int inf=1e9+10;
const ll INF=4e18;
const ll dy[9]={0,1,0,-1,1,1,-1,-1,0};
const ll dx[9]={1,0,-1,0,1,-1,1,-1,0};
template<class T> inline bool chmin(T& a, T b) {
    if (a > b) {
        a = b;
        return true;
    }
    return false;
}
template<class T> inline bool chmax(T& a, T b) {
    if (a < b) {
        a = b;
        return true;
    }
    return false;
}
pl add(pl a,pl b){
  ll f=a.first*b.second+a.second*b.first,g=a.second*b.second;
  ll k=gcd(f,g);
  return make_pair(f/k,g/k);
}
int main(){
  ll n;cin >> n;
  vvpl lin(3);
  rep(i,n){
    ll k,a,b;cin >> k >> a >> b;
    ll g=gcd(a,b);
    lin[k].emplace_back(make_pair(a/g,(a+b)/g));
  }
  rep(i,3)lin[i].emplace_back(make_pair(1,1));
  rep(i,3){
    sort(all(lin[i]),[](pl a,pl b){return a.first*b.second<a.second*b.first;});
  }
  set<pl> st;for(auto p:lin[2])st.insert(p);
  ll ans=0;
  for(auto p:lin[0]){
    for(auto q:lin[1]){
      if(p.first*q.second+q.first*p.second<p.second*q.second)continue;
      ll ok=lin[2].size(),ng=-1;
      while(ok-ng>1){
        ll mid=(ok+ng)>>1;
        auto r=lin[2][mid];
        if(max(q.second*r.second*(p.second-p.first),p.second*r.second*(q.second-q.first))<=p.second*q.second*r.first)ok=mid;
        else ng=mid;
      }
      ans+=lin[2].size()-ok;
      auto x=add(p,q);
      x.first=x.second*2-x.first;
      ll g=gcd(x.first,x.second);
      x.first/=g;x.second/=g;
      if(st.count(x))ans--;
      /*for(auto r:lin[2]){
        if(p.second*q.second*r.second*2-p.first*q.second*r.second-q.first*p.second*r.second==r.first*p.second*q.second)ans--;
      }*/
    }
  }
  cout << ans << endl;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0