数学中国

 找回密码
 注册
搜索
热搜: 活动 交友 discuz
楼主: ysr

几个vb小程序

[复制链接]
发表于 2023-2-7 23:22 | 显示全部楼层
本帖最后由 Treenewbee 于 2023-2-7 23:34 编辑

各类计算,还是建议选择MATHEMATICA, PARI或PYTHON
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-2-8 05:58 | 显示全部楼层
13w+1是完全平方数的情况:
n取20000内的值:

n=1 d=104n^2+65=13^2,  w=0
n=14 d=104n^2+65=143^2,  w=15
n=116 d=104n^2+65=1183^2,  w=1035
n=1429 d=104n^2+65=14573^2,  w=157080
n=11831 d=104n^2+65=120653^2,  w=10767120
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-2-8 06:25 | 显示全部楼层
17w+1为完全平方数的情况:

n取20000内:
n=1  d=136*n^2+153=17^2,  w=0
n=16  d=136*n^2+153=187^2,  w=15
n=86  d=136*n^2+153=1003^2,  w=435
n=1121  d=136*n^2+153=13073^2,  w=73920
n=6019  d=136*n^2+153=70193^2,  w=2131080
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-2-8 14:00 | 显示全部楼层
本帖最后由 ysr 于 2023-2-8 06:59 编辑

while True:
    try:
        a = 1
        b = 0
        c = input('请输入一个数字:')
        c = eval(c)
        
        for i in range(c):
            
            e1 = int(i)
            w = ((15 + 4 * 14**0.5) ** (2 *  e1 + 1) + (15 - 4 * 14**0.5) ** (2 *  e1 + 1)-30)/224
            print(format(w))
            
    except:
        print("错误")
        exit()
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-2-8 14:03 | 显示全部楼层
本帖最后由 ysr 于 2023-2-8 07:12 编辑

python不精确,下面就是程序结果:
请输入一个数字:11
0.0
119.99999999999997
107879.99999999997
96876239.99999996
86994755759.99995
78121193796359.95
7.0152745034375576e+16
6.299708691967547e+19
5.657131390112352e+22
5.0800976886122e+25
4.561922067242364e+28
请输入一个数字:
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-2-8 14:05 | 显示全部楼层
如下是vb程序结果:(输入10)
w=0
w=120
w=107880
w=96876240
w=86994755760
w=78121193796360
w=7.01527450343756E+16
w=6.29970869196755E+19
w=5.65713139011235E+22
w=5.0800976886122E+25
w=4.56192206724236E+28
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-2-8 14:07 | 显示全部楼层
vb代码:(前面有了再发一下)

Private Sub Command1_Click()
Dim a, b, c
a = Val(Text1)
b = 0

Do While b <= a
b1 = Int((b + 1) / 2)
w = Val(((15 + 4 * Sqr(14)) ^ (2 * b + 1) + (15 - 4 * Sqr(14)) ^ (2 * b + 1) - 30) / 224)
s1 = s1 & " w=" & w & vbCrLf
b = Val(b + 1)

Loop
Text2 = s1
End Sub
Private Sub Command2_Click()
Text1 = ""
Text2 = ""
End Sub
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-2-8 15:58 | 显示全部楼层
Private Sub Command1_Click()
Dim a, b, c
a = Val(Text1)
Do While b <= a
d = Abs(44 * b ^ 2 + 77)
If InStr(Sqr(d), ".") = 0 Then
s = s + 1
m = (Sqr(d) - 11) / 22
w = m * (m + 1) / 2
s1 = s1 & "n=" & b & " d=44n^2+77=" & Sqr(d) & "^2,  w=" & w & vbCrLf
Else
s = s
End If

b = b + 1
Loop
If s > 0 Then
Text2 = s1
Else
Text2 = "无解"
End If


End Sub

Private Sub Command2_Click()
Text1 = ""
Text2 = ""

End Sub
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-2-8 15:59 | 显示全部楼层
本帖最后由 ysr 于 2023-2-8 08:01 编辑

w为三角数,22w+1是完全平方数
n=1 d=44n^2+77=11^2,  w=0
n=131 d=44n^2+77=869^2,  w=780
n=529 d=44n^2+77=3509^2,  w=12720
n=52139 d=44n^2+77=345851^2,  w=123567060
n=210541 d=44n^2+77=1396571^2,  w=2014886940
n取500000内的数
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-2-8 23:05 | 显示全部楼层
兔子数列前101项:

w=0
w=1
w=1
w=2
w=3
w=5
w=8
w=13
w=21
w=34
w=55
w=89
w=144
w=233
w=377
w=610
w=987
w=1597
w=2584
w=4181
w=6765
w=10946
w=17711
w=28657
w=46368
w=75025.0000000001
w=121393
w=196418
w=317811
w=514229
w=832040.000000001
w=1346269
w=2178309
w=3524578
w=5702887.00000001
w=9227465.00000001
w=14930352
w=24157817
w=39088169
w=63245986.0000001
w=102334155
w=165580141
w=267914296
w=433494437.000001
w=701408733.000001
w=1134903170
w=1836311903
w=2971215073
w=4807526976.00001
w=7778742049.00001
w=12586269025
w=20365011074
w=32951280099.0001
w=53316291173.0001
w=86267571272.0002
w=139583862445
w=225851433717
w=365435296162.001
w=591286729879.001
w=956722026041.002
w=1548008755920
w=2504730781961
w=4052739537881.01
w=6557470319842.01
w=10610209857723
w=17167680177565
w=27777890035288.1
w=44945570212853.1
w=72723460248141.2
w=117669030460994
w=190392490709135
w=308061521170130
w=498454011879265
w=806515533049395
w=1.30496954492866E+15
w=2.11148507797806E+15
w=3.41645462290672E+15
w=5.52793970088477E+15
w=8.94439432379149E+15
w=1.44723340246763E+16
w=2.34167283484677E+16
w=3.7889062373144E+16
w=6.13057907216118E+16
w=9.91948530947558E+16
w=1.60500643816368E+17
w=2.59695496911123E+17
w=4.20196140727491E+17
w=6.79891637638614E+17
w=1.10008777836611E+18
w=1.77997941600472E+18
w=2.88006719437082E+18
w=4.66004661037554E+18
w=7.54011380474637E+18
w=1.22001604151219E+19
w=1.97402742198683E+19
w=3.19404346349902E+19
w=5.16807088548585E+19
w=8.36211434898487E+19
w=1.35301852344707E+20
w=2.18922995834556E+20
w=3.54224848179263E+20
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|手机版|小黑屋|数学中国 ( 京ICP备05040119号 )

GMT+8, 2024-5-15 21:13 , Processed in 0.087890 second(s), 15 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表