Jzoj4894 SJR的直线

给你n条直线，每条直线的方程为Aix+Biy+Ci=0，求这些直线相交组成了多少个三角形

注意，没有三条直线交于一点

十分好做，我们将所有直线按照斜率排序就好了

对于所有的平行线，我们考虑他们对答案的影响分别是什么

假设有k条直线斜率相同，那么答案就要减去C(k,2)*(n-k)+C(k,3)

注意排序的时候要将斜率算出来而不是用乘法比较

#include<stdio.h>#include<string.h>#include<algorithm>#define L long long#define M 1000000007using namespace std;struct Line{ L A,B,C; } s[300010];inline bool c1(Line a,Line b){ double k1=a.B?(double)-a.A/a.B:1e100;double k2=b.B?(double)-b.A/b.B:1e100;return k1<k2;}inline L c2(L n){ return (n*(n-1)>>1)%M; }inline L c3(L n){ return (n*(n-1)*(n-2)/6)%M; }int n,c[300010],t=0; L A;int main(){freopen("trokuti.in","r",stdin);freopen("trokuti.out","w",stdout);scanf("%d",&n); A=c3(n);for(int i=1;i<=n;++i) scanf("%lld%lld%lld",&s[i].A,&s[i].B,&s[i].C);sort(s+1,s+1+n,c1);for(int i=1,j=1;i<=n;){for(;j<=n&&s[j].A*s[i].B==s[j].B*s[i].A;++j);c[++t]=j-i; i=j;}for(int i=1;i<=t;++i)A=(A+M-(c2(c[i])*(n-c[i])%M)-c3(c[i])+M)%M;printf("%d\n",A);}