login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A331291 a(n) = A001222(n) * A329348(n). 3

%I #9 Jan 15 2020 18:40:27

%S 1,1,1,4,1,4,1,3,2,4,1,12,1,4,12,8,1,6,1,12,12,4,1,12,4,4,3,12,1,15,1,

%T 5,12,4,16,16,1,4,12,4,1,3,1,12,3,4,1,5,2,12,12,12,1,8,8,32,12,4,1,12,

%U 1,4,9,12,26,36,1,12,12,15,1,15,1,4,15,12,4,36,1,10,4,4,1,8,22,4,12,32,1,8,12,12,12,4

%N a(n) = A001222(n) * A329348(n).

%C Records occur at n = 1, 4, 12, 30, 35, 56, 66, 132, 145, 208, 276, 279, 469, ...

%H Antti Karttunen, <a href="/A331291/b331291.txt">Table of n, a(n) for n = 1..12571</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) = A001222(n) * A329348(n) = A001222(n) * (A331188(n) mod A117366(n)).

%F a(p) = 1 for all primes p.

%o (PARI)

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

%o A331291(n) = if(1==n, 1, my(f=factor(n), p=nextprime(1+vecmax(f[, 1]))); ((prod(i=1, #f~, A002110(primepi(f[i, 1]))^(f[i, 2]-(#f~==i)))%p)*bigomega(n)));

%Y Cf. A001222, A117366, A329348, A331188, A331289, A331290.

%K nonn,less

%O 1,4

%A _Antti Karttunen_, Jan 15 2020

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 07:30 EDT 2024. Contains 371235 sequences. (Running on oeis4.)