hdoj 2993 斜率优化

对于这道题没有头绪的，强烈推荐阅读这篇论文：http://wenku.baidu.com/view/b97cd22d0066f5335a8121a3.html

很经典啊！网上有些人觉得论文最后的O（n）的算法不对，其实是正确的，不要质疑论文的正确性，自己认真想一下就知道是对的。

别的不说了，去看论文吧。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAX 100111
int a[MAX];
double sum[MAX];
int q[MAX];//队列
int N,K;
int head,tail;
double ans;
double Max(double a,double b)
{
    if(a>b) return a;
    return b;
}
double getDP(int i,int j)
{
    return (sum[j]-sum[i])/double(j-i);
}
double getUP(int i,int j) //yj-yk部分
{
    return (sum[j]-sum[i]);
}
int main()
{
    int i,j;
    while(scanf("%d%d",&N,&K)!=EOF)
    {
        for(i=1;i<=N;i++)
        {
            scanf("%d",&a[i]);
        }
        sum[0]=0;
        for(i=1;i<=N;i++)
            sum[i]=a[i]+sum[i-1];
        head=tail=0;
        q[tail++]=0;
        ans=0;
        for(i=K;i<=N;i++)
        {
            while(head+1<tail&&getUP(q[head],i)*(i-q[head+1])<=getUP(q[head+1],i)*(i-q[head]))
            {
                head++;
            }
            ans=Max(ans,getDP(q[head],i));
            j=i-K+1;
            while(head+1<tail&&getUP(q[tail-1],j)*(q[tail-1]-q[tail-2])<=getUP(q[tail-2],q[tail-1])*(j-q[tail-1]))
                tail--;
            q[tail++]=j;
        }
        printf("%.2lf\n",ans);
    }
    return 0;
}