OFFSET
1,1
LINKS
EXAMPLE
Case a(1)=2: 13^4 = 119^2 + 120^2, 1^4 = 0^2 + 1^2.
Case a(3)=169: 26^4 = (-67+1)^2 + (-67+2)^2 + ... + (-67+168)^2 + (-67+169)^2.
Case a(5)=177: 295^4 = (6452+1)^2 + (6452+2)^2 + ... + (6452+176)^2 + (6452+177)^2.
...
Case a(10)=163607: 5546^4 = (-22206+1)^2 + (-22206+2)^2 + ... + (-22206+163606)^2 + (-22206+163607)^2.
MATHEMATICA
lst={}; Monitor[Do[mm=6 m^4; div=TakeWhile[Divisors[mm][[2;; -2]], 2mm/#+1>#^2&];
ans=Select[div, IntegerQ[Sqrt[(2mm/#+1-#^2)/3]]&&Mod[#-Sqrt[(2mm/#+1-#^2)/3], 2]==1&];
If[Length[ans]>0, tmp={m, {#, q=Sqrt[(2mm/#+1-#^2)/3], p=(q+1-#)/2}&/@ans}; Print[tmp];
AppendTo[lst, tmp]], {m, 1, 10^4}], m]; lst
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Xianwen Wang, May 04 2025
STATUS
approved
