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.
1<n<231

1、n范围这么大,一定是不能暴力啦
2、写了几项出来,感觉应该是有规律的,但是看不出来>_<
3、解答是这么说的,如果我们把 t 作为n/i得到的值,那么它对应的数的个数是 ntnt+1 ,例:

     n =10时,t=n/i=1有5个数,为6,7,8,8,10,而(10/1)-(10/2)=5
             t=n/i=2有2个数,为4,5,而(10/2)-(10/3)=2
             t=n/i=3有1个数,为3,而(10/3)-(10/4)=1
             ``````

啊啊啊我又想多了,还在想这个式子是怎么推出来的,其实就是简单的思路,像 nn=1 nn1=1 要去找到哪个数t变为2即 n?=2 ,则 ?=n2 ,也就是t=1对应的数的个数为 n1n2 (这么写感觉自己像个笨蛋= =一般人应该一眼就看出来了吧)
4、还有当我最初看到这个递推公式的时候,我感觉没啥用啊。。。又没降低运算的复杂度,后来证明还是我比较迟钝(哭)因为当t也就是n/i变的比较小的时候出现的次数会比较多,就可以用这个公式啦,之前n/1,n/2的时候直接用题目给出的公式就好,就是没搞懂为什么把 n 作为分界点?最后还要判断一下这两个公式在 n 处有没有撞上

#include <cstdio>
#include <algorithm>
#include <iostream>
#include <cstring>
#include <queue>
#include <cmath>
using namespace std;
int main()
{
    int T, t;
    while(~scanf("%d", &T)){
        for(t=1; t<=T; t++){
            int n;
            scanf("%d", &n);
            int m=sqrt(n);
            long long int sum=0;
            for(int i=1; i<=m; i++){
                sum+=n/i;
            }
            for(int t=1; t<=m; t++){
                sum+=t*(n/t-n/(t+1));
            }
            if(n/m==m)sum-=m;

            printf("Case %d: %lld\n", t, sum);
        }
    }
    return 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、付费专栏及课程。

余额充值