OFFSET
1,1
COMMENTS
a(n) is of the form m^2 or 2*m^2.
See the comments in A195268 (numbers such that the sum of the odd divisors is prime).
It is interesting to observe that the intersection of this sequence with A195268 gives {9, 18, 36, 72, 144, 288, 576, 1152, 2304, 4608, 9216, 18432, 36864, 73728, 146689, 147456, 293378, 294912,...} and contains the sequence A005010(n) (numbers of the form 9*2^n), but is not equal to this sequence. For example, up to n = 400000, the numbers 146689 and 293378 are not divisible by 9.
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
The divisors of 18 are { 1, 2, 3, 6, 9, 18}, and the sum of the cube of the odd divisors 1^3 + 3^3 + 9^3 =757 is prime. Hence 18 is in the sequence.
MAPLE
with(numtheory):for n from 1 to 400000 do:x:=divisors(n):n1:=nops(x):s:=0:for m from 1 to n1 do:if irem(x[m], 2)=1 then s:=s+x[m]^3:fi:od:if type(s, prime)=true then printf(`%d, `, n): else fi:od:
MATHEMATICA
Module[{c=Range[800]^2, m}, m=Sort[Join[c, 2c]]; Select[m, PrimeQ[Total[ Select[ Divisors[#], OddQ]^3]]&]](* Harvey P. Dale, Jul 31 2012 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Sep 15 2011
STATUS
approved