#include #pragma GCC optimize("Ofast") #define _GLIBCXX_DEBUG using namespace std; using std::cout; using std::cin; using std::endl; using ll=long long; using ld=long double; ll ILL=1167167167167167167; const int INF=2100000000; const ll mod=1e9+7; #define rep(i,a) for (ll i=0;i using _pq = priority_queue, greater>; template ll LB(vector &v,T a){return lower_bound(v.begin(),v.end(),a)-v.begin();} template ll UB(vector &v,T a){return upper_bound(v.begin(),v.end(),a)-v.begin();} template bool chmin(T &a,const T &b){if(a>b){a=b;return 1;}else return 0;} template bool chmax(T &a,const T &b){if(a void So(vector &v) {sort(v.begin(),v.end());} template void Sore(vector &v) {sort(v.begin(),v.end(),[](T x,T y){return x>y;});} void yneos(bool a){if(a) cout<<"Yes\n"; else cout<<"No\n";} template void vec_out(vector &p){for(int i=0;i<(int)(p.size());i++){if(i) cout<<" ";cout< l,pair r){ if(l.first==0&&r.first==0) return l.second>r.second; if(l.first==0){ return l.second>0||r.first<0; } if(r.first==0){ return r.second>0||l.first<0; } if((l.first<0)^(r.first<0)) return l.first>0; return l.first*r.second>t; rep(i,t) solve(); } void solve(){ int N; cin>>N; vector X(N),Y(N);ll ans=0; rep(i,N) cin>>X[i]>>Y[i]; rep(roop,N){ vector> point; rep(i,N){ if(X[i]==X[roop]&&Y[i]==Y[roop]) continue; point.push_back({X[i]-X[roop],Y[i]-Y[roop]}); } int M=point.size(); vector order(M*2); rep(i,M*2) order[i]=i; sort(all(order),[&](int L,int R){ auto l=point[L%M],r=point[R%M]; if(L>=M) l.first*=-1,l.second*=-1; if(R>=M) r.first*=-1,r.second*=-1; if(l==r) return Lr.second; } if(l.first==0){ return (l.second>0||r.first<0); } if(r.first==0){ return !(r.second>0||l.first<0); } if((l.first<0)^(r.first<0)){ return l.first>0; } return l.first*r.second0||(point[i].first==0&&point[i].second>0)){ x_sum-=point[i].first; y_sum-=point[i].second; }else{ x_sum+=point[i].first; y_sum+=point[i].second; } } rep(i,M*2){ int a=order[i]; if(a>=M){ a-=M; x_sum-=point[a].first*2ll; y_sum-=point[a].second*2ll; continue; } x_sum+=point[a].first*2ll; y_sum+=point[a].second*2ll; x_sum%=mod; y_sum%=mod; ans=(ans+(y_sum*point[a].first)%mod)%mod; ans=(ans-(x_sum*point[a].second)%mod)%mod; //cout<