|
|
|
|
2, 3, 5, 9, 7, 25, 11, 15, 35, 49, 13, 625, 17, 121, 117649, 225, 19, 1225, 23, 2401, 1771561, 169, 29, 875, 717409, 289, 55, 14641, 31, 184877, 37, 21, 4826809, 361, 36226650889, 1500625, 41, 529, 24137569, 77, 43, 143, 47, 28561, 1127357, 841, 53, 1715, 902613283, 514675673281, 47045881, 83521, 59, 3025, 8254129, 214358881, 148035889, 961, 61
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
FORMULA
|
(End)
|
|
MATHEMATICA
|
With[{b = MixedRadix[Reverse@ Prime@ Range@ 120]}, Array[Function[k, Times @@ Power @@@ # &@ Transpose@ {Prime@ Range@ Length@ k, Reverse@ k}]@ IntegerDigits[Apply[Times, Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Times @@ Prime@ Range@ PrimePi@ p, e}]], b] &, 58]] (* Michael De Vlieger, Nov 18 2019 *)
A276086[n0_] := Module[{m = 1, i = 1, n = n0, p}, While[n > 0, p = Prime[i]; m *= p^Mod[n, p]; n = Quotient[n, p]; i++]; m];
b[n_] := b[n] = Module[{pe = FactorInteger[n], p, e}, If[Length[pe] > 1, Times @@ b /@ Power @@@ pe, {{p, e}} = pe; Times @@ (Prime[Range[ PrimePi[p]]]^e)]]; b[1] = 1;
|
|
PROG
|
(PARI)
A034386(n) = prod(i=1, primepi(n), prime(i));
A276086(n) = { my(i=0, m=1, pr=1, nextpr); while((n>0), i=i+1; nextpr = prime(i)*pr; if((n%nextpr), m*=(prime(i)^((n%nextpr)/pr)); n-=(n%nextpr)); pr=nextpr); m; };
|
|
CROSSREFS
|
Cf. A034386, A108951, A117366, A276086, A324887, A324888, A324896, A329047, A342920, A344592, A346106, A346107.
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|