|
|
|
|
1, 1, 1, 1, 1, 3, 1, 1, 2, 9, 1, 5, 1, 27, 1, 1, 1, 1, 1, 25, 3, 81, 1, 7, 4, 243, 2, 125, 1, 5, 1, 1, 9, 729, 2, 5, 1, 2187, 27, 49, 1, 25, 1, 625, 1, 6561, 1, 11, 8, 1, 81, 3125, 1, 3, 1, 343, 243, 19683, 1, 35, 1, 59049, 5, 1, 3, 125, 1, 15625, 729, 5, 1, 7, 1, 177147, 2, 78125, 4, 625, 1, 121, 2, 531441, 1, 245, 9, 1594323, 2187, 2401, 1, 21
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,6
|
|
COMMENTS
|
It appears that these and the terms of A331598 have the same prime signatures, that is, A046523(a(n)) = A046523(A331598(n)) seems to hold for all n.
|
|
LINKS
|
|
|
FORMULA
|
|
|
MATHEMATICA
|
Array[If[# == 1, 1, #2/GCD[#1, #2] & @@ {Block[{k = #, m = 0}, Times @@ Power @@@ Table[k -= m; k = DeleteCases[k, 0]; {Prime@ Length@ k, m = Min@ k}, Length@ Union@ k]] &@ Catenate[ConstantArray[PrimePi[#1], #2] & @@@ #], Function[t, Times @@ Prime@ Accumulate[If[Length@ t < 2, {0}, Join[{1}, ConstantArray[0, Length@ t - 2], {-1}]] + ReplacePart[t, Map[#1 -> #2 & @@ # &, #]]]]@ ConstantArray[0, Transpose[#][[1, -1]]] &[# /. {p_, e_} /; p > 0 :> {PrimePi@ p, e}]} &@ FactorInteger[#]] &, 90] (* Michael De Vlieger, Jan 24 2020, after JungHwan Min at A122111 *)
|
|
PROG
|
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|