poj 2018 Best Cow Fences

题目链接：http://poj.org/problem?id=2018

这道题参考了周源的论文:

最后一个问题就是如何求过pt点，且与折线相切的直线了。一种直接的方法就是二分，每次查找的复杂度是O(log2N)。但是从图形的性质上很容易得到另一种更简便更迅速的方法：由于折线上过每一个点切线的斜率都是一定的3，而且根据下凸函数斜率的单调性，如果在检查点Pt 时找到了折线上的已知一个切点A，那么A以前的所有点都可以删除了：过这些点的切线斜率一定小于已知最优解，不会做出更大的贡献了。于是另外保留一个指针不回溯的向后移动以寻找切线斜率即可，平摊复杂度为为O(1)。

代码如下：

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<string>
#include<queue>
#include<algorithm>
#include<vector>
#include<iostream>
#include<map>
using namespace std;
#define Max 100100
#define inf 0x3f3f3f3f
struct node
{
 double x,y;
}q[Max];
double sum[Max],temp;
int num,max; 
double getk(node &p0,node &p1)
{
 return (p1.y-p0.y)/(p1.x-p0.x);
}
double cross(node &p0,node &p1,node &p2)
{
 return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
int main()
{
 int n,f,i,j;
 scanf("%d%d",&n,&f);
 for(i=1;i<=n;i++)
 {
  scanf("%d",&num);
  sum[i]=sum[i-1]+num;
 } 
 int max=0;
 int top=0,tail=0;
 for(i=f;i<=n;i++)
 {
  node newp;
  newp.x=i-f;
  newp.y=sum[i-f];
  while(tail-top>=2&&cross(q[tail-2],q[tail-1],newp)<0)
   tail--;
  q[tail++]=newp;
  node z;
  z.x=i;z.y=sum[i];
  double k=getk(q[top],z);
  while(tail-top>0&&(temp=getk(q[top+1],z))>k)
  {
   k=temp;
   top++;
  }
  if((1000*k)>max)
   max=1000*k;
 }
 printf("%d\n",max);
 return 0;
}