login
a(n) = A108951(n) * A276086(A108951(n)).
3

%I #5 Mar 31 2019 00:17:39

%S 2,6,30,36,210,300,2310,120,1260,2940,30030,15000,510510,50820,

%T 21176820,3600,9699690,88200,223092870,288120,2232166860,780780,

%U 6469693230,42000,645668100,17357340,11880,12298440,200560490130,66555720,7420738134810,672,66899572740,368588220,228227900600700,216090000,304250263527210

%N a(n) = A108951(n) * A276086(A108951(n)).

%H Antti Karttunen, <a href="/A324887/b324887.txt">Table of n, a(n) for n = 1..210</a>

%H <a href="/index/Pri#prime_indices">Index entries for sequences computed from indices in prime factorization</a>

%H <a href="/index/Pri#primorialbase">Index entries for sequences related to primorial base</a>

%H <a href="/index/Pri#primorial_numbers">Index entries for sequences related to primorial numbers</a>

%F a(n) = A324580(A108951(n)) = A108951(n) * A324886(n).

%o (PARI)

%o A034386(n) = prod(i=1, primepi(n), prime(i));

%o A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A034386(f[i, 1])^f[i, 2]) }; \\ From A108951

%o 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; };

%o A324886(n) = A276086(A108951(n));

%o A324887(n) = (A324886(n)*A108951(n));

%Y Cf. A034386, A108951, A276086, A324580, A324886, A324888.

%Y Permutation of A324577.

%K nonn

%O 1,1

%A _Antti Karttunen_, Mar 30 2019