状压DP
#include<cstdio>
using namespace std;
const int mod=1e8+7;
int F[1000005][25],dis[25][25],lim[1000005];
struct node{
int x,y;
}E[25];
int main(){
int n;
scanf("%d",&n);
for (int i=0; i<n; i++) scanf("%d%d",&E[i].x,&E[i].y);
for (int S=0; S<n; S++)
for (int T=S+1; T<n; T++)
for (int k=0; k<n; k++){
if (k==S || k==T) continue;
int x=E[S].x-E[k].x,y=E[S].y-E[k].y;
int X=E[T].x-E[k].x,Y=E[T].y-E[k].y;
if (X*y==x*Y) {
if (E[k].x>E[S].x && E[k].x>E[T].x) continue;
if (E[k].x<E[S].x && E[k].x<E[T].x) continue;
if (E[k].y>E[S].y && E[k].y>E[T].y) continue;
if (E[k].y<E[S].y && E[k].y<E[T].y) continue;
dis[S][T]|=(1<<k);
dis[T][S]|=(1<<k);
}
}
for (int i=0; i<n; i++) F[1<<i][i]=1;
for (int i=1; i<(1<<n); i++) lim[i]=lim[i>>1]+(i&1);
for (int i=0; i<(1<<n); i++)
for (int S=0; S<n; S++)
if (F[i][S]){
for (int T=0; T<n; T++)
if (!(i&(1<<T)) && (i&dis[S][T])==dis[S][T]) (F[i|(1<<T)][T]+=F[i][S])%=mod;
}
int ans=0;
for (int i=0; i<(1<<n); i++)
for (int S=0; S<n; S++)
if (lim[i]>=4) (ans+=F[i][S])%=mod;
printf("%d\n",ans);
return 0;
}