cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

User: Jean-Luc Chabert

Jean-Luc Chabert's wiki page.

Jean-Luc Chabert has authored 2 sequences.

A212429 a(n) is the LCM of denominators of polynomials of degree n which are integer-valued on primes together with their first divided differences.

Original entry on oeis.org

1, 1, 2, 4, 48, 96, 1152, 2304, 276480, 552960, 6635520, 13271040, 33443020800, 66886041600, 802632499200, 1605264998400, 385263599616000, 770527199232000, 194172854206464000, 388345708412928000, 512616335105064960000, 1025232670210129920000
Offset: 1

Author

Jean-Luc Chabert, Jun 21 2012

Keywords

Comments

a(n) is also the n-th Bhargava's factorial n_P^{{1}} of the set P of primes with respect to the first divided difference.

Examples

			a(5) = 48 because f(x) = (x-1)(x-2)(x-3)(x-5)(x-7)/48 satisfies f(p) and (f(p)-f(q))/(p-q) are integers for all primes p,q.
		

Crossrefs

Cf. A053657.

Programs

  • Maple
    a:= proc(n) local i, p, wp, r;
          r:=1;
          for i do p:= ithprime(i);
            wp:= p^(w(p,n-1));
            if wp=1 then break fi;
            r:= r*wp
          od; r
        end:
    w:= proc(p, n) local d, k, r;
          r:= 0;
          for k from 0 do d:= floor(n/((p-1)*p^k));
            if d=0 then break fi;
            r:= r+d;
          od;
          r -t(n,p)
        end:
    t:= proc(n, p) local h, q;
          q:= n/(p-1);
          for h from 0 while q>= p^h do od; h
        end:
    seq (a(n), n=1..30);  # Alois P. Heinz, Jun 25 2012
  • Mathematica
    a[n_] := Module[{i, p, wp, r}, r = 1; For[i = 1, True, i++, p = Prime[i]; wp = p^w[p, n - 1]; If[wp == 1, Break[]]; r = r*wp]; r];
    w[p_, n_] := Module[{d, k, r}, r = 0; For[k = 0, True, k++, d = Floor[n/((p - 1)*p^k)]; If[d == 0, Break[]]; r = r + d]; r - t[n, p]];
    t[n_, p_] := Module[{h, q}, q = n/(p - 1); For[h = 0, q >= p^h , h++]; h];
    a /@ Range[1, 30] (* Jean-François Alcover, Oct 14 2019, after Alois P. Heinz *)

Formula

a(n) = Prod_{p prime} p^w_p(n-1) where w_p(n) = Sum_{k>=0} floor(n / ((p-1)*p^k)) - t_{p,n} and p^(t_{p,n}-1) <= n/(p-1) < p^t_{p,n}.

A053657 a(n) = Product_{p prime} p^{ Sum_{k>=0} floor[(n-1)/((p-1)p^k)]}.

Original entry on oeis.org

1, 2, 24, 48, 5760, 11520, 2903040, 5806080, 1393459200, 2786918400, 367873228800, 735746457600, 24103053950976000, 48206107901952000, 578473294823424000, 1156946589646848000, 9440684171518279680000, 18881368343036559360000, 271211974879377138647040000
Offset: 1

Author

Jean-Luc Chabert, Feb 16 2000

Keywords

Comments

LCM of denominators of the coefficients of x^n*z^k in {-log(1-x)/x}^z as k=0..n, as described by triangle A075264.
Denominators of integer-valued polynomials on prime numbers (with degree n): 1/a(n) is a generator of the ideal formed by the leading coefficients of integer-valued polynomials on prime numbers with degree less than or equal to n.
Also the least common multiple of the orders of all finite subgroups of GL_n(Q) [Minkowski]. Schur's notation for the sequence is M_n = a(n+1). - Martin Lorenz (lorenz(AT)math.temple.edu), May 18 2005
This sequence also occurs in algebraic topology where it gives the denominators of the Laurent polynomials forming a regular basis for K*K, the hopf algebroid of stable cooperations for complex K-theory. Several different equivalent formulas for the terms of the sequence occur in the literature. An early reference is K. Johnson, Illinois J. Math. 28(1), 1984, pp.57-63 where it occurs in lines 1-5, page 58. A summary of some of the other formulas is given in the appendix to K. Johnson, Jour. of K-theory 2(1), 2008, 123-145. - Keith Johnson (johnson(AT)mscs.dal.ca), Nov 03 2008
a(n) is divisible by n!, by Legendre's formula for the highest power of a prime that divides n!. Also, a(n) is divisible by (n+1)! if and only if n+1 is not prime. - Jonathan Sondow, Jul 23 2009
Triangle A163940 is related to the divergent series 1^m*1! - 2^m*2! + 3^m*3! - 4^m*4! + ... for m =>-1. The left hand columns of this triangle can be generated with the MC polynomials, see A163972. The Minkowski numbers appear in the denominators of these polynomials. - Johannes W. Meijer, Oct 16 2009
Unsigned Stirling numbers of the first kind as [s + k, k] (Karamata's notation) where k = {0, 1, 2, ...} and s is in general complex results in Pochhammer[s,k]*(integer coefficient polynomial of (k-1) degree in s) / M[k], where M[k] is the least common multiple of the orders of all finite groups of n x n-matrices over rational numbers (Minkowiski's theorem) which is sequence A053657. - Lorenz H. Menke, Jr., Feb 02 2010
From Peter Bala, Feb 21 2011: (Start)
Given a subset S of the integers Z, Bhargava has shown how to associate with S a generalized factorial function, denoted n!_S, which shares many properties of the classical factorial function n!.
The present sequence is the generalized factorial function n!S associated with the set of primes S = {2,3,5,7,...}. The associated generalized exponential function E(x) = Sum{n>=1} x^(n-1)/a(n) vanishes at x = -2: i.e. Sum_{n>=1} (-2)^n/a(n) = 0.
For the table of associated generalized binomial coefficients n!_S/(k!_S*(n-k)!_S) see A186430.
This sequence is related to the Bernoulli polynomials in two ways [Chabert and Cahen]:
(1) a(n) = (n-1)!*A001898(n-1).
(2) (t/(exp(t)-1))^x = sum {n = 0..inf} P(n,x)*t^n/a(n+1),
where the P(n,x) are primitive polynomials in the ring Z[x].
If p_1,...,p_n are any n primes then the product of their pairwise differences Product_{i
(End)
LCM of denominators of the coefficients of S(m+n-1,m) as polynomial in m of degree 2*(n-1), as described by triangle A202339. - Vladimir Shevelev, Dec 17 2011
Sometimes called "Minkowski numbers" (e.g., by Guralnick and Lorenz), after the German mathematician Hermann Minkowski (1864-1909). - Amiram Eldar, Aug 24 2024

Examples

			a(7)=24^3*Product_{i=1..3} A202318(i)=24^3*1*10*21=2903040. - _Vladimir Shevelev_, Dec 17 2011
		

References

  • Jean-Luc Chabert, Scott T. Chapman, and William W. Smith, A basis for the ring of polynomials integer-valued on prime numbers, in: Daniel Anderson (ed.), Factorization in integral domains, Lecture Notes in Pure and Appl. Math. 189, Dekker, New York, 1997.

Crossrefs

a(n) = n!*A163176(n). - Jonathan Sondow, Jul 23 2009
Cf. A202318.
Appears in A163972. - Johannes W. Meijer, Oct 16 2009

Programs

  • Maple
    A053657 := proc(n) local P,p,q,s,r;
    P := select(isprime,[$2..n]); r:=1;
    for p in P do s := 0; q := p-1;
    do if q > (n-1) then break fi;
    s := s + iquo(n-1,q); q := q*p; od;
    r := r * p^s; od; r end: # Peter Luschny, Jul 26 2009
    ser := series((y/(exp(y)-1))^x, y, 20): a := n -> denom(coeff(ser, y, n-1)):
    seq(a(n), n=1..19); # Peter Luschny, May 13 2019
  • Mathematica
    m = 16; s = Expand[Normal[Series[(-Log[1-x]/x)^z, {x, 0, m}]]];
    a[n_, k_] := Denominator[ Coefficient[s, x^n*z^k]];
    Prepend[Apply[LCM, Table[a[n,k], {n,m}, {k,n}], {1}], 1]
    (* Jean-François Alcover, May 31 2011 *)
    a[n_] := Product[p^Sum[Floor[(n-1)/((p-1) p^k)], {k, 0, n}], {p, Prime[ Range[n] ]}]; Array[a, 30] (* Jean-François Alcover, Nov 22 2016 *)
  • PARI
    {a(n)=local(X=x+x^2*O(x^n),D);D=1;for(j=0,n-1,D=lcm(D,denominator( polcoeff(polcoeff((-log(1-X)/x)^z+z*O(z^j),j,z),n-1,x))));return(D)} /* Paul D. Hanna, Jun 27 2005 */
    
  • PARI
    {a(n)=prod(i=1,#factor(n!)~,prime(i)^sum(k=0,#binary(n), floor((n-1)/((prime(i)-1)*prime(i)^k))))} /* Paul D. Hanna, Jun 27 2005 */
    
  • PARI
    S(n, p) = {
         my(acc = 0, tmp = p-1);
         while (tmp < n, acc += floor((n-1)/tmp); tmp *= p);
         return(acc);
    };
    a(n) = {
         my(rv = 1);
         forprime(p = 2, n, rv *= p^S(n,p));
         return(rv);
    };
    vector(17, i, a(i))  \\ Gheorghe Coserea, Aug 24 2015

Formula

a(2n) = 2*a(2n-1). - Jonathan Sondow, Jul 23 2009
a(2*n+1) = 24^n * Product_{i=1..n} A202318(i). - Vladimir Shevelev, Dec 17 2011
For n>=0, A007814(a(n+1)) = n+A007814(n!). - Vladimir Shevelev, Dec 28 2011
a(n) = denominator([y^(n-1)] (y/(exp(y)-1))^x). - Peter Luschny, May 13 2019
Sum_{n>=1} 1/a(n) = A346046. - Amiram Eldar, Jul 02 2023

Extensions

More terms from Paul D. Hanna, Jun 27 2005