OFFSET
1,1
LINKS
FORMULA
From Antti Karttunen, Jul 09 2021: (Start)
(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 is A108951 *)
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;
a[n_] := A276086[b[n]];
Array[a, 100] (* Jean-François Alcover, Dec 01 2021, after _Antti Karttunen in A296086 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Mar 30 2019
STATUS
approved