OFFSET
0,2
COMMENTS
Partial sums of A382720. - James C. McMahon, Aug 15 2025
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Hsien-Kuei Hwang, Svante Janson, and Tsung-Hsi Tsai, Periodic minimum in the count of binomial coefficients not divisible by a prime, arXiv:2408.06817 [math.NT], 2024.
MATHEMATICA
a[n_]:=(n^2+3n+2)/2-Count[Mod[Flatten[Table[Binomial[m, k], {m, 0, n}, {k, 0, m}]] , 7], 0]; Array[a, 69, 0] (* James C. McMahon, Aug 15 2025 *)
Table[Sum[Times@@(IntegerDigits[m, 7]+1), {m, 0, n}], {n, 0, 68}] (* Vincenzo Librandi, Feb 17 2026 *)
PROG
(Python)
from math import prod
from gmpy2 import digits
def A382726(n): return sum(prod(int(d)+1 for d in digits(m, 7)) for m in range(n+1)) # Chai Wah Wu, Aug 10 2025
(Python)
from math import prod
from gmpy2 import digits
def A382726(n):
d = list(map(lambda x:int(x)+1, digits(n+1, 7)[::-1]))
return sum((b-1)*prod(d[a:])*28**a for a, b in enumerate(d))>>1 # Chai Wah Wu, Aug 13 2025
(Magma) [&+[ &*[ d+1 : d in (m eq 0 select [0] else IntegerToSequence(m, 7))]: m in [0..n]]: n in [0..68] ]; // Vincenzo Librandi, Feb 17 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Apr 23 2025
STATUS
approved
