【斜率优化dp】HDU 2993

第一道斜率优化dp，看着别人的报告做的，做之前建议看一下04年周源的《浅谈数形结合思想在信息学竞赛中的应用》论文         http://wenku.baidu.com/view/f6421ce8b8f67c1cfad6b8ac.html

斜率优化其实就是把每个状态看作直角坐标系上离散的点抽象出x，y 表示斜率 (y2 - y1) / (x2 - x1) 于一个关系状态i个函数的关系，然后维护点见斜率的上凸性或者下凸性。具体的情况要看于i有关的函数的单调性。
说一些注意事项：
注意斜率尽量中乘法，不要中实数，如有必要全用long long，以免溢出。注意 x2 - x1可能小于0不等式的负号要改变！

这道是论文的例题，求的是大于等于k范围的最大均值，所以需要造出下凹折线......

#define N 100005
double sum[N];
int que[N],head,tail;
int main(){
    int n,k;
    while(scanf("%d%d",&n,&k) != -1){
        int i,j;
        sum[0] = 0;
        for(i=1;i<=n;i++){
            int a;
            scanf("%d",&a);
            sum[i] = sum[i-1]+a*1.0;
        }
        head = tail = 0;
        que[0] = 0;
        double ans = 0.0;
        for(i=k;i<=n;i++){
            int now = i-k;//i-k点跟第i点比较，保证范围>=k
            //去掉上凸点
            while(head<tail){
                double y1 = sum[que[tail]] - sum[que[tail-1]];
                double x1 = que[tail] - que[tail-1];
                double y2 = sum[now] - sum[que[tail]];
                double x2 = now - que[tail];
                if(y1*x2>=y2*x1)tail--;
                else break;
            }
            que[++tail] = now;
            //去掉下凹折线的斜率小的点，因为这些点对后面的最优解没有贡献
            while(head<tail){
                double y1 = sum[que[head]] - sum[i];
                double x1 = que[head] - i;
                double y2 = sum[que[head+1]] - sum[i];
                double x2 = que[head+1] - i;
                if(y1*x2<=y2*x1)head++;
                else break;
            }
            double tmp = (sum[que[head]] - sum[i])/((que[head] - i)*1.0);
            ans = max(ans,tmp);
        }
        printf("%.2lf\n",ans);
    }
    return 0;
}

mark一下这个bloghttp://www.cnblogs.com/ronaflx/archive/2011/02/05/1949278.html，决定这几天每天都做1~2到斜率dp，好好吸收它的思维精髓！