【jzoj4894】【SJR的直线】【数论】

题目大意

给出没有三线交在一点的一些直线，求围成了多少个三角形。

解题思路

求出交点的个数组合数一发，再求出三点共线的情况，减一发即可。

code

#include<cmath>
#include<cstdio>
#include<algorithm>
#define LL long long
#define fo(i,j,k) for(int i=j;i<=k;i++)
#define fd(i,j,k) for(int i=j;i>=k;i--)
using namespace std;
int const maxn=3*1e5;
LL n,cnt[maxn+10],mod=1e9+7;
struct rec{
    LL a,b,c;
};
rec a[maxn+10];
bool cmp(rec x,rec y){
    return -1.0*x.a/x.b+1e-10<-1.0*y.a/y.b;
}
int main(){
    //freopen("trokuti.in","r",stdin);
    //freopen("trokuti.out","w",stdout);
    freopen("d.in","r",stdin);
    freopen("d.out","w",stdout);
    scanf("%lld",&n);int tmp=0;cnt[0]=2;
    fo(i,1,n){
        LL aa,bb,cc;
        scanf("%lld%lld%lld",&aa,&bb,&cc);
        if(!aa)cnt[1]++;
        else if(!bb)cnt[2]++;
        else{
            tmp++;
            a[tmp].a=aa;
            a[tmp].b=bb;
            a[tmp].c=cc;
        }
    }
    sort(a+1,a+tmp+1,cmp);
    cnt[++cnt[0]]=1;LL sum=0;
    fo(i,2,tmp)
        if(-a[i-1].a*a[i].b==-a[i].a*a[i-1].b)cnt[cnt[0]]++;
        else{
            sum+=cnt[cnt[0]]*cnt[cnt[0]];
            cnt[++cnt[0]]=1;
        }
    if(cnt[0]>2)sum+=cnt[cnt[0]]*cnt[cnt[0]];
    sum+=cnt[1]*cnt[1];
    sum+=cnt[2]*cnt[2];
    LL ans=0;
    fo(i,1,cnt[0]){
        ans+=cnt[i]*((n*n-sum)/2-(n-cnt[i])*cnt[i]);
    }
    printf("%lld",ans/3%mod);
    return 0;
}