牛客练习赛11 E、求最值

题目链接：https://www.nowcoder.com/acm/contest/59/E

解析：表达式显然为：求点(i,Si)与点(j,Sj)的距离的平方的最小值（i！=j），即找最近点对；

          分治法处理平面最近点对距离问题；旋转卡壳+凸包算法处理平面最远点对距离问题；

  

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <math.h>
#define llt long long
using namespace std;
const int Size=1e5+100;
llt S[Size];//Si记录a1至ai的和；
llt b[Size];

struct Point{
    llt x;
    llt y;
}p[Size];

llt dist_pow(Point &A,Point &B){ return (A.x-B.x)*(A.x-B.x)+(A.y-B.y)*(A.y-B.y);}

bool cmp(int A,int B){
    return (p[A].y<p[B].y);
}
llt Div_Con(int l,int r){
    //cout<<1<<endl;
    if(r==l) return 1e18;
    if(r-l==1) return dist_pow(p[l],p[r]);
    int mid=(l+r)>>1;
    llt d=min(Div_Con(l,mid),Div_Con(mid+1,r));

    int cnt=0;
    llt D=sqrt(d)+1;
    for(int i=mid;i>=l&&p[mid+1].x-p[i].x<D;--i){
        b[cnt++]=i;
    }
    for(int i=mid+1;i<=r&&p[i].x-p[mid].x<D;++i)
        b[cnt++]=i;

    sort(b,b+cnt,cmp);
    for(int i=0;i<cnt;++i)
        for(int j=i+1;j<cnt&&p[b[j]].y-p[b[i]].y<D;++j)
            d=min(d,dist_pow(p[b[j]],p[b[i]]));

    return d;
}

int main(){
    int n;scanf("%d",&n);
    p[0].y=0;
    llt a;
    for(int i=1;i<=n;++i){
        p[i].x=i;
        scanf("%lld",&a);
        p[i].y=p[i-1].y+a;
    }

    printf("%lld\n",Div_Con(1,n));
    return 0;
}