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

来源:互联网 发布:医院大数据风控 编辑:程序博客网 时间:2024/05/20 05:53

题解:

暴力肯定过不了,

但是我们可以把除数模糊化一下,让一部分数统一用一个除数,这样就可以在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;}



0 0
原创粉丝点击