Harmonic Number (II) [数学]

I was trying to solve problem ‘1234 - Harmonic Number’, I wrote the following code

long long H( int n ) {
    long long res = 0;
    for( int i = 1; i <= n; i++ )
        res = res + n / i;
    return res;
}

Yes, my error was that I was using the integer divisions only. However, you are given n, you have to find H(n) as in my code.

Input

Input starts with an integer T (≤ 1000), denoting the number of test cases.

Each case starts with a line containing an integer n (1 ≤ n < 231).

Output

For each case, print the case number and H(n) calculated by the code.

Sample Input

11
1
2
3
4
5
6
7
8
9
10
2147483647

Sample Output

Case 1: 1
Case 2: 3
Case 3: 5
Case 4: 8
Case 5: 10
Case 6: 14
Case 7: 16
Case 8: 20
Case 9: 23
Case 10: 27
Case 11: 46475828386

题解

发现一下规律就可以写出 o(n) 的代码了

#include<stdio.h>
typedef long long LL;

int main()
{
    int T;LL n;
    scanf("%d",&T);
    for(int t=1;t<=T;t++){
        scanf("%lld",&n);
        LL l=2,r=n,ans=0;
        while(l<=n/l){
            ans+=(r-n/l)*(l-1);
            r=n/l;
            ans+=n/l;
            l++;
        }
        for(LL i=l;i<=r;i++) ans+=n/l;
        printf("Case %d: %lld\n",t,ans+n);
    }
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Duffing谐振子是一种非线性振动系统,具有能量耗散和非线性耦合的特点。在Matlab中,可以通过数值计算来模拟和分析Duffing谐振子的行为。 首先,我们需要定义Duffing谐振子的数学模型。Duffing谐振子的微分方程可以表示为: m*x'' + c*x' + k*x + a*x^3 = F*cos(ω*t) 其中,m是质量,c是阻尼系数,k是刚度系数,a是非线性刚度系数,F是外加力的振幅,ω是外加力的角频率。这个方程描述了Duffing谐振子的运动。 在Matlab中,可以使用ode45函数来数值求解微分方程。我们可以定义一个函数,将微分方程转化为一阶微分方程的形式。然后,可以使用ode45函数来求解这个一阶微分方程,得到系统的时间响应。 除了求解时间响应外,还可以分析Duffing谐振子的振幅-频率响应曲线。通过改变外加力的角频率ω,可以得到不同频率下的振幅。可以构建一个循环,改变外加力的频率,在每个频率下求解微分方程,得到振幅的曲线。这样就可以分析Duffing谐振子的频率特性。 此外,还可以通过绘制Duffing谐振子的相图来分析其运动特性。相图展示了系统状态变量 x 和速度变量 x' 之间的关系。可以通过定义一组初始条件,求解微分方程,并将系统的状态变量绘制在相图上。这样可以观察到系统的稳定点、周期轨道等特性。 总之,通过在Matlab中数值求解微分方程、分析振幅-频率响应曲线和绘制相图,可以对Duffing谐振子的行为进行模拟和分析。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值