(In reply to
A solution by Math Man)
No prime other than the one already mentioned has been found under 18.8 million (that's where my program was stopped). But of course that's not a proof.
Do
p = nxtprm(p)
DoEvents
pp1 = p + 1
sq = pp1 / 2
sr = Int(Sqr(sq) + 0.5)
If sr * sr = sq Then
sr1 = sr
p2p1 = p * p + 1
sq = p2p1 / 2
sr = Int(Sqr(sq) + 0.5)
If sr * sr = sq Then
Text1.Text = Text1.Text & p & Str(p + 1) & Str(p * p + 1) & " "
Text1.Text = Text1.Text & Str((p + 1) / 2) & Str((p * p + 1) / 2) & " "
Text1.Text = Text1.Text & sr1 & Str(sr) & crlf
End If
End If
Loop
|
Posted by Charlie
on 2016-10-22 15:54:41 |