思路：

1.銳角三角形總銳角個數=總銳角數-非銳角三角形提供銳角數。則銳角三角形個數=總銳角數/3（即（銳角數-2*（直角+鈍角數））/3，每鈍角和直角三角形提供兩銳角）。

2.列舉每個點p[i]，以p[i]為原點，求其他n-1個點與原點組成的向量，按極角（小於0時加2*PI）遞增排序。設定三個指標l、r、equ，分別代表第一個直角、第一個平角、第一個非0角的位置。列舉n-1個向量j，由於極角遞增，兩向量角度即為極角l-極角j，l、r、equ遞增則形成的角度也遞增（相當於尺取法）。則以點p[i]為頂點的銳角個數為l-equ（所有小於90度的角減去0度的角），鈍角個數為r-l（所有小於180度的角減去小於等於90度的角）。

3.當指標l、r、equ旋轉一週時（l、r、equ下標超過n-1，但此時所有角度尚未列舉完畢），為了便於求角度（即最後一部分向量與起始部分向量所成角度），列舉之前將所有極角加2*PI複製一遍。

#include<cmath>
#include<cstdio>
#include<vector>
#include<cstring>
#include<iostream>
#include<algorithm>
#define debug
using namespace std;

const double eps=1e-11;
const double PI=acos(-1.0);
const int maxn=4000+50;

int n;

struct Point
{
    double x,y;
    Point() {}
    Point(double x,double y):x(x),y(y) {}
    void read()
    {
        scanf("%lf%lf",&x,&y);
    }
};

typedef Point Vector;

Vector operator - (Vector A,Vector B)
{
    return Vector(A.x-B.x,A.y-B.y);
}

int dcmp(double x)
{
    if(fabs(x)<eps) return 0;
    else return x<0?-1:1;
}

int cnt;
Point p[maxn];
double a[maxn];

int main()
{
#ifdef debu
    freopen("in.txt","r",stdin);
#endif // debug
    while(scanf("%d",&n)!=EOF)
    {
        int acute=0,notacute=0;
        for(int i=1; i<=n; i++) p[i].read();
        for(int i=1; i<=n; i++)
        {
            cnt=0;
            for(int j=1; j<=n; j++)
            {
                if(i==j)continue;
                Vector tmp=p[j]-p[i];
                double ang=atan2(tmp.y,tmp.x);
                if(ang<0) ang+=2*PI;
                a[++cnt]=ang;
            }

            sort(a+1,a+cnt+1);

            for(int j=1;j<=cnt;j++) a[cnt+j]=a[j]+2*PI;

            int l=1,r=1,equ=1;

            for(int j=1; j<=cnt; j++)
            {
                while(l<=2*cnt&&dcmp(a[l]-a[j]-PI/2)<0) l++;
                while(r<=2*cnt&&dcmp(a[r]-a[j]-PI)<0) r++;
                while(equ<=2*cnt&&dcmp(a[equ]-a[j])==0) equ++;


                acute+=l-equ;
                notacute+=r-l;
            }
        }

        printf("%d\n",(acute-2*notacute)/3);
    }
    return 0;
}