Harmonic Number LightOJ - 1234

最新推荐文章于 2020-08-13 00:36:27 发布

klchen0112

最新推荐文章于 2020-08-13 00:36:27 发布

阅读量161

点赞数

CC 4.0 BY-SA版权

分类专栏： number theory

本文链接：https://blog.youkuaiyun.com/weixin_40959045/article/details/88381903

number theory 专栏收录该内容

24 篇文章

订阅专栏

本文介绍了一种高效的算法，用于每隔64个元素进行一次求和操作，并存储中间结果，以加速后续的求和计算。该算法适用于大规模数据集，通过预处理将计算复杂度降低，特别适合于内存和CPU资源有限的场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

每隔64个打个表

 #include<bits/stdc++.h>
 using namespace std;
 #define fst first
 #define sec second
 #define sci(num) scanf("%d",&num)
 #define scl(num) scanf("%lld",&num)
 #define mem(a,b) memset(a,b,sizeof a)
 #define cpy(a,b) memcopy(a,b,sizeof b)
 typedef long long LL;
 typedef pair<int,int> P;
 const int MAX_N = 1e8;
 const int MAX_M = 10000;
 double sum[(MAX_N >> 6) + 10];
 void init() {
     double t = 0;
     for (int i = 1;i <= MAX_N;i++) {
         t += 1.0 / i;
         if ((i & 63) == 0) sum[i >> 6] = t;
     }
 }
 int main() {
     init();
     int T;
     sci(T);
     for (int cs = 1;cs <= T;cs++) {
         int N;
         sci(N);
         double ans = sum[(N >> 6)];
         for (int i = ((N>> 6)<<6) + 1;i <= N;i++) {
             ans += 1.0/i;
         }
         printf("Case %d: %.15lf\n",cs,ans);
     }
     return 0;
 }