結果
問題 | No.1144 Triangles |
ユーザー | leaf_1415 |
提出日時 | 2020-07-31 23:13:44 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 410 ms / 3,000 ms |
コード長 | 1,685 bytes |
コンパイル時間 | 781 ms |
コンパイル使用メモリ | 89,752 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-07-06 21:11:42 |
合計ジャッジ時間 | 6,443 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 25 |
ソースコード
#include <iostream> #include <cstdio> #include <cmath> #include <ctime> #include <cstdlib> #include <cassert> #include <vector> #include <list> #include <stack> #include <queue> #include <deque> #include <map> #include <set> #include <bitset> #include <string> #include <algorithm> #include <utility> #define llint long long #define inf 1e18 #define rep(x, s, t) for(llint (x) = (s); (x) < (t); (x)++) #define Rep(x, s, t) for(llint (x) = (s); (x) <= (t); (x)++) #define chmin(x, y) (x) = min((x), (y)) #define chmax(x, y) (x) = max((x), (y)) #define PI 3.141592653589793238462643383279502884197169393 #define eps 1e-12 #define mod 1000000007 using namespace std; typedef pair<double, llint> P; llint n; llint x[2005], y[2005]; vector<P> vec; llint get(llint x) { return (x%mod+mod)%mod; } int main(void) { ios::sync_with_stdio(0); cin.tie(0); cin >> n; for(int i = 1; i <= n; i++) cin >> x[i] >> y[i]; llint ans = 0; for(int i = 1; i <= n; i++){ vec.clear(); for(int j = 1; j <= n; j++){ if(i == j) continue; if(x[i] == x[j] && y[i] == y[j]) continue; vec.push_back(P(atan2((double)y[j]-y[i], (double)x[j]-x[i]), j)); } sort(vec.begin(), vec.end()); llint m = vec.size(); for(int j = 0; j < m; j++){ vec.push_back(vec[j]); vec.back().first += 2 * PI; } llint l = 0, r = 0; for(int j = 0; j < m; j++){ double ang = vec[j].first; llint id = vec[j].second; while(vec[l+1].first <= ang+eps) l++; while(vec[r+1].first < ang+PI-eps) r++; llint tmp = (get(x[i])*get(y[id])%mod + mod - get(x[id])*get(y[i])%mod) % mod; tmp *= max(0LL, r-l), tmp %= mod; ans += tmp, ans %= mod; } } cout << ans << endl; return 0; }