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 1000000007
using 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);
}

posted @ 2017-11-03 17:10  扩展的灰(Extended_Ash)  阅读(125)  评论(0编辑  收藏  举报