login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A141573
Prime powers of pairs of the infinite list generated by the flat list of the factorizations of all integers.
0
1, 9, 32, 8, 49, 4, 27, 32, 121, 8, 169, 343, 25, 4, 131072, 8, 1162261467, 4, 125, 49, 895430243255237372246531, 4, 8, 3125, 8192, 27, 9, 128, 841, 243, 961, 4, 4, 177147, 131072, 78125, 4, 27, 1369, 6859, 169, 4, 45474735088646411895751953125, 8
OFFSET
1,2
COMMENTS
Factorize n to its full extent into the list 1, 2, 3, 2*2, 5, 2*3, 7, 2*2*2, 3*3, 2*5, 11, 2*2*3,..
Remove delimiting commas and the multiplication signs and build new powers by pairing consecutive pairs of numbers of this list, 1^2, 3^2, 2^5, 2^3, 7^2, 2^2, 3^3, 2^5, 11^2, 2^3.. to construct the sequence.
MAPLE
pflat := proc(nmax) local a, ifs, n, p, c ; a := [1] ; for n from 2 to nmax do ifs := ifactors(n)[2] ; for p in ifs do q := op(1, p) ; for c from 1 to op(2, p) do a := [op(a), q] ; od: od: od: a ; end: pL := pflat(90) : for n from 1 to nops(pL)-2 by 2 do printf("%d, ", op(n, pL)^op(n+1, pL)) ; od: # R. J. Mathar, Aug 21 2008
CROSSREFS
Sequence in context: A159279 A048547 A351713 * A075433 A018833 A130510
KEYWORD
nonn,less
AUTHOR
EXTENSIONS
Edited, corrected and extended by R. J. Mathar, Aug 21 2008
STATUS
approved