poj 2018 数形结合 斜率 下凸折线

http://poj.org/problem?id=2018

题意：给定n个数，让你找一段区间长度大于等于F的平均值最大的区间

这道题这篇论文里面有详细的论述，具体的看论文吧，我讲一下关键的步骤。

1：一段区间的和可以表示为(s[b]-s[a-1])/(b-a+1),类似于两个点的斜率，即(a-1,s[a-1])  (b,s[b]) 两个点的斜率可以表示[a,b]区间和，s为前缀和

2：因此问题就转换为给你n个坐标点，找出x距离大于F的连线之后斜率最大的两个点

3：用类似维护凸包的方法维护一条下凸折线，这条折线具有单调性，即当前折线最优值左边的点在以后都不可能比它更优，都可以删掉，所以维护两个指针表示折线的左边右边两个点，每个点最多被加入 和 删除一次，复杂度为O(n);

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn = 100010;
struct Point{
	double x,y;
}q[maxn];
double cross(Point a,Point b,Point c){
	return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);
}
double sum[maxn];
double getk(Point a,Point b){
	 return (b.y-a.y)/(b.x-a.x);
}
int main()
{
    int N,F;
	while(scanf("%d%d",&N,&F)!=EOF)
	{
		sum[0]=0;
		for(int i=1;i<=N;i++)
		{
			scanf("%lf",&sum[i]);
			sum[i]+=sum[i-1];
		}
        int front=0,tail=0;
		int ans=0;
		for(int i=F;i<=N;i++)
		{
			Point tmp;
			tmp.x=i-F;tmp.y=sum[i-F];
			while(tail-front>=2 && cross(q[tail-2],q[tail-1],tmp) < 0)
				tail--;
			q[tail++]=tmp;
			Point now;
			now.x=i,now.y=sum[i];
			double k=getk(q[front],now);
			while(tail-front>0 && getk(q[front+1],now)>k)	
			{
				k=getk(q[front+1],now);
				front++;
			}
			if(1000*k>ans) ans=1000*k;
		}
		printf("%d\n",ans);
	}
	return 0;
}