Cod sursa(job #73039)

Utilizator Binary_FireFlorin Pogocsan Binary_Fire Data 16 iulie 2007 14:06:55
Problema Secventa 3 Scor 10
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.32 kb
#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;
	}
	//tmp=go(0.83);
	//printf("%lf\n",tmp);
	//int aux;
	//m*=100;
	//aux=(int)m;
	//printf("%d.%d%d\n",aux/100,aux%100/10,aux%10);
	//printf("%lf %lf %lf\n",m,lo,hi);
	printf("%.10lf\n",m);
	return 0;
}