Pagini recente » Istoria paginii utilizator/nimeni_altu | Rating Paul Tirlisan (PaulT) | Cod sursa (job #3280464) | Cod sursa (job #454758) | Cod sursa (job #73040)
Cod sursa(job #73040)
#include <cstdio>
#define fin "secv3.in"
#define fout "secv3.out"
#define Nmax 30001
#define EPS 0.000001
#define DBG
#define FL
int N,L,U,a[Nmax],b[Nmax];
//dq
double dq[Nmax];
int pz[Nmax],vf,pr;
//end_dq
double go(double X)
{
int i;
double bst,tmp;
vf=0; pr=1; bst=(double)(-Nmax*1001);
pz[0]=1; pz[1]=1; tmp=0;
dq[1]=0;
for (i=L;i<=N;++i)
{
if (pz[pr] < i-U) ++pr;
tmp=tmp+(double)a[i-L+1]-(double)b[i-L+1]*X;
if (tmp-dq[pr]>bst)
bst=tmp-dq[pr];
while ( vf>0 && dq[vf]>tmp)
--vf;
++vf;
dq[vf]=tmp;
pz[vf]=i-L+1;
#ifdef DBGG
int j;
printf("%lf\n",bst);
for (j=pr;j<=vf;++j)
printf("%lf ",dq[j]);
printf("\n\n");
#endif
}
return bst;
}
int main()
{
int i;
double hi,lo,m,tmp;
freopen(fin,"r",stdin);
#ifdef FL
freopen(fout,"w",stdout);
#endif
scanf("%d%d%d",&N,&L,&U);
for (i=1;i<=N;++i)
scanf("%d",&a[i]);
for (i=1;i<=N;++i)
scanf("%d",&b[i]);
lo=0; hi=(double)Nmax;
while (hi-lo > EPS)
{
m=(lo+hi)/2;
tmp=go(m);
#ifdef DBGg
printf("%lf %lf %lf %lf\n",m,lo,hi,tmp);
#endif
if (tmp > 0 && tmp < EPS)
break;
if (tmp>EPS)
lo=m+EPS;
if (tmp<0)
hi=m-EPS;
}
int aux;
m*=100;
aux=(int)m;
printf("%d.%d%d\n",aux/100,aux%100/10,aux%10);
return 0;
}