SJR的直线

题目大意

给定n条直线，保证任意三条不共点，求形成了多少个三角形。

组合数算算

至多能形成$C_n^3$个三角形。
我们看看什么时候会不合法。
1、两条平行第三条不与它们平行
2、三条都平行
一组平行线假如有t条，那么非法贡献是$C_t^2*(n-t)+C_t^3$
用map统计每种斜率的直线个数然后按上面算即可。

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<map>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fd(i,a,b) for(i=a;i>=b;i--)
using namespace std;
const int maxn=300000+10,mo=1000000007;
typedef long long ll;
typedef double db;
map<db,int> cnt;
int a[maxn][3],fac[maxn],inv[maxn];
int i,j,k,l,t,n,m,top,ans;
db r,sta[maxn];
int read(){
    int x=0,f=1;
    char ch=getchar();
    while (ch<'0'||ch>'9'){
        if (ch=='-') f=-1;
        ch=getchar();
    }
    while (ch>='0'&&ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}
int quicksortmi(int x,int y){
    if (!y) return 1;
    int t=quicksortmi(x,y/2);
    t=(ll)t*t%mo;
    if (y%2) t=(ll)t*x%mo;
    return t;
}
int C(int n,int m){
    if (n<m) return 0;
    return (int)((ll)fac[n]*inv[m]%mo*inv[n-m]%mo);
}
int main(){
    freopen("trokuti.in","r",stdin);freopen("trokuti.out","w",stdout);
    n=read();
    fo(i,1,n){
        j=a[i][0]=read();k=a[i][1]=read();a[i][2]=read();
        if (j&&k){
            r=-(db)j/k;
            cnt[r]++;
            if (cnt[r]==2) sta[++top]=r;
        }
    }
    fac[0]=1;
    fo(i,1,n) fac[i]=(ll)fac[i-1]*i%mo;
    inv[n]=quicksortmi(fac[n],mo-2);
    fd(i,n-1,0) inv[i]=(ll)inv[i+1]*(i+1)%mo;
    ans=C(n,3);
    fo(i,1,top){
        t=cnt[sta[i]];
        l=(C(t,3)+(ll)C(t,2)*(n-t)%mo)%mo;
        (ans-=l)%=mo;
    }
    t=0;
    fo(i,1,n) t+=(a[i][0]==0);
    l=(C(t,3)+(ll)C(t,2)*(n-t)%mo)%mo;
    (ans-=l)%=mo;
    t=0;
    fo(i,1,n) t+=(a[i][1]==0);
    l=(C(t,3)+(ll)C(t,2)*(n-t)%mo)%mo;
    (ans-=l)%=mo;
    (ans+=mo)%=mo;
    printf("%d\n",ans);
}