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.

Showing 1-10 of 10 results.

A003961 Completely multiplicative with a(prime(k)) = prime(k+1).

Original entry on oeis.org

1, 3, 5, 9, 7, 15, 11, 27, 25, 21, 13, 45, 17, 33, 35, 81, 19, 75, 23, 63, 55, 39, 29, 135, 49, 51, 125, 99, 31, 105, 37, 243, 65, 57, 77, 225, 41, 69, 85, 189, 43, 165, 47, 117, 175, 87, 53, 405, 121, 147, 95, 153, 59, 375, 91, 297, 115, 93, 61, 315, 67, 111, 275, 729, 119
Offset: 1

Views

Author

Keywords

Comments

Meyers (see Guy reference) conjectures that for all r >= 1, the least odd number not in the set {a(i): i < prime(r)} is prime(r+1). - N. J. A. Sloane, Jan 08 2021
Meyers' conjecture would be refuted if and only if for some r there were such a large gap between prime(r) and prime(r+1) that there existed a composite c for which prime(r) < c < a(c) < prime(r+1), in which case (by Bertrand's postulate) c would necessarily be a term of A246281. - Antti Karttunen, Mar 29 2021
a(n) is odd for all n and for each odd m there exists a k with a(k) = m (see A064216). a(n) > n for n > 1: bijection between the odd and all numbers. - Reinhard Zumkeller, Sep 26 2001
a(n) and n have the same number of distinct primes with (A001222) and without multiplicity (A001221). - Michel Marcus, Jun 13 2014
From Antti Karttunen, Nov 01 2019: (Start)
More generally, a(n) has the same prime signature as n, A046523(a(n)) = A046523(n). Also A246277(a(n)) = A246277(n) and A287170(a(n)) = A287170(n).
Many permutations and other sequences that employ prime factorization of n to encode either polynomials, partitions (via Heinz numbers) or multisets in general can be easily defined by using this sequence as one of their constituent functions. See the last line in the Crossrefs section for examples.
(End)

Examples

			a(12) = a(2^2 * 3) = a(prime(1)^2 * prime(2)) = prime(2)^2 * prime(3) = 3^2 * 5 = 45.
a(A002110(n)) = A002110(n + 1) / 2.
		

References

  • Richard K. Guy, editor, Problems From Western Number Theory Conferences, Labor Day, 1983, Problem 367 (Proposed by Leroy F. Meyers, The Ohio State U.).

Crossrefs

See A045965 for another version.
Row 1 of table A242378 (which gives the "k-th powers" of this sequence), row 3 of A297845 and of A306697. See also arrays A066117, A246278, A255483, A308503, A329050.
Cf. A064989 (a left inverse), A064216, A000040, A002110, A000265, A027746, A046523, A048673 (= (a(n)+1)/2), A108228 (= (a(n)-1)/2), A191002 (= a(n)*n), A252748 (= a(n)-2n), A286385 (= a(n)-sigma(n)), A283980 (= a(n)*A006519(n)), A341529 (= a(n)*sigma(n)), A326042, A049084, A001221, A001222, A122111, A225546, A260443, A245606, A244319, A246269 (= A065338(a(n))), A322361 (= gcd(n, a(n))), A305293.
Cf. A249734, A249735 (bisections).
Cf. A246261 (a(n) is of the form 4k+1), A246263 (of the form 4k+3), A246271, A246272, A246259, A246281 (n such that a(n) < 2n), A246282 (n such that a(n) > 2n), A252742.
Cf. A275717 (a(n) > a(n-1)), A275718 (a(n) < a(n-1)).
Cf. A003972 (Möbius transform), A003973 (Inverse Möbius transform), A318321.
Cf. A300841, A305421, A322991, A250469, A269379 for analogous shift-operators in other factorization and quasi-factorization systems.
Cf. also following permutations and other sequences that can be defined with the help of this sequence: A005940, A163511, A122111, A260443, A206296, A265408, A265750, A275733, A275735, A297845, A091202 & A091203, A250245 & A250246, A302023 & A302024, A302025 & A302026.
A version for partition numbers is A003964, strict A357853.
A permutation of A005408.
Applying the same transformation again gives A357852.
Other multiplicative sequences: A064988, A357977, A357978, A357980, A357983.
A056239 adds up prime indices, row-sums of A112798.

Programs

  • Haskell
    a003961 1 = 1
    a003961 n = product $ map (a000040 . (+ 1) . a049084) $ a027746_row n
    -- Reinhard Zumkeller, Apr 09 2012, Oct 09 2011
    (MIT/GNU Scheme, with Aubrey Jaffer's SLIB Scheme library)
    (require 'factor)
    (define (A003961 n) (apply * (map A000040 (map 1+ (map A049084 (factor n))))))
    ;; Antti Karttunen, May 20 2014
    
  • Maple
    a:= n-> mul(nextprime(i[1])^i[2], i=ifactors(n)[2]):
    seq(a(n), n=1..80);  # Alois P. Heinz, Sep 13 2017
  • Mathematica
    a[p_?PrimeQ] := a[p] = Prime[ PrimePi[p] + 1]; a[1] = 1; a[n_] := a[n] = Times @@ (a[#1]^#2& @@@ FactorInteger[n]); Table[a[n], {n, 1, 65}] (* Jean-François Alcover, Dec 01 2011, updated Sep 20 2019 *)
    Table[Times @@ Map[#1^#2 & @@ # &, FactorInteger[n] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}] - Boole[n == 1], {n, 65}] (* Michael De Vlieger, Mar 24 2017 *)
  • PARI
    a(n)=local(f); if(n<1,0,f=factor(n); prod(k=1,matsize(f)[1],nextprime(1+f[k,1])^f[k,2]))
    
  • PARI
    a(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ Michel Marcus, May 17 2014
    
  • Perl
    use ntheory ":all";  sub a003961 { vecprod(map { next_prime($) } factor(shift)); }  # _Dana Jacobsen, Mar 06 2016
    
  • Python
    from sympy import factorint, prime, primepi, prod
    def a(n):
        f=factorint(n)
        return 1 if n==1 else prod(prime(primepi(i) + 1)**f[i] for i in f)
    [a(n) for n in range(1, 11)] # Indranil Ghosh, May 13 2017

Formula

If n = Product p(k)^e(k) then a(n) = Product p(k+1)^e(k).
Multiplicative with a(p^e) = A000040(A000720(p)+1)^e. - David W. Wilson, Aug 01 2001
a(n) = Product_{k=1..A001221(n)} A000040(A049084(A027748(n,k))+1)^A124010(n,k). - Reinhard Zumkeller, Oct 09 2011 [Corrected by Peter Munn, Nov 11 2019]
A064989(a(n)) = n and a(A064989(n)) = A000265(n). - Antti Karttunen, May 20 2014 & Nov 01 2019
A001221(a(n)) = A001221(n) and A001222(a(n)) = A001222(n). - Michel Marcus, Jun 13 2014
From Peter Munn, Oct 31 2019: (Start)
a(n) = A225546((A225546(n))^2).
a(A225546(n)) = A225546(n^2).
(End)
Sum_{k=1..n} a(k) ~ c * n^2, where c = (1/2) * Product_{p prime} ((p^2-p)/(p^2-nextprime(p))) = 2.06399637... . - Amiram Eldar, Nov 18 2022

A099884 XOR difference triangle of the powers of 2, read by rows; Square array A(row,col): A(0,col) = 2^col, A(row,col) = A048724(A(row-1, col)) for row > 0, read by descending antidiagonals.

Original entry on oeis.org

1, 2, 3, 4, 6, 5, 8, 12, 10, 15, 16, 24, 20, 30, 17, 32, 48, 40, 60, 34, 51, 64, 96, 80, 120, 68, 102, 85, 128, 192, 160, 240, 136, 204, 170, 255, 256, 384, 320, 480, 272, 408, 340, 510, 257, 512, 768, 640, 960, 544, 816, 680, 1020, 514, 771, 1024, 1536, 1280, 1920
Offset: 0

Views

Author

Paul D. Hanna, Oct 28 2004

Keywords

Comments

Define an "XOR difference triangle" for a sequence A by the following process. Start with A in the leftmost column. Generate the next column by performing the XOR operation between adjacent terms of the prior column. Repeat this process to generate the XOR difference triangle for A. Further, we define the "XOR BINOMIAL transform" of A as the main diagonal in the XOR difference triangle for A. The XOR BINOMIAL transform is its self-inverse. Let a sequence B be the XOR BINOMIAL transform of A, then we may express B by: B(n) = SumXOR_{k=0..n} A047999(n,k)*A(k), which is equivalent to: B(n) = (C(n,0)mod 2)*A(0) XOR (C(n,1)mod 2)*A(1) XOR (C(n,2)mod 2)*A(2) XOR ... XOR (X(n,n)mod 2)*A(n), where the coefficients are C(n,k)(mod 2) = A047999(n,k).
This sequence is a rearrangement of the numbers which are 2^k times distinct Fermat numbers (numbers of the form 2^(2^m) + 1). This matches the sizes of polygons constructible with compass and straightedge (A003401) up to 2^32+1, which is the first nonprime Fermat number. - Franklin T. Adams-Watters, Jun 16 2006

Examples

			The main diagonal equals A001317 (Pascal's triangle mod 2 in decimal):
{1,3,5,15,17,51,85,255,257,771,1285,3855,...}, and defines the XOR BINOMIAL transform of the powers of 2.
Rows begin:
  1;
  2, 3;
  4, 6, 5;
  8, 12, 10, 15;
  16, 24, 20, 30, 17;
  32, 48, 40, 60, 34, 51;
  64, 96, 80, 120, 68, 102, 85;
  128, 192, 160, 240, 136, 204, 170, 255;
  256, 384, 320, 480, 272, 408, 340, 510, 257;
  512, 768, 640, 960, 544, 816, 680, 1020, 514, 771;
  1024, 1536, 1280, 1920, 1088, 1632, 1360, 2040, 1028, 1542, 1285;
  2048, 3072, 2560, 3840, 2176, 3264, 2720, 4080, 2056, 3084, 2570, 3855;
  ...
From _Antti Karttunen_, Sep 19 2016: (Start)
Viewed as a square array, the top left corner looks like this:
     1,    2,     4,     8,    16,     32,     64,    128
     3,    6,    12,    24,    48,     96,    192,    384
     5,   10,    20,    40,    80,    160,    320,    640
    15,   30,    60,   120,   240,    480,    960,   1920
    17,   34,    68,   136,   272,    544,   1088,   2176
    51,  102,   204,   408,   816,   1632,   3264,   6528
    85,  170,   340,   680,  1360,   2720,   5440,  10880
   255,  510,  1020,  2040,  4080,   8160,  16320,  32640
   257,  514,  1028,  2056,  4112,   8224,  16448,  32896
   771, 1542,  3084,  6168, 12336,  24672,  49344,  98688
  1285, 2570,  5140, 10280, 20560,  41120,  82240, 164480
  3855, 7710, 15420, 30840, 61680, 123360, 246720, 493440
  4369, 8738, 17476, 34952, 69904, 139808, 279616, 559232
  ...
(End)
The square array shown above can be viewed as a subtable of a multiplication table with particular relevance to the carryless multiplication defined by A048720, as the first column gives the A048720 powers of 3 (and the first row gives powers of 2, which are the same as in standard arithmetic). - _Peter Munn_, Jan 13 2020
		

Crossrefs

Essentially GF(2)[X] analog of table A036561. - Antti Karttunen, Jan 18 2020
Cf. A047999, A158875 (row sums).
Cf. A000079 (first column of triangular table, the topmost row of square array).
Cf. A001317 (the rightmost diagonal of triangular table, the leftmost column of square array).
Cf. A099885, A117998 (central diagonals).
Cf. A276618 (transpose), A091202, A193231.

Programs

  • Mathematica
    a[n_]:= Sum[Mod[Binomial[n, i], 2]*2^i, {i, 0, n}]; T[n_, k_]:=2^(n - k)a[k]; Table[T[n, k], {n, 0, 20}, {k, 0, n}] // Flatten (* Indranil Ghosh, Apr 11 2017 *)
  • PARI
    {T(n,k)=local(B);B=0;for(i=0,k,B=bitxor(B,binomial(k,i)%2*2^(n-i)));B}
    for(n=0,10,for(k=0,n,print1(T(n,k),", "));print(""))
    
  • Python
    from sympy import binomial
    def a(n):
        return sum((binomial(n, i)%2)*2**i for i in range(n + 1))
    def T(n, k): return 2**(n - k)*a(k)
    for n in range(21): print([T(n, k) for k in range(n + 1)]) # Indranil Ghosh, Apr 11 2017
  • Scheme
    (define (A099884 n) (A099884bi (A002262 n) (A025581 n)))
    ;; Then use either this recurrence:
    (define (A099884bi row col) (if (zero? row) (A000079 col) (A048724 (A099884bi (- row 1) col))))
    ;; or this one:
    (define (A099884bi row col) (if (zero? col) (A001317 row) (* 2 (A099884bi row (- col 1)))))
    ;; Antti Karttunen, Sep 19 2016
    

Formula

T(n, k) = 2^(n-k)*A001317(k). T(n, n) = A001317(n) = SumXOR_{k=0..n} A047999(n, k)*2^k, where SumXOR is the analog of summation under the binary XOR operation.
From Antti Karttunen, Sep 19 2016: (Start)
When viewed as a square array A(row,col), with row >= 0, col >= 0, the following recurrences and formulas are valid:
A(0,col) = A000079(col), for row > 0, A(row,col) = A048724(A(row-1, col)).
A(row,0) = A001317(row), for col > 0, A(row,col) = 2*A(row,col-1).
A(row,col) = A248663(A066117(row+1,col+1)) = A048675(A255483(row,col+1)).
(End)
With the definitions from Antti Karttunen above, A(row+1, col) = A048720(3, A(row, col)). - Peter Munn, Jan 13 2020
A(n,k) = A193231(A(k,n)) = A091202(A036561(n,k)). - Antti Karttunen, Jan 18 2020

Extensions

Square array interpretation added as a second, alternative description by Antti Karttunen, Sep 19 2016

A297845 Encoded multiplication table for polynomials in one indeterminate with nonnegative integer coefficients. Symmetric square array T(n, k) read by antidiagonals, n > 0 and k > 0. See comment for details.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 5, 4, 1, 1, 5, 9, 9, 5, 1, 1, 6, 7, 16, 7, 6, 1, 1, 7, 15, 25, 25, 15, 7, 1, 1, 8, 11, 36, 11, 36, 11, 8, 1, 1, 9, 27, 49, 35, 35, 49, 27, 9, 1, 1, 10, 25, 64, 13, 90, 13, 64, 25, 10, 1, 1, 11, 21, 81, 125, 77, 77, 125, 81
Offset: 1

Views

Author

Rémy Sigrist, Jan 10 2018

Keywords

Comments

For any number n > 0, let f(n) be the polynomial in a single indeterminate x where the coefficient of x^e is the prime(1+e)-adic valuation of n (where prime(k) denotes the k-th prime); f establishes a bijection between the positive numbers and the polynomials in a single indeterminate x with nonnegative integer coefficients; let g be the inverse of f; T(n, k) = g(f(n) * f(k)).
This table has many similarities with A248601.
For any n > 0 and m > 0, f(n * m) = f(n) + f(m).
Also, f(1) = 0 and f(2) = 1.
The function f can be naturally extended to the set of positive rational numbers: if r = u/v (not necessarily in reduced form), then f(r) = f(u) - f(v); as such, f is a homomorphism from the multiplicative group of positive rational numbers to the additive group of polynomials of a single indeterminate x with integer coefficients.
See A297473 for the main diagonal of T.
As a binary operation, T(.,.) is related to A306697(.,.) and A329329(.,.). When their operands are terms of A050376 (sometimes called Fermi-Dirac primes) the three operations give the same result. However the rest of the multiplication table for T(.,.) can be derived from these results because T(.,.) distributes over integer multiplication (A003991), whereas for A306697 and A329329, the equivalent derivation uses distribution over A059896(.,.) and A059897(.,.) respectively. - Peter Munn, Mar 25 2020
From Peter Munn, Jun 16 2021: (Start)
The operation defined by this sequence can be extended to be the multiplicative operator of a ring over the positive rationals that is isomorphic to the polynomial ring Z[x]. The extended function f (described in the author's original comments) is the isomorphism we use, and it has the same relationship with the extended operation that exists between their unextended equivalents.
Denoting this extension of T(.,.) as t_Q(.,.), we get t_Q(n, 1/k) = t_Q(1/n, k) = 1/T(n, k) and t_Q(1/n, 1/k) = T(n, k) for positive integers n and k. The result for other rationals is derived from the distributive property: t_Q(q, r*s) = t_Q(q, r) * t_Q(q, s); t_Q(q*r, s) = t_Q(q, s) * t_Q(r, s). This may look unusual because standard multiplication of rational numbers takes on the role of the ring's additive group.
There are many OEIS sequences that can be shown to be a list of the integers in an ideal of this ring. See the cross-references.
There are some completely additive sequences that similarly define by extension completely additive functions on the positive rationals that can be shown to be homomorphisms from this ring onto the integer ring Z, and these functions relate to some of the ideals. For example, the extended function of A048675, denoted A048675_Q, maps i/j to A048675(i) - A048675(j) for positive integers i and j. For any positive integer k, the set {r rational > 0 : k divides A048675_Q(r)} forms an ideal of the ring; for k=2 and k=3 the integers in this ideal are listed in A003159 and A332820 respectively.
(End)

Examples

			Array T(n, k) begins:
  n\k|  1   2   3    4    5    6    7     8    9    10
  ---+------------------------------------------------
    1|  1   1   1    1    1    1    1     1    1     1  -> A000012
    2|  1   2   3    4    5    6    7     8    9    10  -> A000027
    3|  1   3   5    9    7   15   11    27   25    21  -> A003961
    4|  1   4   9   16   25   36   49    64   81   100  -> A000290
    5|  1   5   7   25   11   35   13   125   49    55  -> A357852
    6|  1   6  15   36   35   90   77   216  225   210  -> A191002
    7|  1   7  11   49   13   77   17   343  121    91
    8|  1   8  27   64  125  216  343   512  729  1000  -> A000578
    9|  1   9  25   81   49  225  121   729  625   441
   10|  1  10  21  100   55  210   91  1000  441   550
From _Peter Munn_, Jun 24 2021: (Start)
The encoding, n, of polynomials, f(n), that is used for the table is further described in A206284. Examples of encoded polynomials:
   n      f(n)        n           f(n)
   1         0       16              4
   2         1       17            x^6
   3         x       21        x^3 + x
   4         2       25           2x^2
   5       x^2       27             3x
   6     x + 1       35      x^3 + x^2
   7       x^3       36         2x + 2
   8         3       49           2x^3
   9        2x       55      x^4 + x^2
  10   x^2 + 1       64              6
  11       x^4       77      x^4 + x^3
  12     x + 2       81             4x
  13       x^5       90   x^2 + 2x + 1
  15   x^2 + x       91      x^5 + x^3
(End)
		

Crossrefs

Row n: n=1: A000012, n=2: A000027, n=3: A003961, n=4: A000290, n=5: A357852, n=6: A191002, n=8: A000578.
Main diagonal: A297473.
Functions f satisfying f(T(n,k)) = f(n) * f(k): A001222, A048675 (and similarly, other rows of A104244), A195017.
Powers of k: k=3: A000040, k=4: A001146, k=5: A031368, k=6: A007188 (see also A066117), k=7: A031377, k=8: A023365, k=9: main diagonal of A329050.
Integers in the ideal of the related ring (see Jun 2021 comment) generated by S: S={3}: A005408, S={4}: A000290\{0}, S={4,3}: A003159, S={5}: A007310, S={5,4}: A339690, S={6}: A325698, S={6,4}: A028260, S={7}: A007775, S={8}: A000578\{0}, S={8,3}: A191257, S={8,6}: A332820, S={9}: A016754, S={10,4}: A340784, S={11}: A008364, S={12,8}: A145784, S={13}: A008365, S={15,4}: A345452, S={15,9}: A046337, S={16}: A000583\{0}, S={17}: A008366.
Equivalent sequence for polynomial composition: A326376.
Related binary operations: A003991, A306697/A059896, A329329/A059897.

Programs

  • PARI
    T(n,k) = my (f=factor(n), p=apply(primepi, f[, 1]~), g=factor(k), q=apply(primepi, g[, 1]~)); prod (i=1, #p, prod(j=1, #q, prime(p[i]+q[j]-1)^(f[i, 2]*g[j, 2])))

Formula

T is completely multiplicative in both parameters:
- for any n > 0
- and k > 0 with prime factorization Prod_{i > 0} prime(i)^e_i:
- T(prime(n), k) = T(k, prime(n)) = Prod_{i > 0} prime(n + i - 1)^e_i.
For any m > 0, n > 0 and k > 0:
- T(n, k) = T(k, n) (T is commutative),
- T(m, T(n, k)) = T(T(m, n), k) (T is associative),
- T(n, 1) = 1 (1 is an absorbing element for T),
- T(n, 2) = n (2 is an identity element for T),
- T(n, 2^i) = n^i for any i >= 0,
- T(n, 4) = n^2 (A000290),
- T(n, 8) = n^3 (A000578),
- T(n, 3) = A003961(n),
- T(n, 3^i) = A003961(n)^i for any i >= 0,
- T(n, 6) = A191002(n),
- A001221(T(n, k)) <= A001221(n) * A001221(k),
- A001222(T(n, k)) = A001222(n) * A001222(k),
- A055396(T(n, k)) = A055396(n) + A055396(k) - 1 when n > 1 and k > 1,
- A061395(T(n, k)) = A061395(n) + A061395(k) - 1 when n > 1 and k > 1,
- T(A000040(n), A000040(k)) = A000040(n + k - 1),
- T(A000040(n)^i, A000040(k)^j) = A000040(n + k - 1)^(i * j) for any i >= 0 and j >= 0.
From Peter Munn, Mar 13 2020 and Apr 20 2021: (Start)
T(A329050(i_1, j_1), A329050(i_2, j_2)) = A329050(i_1+i_2, j_1+j_2).
T(n, m*k) = T(n, m) * T(n, k); T(n*m, k) = T(n, k) * T(m, k) (T distributes over multiplication).
A104244(m, T(n, k)) = A104244(m, n) * A104244(m, k).
For example, for m = 2, the above formula is equivalent to A048675(T(n, k)) = A048675(n) * A048675(k).
A195017(T(n, k)) = A195017(n) * A195017(k).
A248663(T(n, k)) = A048720(A248663(n), A248663(k)).
T(n, k) = A306697(n, k) if and only if T(n, k) = A329329(n, k).
A007913(T(n, k)) = A007913(T(A007913(n), A007913(k))) = A007913(A329329(n, k)).
(End)

Extensions

New name from Peter Munn, Jul 17 2021

A007188 Multiplicative encoding of Pascal triangle: Product p(i+1)^C(n,i).

Original entry on oeis.org

2, 6, 90, 47250, 66852843750, 2806877704512541816406250, 1216935896582703898519354781702537118597533386230468750
Offset: 0

Views

Author

Keywords

Comments

n-th power of x+1 using the encoding of polynomials defined in A206284 and A297845. - Peter Munn, Jul 20 2022

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 828.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Leftmost column of square array A066117.

Programs

  • Mathematica
    c[n_] := CoefficientList[(1 + x)^n, x]; f[n_] := Product[Prime[k]^c[n][[k]], {k, 1, Length[c[n]]}]; Table[f[n], {n, 1, 7}] (* Clark Kimberling, Feb 05 2012 *)

Formula

a(0) = 2; for n > 0, a(n) = A297845(a(n-1), 6). - Peter Munn, Jul 20 2022

A255483 Infinite square array read by antidiagonals downwards: T(0,m) = prime(m), m >= 1; for n >= 1, T(n,m) = T(n-1,m)*T(n-1,m+1)/gcd(T(n-1,m), T(n-1,m+1))^2, m >= 1.

Original entry on oeis.org

2, 3, 6, 5, 15, 10, 7, 35, 21, 210, 11, 77, 55, 1155, 22, 13, 143, 91, 5005, 39, 858, 17, 221, 187, 17017, 85, 3315, 1870, 19, 323, 247, 46189, 133, 11305, 5187, 9699690, 23, 437, 391, 96577, 253, 33649, 21505, 111546435, 46
Offset: 0

Views

Author

N. J. A. Sloane, Feb 28 2015

Keywords

Comments

The first column of the array is given by A123098; subsequent columns are obtained by applying the function A003961, i.e., replacing each prime factor by the next larger prime. - M. F. Hasler, Sep 17 2016
Interpretation with respect to A329329 from Peter Munn, Feb 08 2020: (Start)
With respect to the ring defined by A329329 and A059897, the first row gives powers of 3, the first column gives powers of 6, both in order of increasing exponent, and the body of the table gives their products. A329049 is the equivalent table in which the first column gives powers of 4.
A099884 is the equivalent table for the ring defined by A048720 and A003987. That ring is an image of the polynomial ring GF(2)[x] using a standard representation of the polynomials as integers. A329329 describes a comparable mapping to integers from the related polynomial ring GF(2)[x,y].
Using these mappings, the tables here and in A099884 are matching images: the first row represents powers of x, the first column represents powers of (x+1) and the body of the table gives their products.
Hugo van der Sanden's formula (see formula section) indicates that A019565 provides a mapping from A099884. In the wider terms described above, A019565 is an injective homomorphism between images of the 2 polynomial rings, and maps the image of each GF(2)[x] polynomial to the image of the equivalent GF(2)[x,y] polynomial.
(End)

Examples

			The top left corner of the array, row index 0..5, column index 1..10:
    2,    3,     5,     7,    11,     13,     17,     19,      23,      29
    6,   15,    35,    77,   143,    221,    323,    437,     667,     899
   10,   21,    55,    91,   187,    247,    391,    551,     713,    1073
  210, 1155,  5005, 17017, 46189,  96577, 215441, 392863,  765049, 1363783
   22,   39,    85,   133,   253,    377,    527,    703,     943,    1247
  858, 3315, 11305, 33649, 95381, 198679, 370481, 662929, 1175921, 1816879
		

Crossrefs

First two columns = A123098, A276804.
A kind of generalization of A036262.
Transpose: A276578, terms sorted into ascending order: A276579.
A003987, A048720, A059897, A329049 relate to the A329329 polynomial ring interpretation.

Programs

  • Maple
    T:= proc(n, m) option remember; `if`(n=0, ithprime(m),
          T(n-1, m)*T(n-1, m+1)/igcd(T(n-1, m), T(n-1, m+1))^2)
        end:
    seq(seq(T(n, 1+d-n), n=0..d), d=0..10);  # Alois P. Heinz, Feb 28 2015
  • Mathematica
    T[n_, m_] := T[n, m] = If[n == 0, Prime[m], T[n-1, m]*T[n-1, m+1]/GCD[T[n-1, m], T[n-1, m+1]]^2]; Table[Table[T[n, 1+d-n], {n, 0, d}], {d, 0, 10}] // Flatten (* Jean-François Alcover, Mar 09 2015, after Alois P. Heinz *)
  • PARI
    T=matrix(N=15,N);for(j=1,N,T[1,j]=prime(j));(f(x,y)=x*y/gcd(x,y)^2);for(k=1,N-1,for(j=1,N-k,T[k+1,j]=f(T[k,j],T[k,j+1])));A255483=concat(vector(N,i,vector(i,j,T[j,1+i-j]))) \\ M. F. Hasler, Sep 17 2016
    
  • PARI
    A255483(n,k)=prod(j=0,n,if(bitand(n-j,j),1,prime(j+k))) \\ M. F. Hasler, Sep 18 2016
    
  • Scheme
    (define (A255483 n) (A255483bi (A002262 n) (+ 1 (A025581 n))))
    ;; Then use either an almost standalone version (requiring only A000040):
    (define (A255483bi row col) (if (zero? row) (A000040 col) (let ((a (A255483bi (- row 1) col)) (b (A255483bi (- row 1) (+ col 1)))) (/ (lcm a b) (gcd a b)))))
    ;; Or one based on M. F. Hasler's new recurrence:
    (define (A255483bi row col) (if (= 1 col) (A123098 row) (A003961 (A255483bi row (- col 1)))))
    ;; Antti Karttunen, Sep 18 2016

Formula

T(n,1) = A123098(n), T(n,m+1) = A003961(T(n,m)), for all n >= 0, m >= 1. - M. F. Hasler, Sep 17 2016
T(n,m) = Prod_{k=0..n} prime(k+m)^(!(n-k & k)) where !x is 1 if x=0 and 0 else, and & is binary AND. - M. F. Hasler, Sep 18 2016
From Antti Karttunen, Sep 18 2016: (Start)
For n >= 1, m >= 1, T(n,m) = lcm(T(n-1,m),T(n-1,m+1)) / gcd(T(n-1,m),T(n-1,m+1)).
T(n,k) = A007913(A066117(n+1,k)).
T(n,k) = A019565(A099884(n,k-1)) [After Hugo van der Sanden's observations on SeqFan-list].
(End)
From Peter Munn, Jan 08 2020: (Start)
T(0,1) = 2, and for n >= 0, k >= 1, T(n+1,k) = A329329(T(n,k), 6), T(n,k+1) = A329329(T(n,k), 3).
T(n,k) = A329329(T(n,1), T(0,k)).
(End)

A276588 Square array A(row,col) = Sum_{k=0..row} binomial(row,k)*(1+col+k)!, read by descending antidiagonals as A(0,0), A(0,1), A(1,0), A(0,2), A(1,1), A(2,0), ...

Original entry on oeis.org

1, 2, 3, 6, 8, 11, 24, 30, 38, 49, 120, 144, 174, 212, 261, 720, 840, 984, 1158, 1370, 1631, 5040, 5760, 6600, 7584, 8742, 10112, 11743, 40320, 45360, 51120, 57720, 65304, 74046, 84158, 95901, 362880, 403200, 448560, 499680, 557400, 622704, 696750, 780908, 876809, 3628800, 3991680, 4394880, 4843440, 5343120, 5900520, 6523224, 7219974, 8000882, 8877691
Offset: 0

Views

Author

Antti Karttunen, Sep 19 2016

Keywords

Examples

			The top left corner of the array:
     1,     2,     6,     24,     120,      720,      5040,      40320
     3,     8,    30,    144,     840,     5760,     45360,     403200
    11,    38,   174,    984,    6600,    51120,    448560,    4394880
    49,   212,  1158,   7584,   57720,   499680,   4843440,   51932160
   261,  1370,  8742,  65304,  557400,  5343120,  56775600,  661933440
  1631, 10112, 74046, 622704, 5900520, 62118720, 718709040, 9059339520
		

Crossrefs

Transpose: A276589.
Topmost row (row 0): A000142, Row 1: A001048 (without its initial 2), Row 2: A001344 (from a(1) = 11 onward), Row 3: A001345 (from a(1) = 49 onward), Row 4: A001346 (from a(1) = 261 onward), Row 5: A001347 (from a(1) = 1631 onward).
Leftmost column (column 0): A001339, Column 1: A001340, Columns 2-3: A001341 & A001342 (apparently).
Cf. A276075.
Cf. also arrays A066117, A276586, A099884, A255483.

Programs

  • Mathematica
    T[r_, c_]:=Sum[Binomial[r, k](1 + c + k)!, {k, 0, r}]; Table[T[c, r - c], {r, 0, 10}, {c, 0, r}] // Flatten (* Indranil Ghosh, Apr 11 2017 *)
  • PARI
    T(r, c) = sum(k=0, r, binomial(r, k)*(1 + c + k)!);
    for(r=0, 10, for(c=0, r, print1(T(c, r - c),", ");); print();) \\ Indranil Ghosh, Apr 11 2017
    
  • Python
    from sympy import binomial, factorial
    def T(r, c): return sum([binomial(r, k) * factorial(1 + c + k) for k in range(r + 1)])
    for r in range(11): print([T(c, r - c) for c in range(r + 1)]) # Indranil Ghosh, Apr 11 2017
  • Scheme
    (define (A276588 n) (A276588bi (A002262 n) (A025581 n)))
    (define (A276588bi row col) (A276075 (A066117bi (+ 1 row) (+ 1 col)))) ;; Code for A066117bi given in A066117, and for A276075 under the respective entry.
    

Formula

A(row,col) = Sum_{k=0..row} binomial(row,k)*A000142(1+col+k).
A(row,col) = A276075(A066117(row+1,col+1)).

A267096 a(n) = Product_{i=0..n} prime(i+2)^binomial(n,i).

Original entry on oeis.org

3, 15, 525, 1414875, 41985913344375, 433555011900329243987584396875, 3514495551481947615680580256869117013417604971088496013610671875
Offset: 0

Views

Author

Antti Karttunen, Feb 06 2016

Keywords

Examples

			Terms are obtained by exponentiating the odd primes in range [3 .. prime(2+n)] with the binomial coefficients obtained from row n of Pascal's triangle (A007318) and then multiplying the factors together:
            3^1
         3^1 * 5^1
      3^1 * 5^2 * 7^1
   3^1 * 5^3 * 7^3 * 11^1
3^1 * 5^4 * 7^6 * 11^4 * 13^1
etc.
		

Crossrefs

Second column (or diagonal from right) in A066117.

Programs

  • Scheme
    (define (A267096 n) (mul (lambda (k) (expt (A000040 (+ 2 k)) (A007318tr n k))) 0 n)) ;; Where A007318tr gives binomial coefficients, as in A007318.
    (define (mul intfun lowlim uplim) (let multloop ((i lowlim) (res 1)) (cond ((> i uplim) res) (else (multloop (1+ i) (* res (intfun i)))))))

Formula

a(n) = Product_{i=0..n} prime(i+2)^C(n,i).
a(n) = A003961(A007188(n)).

A276586 Square array A(row,col) = Sum_{k=0..row} binomial(row,k)*A002110(col+k), read by descending antidiagonals as A(0,0), A(0,1), A(1,0), A(0,2), A(1,1), A(2,0), ...

Original entry on oeis.org

1, 2, 3, 6, 8, 11, 30, 36, 44, 55, 210, 240, 276, 320, 375, 2310, 2520, 2760, 3036, 3356, 3731, 30030, 32340, 34860, 37620, 40656, 44012, 47743, 510510, 540540, 572880, 607740, 645360, 686016, 730028, 777771, 9699690, 10210200, 10750740, 11323620, 11931360, 12576720, 13262736, 13992764, 14770535
Offset: 0

Views

Author

Antti Karttunen, Sep 18 2016

Keywords

Examples

			The top left corner of the array:
     1,     2,      6,       30,       210,       2310,        30030
     3,     8,     36,      240,      2520,      32340,       540540
    11,    44,    276,     2760,     34860,     572880,     10750740
    55,   320,   3036,    37620,    607740,   11323620,    253753500
   375,  3356,  40656,   645360,  11931360,  265077120,   7422334920
  3731, 44012, 686016, 12576720, 277008480, 7687412040, 235239464460
		

Crossrefs

Transpose: A276587.
Topmost row: A002110, Leftmost column: A136104.
Cf. also arrays A066117, A276588, A099884, A255483.

Programs

  • Mathematica
    primorial[n_] := Product[Prime[k], {k, 1, n}]; A[n_, k_] := Sum[Binomial[n, j]*primorial[k+j], {j, 0, n}]; Table[A[n-k, k], {n, 0, 9}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Jan 22 2017 *)
  • PARI
    P(n)=prod(i=1, n, prime(i));
    T(n, k) = sum(j=0, n, binomial(n, j)*P(k + j));
    for(n=0, 10, for(k=0, n, print1(T(k, n - k),", ");); print();) \\ Indranil Ghosh, Apr 11 2017
  • Scheme
    (define (A276586 n) (A276586bi (A002262 n) (A025581 n)))
    (define (A276586bi row col) (A276085 (A066117bi (+ 1 row) (+ 1 col))))
    

Formula

A(row,col) = Sum_{k=0..row} binomial(row,k)*A002110(col+k).
A(row,col) = A276085(A066117(row+1,col+1)).

A066116 a(n) = prime(n-2)*prime(n-1)^2*prime(n).

Original entry on oeis.org

90, 525, 2695, 11011, 31603, 71383, 141151, 291479, 599633, 1031153, 1739999, 2674471, 3563023, 5034311, 7789357, 11254073, 14709113, 19441859, 24655531, 29890361, 37814219, 48436559, 63771971, 84577501, 101918191, 114651463
Offset: 3

Views

Author

Henry Bottomley, Dec 05 2001

Keywords

Examples

			a(3) = 2*3^2*5 = 90.
		

Crossrefs

Cf. A047928.

Programs

  • PARI
    { for (n=3, 1000, a=prime(n-2)*prime(n-1)^2*prime(n); write("b066116.txt", n, " ", a) ) } \\ Harry J. Smith, Feb 01 2010

Formula

a(n) = A000040(n-2)*A000040(n-1)^2*A000040(n) = A006094(n-1)*A006094(n) = A066117(3, n).

A375761 Infinite triangle T(n, k), n > 0, k = 1..n, read and filled in the greedy way by rows with distinct positive integers such that for any n > 1, k = 2..n, T(n, k) = T(n-1, k-1) * T(n, k-1).

Original entry on oeis.org

2, 3, 6, 4, 12, 72, 5, 20, 240, 17280, 7, 35, 700, 168000, 2903040000, 8, 56, 1960, 1372000, 230496000000, 669139107840000000000, 10, 80, 4480, 8780800, 12047257600000, 2776844687769600000000000, 1858095376984393503473664000000000000000000000
Offset: 1

Views

Author

Rémy Sigrist, Aug 27 2024

Keywords

Comments

This sequence is a variant of A035312, based on products instead of sums.
The value 1 cannot appear (as it would imply a duplicate term).
We can always start a new row with a prime number that does not divide prior terms, so the sequence is well defined.
Will every integer > 1 appear in the sequence?

Examples

			Triangle T(n, k) begins:
               2
           3      6
         4   12      72
       5  20    240      17280
     7  35  700    168000    2903040000
    8 56 1960 1372000 230496000000 669139107840000000000
		

Crossrefs

Cf. A035312, A064319, A066117, A375846 (right border), A375857 (left border).

Programs

  • PARI
    \\ See Links section.
Showing 1-10 of 10 results.