bzoj1007[HNOI2008]水平可见直线

还记得做过一道简单的题：

给出很多直线(给定斜率和在y轴上的截距)，求可见直线的数量。

那道题比较简单，n^2的枚举就可以过，就是对每条直线维护一个区间[l,r]，只要[l,r]不为空，这条直线就能被看见。

这里差不多也是求可见直线

这里数据是50000条，n^2是肯定不能过的.

考虑把直线按k升序排列

对于一条直线i-1和一条直线i,因为k是升序排列，它们交点的左边一定是i-1被看见，右边是i被看见。

我们再加入一条直线i+1,i+1的斜率更高，此时我们考虑3种情况：

①i+1在i和i-1交点的右边，这时候3条直线都会被看见

②三线共点，交点左边仍然是i-1被看见，此时i不能被看见了。

③i+1在交点左边,i-1仍然会被看见，i被覆盖

总结一下：

①按斜率排序后，斜率最小的一定会被看见（与其它直线最左边的交点以左不能被覆盖）

②对于之前维护了的可以被看见的直线，当前直线与维护的最后两条直线有这3种情况，如果在交点右边无影响，否则就覆盖了上一条。

③直线与交点的关系不好判断，但观察一下我们按k排了序，如果i+1在交点左边，它和任意1条直线交点都会在左边，否则在右边。所以转换一下就成了直线求交。

然后用栈维护：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<queue>
#include<cmath>
using namespace std;
const int maxn=50000+30;
int n;
struct line
{
	int a,b,id;
}L[maxn],stack[maxn];
bool cmp(line x,line y)
{
	if(x.a!=y.a)return x.a<y.a;
	return x.b>y.b;
}
int top;
double T(line x,line y)     
//y1=a1x+b1  y2=a2x+b2               (a1-a2)x=b2-b1
{
	return 1.0*(y.b-x.b)/(x.a-y.a);
}
int ans[maxn];
int cnt;
int main()
{
	top=0;
	cnt=0;
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
		L[i].id=i;
		scanf("%d%d",&L[i].a,&L[i].b);
	}
	sort(L+1,L+n+1,cmp);
	L[0].a=-600000;
	for(int i=1;i<=n;i++)
	{
		if(L[i].a!=L[i-1].a)
		{
			while(top>=2&&T(stack[top],stack[top-1])>=T(stack[top],L[i]))top--;
			stack[++top]=L[i];
		}
	}
	for(int i=1;i<=top;i++)
	{
		ans[++cnt]=stack[i].id;
	}
	sort(ans+1,ans+cnt+1);
	for(int i=1;i<=cnt;i++)printf("%d ",ans[i]);
	return 0;
}