BZOJ3156: 防御准备

3156: 防御准备

Time Limit: 10 Sec  Memory Limit: 512 MB
Submit: 442  Solved: 210
[Submit][Status]

Description

 

Input

第一行为一个整数N表示战线的总长度。

第二行N个整数，第i个整数表示在位置i放置守卫塔的花费Ai。

Output

共一个整数，表示最小的战线花费值。

Sample Input

10
2 3 1 5 4 5 6 3 1 2

Sample Output

18

HINT

1<=N<=10^6,1<=Ai<=10^9

Source

Katharon+#1

题解：

裸的斜率优化，比云神的题还简单吧

爆intWA了好久
代码：

 #include<cstdio>
 
 #include<cstdlib>
 
 #include<cmath>
 
 #include<cstring>
 
 #include<algorithm>
 
 #include<iostream>
 
 #include<vector>
 
 #include<map>
 
 #include<set>
 
 #include<queue>
 
 #include<string>
 
 #define inf 1000000000
 
 #define maxn 1000000+5
 
 #define maxm 500+100
 
 #define eps 1e-10
 
 #define ll long long
 
 #define pa pair<int,int>
 
 #define for0(i,n) for(int i=0;i<=(n);i++)
 
 #define for1(i,n) for(int i=1;i<=(n);i++)
 
 #define for2(i,x,y) for(int i=(x);i<=(y);i++)
 
 #define for3(i,x,y) for(int i=(x);i>=(y);i--)
 
 #define mod 1000000007
 
 using namespace std;
 
 inline int read()
 
 {
 
     int x=,f=;char ch=getchar();
 
     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
 
     while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}
 
     return x*f;
 
 }
 int n,q[maxn];
 ll a[maxn],f[maxn];
 inline double k(ll i,ll j)
 {
 
     return ((double)(f[i]+(i*i+i)/-f[j]-(j*j+j)/))/(double)(i-j);
 }
 
 int main()
 
 {
 
     freopen("input.txt","r",stdin);
 
     freopen("output.txt","w",stdout);
 
     n=read();
     for1(i,n)a[i]=read();
     int l=,r=;
     for1(i,n)
     {
         while(l<r&&k(q[l+],q[l])<i)l++;
         f[i]=f[q[l]]+(ll)(i-q[l])*(ll)(i-q[l]-)/+a[i];
         while(l<r&&k(i,q[r])<k(q[r],q[r-]))r--;
         q[++r]=i;
     }
     cout<<f[n]<<endl;
 
     return ;
 
 }