Harmonic Number (II)(模拟找规律)

题目来源https://cn.vjudge.net/problem/LightOJ-1245
【题意】
求前n项和,每一项是n/i;
【思路】
分解单个数,可以找到规律,举例说明。
当n为10时:
i ___1 2 3 4 5 6 7 8 9 10
n/i 10 5 3 2 2 1 1 1 1 1
当n为12时:
i ___1 2 3 4 5 6 7 8 9 10
n/i 12 6 4 3 2 2 1 1 1 1
此时就可以发现n/i后面的几个数是连在一起的,3,2,1而这些数完全可以按乘法去计算,没必要逐个相加,再看,第一个数10,i是从3开始进入循环的,第二个数也是从三,而3刚好是sqrt(n)(有时候需要大胆去猜),然后,就是计数,会发现,比如n==10时,有6个1,2个2,10/1-10/2==6,10/2-10/3==2,这个时候,看,答案就出来了。而前面的数,就是2^31,开方之后也会变得很小,for循环枚举就可以了。
注意一点哦,关键是我们枚举到哪呢,若是枚举到sqrt(n),要考虑会不会有重复,什么时候才会有重复呢?10/3=3,12/3=4,15/3=5;而他们的开方之后的结果都是3,但是如果模拟的话,会发先10的3算多了一次,也就是说,如果n/sqrt(n)==sqrt(n),那就减去sqrt(n)。
【代码】

#include<set>
#include<map>
#include<stack>
#include<cmath>
#include<queue>
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
const int INF=1e9;
typedef unsigned long long ll;
typedef long long LL;
int main()
{
    int T,cases=1;
    scanf("%d",&T);
    while(T--)
    {
        int n;
        scanf("%d",&n);
        int p=sqrt(n);
        LL sum=0;
        for(int i=1;i<=p;i++)
            sum+=((n/i)-(n/(i+1)))*i;
//        printf("%lld\n",sum);
        for(int i=1;i<=p;i++)
            sum+=n/i;
        if(p==n/p)
        sum-=p*1;
        printf("Case %d: %lld\n",cases++,sum);
    }
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
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谐振子的行为进行模拟和分析。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值