Integral Points on Conic Curves: x^2-Dy^2=±1(D \in [30000,39999])

[2016.10.22]x²-Dy²=±1 (D \in [30000,39999])の整点

• x^2-Dy^2=±1の基本解を計算するプログラム(ruby)

• D=39901
  x=11014531262071775139254857546254470587809492269499820687249462076724901879133968567573112006797369650685219620230182197005239472938348898948646187246351773125262466634425365325961551264697038934743273358781893695496147313054339601531409703779999804490704281688605707014759622883606167011136045065471332108397000481701104735716607340405235792954772896839522157297234609767206096169312297932192591384646779247672017914511529776087381817374768317586249,
  y=55140935491946958397054994797773126053184793124826483006624667226022477045736347117247141822513682441399848213782507999727160237289450363518450454070599207426051338569175139097142455575521497261368516637474292892284365576711783366721156760439812205173247275239904480997663437533110540333689646491988339876823111896306827473281260798739001091591800954772821622444796202370404268661633136005310935075991846904593090831907244580561566624216065476700;
  
• D=36061
  x=231446625469330547063653466380341650884768415567920847120798664518358126368924132836346754515774790386772530815313607071647252706791491970256814437146455834947553996272180066550370609346023922429497967074555308419257185733092605620947616443018734049595421511838831572971889789472830615589210733725955297122467621706874211608495667020868356856433380951838566277999358873515617529778484190283045449,
  y=1218798665805420760760993515869111494876030237401912760054608339994325362300034355436404461428880277087165313111628583207028913304790205289657805445957815141497046366560804440969175747705628824296206761455145989381080552138151595439261491097143632851339717135417771494953551668090322643397287222460679223936916403146798434665640970802981350640027570183733591802709127746607994166234064279677340.