【BZOJ1011】【HNOI2008】遥远的行星 根据所允许的误差进行模糊DP

3 篇文章 0 订阅
1 篇文章 0 订阅

题解:

暴力肯定过不了,

但是我们可以把除数模糊化一下,让一部分数统一用一个除数,这样就可以在5%误差以内快速出解了。

额,但是前面一大段的动规值需要暴力求,要不就很难满足5%了。。


给个神博客http://hi.baidu.com/zeonsgtr/item/789da6f2838a3dc742c36ab7

再贴个代码

#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 101000
#define eps 1e-7
#define LL long long 
using namespace std;
int n;
LL m[N],sum[N];
double a,F[N];
int main()
{
//	freopen("test.in","r",stdin);
	scanf("%d%lf",&n,&a);
	for(int i=1;i<=n;i++)scanf("%lld",&m[i]),sum[i]=sum[i-1]+m[i];
	for(int i=min(n,2000);i>=1;i--)// 暴力求出前面的f值
	{
		int k=(int)(a*i);
		if(fabs((double)(k+1)/a-i)<eps)k++;
		for(int j=1;j<=k;j++)F[i]=F[i]+m[i]*m[j]/(double)(i-j);
	}
	for(int i=2001;i<=n;i++) // 因为题目允许有小误差,所以我们可以把一些除数来近似一下。
	{
		int k=(int)(a*i),mid=(1+k)>>1;
		F[i]=sum[k]*m[i]/(double)(i-mid);
	}
	for(int i=1;i<=n;i++)printf("%lf\n",F[i]);
	return 0;
}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值