MZ Training 2014 #3 C题

看见这个题我已开始就想用O(n*t)做，结果WA了。题解居然是o(n^2*t)的！

后来才发现是我的想法错了。

先上题：

Evelyn 非常喜欢画画。今天，她在白纸上画了很多不同颜色的彩虹。由于现在已经有很多
彩虹了，她想知道有多少能被看见？
简单地说，每个彩虹都被看做一条直线 y = aix + bi。
如果在 x 轴上存在一点 x0，使得其 y 值严格大于其他所有彩虹线段在 x0 处的 y 值（即
aix0+bi > ajx0+bj for all j != i），那么这条彩虹就能被看见。
现在你的任务是给出一些彩虹，算出能被看见的彩虹的数量
[Input]
第一行一个整数 T，表示 T 组数据。每组数据前面都有一个空行。
对于每组数据，第一行一个数 n，表示彩虹的数量。然后 n 行每行两个实数 ai, bi，表示这
条彩虹为  y = aix + bi。没有两条彩虹是一样的，即没有 a 和 b 都相同的两条直线。
[Output]
每组数据一行一个整数，能看见的彩虹的数量。
[Sample Input]
2
1
1 1
3
1 0
2 0
3 0
[Sample Output]
1
2
[Hint]
1 <= T <= 60, 1 <= n <= 5000

其实枚举每条边，并在计算与每条边相交时更新能被露出来的线段的左端点横坐标l和右端点横坐标r。if(l<r)ans++;就完了。

方程y=a1x+b1=a2x+b2解：

x=(b1-b2)/(a2-a1);

y=(a2*b1-a1*b2)/(a2-a1);

<strong><span style="font-size:18px;color:#3366ff;">#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
int t,n;
double a[5010],b[5010];
int k=1;
int ans=0;
int main()
{
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d",&n);
		ans=0;
		for(int i=1;i<=n;i++)
		{
			scanf("%lf%lf",&a[i],&b[i]);
		}
		double l,r,x;
		for(int i=1;i<=n;i++)
		{
			l=-10000000000.0;
			r=+10000000000.0;
			for(int j=1;j<=n;j++)
			{
				if(i==j) continue;
				if(a[i]==a[j])
				{
					if(b[i]<b[j])
					{
						l=1;r=0;break;
					}
				}
				else
				{
					if(a[i]>a[j])
					{
						x=(b[j]-b[i])/(a[i]-a[j]);
						l=max(l,x);
					}
					else
					{
						x=(b[j]-b[i])/(a[i]-a[j]);
						r=min(r,x);
					}
				}
				if(l>=r) break;
			}
			if(l<r) ans++;
		}
		printf("%d\n",ans);
	}
	return 0; 
}</span></strong>

惨金坷垃，吸氮磷钾，不会蒸发，亩产千八，日本粮食，再也不用，向美国进口啦~~