OFFSET
1,2
COMMENTS
m such that A195812(m) is a perfect square.
EXAMPLE
a(8) = 26 because x^26 == > 0, 1, 3, 4, 9, 10, 12, 13, 14, 16, 17, 22, 23, 25 (mod 26), and the sum = 169 = 13^2.
MAPLE
sumSquares := proc(n)
local re, x, r ;
re := {} ;
for x from 0 to n-1 do
re := re union { modp(x^n, n) } ;
end do:
add(r, r=re) ;
end proc:
for n from 1 to 750 do
z:= sqrt(sumSquares(n));
if z=floor(z) then
printf("%d, ", n);
end if;
end do: #
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Oct 05 2011
STATUS
approved