Wannafly挑战赛15-E-小W的斜率（角度排序）

链接： https://www.nowcoder.com/acm/contest/112/E
来源：牛客网

题目描述

小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的斜率

#include<math.h>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define ll long long 
#define maxn 1000050
ll p,q,n;
struct node
{
	ll x,y;
}a[maxn];
ll gcd(ll a,ll b)
{
	if(a%b==0)
		return b;
	return gcd(b,a%b);
}
bool comp(node a,node b)
{
	return q*a.y-p*a.x<q*b.y-p*b.x;
}
double work(int x)
{
	double k1=a[x+1].x==a[x].x?1e18:(a[x+1].y-a[x].y)*1.0/(a[x+1].x-a[x].x);
	double k2=p*1.0/q;
	return fabs(k1-k2);
}
int main(void)
{
	scanf("%lld%lld%lld",&n,&p,&q);
	for(int i=1;i<=n;i++)
		scanf("%lld%lld",&a[i].x,&a[i].y);
	sort(a+1,a+n+1,comp);
	int id=1;double mx=work(1);
	for(int i=2;i<n;i++)
	{
		double tmp=work(i);
		if(mx>tmp) id=i,mx=tmp;
	}
	ll ansx=abs(a[id+1].x-a[id].x);
	ll ansy=abs(a[id+1].y-a[id].y);
	ll ans_gcd=gcd(ansx,ansy);
	printf("%lld/%lld\n",ansy/ans_gcd,ansx/ans_gcd);
	return 0;
}