BZOJ1011
2015-07-27 16:46
302 查看
传送门:BZOJ1011
看完这道题的题解后,我的内心就是崩溃的。
又TM是乱搞题。
我们暴力计算出前2000个星球的受力,而2000个星球之后就被我们直接用1到最后一个影响它的星球的质量的和除以他们编号的平均数……的平均数……的平均数。
我觉得我要学一下浮点误差分析了。
我的内心几乎是崩溃的。
我的内心几乎是崩溃的。
我的内心几乎是崩溃的。
重要的事情说三遍。
看完这道题的题解后,我的内心就是崩溃的。
又TM是乱搞题。
我们暴力计算出前2000个星球的受力,而2000个星球之后就被我们直接用1到最后一个影响它的星球的质量的和除以他们编号的平均数……的平均数……的平均数。
我觉得我要学一下浮点误差分析了。
我的内心几乎是崩溃的。
我的内心几乎是崩溃的。
我的内心几乎是崩溃的。
重要的事情说三遍。
[code]#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; }
相关文章推荐
- [Android]ToolBar使用详解(一)——项目配置
- JSON数据的生成
- easyUI的datagrid 分页保存选中状态
- 浅谈程序员接私单那点事及接私单需要注意的问题
- IIC 概述之1
- base64的js实现
- Hama查询zookeerp节点
- 两种cookie操作方式代码
- java性能调优(转载)
- Keepalived+lvs+apache
- fragment的一些探讨
- java-web 过滤器 & 监听器 & 拦截器
- RabbitMQ安装不当出现的问题:TCP connection succeeded but Erlang distribution failed
- Linux之文件操作(1)
- Mybatis的小错误--Could not find resource
- 网络编程 TCP
- 17个新手常见Python运行时错误
- BJFU1274 最大子序列和(分治)
- 如何成为一名软件架构师
- [转]Centos 6.5 优化 一些基础优化和安全设置