%I #46 May 30 2022 13:00:42

%S 1,2,6,30,6,210,6,2310,2730,30,39270,7410,42,7590,46410,1272810,30,

%T 930930,82110,6,21111090,1230,48969690,1738215570,2310,102,144690,6,

%U 85470,29594505363092670,16770,49990710,138,7849357706190,30,300690390,20223210,1122990,37916970

%N Product of the distinct primes dividing the product of composite numbers between consecutive primes.

%C Equivalently, the largest squarefree number that divides the product of composite numbers between successive primes.

%C From _Robert G. Wilson v_, Dec 02 2020: (Start)

%C All terms greater than one are even.

%C Omega(a(n)): 0, 1, 2, 3, 2, 4, 2, 5, 5, 3, 6, 5, 3, 5, 6, 7, 3, 7, 6, 2, 8, 4, 8, 9, 5, ..., .

%C Records: 1, 2, 6, 30, 210, 2310, 2730, 39270, 46410, 1272810, 21111090, ..., (2*A354218).

%C Factored: 1, 2, 2*3, 2*3*5, 2*3*5*7, 2*3*5*7*11, 2*3*5*7*13, 2*3*5*7*11*17, 2*3*5*7*13*17, 2*3*5*7*11*19*29, ..., .

%C (End)

%H Robert G. Wilson v, <a href="/A076978/b076978.txt">Table of n, a(n) for n = 1..10000</a>

%F From _Michel Marcus_, May 29 2022: (Start)

%F a(n) = A007947(A074167(n)).

%F a(n) = A007947(A061214(n)). (End)

%e a(4) = product of prime divisors of the product of composite numbers between 7 and 11 = 2 * 3 * 5 = 30.

%e a(5)=6 because 12 is the only composite number between the 5th and the 6th primes (11 and 13) and largest squarefree divisor of 12 is 6.

%p with(numtheory): b:=proc(j) if issqrfree(j) then j else fi end: a:=proc(n) local B,BB: B:=divisors(product(i,i=ithprime(n)+1..ithprime(n+1)-1)): BB:=(seq(b(B[j]),j=1..nops(B))): max(BB); end: seq(a(n),n=1..33); # _Emeric Deutsch_, Jul 28 2006

%t f[n_] := Times @@ (First@# & /@ FactorInteger[Times @@ Range[Prime[n] + 1, Prime[n + 1] - 1]]); Array[f, 50] (* _Robert G. Wilson v_, Dec 02 2020 *)

%o (PARI) a(n) = my(p=1); forcomposite(c=prime(n), prime(n+1), p*=c); factorback(factorint(p)[, 1]); \\ _Michel Marcus_, May 29 2022

%o (Python) from sympy import sieve as p, primefactors

%o def A076978(n):

%o result = 1

%o for composites in range(p[n]+1, p[n+1]):

%o for primefactor in primefactors(composites):

%o if result % primefactor != 0: result *= primefactor

%o return result # _Karl-Heinz Hofmann_, May 30 2022

%Y Cf. A007947, A061214, A074167, A354217, A354218.

%K nonn

%O 1,2

%A _Amarnath Murthy_, Oct 23 2002

%E More terms from _Emeric Deutsch_, Jul 28 2006

%E More terms from _Robert G. Wilson v_, Dec 02 2020

%E Entry revised by _N. J. A. Sloane_, Dec 02 2020