bzoj1011: [HNOI2008]遥远的行星
题目链接
题解
emmm玄学误差题
之前有做过一道min_hash的trick题,有空整理一下 (坑
a的上限0.35
\(f_i=\sum_{j=1}^{a\times{i}}\frac{m_i\times{m_j}}{i-j}\)
暴力复杂的(n*0.35n)
过不掉
对于一个i来说若可以 O1 求出\(\sum_{j = 1}^{a*i} \frac{1}{i-j} / (i * a)\)的近似值
就可利用m的前缀和O1计算了
发现j的范围是\([1,a \times{i})\)
令-\(a\times{i}\over 2\)这个定值替换到上式分母
a的上限0.35,当i增大的时,分母越来越大,误差是越来越小的
然后你可算一下i增大到何时时误差小于5%,比这个i小是暴力,大时....
代码
#include<cstdio>
#include<algorithm>
int n;
#define eps 1e-8
const int maxn = 100005;
double a,f[maxn],m[maxn];
double sum[maxn];
int main() {
scanf("%d%lf",&n,&a);
for(int i = 1;i <= n;++ i) scanf("%lf",m + i),sum[i] =sum[i - 1] + m[i];
for(int i = 1;i <= n;++ i) {
int t = int(a * i + eps);
if(i <= 500)
for(int j = 1;j <= t;++ j)
f[i] += m[i] * m[j] / (i - j);
else f[i] += sum[t] * m[i] / double ((i - t / 2 )) ;
}
for(int i = 1;i <= n;++ i) {
printf("%lf\n",f[i]);
}
return 0;
}
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· .NET 9 new features-C#13新的锁类型和语义
· Linux系统下SQL Server数据库镜像配置全流程详解
· 现代计算机视觉入门之:什么是视频
· 你所不知道的 C/C++ 宏知识
· 聊一聊 操作系统蓝屏 c0000102 的故障分析
· DeepSeek V3 两周使用总结
· 回顾我的软件开发经历(1)
· C#使用yield关键字提升迭代性能与效率
· 低成本高可用方案!Linux系统下SQL Server数据库镜像配置全流程详解
· 4. 使用sql查询excel内容
2017-05-08 luogu P3374 【模板】树状数组 1
2017-05-08 线段树标版
2017-05-08 codevs 1082 线段树练习 3
2017-05-08 codevs1081 线段树练习 2<区间修改>