计算标准佩尔方程x^2-Dy^2=1的最小解的vb程序

ysr

while True:
    try:
        a = 1
        b = 0
        c = input('请输入一个数字：')
        c = eval(c)
        i=42905136518303
        while i<=c:
            y=abs(i*i+1)/19633
            y=y**0.5
            s=i*i-19633*y*y
            if int(y)==y:
                if y<c and i*i-19633*y*y>=-1 :
                        print("x=",i,"s=",s)
            i+=1
    except:exit()

   
请输入一个数字：42905136518309
x= 42905136518304 s= 0.0
请输入一个数字：

ysr

while True:
    try:
        a = 1
        b = 0
        c = input('请输入一个数字：')
        c = eval(c)
        i=100000000
        while i<=c:
            y=abs(i*i+1)/634
            y=y**0.5
            s=634*int(y)*int(y)
            s1=i*i
            if int(y)==y:
                if y<c and s1-s==-1 :
                        print("x=",i,"y=",int(y))
            i+=1
    except:exit()

ysr

while True:
    try:
        a = 1
        b = 0
        c = input('请输入一个数字：')
        c = eval(c)
        c=c*10**200
        d=c**0.5
        d=int(d)
        d=str(d)
        a0=int(d[:len(d)-100])
        a1=a0
        r=d
        p0=a0
        q0=1
        r=10**200/int(r[len(r)-100:])
        r=int(r)
        a1=int(r/10**100)
        r=str(r)
        p1=a0*a1+1
        q1=a1
        s=1
        while int(r[:len(r)-100])!=2*a0 :
            s+=1
            r=10**200/int(r[len(r)-100:])
            r=int(r)
            a1=int(r/10**100)
            r=str(r)
            p2=a1*p1+p0
            q2=a1*q1+q0
            p0=p1
            q0=q1
            p1=p2
            q1=q2
        if s % 2 == 0 :
            x0=p0
            y0=q0
        else :
            x0=2*p0*p0+1
            y0=2*p0*q0
        print(x0,y0)
        print("循环节长度",s)
    except:exit()

ysr

Private Sub Command1_Click()
'验证佩尔方程解的程序vb程序
Dim A, B
A = Trim(Text1): B = Trim(Text2)
B = MbC(Trim(B), Trim(B))
'c = MPC(MbC(Trim(B), 1018), MbC(Trim(A), Trim(A)))
c = MPC(MbC(Trim(A), Trim(A)), MbC(Trim(B), 922))
Text3 = "左边=" & c

End Sub

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

End Sub

ysr

while True:
    try:
        a = 1
        b = 0
        c = input('请输入一个数字：')
        c = eval(c)
        i=1
        while i<=c:
            y=abs(i*i+1)/1114
            y=y**0.5
            s=1114*int(y)*int(y)
            s1=i*i
            if int(y)==y:
                if y<c and s1-s==-1 :
                        print("x=",i,"y=",int(y))
            i+=1
    except:exit()

ysr

while True:
    try:
        a = 1
        b = 0
        c = input('请输入一个数字：')
        c = eval(c)
        i=100000000000
        while i<=c:
            y=abs(i*i+1)/1594
            y=y**0.5
            s=1594*int(y)*int(y)
            s1=i*i
            if int(y)==y:
                if y<c and s1-s==-1 :
                        print("x=",i,"y=",int(y))
            i+=1
    except:exit()

   
请输入一个数字：1000000000000
请输入一个数字：

ysr

while True:
    try:
        a = 1
        b = 0
        c = input('请输入一个数字：')
        c = eval(c)
        i=1
        while i<=c:
            y=abs(i*i+7)/26
            y=y**0.5
            s=26*int(y)*int(y)
            s1=i*i
            if int(y)==y:
                if y<c and s1-s==-7 :
                        print("x=",i,"y=",int(y))
            i+=1
    except:exit()

   
请输入一个数字：100000000
请输入一个数字：

ysr

while True:
    try:
        a = 1
        b = 0
        c = input('请输入一个数字：')
        c = eval(c)
        i=100000000
        while i<=c:
            y=abs(i*i-7)/26
            y=y**0.5
            s=26*int(y)*int(y)
            s1=i*i
            if int(y)==y:
                if y<c and s1-s==7 :
                        print("x=",i,"y=",int(y))
            i+=1
    except:exit()

   
请输入一个数字：10000000000
请输入一个数字：

ysr

最新版数论探秘的目录：
第一章  孪生素数猜想和哥德巴赫猜想的初等证明
第一节  几个概念......................................................1
第二节  孪生素数猜想的证明和哥德巴赫猜想的证明..........................3
第三节  抛物线数列中的孪生素数对和相邻素数对的差的定理..................6
第四节  孪生素数对总个数及其分布规律....................................10
第五节  差为2,4,6,8，……的相邻素数对是无穷多的........................12
第六节  抛物线数列中素因子的周期性和同一周期中的对称性.................12
第二章  哥德巴赫猜想成立的必要条件和充分条件
第一节  哥德巴赫猜想成立的条件.........................................13
第二节  哥德巴赫猜想解的个数的绝对下限.................................14
第三节  偶数哥德巴赫猜想解中的最小素数的求证...........................24
第四节  偶数的哥猜拆分素数和对的下限公式及程序等.......................29
第三章  素数分布规律和哥德巴赫猜想的验证
第一节  素数的分布规律.................................................30
第二节  哥德巴赫猜想的验证.............................................32
第三节  某数内相邻素数的最大间距的公式及推导...........................33
第四章  研究素数的几个常用公式
第一节  几个常用公式...................................................37
第二节  我证明的定理...................................................39
第三节  关于素数对个数的几个命题.......................................40
第四节  关于精确的素数个数公式和素数对个数公式及哥德巴赫猜想解个数公式的
推导和探索.....................................................41
第五节  差为2m的素数对个数的比例以及特殊K生素数探索……………………48
第五章  费尔马大定理的初等证明
第一节  费尔马大定理的初等证明.........................................49
第二节  证明a^(2/3),b^(2/3),c^(2/3)之中(abc为勾股数)必有1个无理数....69
第三节  勾股小题（1）..................................................72
第四节  勾股小题（2）..................................................72
第六章  知识储备
第一节  费尔马小定理.........................................74
第二节  欧拉原理等...........................................74
第三节  中国剩余定理和求乘法的逆元...........................75
第七章  知识扩展
第一节  傅立叶变换与大整数的快速计算.........................81
第二节  朋友的一元三次方程根式解的研究.......................88
第三节  RSA密码体制及大整数的快速分解和快速素性测试.........90
第四节  梅森素数和费马数的密码特性等.........................95
第五节  李明波孪中猜想的证明.................................97
第八章  几个趣味问题
第一节  素数小题.............................................108
第二节  电话号码问题.........................................109
第三节  传令兵走多远等....................................... 109
第四节  勾股定理的平民证法....................................114
后记.............................................................117
附录1，素数表.....................................................119
附录2，两个可调用程序............................................127
附录3，李明波给美国人的挑战书.....................................132
个人简介.........................................................133