【BZOJ 1011】 [HNOI2008]遥远的行星
1011: [HNOI2008]遥远的行星
Time Limit: 10 Sec Memory Limit: 162 MBSec Special Judge
Submit: 2131 Solved: 754
[ Submit][ Status]
Description
直线上N颗行星,X=i处有行星i,行星J受到行星I的作用力,当且仅当i<=AJ.此时J受到作用力的大小为 Fi->j=Mi*Mj/(j-i) 其中A为很小的常量,故直观上说每颗行星都只受到距离遥远的行星的作用。请计算每颗行星的受力,只要结果的相对误差不超过5%即可.
Input
第一行两个整数N和A. 1<=N<=10^5.0.01< a < =0.35
接下来N行输入N个行星的质量Mi,保证0<=Mi<=10^7
Output
N行,依次输出各行星的受力情况
Sample Input
3
5
6
2
4
Sample Output
0.000000
0.000000
1.968750
2.976000
HINT
精确结果应该为0 0 0 2 3,但样例输出的结果误差不超过5%,也算对
用近似估计来做。
很神的一道题,读题的时候就觉得 结果误差不超过5%,也算对这句话非常奇怪。
然后看了题解才明白这句话是这道题的关键啊!
因为当行星的相距很远时,作用力会非常小,完全可以估算,即把前1-now的行星到当前行星的距离都算做到now/2的距离,O(1)求解。
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstdio>
#define eps 1e-8
using namespace std;
int n;
double a,ans[100005],m[100005],pre[100005];
int main()
{scanf("%d %lf",&n,&a);for (int i=1;i<=n;i++)scanf("%lf",&m[i]);int now=0;ans[1]=0.0;for (int i=2;i<=min(2000,n);i++){if ((double)i*a+eps>=(double)(now+1)) now++;for (int j=1;j<=now;j++)ans[i]=ans[i]+m[i]*m[j]/(double)(i-j);}pre[0]=0.0;for (int i=1;i<=n;i++)pre[i]=pre[i-1]+m[i];for (int i=2001;i<=n;i++){if ((double)i*a+eps>=(double)(now+1)) now++;ans[i]=pre[now]*m[i]/(double)(i-(int)now/2);}for (int i=1;i<=n;i++)printf("%.6lf\n",ans[i]);return 0;
}
感悟:
1.近似估计这种做法感觉非常神~
2.一开始WA了,加上eps之后就过了。
第一次遇到这种问题,以后做小数的题要加上浮点误差,尤其是数据精度非常大的时候。
【BZOJ 1011】 [HNOI2008]遥远的行星相关推荐
- BZOJ 1011: [HNOI2008]遥远的行星
1011: [HNOI2008]遥远的行星 Time Limit: 10 Sec Memory Limit: 162 MBSec Special Judge Submit: 4814 Solve ...
- [BZOJ 1011][HNOI2008]遥远的行星(奇技淫巧)
Description 直线上N颗行星,X=i处有行星i,行星J受到行星I的作用力,当且仅当i<=AJ.此时J受到作用力的大小为 Fi->j= Mi*Mj/(j-i) 其中A为很小的常量, ...
- 1011: [HNOI2008]遥远的行星
1011: [HNOI2008]遥远的行星 Time Limit: 10 Sec Memory Limit: 162 MBSec Special Judge Submit: 2241 Solve ...
- [HNOI2008]遥远的行星
https://www.luogu.org/problem/show?pid=3198 小于500的挨着算 大于的分母取平均值 #include<cstdio> using namespa ...
- [bzoj1011] [HNOI2008]遥远的行星
Description 直线上N颗行星,X=i处有行星i,行星J受到行星I的作用力,当且仅当i<=AJ.此时J受到作用力的大小为 Fi->j=Mi*Mj/(j-i) 其中A为很小的常量,故 ...
- BZOJ 1005: [HNOI2008]明明的烦恼
BZOJ 1005: [HNOI2008]明明的烦恼 Description 自从明明学了树的结构,就对奇怪的树产生了兴趣......给出标号为1到N的点,以及某些点最终的度数,允许在 任意两点间连线 ...
- BZOJ 1009 [HNOI2008]GT考试
1009: [HNOI2008]GT考试 Time Limit: 1 Sec Memory Limit: 162 MB Submit: 2154 Solved: 1327 [Submit][Sta ...
- bzoj 1010: [HNOI2008]玩具装箱toy 2011-12-27
1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec Memory Limit: 162 MB Submit: 2254 Solved: 744 [Submit][S ...
- BZOJ 1008 [HNOI2008]越狱
1008: [HNOI2008]越狱 Time Limit: 1 Sec Memory Limit: 162 MB Submit: 5166 Solved: 2242 [Submit][Statu ...
最新文章
- [字符编码]Invalid byte 1 of 1-byte UTF-8 sequence终极解决方案
- 从离散值中把值相近的放在一起
- 十家全国学会就IEEE“审稿门”事件发表联合声明
- MFC中的几个常用类——CFileDialog
- Docker容器中的Linux机器快速设置国内源
- VMware 禁用虚拟内存文件,提升虚拟机响应速度 .
- 最简单的基于FFmpeg的移动端例子附件:IOS自带播放器
- Mybatis学习(2)—— 对象关系映射
- HC3i论坛电子病历系统热门资料30篇(二)(
- php led显示屏控制软件下载,中航led控制软件
- steamcom启动服务:443端口被占用,请关闭占用该端口的进程后再点击启动服务!
- suse linux 10 下载,SUSE Linux 10下载
- 台式电脑锁定计算机,电脑锁住了怎么解锁
- Java解析富文本rtf中文乱码
- 浅谈矩阵加速——以时间复杂度为O(log n)的算法实现裴波那契数列第n项及前n之和使用矩阵加速法的优化求法
- 技术人生:故事之八	OFFICE是软件打字机?
- Mac安装与使用MacTeX
- 【java毕业设计】基于javaEE+原生Servlet+MySql的企业员工信息管理系统设计与实现(毕业论文+程序源码)——企业员工信息管理系统
- 全球与中国回音壁市场深度研究分析报告
- 数字证书连接服务器异常,连接时Socket.io + SSL +自签名CA证书出现错误