|
%I #7 Oct 16 2019 13:18:58
%S 1,1,1,2,1,4,3,2,5,1,5,4,6,3,3,7,7,4,7,5,16,6,12,27,21,35,28,23,31,28,
%T 2,2,4,5,4,5,4,10,9,2,11,6,7,10,12,7,30,6,10,15,14,7,23,37,26,32,28,
%U 33,24,28,8,3,17,11,3,5,6,11,7,12,30,21,28,15,28,11,24,30,14,16,43,17,52,26,19,29,27,33,46,27,12,15,28,28,24,27,11,20,16,20
%N Maximal digit value in primorial base expansion of A276087(n): a(n) = A328114(A276087(n)).
%H Antti Karttunen, <a href="/A328394/b328394.txt">Table of n, a(n) for n = 0..32768</a>
%H <a href="/index/Pri#primorialbase">Index entries for sequences related to primorial base</a>
%F a(n) = A328389(A276086(n)) = A328114(A276087(n)) = A051903(A328403(n)).
%t Block[{b = MixedRadix[Reverse@ Prime@ Range@ 120], f}, f[n_] := Times @@ Power @@@ # &@ Transpose@ {Prime@ Range@ Length@ #, Reverse@ #} &@ IntegerDigits[n, b]; Array[Max@ IntegerDigits[Nest[f, #, 2], b] &, 100, 0]] (* _Michael De Vlieger_, Oct 15 2019 *)
%o (PARI)
%o A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
%o A276087(n) = A276086(A276086(n));
%o A328114(n) = { my(s=0, p=2); while(n, s = max(s,(n%p)); n = n\p; p = nextprime(1+p)); (s); };
%o A328394(n) = A328114(A276087(n));
%Y Cf. A051903, A276086, A276087, A328114, A328322, A328389, A328403.
%K nonn
%O 0,4
%A _Antti Karttunen_, Oct 15 2019
|
