试题B：直线

答案：40257

思路一：
刚开始想的是求k和b，用pair<double,double>存储，然后再用set装数，再用set.size()求，但是发现一件事，double类型好像用set就不怎么容易去重了。（新手不推荐）
需要注意！！！
计算b时候用b=y-kx，k是浮点数会造成二次误差。
因C++无控制double变量小数位数的函数，所以要注意截距b在计算时不可直接用上述算出的k来直接带入计算，因为double本身取值会造成误差，带入计算则会造成二次误差，影响结果。在这里对b的计算进行化简合并：
b=(y2x1-y1x2)/(x2-x1)

int main(){
	set<pair<double,double> > s;
	pair<int,int> p[430];
	int cnt=0;
	for(int i=0;i<20;i++)
		for(int j=0;j<21;j++){
			p[cnt].first=i,p[cnt].second=j;
			cnt++;
		}
	for(int i=0;i<cnt;i++)
		for(int j=i+1;j<cnt;j++){
			pair<int,int> a=p[i],b=p[j];
			pair<double,double> line;
			double x=a.first-b.first,y=a.second-b.second;
			if(x==0) {
				line.first=0;
				line.second=a.second;
			}
			else{
				line.first=y/x;
				line.second=(a.first*b.second-a.second*b.first)/x;
			}
			s.insert(line);
		}
	cout<<s.size();
}

思路二：
所以，改用直线：ax+by+c的形式。求a,b,c，所以-a/b=k，且c/b=b，所以令a=y2-y1,b=x2-x1,c=y2x1-y1x2。a是[-20,20],b是[-19,19],c是[-380,380]，a,b,c都是整数，所以可以用set进行求解，再把set的cmp进行重写。
直线的确立：2个点确立一个直线，4个点则对应3，2，1，所以for循环时候要注意去重。
注意：由于ax+by+c=0，可以约分，所以求a,b,c最大公因数，再/最大公因数，然后再存

struct node{
	int a,b,c;
	bool operator <(const node &n) const{
		if(a==n.a) {
			if(b==n.b) return c<n.c;
			return b<n.b;
		}
		return a<n.a;
	}
}nod[500];

int gcd(int a,int b){
	return b==0?a:gcd(b,a%b);
}

int main(){
	set<node> s;
	int cnt=0;
	for(int i=0;i<20;i++)
		for(int j=0;j<21;j++){
			nod[cnt].a=i,nod[cnt].b=j,nod[cnt].c=0;
			cnt++;
		}
	
	for(int i=0;i<cnt;i++)
		for(int j=i+1;j<cnt;j++){
			node a=nod[i],b=nod[j],line;
			line.a=b.b-a.b;
			line.b=b.a-a.a;
			line.c=b.b*a.a-b.a*a.b;
			//由于ax+by+c，可以约分，所以求a,b,c最大公因数，再/最大公因数，然后再存 
			int tmp=gcd(line.a,gcd(line.b,line.c));
			line.a/=tmp;
			line.b/=tmp;
			line.c/=tmp;
			s.insert(line);
		}
	cout<<s.size();
}