A386855 Positive nonsquare integers of the form (r^2+s^2) / (1+r*s) for rational numbers r and s.
10, 20, 34, 52, 65, 73, 74, 130, 148, 160, 164, 202, 226, 241, 244, 265, 281, 290, 340, 394, 416, 436, 450, 452, 505, 514, 569, 577, 580, 586, 601, 641, 650, 720, 724, 745, 801, 802, 820, 848, 865, 884, 898, 916, 929, 970, 976, 1044, 1060, 1073, 1098, 1105, 1152, 1154, 1226, 1252, 1280, 1305, 1321, 1345
Offset: 1
Keywords
Links
- Xianwen Wang, Table of n, a(n) for n = 1..10000
- Encyclopedia of Mathematics, Legendre theorem
Programs
-
Mathematica
pool=Association[];mSize=100;Block[{bc,y},Monitor[Do[bc=Table[Times@@(Select[FactorInteger[d],Mod[#[[2]],2]==1&][[All,1]]),{d,{y^2-4,y}}];If[AllTrue[bc,#>1&],If[AllTrue[{Length@Solve[x^2==bc[[1]],x,Modulus->bc[[2]]],Length@Solve[x^2==bc[[2]],x,Modulus->bc[[1]]]},#>0&],pool[y]=Lookup[pool,y,0]+1];If[Length[pool]==mSize,Break[]]],{y,1,10^10}],{y}]];Keys[pool]
Comments