The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
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!)
A352447 Numbers k such that BarnesG(k) is divisible by Gamma(k). 0
1, 2, 7, 9, 10, 11, 13, 15, 16, 17, 19, 21, 22, 23, 25, 26, 27, 28, 29, 31, 33, 34, 35, 36, 37, 39, 40, 41, 43, 45, 46, 47, 49, 50, 51, 52, 53, 55, 56, 57, 58, 59, 61, 63, 64, 65, 66, 67, 69, 70, 71, 73, 75, 76, 77, 78, 79, 81, 82, 83, 85, 86, 87, 88, 89, 91, 92, 93, 94, 95 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
These are k such that G(k)/Gamma(k) = 1!*2!*3!*...*(k-2)!/(k-1)! = 1!*2!*3!*...*(k-3)!/(k-1) are integer. Let k=1+c, so require 1!*2!*3!*...*(c-2)!/c to be integer. If c is composite, take any factorization of c=c_1*c_2 with 2<=c_1<=c_2<=c/2; then matching factors exist in the product 1!*2!*3!*...*(c-2)! that cancel this factor [either c_1! and c_2! if c_1 <> c_2, or c_1! and (c_1+1)! if c_1=c_2 and c-2 >= c_1+1]. If c is prime, none of the 1!*2!*..(c-2)! contains a factor matching that prime. So this shows that the sequence is (apart from offset and at c=4) the same as A079696. - R. J. Mathar, Mar 25 2022
LINKS
Eric Weisstein's World of Mathematics, Barnes G-Function.
Eric Weisstein's World of Mathematics, Divisible.
Eric Weisstein's World of Mathematics, Gamma Function.
Eric Weisstein's World of Mathematics, Superfactorial.
FORMULA
Conjecture: a(n) = A079696(n-1), n>1. - R. J. Mathar, Mar 20 2022
EXAMPLE
BarnesG(7) = 34560, Gamma(7) = 720, 34560 is divisible by 720, so 7 is in this sequence.
MATHEMATICA
Table[If[Divisible[BarnesG[k], Gamma[k]], k, Nothing], {k, 115}]
PROG
(Python)
from itertools import count, islice
from collections import Counter
from sympy import factorint
def A352447_gen(): # generator of terms
yield 1
a = Counter()
for k in count(2):
b = Counter(factorint(k-1))
if all(b[p] <= a[p] for p in b):
yield k
a += b
A352447_list = list(islice(A352447_gen(), 100)) # Chai Wah Wu, Mar 17 2022
CROSSREFS
Sequence in context: A086725 A019950 A201892 * A323528 A073074 A034796
KEYWORD
nonn
AUTHOR
STATUS
approved

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 May 19 07:05 EDT 2024. Contains 372666 sequences. (Running on oeis4.)