计算标准佩尔方程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()