login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A125022
Numbers with a unique partition as the sum of 2 squares x^2 + y^2.
8
0, 1, 2, 4, 5, 8, 9, 10, 13, 16, 17, 18, 20, 26, 29, 32, 34, 36, 37, 40, 41, 45, 49, 52, 53, 58, 61, 64, 68, 72, 73, 74, 80, 81, 82, 89, 90, 97, 98, 101, 104, 106, 109, 113, 116, 117, 121, 122, 128, 136, 137, 144, 146, 148, 149, 153, 157, 160, 162, 164, 173, 178, 180, 181
OFFSET
1,3
COMMENTS
A000161(a(n)) = 1. [Reinhard Zumkeller, Aug 16 2011]
LINKS
FORMULA
a(n) = A125021(n)/2.
MATHEMATICA
Select[Range[0, 200], Length@PowersRepresentations[#, 2, 2]==1&] (* Giorgos Kalogeropoulos, Mar 21 2021 *)
PROG
(Haskell)
import Data.List (elemIndices)
a125022 n = a125022_list !! (n-1)
a125022_list = elemIndices 1 a000161_list
-- Reinhard Zumkeller, Aug 16 2011
KEYWORD
nonn
AUTHOR
Artur Jasinski, Nov 16 2006
EXTENSIONS
Name edited by Giorgos Kalogeropoulos, Mar 21 2021
STATUS
approved