結果

問題 No.947 ABC包囲網
ユーザー face4
提出日時 2019-12-10 22:11:54
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,744 bytes
コンパイル時間 875 ms
コンパイル使用メモリ 88,272 KB
実行使用メモリ 10,524 KB
最終ジャッジ日時 2024-06-24 02:26:04
合計ジャッジ時間 7,860 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 12 WA * 16 TLE * 1 -- * 31
権限があれば一括ダウンロードができます

ソースコード

diff # 

// ?
#include<iostream>
#include<vector>
#include<algorithm>
#include<cmath>
#include<iomanip>
using namespace std;

#define EPS 0.00000001
#define equals(a, b) (fabs(a-b) < EPS)

const long double PI = 2*acos(0);

long double arg(long double x, long double y){
    return (atan2(y, x)>0 ? atan2(y,x) : 2*PI+atan2(y,x)) / PI * 180.0;
}

typedef long long ll;

int main(){
    int n;
    cin >> n;
    vector<long double> v, x(n), y(n);
    for(int i = 0; i < n; i++){
        cin >> x[i] >> y[i];
        v.push_back(arg(x[i], y[i]));
    }
    sort(v.begin(), v.end());
    ll ans = 0;
    for(int i = 0; i < n; i++){
        for(int j = i+1; j < n; j++){
            if(equals(v[i],v[j]) || equals(fabs(v[i]-v[j]),180))    continue;
            long double f = v[i]+180, g = v[j]+180;
            if(f < 360 && g < 360){
                ans += lower_bound(v.begin(),v.end(),g)-upper_bound(v.begin(),v.end(),f);
            }else if(f < 360 && g >= 360){
                long double diff = min(v[j]-v[i], 360-(v[j]-v[i]));
                g -= 360;
                if(equals(g+360-f, diff)){
                    ans += v.end() - upper_bound(v.begin(),v.end(),f);
                    ans += lower_bound(v.begin(),v.end(),g) - v.begin();
                }else{
                    swap(f, g);
                    ans += lower_bound(v.begin(),v.end(),g)-upper_bound(v.begin(),v.end(),f);
                }
            }else if(f >= 360 && g >= 360){
                f -= 360;
                g -= 360;
                ans += lower_bound(v.begin(),v.end(),g)-upper_bound(v.begin(),v.end(),f);
            }
            cerr << i << " " << j << " " << ans << endl;
        }
    }
    cerr << ans << endl;
    cout << ans/3 << endl;
    return 0;
}
0