小W的斜率

题目描述

小W非常喜欢某个有理数P/Q，而且他非常喜欢用它来进行一些玄学操作。他在二维平面上撒下了n个点，这些点相互不同，但一番观察之后，他失望的发现并没有任何两个点连成的直线的斜率是P/Q。你能不能告诉他在这些斜率中最接近P/Q的是多少。

输入描述:

第一行包括三个正整数n,P,Q(5<=n<=10⁶,1<=P,Q<=10⁵)

接下来n行，每行两个正整数x,y表示点的坐标(0<=x,y<=10⁹)

输出描述:

一个有理数P’/Q’表示最接近P/Q的斜率

题解：

1. 求(y2-y1)/(x2-x1)最接近P/Q,移项得Q*y2-P*x2最近接Q*y1-P*x1  
2. 排序，判断相邻的两项即可,因为相邻的两项最满足条件1，那么就将相邻的两点与p/q的值比较即可。
3. 这里要约分，并且注意正负 
   

#include<bits/stdc++.h>
using namespace std;

const int MAX = 1e6+10; 
struct node{
	long long int x,y;
	long long int s;
}inf[MAX];

int cmp(node a,node b)
{
	return a.s < b.s;	
} 

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

int main()
{
	long long int n,p,q;
	long long ax,ay;
	
	cin >> n >> p >> q;
	 
	for(int i = 0; i < n; i++){
		cin >> inf[i].x >> inf[i].y;
		inf[i].s = q*inf[i].y - p*inf[i].x;
	}
	
	sort(inf,inf+n,cmp);
	
	double minn = 0x9f9f9f;
	for(int i = 1; i < n; i++){
		long long x = inf[i].x - inf[i-1].x;
		long long y = inf[i].y - inf[i-1].y;
		
		double s = fabs(1.0*y/x-1.0*p/q);
		if (minn > s){
			minn = s;
			ax = x;
			ay = y;
		}
	}
	
	long long g = gcd(ax,ay);
	printf("%lld/%lld\n",ay/g,ax/g);
	
	return 0;
}