Cod sursa(job #73031)

Utilizator Binary_FireFlorin Pogocsan Binary_Fire Data 16 iulie 2007 13:24:16
Problema Secventa 3 Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.22 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;
	#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 (lo<hi)
	{
		m=(lo+hi)/2;
		tmp=go(m);
		
	#ifdef DBGg
		printf("%lf %lf %lf %lf\n",m,lo,hi,tmp);
	#endif
		
		if (tmp<EPS && tmp>0) break;
		if (tmp>EPS) 
			lo=m+1;
		if (tmp<0)
			hi=m-1;
	}

	//int aux;
	//m*=100;
	//aux=(int)m;
	//printf("%d.%d%d\n",aux/100,aux%100/10,aux%10);
	printf("%.2lf\n",m);
	return 0;
}