cp's OEIS Frontend

A348500 a(n) = A348494(A276086(n)).

%I A348500 #12 Nov 13 2021 22:33:43
%S A348500 1,1,1,5,1,7,1,1,1,1,1,1,1,3,1,1,1,1,1,1,1,1,1,1,1,1,1,7,1,1,1,3,5,1,
%T A348500 1,1,3,1,1,13,1,1,1,1,1,1,1,1,13,3,1,1,1,1,3,1,1,1,13,1,1,1,1,1,1,1,1,
%U A348500 3,1,1,1,1,1,1,1,1,1,1,1,1,1,19,1,1,19,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,1,1
%N A348500 a(n) = A348494(A276086(n)).
%H A348500 Antti Karttunen, <a href="/A348500/b348500.txt">Table of n, a(n) for n = 0..11550</a>
%H A348500 <a href="/index/Pri#primorialbase">Index entries for sequences related to primorial base</a>
%F A348500 a(n) = A348494(A276086(n)), where A348494(n) = A348492(n) / A003557(n).
%t A348500 Block[{b = MixedRadix[Reverse@ Prime@ Range@ 12]}, Map[GCD[Total@GCD[#1, Range[#1]], #1 Total[#2/#1 & @@@ #2]]/Apply[Times, Map[#1^(#2 - 1) & @@ # &, #2]] & @@ {#, FactorInteger[#]} &, Table[Function[k, Times @@ Power @@@ # &@ Transpose@ {Prime@ Range@ Length@ k, Reverse@ k}]@ IntegerDigits[n, b], {n, 0, 105}] ] ] (* _Michael De Vlieger_, Oct 21 2021 *)
%o A348500 (PARI)
%o A348500 \\ Needs also code from A348494:
%o A348500 A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
%o A348500 A348500(n) = A348494(A276086(n));
%Y A348500 Cf. A003557, A276086, A348492, A348494.
%Y A348500 Cf. also A348501.
%K A348500 nonn,base
%O A348500 0,4
%A A348500 _Antti Karttunen_, Oct 21 2021