OFFSET
1,2
EXAMPLE
159 is in the sequence because L((1+5+9)/159) = L(15/159) = 1 where L(a/b) is the Legendre symbol of a and b, which is defined to be 1 if a is a quadratic residue (mod b) and -1 if a is a quadratic non-residue (mod b).
MAPLE
isA177053 := proc(n) is(numtheory[quadres](A007953(n), n) =1 ); end proc:
for n from 1 to 140 do if isA177053(n) then printf("%d, ", n); end if; end do: # R. J. Mathar, Dec 12 2010
PROG
(PARI) isok(n) = issquare(Mod(sumdigits(n), n)); \\ Michel Marcus, Feb 01 2016
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, Dec 09 2010
STATUS
approved