OFFSET
1,1
COMMENTS
a(17) <= 610351562499. - David A. Corneth, Jul 23 2020
EXAMPLE
a(2) = 81 because 81^2 - 1 = 28^2 + 76^2 = 44^2 + 68^2.
PROG
(PARI) A025426(n) = my(v=valuation(n, 2), f=factor(n>>v), t=1); for(i=1, #f~, if(f[i, 1]%4>1, if(f[i, 2]%2, return(0)), t*=f[i, 2]+1)); if(t%2, t-(-1)^v, t)/2;
a(n) = my(k=2); while(A025426(k^2-1)!=n, k++); k; \\ Michel Marcus, Sep 22 2025
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Altug Alkan, Jun 28 2016
EXTENSIONS
a(10) from Chai Wah Wu, Jul 22 2020
a(13) from David A. Corneth confirmed and a(14)-a(16) added by Max Alekseyev, Sep 21 2025
STATUS
approved
