（计算几何）2020牛客寒假算法基础集训营1C. umi和弓道

2020牛客寒假算法基础集训营1C. umi和弓道

思路：

被传说中的计算几何吓到了，一直没有看，等到最后五分钟才看，一看好像不是很难欸，然后和队友口胡了算法，最后看题解的时候发现是对的QAQ。
判断$(x,y)$于$(x_0,y_0)$是不是在同一个象限，在的话就不可能挡住，不在的话，求其与x，y轴的交点的坐标，最后排序，区间判断一下最小值就可以了。

代码：

#include<bits/stdc++.h>
#define pii pair<int,int>
#define ll long long
#define cl(x,y) memset(x,y,sizeof(x))
#define ct cerr<<"Time elapsed:"<<1.0*clock()/CLOCKS_PER_SEC<<"s.\n";
const int N=1e6+10;
const int mod=1e7+9;
const int maxn=0x3f3f3f3f;
const int minn=0xc0c0c0c0;
const int inf=99999999;
using namespace std;
vector<double> disx,disy;
int main()
{
	int n,k,i,j;
	double x0,y0;
	scanf("%lf%lf%d%d",&x0,&y0,&n,&k);
	k=n-k;
	for(i=1;i<=n;i++)
	{
		double x,y;
		scanf("%lf%lf",&x,&y);
		if(x*x0<0)
			disy.push_back(y0-x0*(y0-y)/(x0-x));
		if(y*y0<0)
			disx.push_back(x0-y0*(x0-x)/(y0-y));
	}
	sort(disx.begin(),disx.end());
	sort(disy.begin(),disy.end());
	double ans=1e10;
	if(disx.size()>=k)
	{
		i=0;j=k-1;
		for(;j<disx.size();i++,j++)
			ans=min(ans,disx[j]-disx[i]);
	}
	if(disy.size()>=k)
	{
		i=0;j=k-1;
		for(;j<disy.size();i++,j++)
			ans=min(ans,disy[j]-disy[i]);
	}
	if(ans==1e10)
		printf("-1\n");
	else
		printf("%.7lf\n",ans); 
	return 0;
}