Pagini recente » Cod sursa (job #1891269) | Cod sursa (job #1571053) | Cod sursa (job #373091) | Cod sursa (job #1034958) | Cod sursa (job #485958)
Cod sursa(job #485958)
#include<fstream>
#include<set>
using namespace std;
const char iname[]="euro.in";
const char oname[]="euro.out";
const int maxn=36000;
ifstream f(iname);
ofstream g(oname);
long long a[maxn],s[maxn],maxt,rez,dp[maxn],zi[maxn];
int i,n,t,x;
struct maxv
{
bool operator()(const int &x,const int &y)
{
return dp[x]-s[x]*i<dp[y]-s[y]*i;
}
};
struct zile
{
bool operator()(const int &x,const int &y)
{
return zi[x]<zi[y];
}
};
set<int,zile> R;
set<int,maxv> S;
int main()
{
f>>n>>t;
for(i=1;i<=n;++i)
f>>a[i],s[i]=s[i-1]+a[i];
maxt=rez=0;
S.insert(0);
for(i=1;i<=n;++i)
{
set<int,maxv>::iterator it,jt;
x=*(--S.end());
dp[i]=dp[x]+(s[i]-s[x])*i-t;
it=S.insert(i).first;
if(it!=S.end())
{
jt=it;
jt++;
if(jt!=S.end())
if(s[*it]>=s[*jt])
{
S.erase(it);
continue;
}
}
if(it!=S.begin())
{
jt=it;
jt--;
R.erase(*jt);
}
while(it!=S.begin())
{
jt=it;
jt--;
if(s[*jt]>=s[i])
R.erase(*jt),S.erase(jt),it=S.find(i);
else
break;
}
if(it!=S.begin())
{
x=dp[*it]-dp[*jt];
x=x/(s[*it]-s[*jt])+(x%(s[*it]-s[*jt])>0);
zi[*jt]=x;
R.insert(*jt);
}
if(it!=S.end())
{
jt=it;
++it;
if(jt!=S.end())
{
x=dp[*it]-dp[*jt];
x=x/(s[*it]-s[*jt])+(x%(s[*it]-s[*jt])>0);
zi[*jt]=x;
R.insert(*jt);}
}
// for(set<int,maxv>::iterator it=S.begin();it!=S.end();++it)
// fprintf(stderr,"%d ",*it);
// printf("\n");
while(R.size()>0&&zi[*R.begin()]<=i+1)
{
set<int,maxv>::iterator it,jt;
x=*R.begin();
R.erase(*R.begin());
it=S.find(x);
while(it!=S.end())
{
jt=it;
++jt;
if(jt!=S.end()&&dp[*it]-s[*it]*(i+1)>=dp[*jt]-s[*jt]*(i+1))
R.erase(*jt),S.erase(jt),it=S.find(x);
else
break;
}
if(it!=S.end())
{
jt=it;
++it;
if(jt!=S.end())
continue;
x=dp[*it]-dp[*jt];
x=x/(s[*it]-s[*jt])+(x%(s[*it]-s[*jt])>0);
zi[*jt]=x;
R.insert(*jt);
}
}
}
g<<dp[n]<<"\n";
}