BZOJ1011

来源:互联网 发布:联通网络ip地址丢失 编辑:程序博客网 时间:2024/06/05 11:39

传送门:BZOJ1011

看完这道题的题解后,我的内心就是崩溃的。
又TM是乱搞题。
我们暴力计算出前2000个星球的受力,而2000个星球之后就被我们直接用1到最后一个影响它的星球的质量的和除以他们编号的平均数……的平均数……的平均数。
我觉得我要学一下浮点误差分析了。
我的内心几乎是崩溃的。
我的内心几乎是崩溃的。
我的内心几乎是崩溃的。
重要的事情说三遍。

#include <cstdio>#include <cstdlib>#include <cmath>#include <cstring>#include <algorithm>#include <iostream>using namespace std;const double eps=1e-6;int n,m;double A;double a[100005];double s[100005];double f[100005];void Readdata(){    freopen("loli.in","r",stdin);    scanf("%d%lf",&n,&A);    for(int i=1;i<=n;i++){        scanf("%lf",&a[i]);        s[i]=s[i-1]+a[i];    }}void Solve(){    int Q=min(n,2000);    for(int i=1;i<=Q;i++){        int p=floor(A*i+eps);        for(int j=1;j<=p;j++)            f[i]+=(a[j]*a[i]/(i-j)+eps);    }    for(int i=Q+1;i<=n;i++){        int p=floor(A*i+eps);        f[i]=s[p]*a[i]/(i-(1+p)/2)+eps;    }    for(int i=1;i<=n;i++)        printf("%.6lf\n",f[i]);}void Close(){    fclose(stdin);    fclose(stdout);}int main(){    Readdata();    Solve();    Close();    return 0;}
0 0
原创粉丝点击