数学中国

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

求无理平方根连分数表示法的循环节的程序

[复制链接]
 楼主| 发表于 2025-8-20 07:55 | 显示全部楼层
蔡家雄 发表于 2025-8-19 13:42
模 257 的平方剩余奇质数之一是 241,

求 \(x^2 - 257*y^2= 241\) 的正整数解,,

请输入一个数字:100000000
x= 193 y= 12
x= 321 y= 20
x= 197697 y= 12332
x= 329153 y= 20532
请输入一个数字:
回复 支持 反对

使用道具 举报

 楼主| 发表于 2025-8-20 07:57 | 显示全部楼层
蔡家雄 发表于 2025-8-19 13:44
模 257 的平方剩余奇质数之一是 241,

求 \(x^2 - 257*y^2= - 241\) 的正整数解,,

请输入一个数字:100000000
x= 4 y= 1
x= 6172 y= 385
x= 10276 y= 641
x= 6332476 y= 395009
x= 10543172 y= 657665
请输入一个数字:
回复 支持 反对

使用道具 举报

 楼主| 发表于 2025-8-20 22:41 | 显示全部楼层
蔡家雄 发表于 2025-8-20 12:27
模 257 的平方剩余奇质数之一是 157,

求 \(x^2 - 257*y^2= 157\) 的正整数解,,

请输入一个数字:100000000
x= 97 y= 6
x= 417 y= 26
x= 99105 y= 6182
x= 427745 y= 26682
请输入一个数字:
回复 支持 反对

使用道具 举报

 楼主| 发表于 2025-8-20 22:42 | 显示全部楼层
蔡家雄 发表于 2025-8-20 12:27
模 257 的平方剩余奇质数之一是 157,

求 \(x^2 - 257*y^2= - 157\) 的正整数解,,

请输入一个数字:100000000
x= 10 y= 1
x= 3094 y= 193
x= 13354 y= 833
x= 3174454 y= 198017
x= 13701194 y= 854657
请输入一个数字:
回复 支持 反对

使用道具 举报

 楼主| 发表于 2025-8-22 08:25 | 显示全部楼层
蔡家雄 发表于 2025-8-22 00:20
模 37 的平方剩余奇质数之一是 7,

求 \(x^2 - 37*y^2= 67\) 的正整数解,,

请输入一个数字:100000000
x= 20 y= 3
x= 128 y= 21
x= 2792 y= 459
x= 18668 y= 3069
x= 407612 y= 67011
x= 2725400 y= 448053
x= 59508560 y= 9783147
回复 支持 反对

使用道具 举报

 楼主| 发表于 2025-8-22 08:25 | 显示全部楼层
蔡家雄 发表于 2025-8-22 00:21
模 37 的平方剩余奇质数之一是 7,

求 \(x^2 - 37*y^2= - 67\) 的正整数解,,

请输入一个数字:100000000
x= 9 y= 2
x= 231 y= 38
x= 1545 y= 254
x= 33735 y= 5546
x= 225561 y= 37082
x= 4925079 y= 809678
x= 32930361 y= 5413718
回复 支持 反对

使用道具 举报

 楼主| 发表于 2025-8-22 09:01 | 显示全部楼层
蔡家雄 发表于 2025-8-22 00:58
模 37 的平方剩余奇质数之一是 7,

求 \(x^2 - 37*y^2= 881\) 的正整数解,,

请输入一个数字:100000000
x= 57 y= 8
x= 609 y= 100
x= 7713 y= 1268
x= 88857 y= 14608
x= 1126041 y= 185120
x= 12972513 y= 2132668
请输入一个数字:
回复 支持 反对

使用道具 举报

 楼主| 发表于 2025-8-22 09:02 | 显示全部楼层
蔡家雄 发表于 2025-8-22 00:59
模 37 的平方剩余奇质数之一是 7,

求 \(x^2 - 37*y^2= - 881\) 的正整数解,,

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

   
请输入一个数字:100000000
x= 46 y= 9
x= 638 y= 105
x= 7354 y= 1209
x= 93194 y= 15321
x= 1073638 y= 176505
x= 13605686 y= 2236761
回复 支持 反对

使用道具 举报

 楼主| 发表于 2025-8-22 09:38 | 显示全部楼层
蔡家雄 发表于 2025-8-22 01:31
求 \(x^2 - 37*y^2= 137\) 的正整数解,,

请输入一个数字:100000000
x= 27 y= 4
x= 195 y= 32
x= 3747 y= 616
x= 28443 y= 4676
x= 547035 y= 89932
x= 4152483 y= 682664
x= 79863363 y= 13129456
回复 支持 反对

使用道具 举报

 楼主| 发表于 2025-8-22 09:39 | 显示全部楼层
蔡家雄 发表于 2025-8-22 01:32
求 \(x^2 - 37*y^2= - 137\) 的正整数解,,

请输入一个数字:100000000
x= 14 y= 3
x= 310 y= 51
x= 2354 y= 387
x= 45274 y= 7443
x= 343670 y= 56499
x= 6609694 y= 1086627
x= 50173466 y= 8248467
请输入一个数字:
回复 支持 反对

使用道具 举报

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

本版积分规则

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

GMT+8, 2026-7-25 04:56 , Processed in 0.135453 second(s), 13 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

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