login
Odd refactorable numbers.
14

%I #36 Oct 04 2024 20:48:46

%S 1,9,225,441,625,1089,1521,2025,2601,3249,4761,5625,6561,7569,8649,

%T 12321,15129,16641,19881,25281,31329,33489,35721,40401,45369,47961,

%U 50625,56169,62001,71289,84681,91809,95481,99225,103041,106929,114921

%N Odd refactorable numbers.

%C Odd refactorable numbers are always squares.

%C All terms = 1 (mod 8). [_Zak Seidov_, May 25 2010]

%H Giovanni Resta, <a href="/A036896/b036896.txt">Table of n, a(n) for n = 1..10000</a> (first 1001 terms from Harvey P. Dale)

%H S. Colton, <a href="http://www.cs.uwaterloo.ca/journals/JIS/colton/joisol.html">Refactorable Numbers - A Machine Invention</a>, J. Integer Sequences, Vol. 2, 1999, #2.

%H S. Colton, <a href="http://web.archive.org/web/20070831060523/http://www.dai.ed.ac.uk/homes/simonco/research/hr/">HR - Automatic Theory Formation in Pure Mathematics</a>

%e 9 is refactorable because tau(9)=3 and 3 divides 9.

%t Do[If[IntegerQ[n/DivisorSigma[0, n]], Print[n]], {n, 1, 100000, 2}]

%t Select[Range[1,1001,2]^2,Divisible[#,DivisorSigma[0,#]]&] (* _Harvey P. Dale_, Jan 22 2012 *)

%o (PARI) is(n)=n%2&&issquare(n)&&n%numdiv(n)==0 \\ _Charles R Greathouse IV_, Apr 23 2013

%o (PARI) list(lim)=my(v=List(),f);forstep(n=1,sqrtint(lim\1),2,f=factor(n)[,2];if(n^2%prod(i=1,#f,2*f[i]+1)==0,listput(v,n^2))); Vec(v) \\ _Charles R Greathouse IV_, Apr 23 2013

%o (Python)

%o from itertools import count, islice

%o from math import prod

%o from sympy import factorint

%o def A036896_gen(): # generator of terms

%o for n in count(1,2):

%o if not (m:=n**2)%prod(e<<1|1 for e in factorint(n).values()): yield m

%o A036896_list = list(islice(A036896_gen(),40)) # _Chai Wah Wu_, Oct 04 2024

%Y Subsequence of A033950 and of A016754.

%K nonn,easy

%O 1,2

%A Simon Colton (simonco(AT)cs.york.ac.uk)