|
%I #5 Jul 09 2021 22:16:37
%S 1,2,6,2,30,6,210,8,36,30,2310,6,30030,210,30,8,510510,36,9699690,30,
%T 210,2310,223092870,36,240,30030,216,210,6469693230,240,200560490130,
%U 32,2310,510510,2520,36,7420738134810,9699690,30030,240,304250263527210,2520,13082761331670030,2310,240,223092870,614889782588491410,36,32550
%N a(n) = A276085(A328571(A108951(n))).
%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) = A276085(A346091(n)) = A276085(A328571(A108951(n))).
%F a(n) = A108951(n) - A346092(n).
%o (PARI)
%o A002110(n) = prod(i=1,n,prime(i));
%o A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A002110(primepi(f[i, 1]))^f[i, 2]) };
%o A276085(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*A002110(primepi(f[k, 1])-1)); };
%o A328571(n) = { my(m=1, p=2); while(n, m *= (p^!!(n%p)); n = n\p; p = nextprime(1+p)); (m); };
%o A346093(n) = A276085(A328571(A108951(n)));
%Y Cf. A002110, A108951, A276085, A328571, A346091, A346092.
%K nonn
%O 1,2
%A _Antti Karttunen_, Jul 09 2021
|
