蓝桥杯2021年第十二届C++省赛第三题-直线

最新推荐文章于 2023-03-30 21:48:08 发布

hero_th

最新推荐文章于 2023-03-30 21:48:08 发布

阅读量461

点赞数 1

分类专栏：蓝桥杯 # 2021年第十二届省赛文章标签：蓝桥杯

本文链接：https://blog.csdn.net/weixin_45774972/article/details/123896626

版权

蓝桥杯同时被 2 个专栏收录

103 篇文章 0 订阅

订阅专栏

2021年第十二届省赛

6 篇文章 0 订阅

订阅专栏

在这里插入图片描述

Idea

两点式直线方程：Ax + By + C = 0
如已知 a (x1,y1),b(x2,y2),则
(y1-y2)x + (x2-x1)y + x1y2-x2y1 = 0
对于这个公式的推导也比较简单：
我们已知直线的斜率K = (y1-y2) / (x1-x2)
故我们直线方程现在可以写成
y = kx+b

我们先 a,b 两点代入该方程，我们得到两个方程

y1 = k x1 + b
y2 = k x2 + b

然后我们两个方程相加，代入斜率k的表达式，然后将新的等式进行通分，展开之后，我们可以很轻易的得到 b = (x1y2-x2y1) / (x1-x2)

得到k，b的表达式之后，我们就可以用两点来确定一条唯一的直线。
但是这里我们需要注意的是，我们这样得到的斜率k和截距b有可能得到的不是整数类型，这个时候呢，对于我们计算机代码的处理就会产生误差，数据不好处理。
因此我们将y = kx + b ，分别把k,b的表达式代入之后，我们依然进行通分，移项。我们就可以得到一个一般式:(y1-y2)x + (x2-x1)y + x1y2-x2y1 = 0
最后，我们还需要对A,B,C进行一个处理就是，化简最简型，分别除以三个数的最大公因数，存储到set集合中，进行去重。

Code

#include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
#include<queue>
#include<string>
#include<unordered_set> 
#include<set>
#include<unordered_map> 
using namespace std;
#define MAX_SIZE 100005
#define MOD 1000000009
typedef long long ll;
typedef pair<int,int> PII;

set<pair<PII,int>> map;
vector<PII> vec;


int gcd(int a,int b){
	if(!b) return a;
	return gcd(b,a % b);
}

int main(){
	int m,n;
	cin>>m>>n;
	for(int i = 0 ;i <=m ;i++){
		for(int j = 0 ; j <= n;j++){
			vec.push_back({i,j});
		}
	}
	for(int i = 0 ; i < vec.size() ;i++){
		for(int j = i + 1 ; j < vec.size() ;j++){
			int x1 = vec[i].first,y1 = vec[i].second;
			int x2 = vec[j].first,y2 = vec[j].second;
			
			int A = y1 - y2 ;
			int B = x2 - x1 ;
			int C = x1 * y2 - x2 * y1;
			int fact = gcd(gcd(A,B),C);
			map.insert({{A / fact , B / fact} , C / fact});
		} 
	}
	cout<< map.size();
	return 0;
}