%I #29 Dec 01 2014 00:23:43
%S 1,3,5,7,11,13,14,15,17,19,21,23,24,29,31,33,35,37,39,41,42,43,46,47,
%T 51,52,53,55,56,57,59,61,62,65,66,67,69,70,71,73,77,79,80,83,85,87,89,
%U 91,93,94,95,97,101,103,105,107,109,111,113,114,115,117,119,120,123,127,129,131,133,137,138,139
%N Numbers n such that n^3 is an arithmetic number.
%C A065091 is a subsequence.
%H Reinhard Zumkeller and Jens Kruse Andersen, <a href="/A245644/b245644.txt">Table of n, a(n) for n = 1..10000</a> (first 147 terms from Reinhard Zumkeller)
%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Arithmetic_number">Arithmetic number</a>
%F A245656(a(n)^3) = 1. - _Reinhard Zumkeller_, Jul 28 2014
%p isArithPow := proc(n,e)
%p local dvs,d ;
%p dvs := numtheory[divisors](n^e) ;
%p add(d,d=dvs)/nops(dvs) ;
%p if type(%,'integer') then
%p true;
%p else
%p false;
%p end if;
%p end proc:
%p for n to 300 do
%p if isArithPow(n,3) then
%p printf("%d,",n) ;
%p end if;
%p end do:
%t Select[Range[120], IntegerQ[DivisorSigma[1, #^3 ]/DivisorSigma[0, #^3 ]] &] (* _Michael De Vlieger_, Aug 05 2014 after _Stefan Steinerberger_ at A003601 *)
%o (Haskell)
%o a245644 n = a245644_list !! (n-1)
%o a245644_list = filter ((== 1) . a245656 . (^ 3)) [1..]
%o -- _Reinhard Zumkeller_, Jul 28 2014
%o (Python) from sympy import divisors, divisor_count
%o [n for n in range(1,10**3) if not sum(divisors(n**3)) % divisor_count(n**3)] # _Chai Wah Wu_, Aug 04 2014
%Y Cf. A003601, A107924, A107925.
%Y Cf. A245656.
%K nonn
%O 1,2
%A _R. J. Mathar_, Jul 28 2014