A318256 -id:A318256 - cp's OEIS frontend

A306745 Even bisection of A318256.

1, 2, 6, 6, 10, 6, 210, 2, 30, 210, 110, 30, 546, 14, 30, 462, 1190, 6, 51870, 70, 2310, 2310, 322, 210, 6630, 286, 330, 798, 290, 30, 930930, 1430, 1122, 82110, 70, 330, 21111090, 38038, 390, 210, 536690, 39270, 7489482, 10010, 62790, 4470510, 687610, 462

Offset: 0

Jean-François Alcover and Paul Curtz, Mar 07 2019

nonn

The ratio sequence of a(n)/A006955(n) begins: 1,1,1,1,1,1,1,1,1,5,1,5,1,7,1,1,7,1,1,35, ...

This divisibility by A006955 was confirmed by Peter Luschny in the SeqFan mailing list.

Jean-François Alcover, Table of n, a(n) for n = 0..1000

Cf. A006955, A318256.

sfk[n_] := Times @@ FactorInteger[n][[All, 1]];
a[n_] := (BernoulliB[2 n, x] // Together // Denominator)/sfk[2 n + 1];
Table[a[n], {n, 0, 60}]

a(n) = (denominator of B(2n,x)) / (the squarefree kernel of 2n+1), where B(2n,x) is the 2n-th Bernoulli polynomial.

A324370 Product of all primes p not dividing n such that the sum of the base-p digits of n is at least p, or 1 if no such prime exists.

Original entry on oeis.org

1, 1, 2, 1, 6, 1, 6, 3, 10, 1, 6, 1, 210, 15, 2, 3, 30, 5, 210, 21, 110, 15, 30, 5, 546, 21, 14, 1, 30, 1, 462, 231, 1190, 105, 6, 1, 51870, 1365, 70, 21, 2310, 55, 2310, 105, 322, 105, 210, 35, 6630, 663, 286, 33, 330, 55, 798, 57, 290, 15, 30, 1, 930930, 15015, 1430, 2145, 1122, 85, 82110, 2415, 70, 3, 330, 55, 21111090, 285285

Offset: 1

Bernd C. Kellner and Jonathan Sondow, Feb 24 2019

The product is finite, as the sum of the base-p digits of n is n if p > n.
a(198) = 2465 is the only term below 10^6 that is a Carmichael number (A002997).
It appears that a(n)=1 if and only if n is in A094960. - Robert Israel, Mar 30 2020
It turns out that a(n) equals the denominator of the first derivative of the Bernoulli polynomial B(n,x). So a(n)=1 if and only if n is in A094960, also impyling that n+1 is prime. A324370 is also involved in such formulas regarding higher derivatives. See Kellner 2023. - Bernd C. Kellner, Oct 12 2023

			For p = 2, 3, and 5, the sum of the base p digits of 7 is 1+1+1 = 3 >= 2, 2+1 = 3 >= 3, and 1+2 = 3 < 5, respectively, so a(7) = 2*3 = 6.

Robert Israel, Table of n, a(n) for n = 1..5000
Bernd C. Kellner, On a product of certain primes, J. Number Theory, 179 (2017), 126-141; arXiv:1705.04303 [math.NT], 2017.
Bernd C. Kellner, On the finiteness of Bernoulli polynomials whose derivative has only integral coefficients, J. Integer Seq. 27 (2024), Article 24.2.8, 11 pp.; arXiv:2310.01325 [math.NT], 2023.
Bernd C. Kellner and Jonathan Sondow, Power-Sum Denominators, Amer. Math. Monthly, 124 (2017), 695-709; arXiv:1705.03857 [math.NT], 2017.
Bernd C. Kellner and Jonathan Sondow, On Carmichael and polygonal numbers, Bernoulli polynomials, and sums of base-p digits, Integers 21 (2021), #A52, 21 pp.; arXiv:1902.10672 [math.NT], 2019.

Cf. A002997, A094960, A144845, A195441, A324315, A324316, A324317, A324318, A324319, A324320, A324369, A324371, A324404, A324405, A366168, A366169, A366186, A366187, A366188.

N:= 100: # for a(1)..a(N)
V:= Vector(N,1):
p:= 1:
for iter from 1 do
   p:= nextprime(p);
   if p >= N then break fi;
   for n from p+1 to N do
     if n mod p <> 0 and convert(convert(n,base,p),`+`)>= p then
       V[n]:= V[n]*p
     fi
od od:
convert(V,list); # Robert Israel, Mar 30 2020
# Alternatively, note that this formula is suggesting offset 0 and a(0) = 1:
seq(denom(diff(bernoulli(n, x), x)), n = 1..51); # Peter Luschny, Oct 13 2023

SD[n_, p_] := If[n < 1 || p < 2, 0, Plus @@ IntegerDigits[n, p]];
DD2[n_] := Times @@ Select[Prime[Range[PrimePi[(n+1)/(2+Mod[n+1, 2])]]], !Divisible[n, #] && SD[n, #] >= # &];
Table[DD2[n], {n, 1, 100}]
(* From Bernd C. Kellner, Oct 12 2023 (Start) *)
(* Denominator of first derivative of BP *)
k = 1; Table[Denominator[Together[D[BernoulliB[n, x], {x, k}]]], {n, 1, 100}]
(* End *)

from math import prod
from sympy.ntheory import digits
from sympy import primefactors, primerange
def a(n):
    nonpf = set(primerange(1, n+1)) - set(primefactors(n))
    return prod(p for p in nonpf if sum(digits(n, p)[1:]) >= p)
print([a(n) for n in range(1, 75)]) # Michael S. Branicky, Jul 03 2022

a(n) * A324369(n) = A195441(n-1) = denominator(Bernoulli_n(x) - Bernoulli_n).
a(n) * A324369(n) * A324371(n) = A144845(n-1) = denominator(Bernoulli_{n-1}(x)).
a(n+1) = A195441(n)/A324369(n+1) = A144845(n)/A007947(n+1) = A318256(n). Essentially the same as A318256. - Peter Luschny, Mar 05 2019
From Bernd C. Kellner, Oct 12 2023: (Start)
a(n) = denominator(Bernoulli_n(x)').
k-th derivative: let (n)_m be the falling factorial.
For n > k, a(n-k+1)/gcd(a(n-k+1), (n)_{k-1}) = denominator(Bernoulli_n(x)^(k)). Otherwise, the denominator equals 1. (End)

A195441 a(n) = denominator(Bernoulli_{n+1}(x) - Bernoulli_{n+1}).

Original entry on oeis.org

1, 1, 2, 1, 6, 2, 6, 3, 10, 2, 6, 2, 210, 30, 6, 3, 30, 10, 210, 42, 330, 30, 30, 30, 546, 42, 14, 2, 30, 2, 462, 231, 3570, 210, 6, 2, 51870, 2730, 210, 42, 2310, 330, 2310, 210, 4830, 210, 210, 210, 6630, 1326, 858, 66, 330, 110, 798, 114, 870, 30, 30, 6

Offset: 0

Peter Luschny, Sep 18 2011

nonn

If s(n) is the smallest number such that s(n)*(1^n + 2^n + ... + x^n) is a polynomial in x with integer coefficients then a(n)=s(n)/(n+1) (see A064538).
a(n) is squarefree, by the von Staudt-Clausen theorem on the denominators of Bernoulli numbers. - Kieren MacMillan and Jonathan Sondow, Nov 20 2015
Kellner and Sondow give a detailed analysis of this sequence and provide a simple way to compute the terms without using Bernoulli polynomials and numbers. They prove that a(n) is the product of the primes less than or equal to (n+2)/(2+(n mod 2)) such that the sum of digits of n+1 in base p is at least p. - Peter Luschny, May 14 2017
The equation a(n-1) = denominator(Bernoulli_n(x) - Bernoulli_n) = rad(n+1) has only finitely many solutions, where rad(n) = A007947(n) is the radical of n. It is conjectured that S = {3, 5, 8, 9, 11, 27, 29, 35, 59} is the full set of all such solutions. Note that (S\{8})+1 joined with {1,2} equals A094960. More precisely, the set S implies the finite sequence of A094960. See Kellner 2023. - Bernd C. Kellner, Oct 18 2023
As was observed in the example section of A318256: denominator(B_n(x)) = rad(n+1) if n is in {0, 1, 3, 5, 9, 11, 27, 29, 35, 59} = {A094960(n) - 1: 1 <= n <= 10}. - Peter Luschny, Oct 18 2023

Peter Luschny, Table of n, a(n) for n = 0..10000 (terms 0..1000 from G. C. Greubel)
Olivier Bordellès, Florian Luca, Pieter Moree, and Igor E. Shparlinski, Denominators of Bernoulli polynomials, Mathematika 64 (2018), 519-541.
Harald Hofstätter, Denominators of coefficients of the Baker-Campbell-Hausdorff series, arXiv:2010.03440 [math.NT], 2020. Mentions this sequence.
Bernd C. Kellner, On a product of certain primes, J. Number Theory, 179 (2017), 126-141; arXiv:1705.04303 [math.NT], 2017.
Bernd C. Kellner, On the finiteness of Bernoulli polynomials whose derivative has only integral coefficients, J. Integer Seq. 27 (2024), Article 24.2.8, 11 pp.; arXiv:2310.01325 [math.NT], 2023.
Bernd C. Kellner and Jonathan Sondow, Power-Sum Denominators, Amer. Math. Monthly, 124 (2017), 695-709; arXiv:1705.03857 [math.NT], 2017.
Bernd C. Kellner and Jonathan Sondow, The denominators of power sums of arithmetic progressions, Integers 18 (2018), #A95, 17 pp.; arXiv:1705.05331 [math.NT], 2017.
Bernd C. Kellner and Jonathan Sondow, On Carmichael and polygonal numbers, Bernoulli polynomials, and sums of base-p digits, Integers 21 (2021), #A52, 21 pp.; arXiv:1902.10672 [math.NT], 2019.

Cf. A002997, A064538, A094960, A144845, A286516, A286762, A286763, A318256, A324369, A324370, A324371.

using Nemo, Primes
function A195441(n::Int)
    n < 4 && return ZZ([1,1,2,1][n+1])
    P = primes(2, div(n+2, 2+n%2))
    prod([ZZ(p) for p in P if p <= sum(digits(n+1, base=p))])
end
println([A195441(n) for n in 0:59]) # Peter Luschny, May 14 2017

A195441 := n -> denom(bernoulli(n+1, x)-bernoulli(n+1)):
seq(A195441(i),i=0..59);
# Formula of Kellner and Sondow:
a := proc(n) local s; s := (p,n) -> add(i,i=convert(n,base,p));
select(isprime,[$2..(n+2)/(2+irem(n,2))]); mul(i,i=select(p->s(p,n+1)>=p,%)) end: seq(a(n), n=0..59); # Peter Luschny, May 14 2017

a[n_] := Denominator[Together[(BernoulliB[n + 1, x] - BernoulliB[n + 1])]]; Table[a[n], {n, 0, 59}] (* Jonathan Sondow, Nov 20 2015 *)
SD[n_, p_] := If[n < 1 || p < 2, 0, Plus @@ IntegerDigits[n, p]]; DD[n_] := Times @@ Select[Prime[Range[PrimePi[(n+2)/(2+Mod[n, 2])]]], SD[n+1, #] >= # &]; Table[DD[n], {n, 0, 59}] (* Bernd C. Kellner, Oct 18 2023 *)

a(n) = {my(vp = Vec(bernpol(n+1, x)-bernfrac(n+1))); lcm(vector(#vp, k, denominator(vp[k])));} \\ Michel Marcus, Feb 08 2016

from math import prod
from sympy.ntheory.factor_ import primerange, digits
def A195441(n): return prod(p for p in primerange((n+2)//(2|n&1)+1) if sum(digits(n+1,p)[1:])>=p) # Chai Wah Wu, Oct 04 2023

A195441 = lambda n: mul([p for p in (2..(n+2)//(2+n%2)) if is_prime(p) and sum((n+1).digits(base=p))>=p])
print([A195441(n) for n in (0..59)]) # Peter Luschny, May 14 2017

a(n) = A064538(n)/(n+1). - Jonathan Sondow, Nov 12 2015
A001221(a(n)) = A001222(a(n)). - Kieren MacMillan and Jonathan Sondow, Nov 20 2015
a(2*n)/a(2*n+1) = A286516(n+1). - Bernd C. Kellner and Jonathan Sondow, May 24 2017
a(n) = A007947(A338025(n+1)). - Harald Hofstätter, Oct 10 2020
From Bernd C. Kellner, Oct 18 2023: (Start)
Note that the formulas here are shifted in index by 1 due to the definition of a(n) using index n+1!
a(n) = A324369(n+1) * A324370(n+1).
a(n) = A144845(n) / A324371(n+1).
a(n-1) = lcm(a(n), rad(n+1)), if n >= 3 is odd.
If n+1 is composite, then rad(n+1) divides a(n-1).
If m is a Carmichael number (A002997), then m divides both a(m-1) and a(m-2).
See papers of Kellner and Kellner & Sondow. (End)

Definition simplified by Jonathan Sondow, Nov 20 2015

A341109 a(n) = denominator(p(n, x)) / (n!denominator(bernoulli(n, x))), where p(n, x) = Sum_{k=0..n} E2(n, k)binomial(x + k, 2*n) / Product_{j=1..n} (j - x) and E2(n, k) are the second-order Eulerian numbers A201637.

Original entry on oeis.org

1, 1, 2, 4, 8, 16, 96, 192, 1152, 768, 1536, 3072, 18432, 36864, 221184, 147456, 884736, 1769472, 10616832, 21233664, 637009920, 424673280, 2548039680, 5096079360, 152882380800, 61152952320, 366917713920, 81537269760, 163074539520, 326149079040, 1956894474240

Offset: 0

Peter Luschny, Feb 06 2021

nonn

The challenge is to characterize the sequence purely arithmetically, i.e., without reference to the Eulerian numbers or the Bernoulli polynomials.

András Bazsó and István Mező, On the coefficients of power sums of arithmetic progressions, J. Number Th., 153 (2015), 117-123.
J.-L. Chabert and P.-J. Cahen, Old problems and new questions around integer-valued polynomials and factorial sequences In: J. W. Brewer, S. Glaz, W. J. Heinzer, B. M. Olberding (eds), Multiplicative Ideal Theory in Commutative Algebra. Springer, Boston, MA., 2006.
Bernd C. Kellner and Jonathan Sondow, Power-Sum Denominators, arXiv:1705.03857 [math.NT] 2017, Amer. Math. Monthly.

Cf. A100655, A053657 (Minkowski), A341107, A341108, A318256, A144845, A163176, A201637 (Eulerian2), A036689, A324370, A007947, A324369, A195441.

Epoly := proc(n, x) add(combinat:-eulerian2(n, k)*binomial(x+k, 2*n), k = 0..n) / mul(j-x, j = 1..n): simplify(expand(%)) end:
seq(denom(Epoly(n, x)) / (n!*denom(bernoulli(n, x))), n = 0..30);

A053657[n_] := Product[p^Sum[Floor[(n-1)/((p-1) p^k)], {k,0,n}],{p, Prime[Range[n]]}];
A144845[n_] := Denominator[Together[BernoulliB[n, x]]];
A163176[n_] := A053657[n] / n!;
Table[(n + 1) A163176[n + 1] / A144845[n], {n, 0, 30}]

def A341109(n): # uses[A341108, A318256]
    return A341108(n)//A318256(n)
print([A341109(n) for n in (0..30)])

a(n) = A053657(n+1)/(n!*A144845(n)).
a(n) = (n+1)*A163176(n+1)/A144845(n).
a(n) = A341108(n)/A318256(n).
a(n) = A341107(n)*A324369(n+1).
a(n) = A341108(n)/A324370(n+1).
a(n) = A341108(n)*A007947(n+1)/A144845(n).
a(n) = A341108(n)*A324369(n+1)/A195441(n).
prime(n) divides a(k) for k >= A036689(n).
2^(n-1) divides exactly a(n) for n >= 2.

cp's OEIS Frontend

A306745 Even bisection of A318256.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A324370 Product of all primes p not dividing n such that the sum of the base-p digits of n is at least p, or 1 if no such prime exists.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Python

Formula

A195441 a(n) = denominator(Bernoulli_{n+1}(x) - Bernoulli_{n+1}).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Julia

Maple

Mathematica

PARI

Python

Sage

Formula

Extensions

A341109 a(n) = denominator(p(n, x)) / (n!*denominator(bernoulli(n, x))), where p(n, x) = Sum_{k=0..n} E2(n, k)*binomial(x + k, 2*n) / Product_{j=1..n} (j - x) and E2(n, k) are the second-order Eulerian numbers A201637.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Sage

Formula

A319084 Numbers k such that the denominator of the Bernoulli polynomial B_k(x) is the squarefree kernel of k+1.

Views

Author

Keywords

Comments

Crossrefs

Formula

A341109 a(n) = denominator(p(n, x)) / (n!denominator(bernoulli(n, x))), where p(n, x) = Sum_{k=0..n} E2(n, k)binomial(x + k, 2*n) / Product_{j=1..n} (j - x) and E2(n, k) are the second-order Eulerian numbers A201637.