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.

Previous Showing 11-20 of 116 results. Next

A011764 a(n) = 3^(2^n) (or: write in base 3, read in base 9).

Original entry on oeis.org

3, 9, 81, 6561, 43046721, 1853020188851841, 3433683820292512484657849089281, 11790184577738583171520872861412518665678211592275841109096961
Offset: 0

Views

Author

Stephan Y Solomon (ilans(AT)way.com)

Keywords

Comments

a(n) is the second-highest value k such that A173419(k) = n+2. - Charles R Greathouse IV, Oct 03 2012
Let b(0) = 6; b(n+1) = smallest number such that b(n+1) + Product_{i=0..n} b(i) divides b(n+1)*Product_{i=0..n} b(i). Then b(n+1) = a(n) for n >= 0. - Derek Orr, Dec 13 2014
Changing "+" to "-": Let b(0) = 6; b(n+1) = smallest number such that b(n+1) - Product_{i=0..n} b(i) divides b(n+1)*Product_{i=0..n} b(i). Then b(n+2) = a(n) for n >= 0. - Derek Orr, Jan 04 2015
With offset = 1, a(n) is the number of collections C of subsets of {1,2,...,n} such that if S is in C then the complement of S is not in C. - Geoffrey Critzer, Feb 06 2017

Crossrefs

Subsequence of A000244 (powers of 3).

Programs

Formula

a(0) = 3 and a(n+1) = a(n)^2. - Benoit Jubin, Jun 27 2009
Sum_{n>=0} 1/a(n) = A078885. - Amiram Eldar, Nov 09 2020
Product_{n>=0} (1 + 1/a(n)) = 3/2. - Amiram Eldar, Jan 29 2021
a(n) = A000244(A000079(n)), or A011764 = A000244 o A000079. - M. F. Hasler, Jul 20 2023

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

A331590 Square array A(n,k) = A225546(A225546(n) * A225546(k)), n >= 1, k >= 1, read by descending antidiagonals.

Original entry on oeis.org

1, 2, 2, 3, 3, 3, 4, 6, 6, 4, 5, 8, 5, 8, 5, 6, 10, 12, 12, 10, 6, 7, 5, 15, 9, 15, 5, 7, 8, 14, 10, 20, 20, 10, 14, 8, 9, 12, 21, 24, 7, 24, 21, 12, 9, 10, 18, 24, 28, 30, 30, 28, 24, 18, 10, 11, 15, 27, 18, 35, 15, 35, 18, 27, 15, 11, 12, 22, 30, 36, 40, 42, 42, 40, 36, 30, 22, 12, 13, 24, 33, 40, 45, 20, 11, 20, 45, 40, 33, 24, 13
Offset: 1

Views

Author

Peter Munn, Jan 21 2020

Keywords

Comments

As a binary operation, this sequence defines a commutative monoid over the positive integers that is isomorphic to multiplication. The self-inverse permutation A225546(.) provides an isomorphism. This monoid therefore has unique factorization. Its primes are the even terms of A050376: 2, 4, 16, 256, ..., which in standard integer multiplication are the powers of 2 with powers of 2 as exponents.
In this monoid, in contrast, the powers of 2 run through the squarefree numbers, the k-th power of 2 being A019565(k). 4 is irreducible and its powers are the squares of the squarefree numbers, the k-th power of 4 being A019565(k)^2 (where "^2" denotes standard integer squaring); and so on with powers of 16, 256, ...
In many cases the product of two numbers is the same here as in standard integer multiplication. See the formula section for details.

Examples

			From _Antti Karttunen_, Feb 02 2020: (Start)
The top left 16 X 16 corner of the array:
   1,  2,  3,  4,  5,  6,   7,   8,   9,  10,  11,  12,  13,  14,  15,  16, ...
   2,  3,  6,  8, 10,  5,  14,  12,  18,  15,  22,  24,  26,  21,  30,  32, ...
   3,  6,  5, 12, 15, 10,  21,  24,  27,  30,  33,  20,  39,  42,   7,  48, ...
   4,  8, 12,  9, 20, 24,  28,  18,  36,  40,  44,  27,  52,  56,  60,  64, ...
   5, 10, 15, 20,  7, 30,  35,  40,  45,  14,  55,  60,  65,  70,  21,  80, ...
   6,  5, 10, 24, 30, 15,  42,  20,  54,   7,  66,  40,  78,  35,  14,  96, ...
   7, 14, 21, 28, 35, 42,  11,  56,  63,  70,  77,  84,  91,  22, 105, 112, ...
   8, 12, 24, 18, 40, 20,  56,  27,  72,  60,  88,  54, 104,  84, 120, 128, ...
   9, 18, 27, 36, 45, 54,  63,  72,  25,  90,  99, 108, 117, 126, 135, 144, ...
  10, 15, 30, 40, 14,  7,  70,  60,  90,  21, 110, 120, 130, 105,  42, 160, ...
  11, 22, 33, 44, 55, 66,  77,  88,  99, 110,  13, 132, 143, 154, 165, 176, ...
  12, 24, 20, 27, 60, 40,  84,  54, 108, 120, 132,  45, 156, 168,  28, 192, ...
  13, 26, 39, 52, 65, 78,  91, 104, 117, 130, 143, 156,  17, 182, 195, 208, ...
  14, 21, 42, 56, 70, 35,  22,  84, 126, 105, 154, 168, 182,  33, 210, 224, ...
  15, 30,  7, 60, 21, 14, 105, 120, 135,  42, 165,  28, 195, 210,  35, 240, ...
  16, 32, 48, 64, 80, 96, 112, 128, 144, 160, 176, 192, 208, 224, 240,  81, ...
(End)
		

Crossrefs

Isomorphic to A003991 with A225546 as isomorphism.
Cf. A003961(main diagonal), A048675, A059895, A059896, A059897.
Rows/columns, sorted in ascending order: 2: A000037, 3: A028983, 4: A252849.
A019565 lists powers of 2 in order of increasing exponent.
Powers of k, sorted in ascending order: k=2: A005117, k=3: A056911, k=4: A062503, k=5: A276378, k=6: intersection of A325698 and A005117, k=7: intersection of A007775 and A005117, k=8: A062838.
Irreducibles are A001146 (even terms of A050376).

Programs

  • PARI
    up_to = 1275;
    A019565(n) = factorback(vecextract(primes(logint(n+!n, 2)+1), n));
    A048675(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; };
    A331590sq(x,y) = if(1==x,y,if(1==y,x, my(fx=factor(x),fy=factor(y),u=max(#binary(vecmax(fx[, 2])),#binary(vecmax(fy[, 2]))),prodsx=vector(u,x,1),m=1); for(i=1,u,for(k=1,#fx~, if(bitand(fx[k,2],m),prodsx[i] *= fx[k,1])); for(k=1,#fy~, if(bitand(fy[k,2],m),prodsx[i] *= fy[k,1])); m<<=1); prod(i=1,u,A019565(A048675(prodsx[i]))^(1<<(i-1)))));
    A331590list(up_to) = { my(v = vector(up_to), i=0); for(a=1,oo, for(col=1,a, i++; if(i > up_to, return(v)); v[i] = A331590sq(col,(a-(col-1))))); (v); };
    v331590 = A331590list(up_to);
    A331590(n) = v331590[n]; \\ Antti Karttunen, Feb 02 2020

Formula

Alternative definition: A(n,1) = n; A(n,k) = A(A059897(n,k), A003961(A059895(n,k))).
Main derived identities: (Start)
A(n,k) = A(k,n).
A(1,n) = n.
A(n, A(m,k)) = A(A(n,m), k).
A(m,m) = A003961(m).
A(n^2, k^2) = A(n,k)^2.
A(A003961(n), A003961(k)) = A003961(A(n,k)).
A(A019565(n), A019565(k)) = A019565(n+k).
(End)
Characterization of conditions for A(n,k) = n * k: (Start)
The following 4 conditions are equivalent:
(1) A(n,k) = n * k;
(2) A(n,k) = A059897(n,k);
(3) A(n,k) = A059896(n,k);
(4) A059895(n,k) = 1.
If gcd(n,k) = 1, A(n,k) = n * k.
If gcd(n,k) = 1, A(A225546(n), A225546(k)) = A225546(n) * A225546(k).
The previous formula implies A(n,k) = n * k in the following cases:
(1) for n = A005117(m), k = j^2;
(2) more generally for n = A005117(m_1)^(2^i_1), k = A005117(m_2)^(2^i_2), with A004198(i_1, i_2) = 0.
(End)

A135764 Distribute the natural numbers in columns based on the occurrence of "2" in each prime factorization; square array A(row,col) = 2^(row-1) * ((2*col)-1), read by descending antidiagonals.

Original entry on oeis.org

1, 3, 2, 5, 6, 4, 7, 10, 12, 8, 9, 14, 20, 24, 16, 11, 18, 28, 40, 48, 32, 13, 22, 36, 56, 80, 96, 64, 15, 26, 44, 72, 112, 160, 192, 128, 17, 30, 52, 88, 144, 224, 320, 384, 256, 19, 34, 60, 104, 176, 288, 448, 640, 768, 512, 21, 38, 68, 120, 208, 352, 576, 896, 1280, 1536, 1024, 23, 42, 76, 136, 240, 416, 704, 1152, 1792, 2560, 3072, 2048, 25, 46, 84, 152, 272, 480, 832, 1408, 2304, 3584, 5120, 6144, 4096, 27, 50, 92, 168, 304, 544, 960, 1664, 2816
Offset: 1

Views

Author

Alford Arnold, Nov 29 2007

Keywords

Comments

The array in A135764 is identical to the array in A054582 [up to the transposition and different indexing. - Clark Kimberling, Dec 03 2010; comment amended by Antti Karttunen, Feb 03 2015; please see the illustration in Example section].
The array gives a bijection between the natural numbers N and N^2. A more usual bijection is to take the natural numbers A000027 and write them in the usual OEIS square array format. However this bijection has the advantage that it can be formed by iterating the usual bijection between N and 2N. - Joshua Zucker, Nov 04 2011
The array can be used to determine the configurations of k-th Towers of Hanoi moves, by labeling odd row terms C,B,A,C,B,A,... and even row terms B,C,A,B,C,A,.... Then given k equal to or greater than term "a" in each n-th row, but less than the next row term, record the label A, B, or C for term "a". This denotes the peg position for the disc corresponding to the n-th row. For example, with k = 25, five discs are in motion since the binary for 25 = 11001, five bits. We find that 25 in row 5 is greater than 16 labeled C, but less than 48. Thus, disc 5 is on peg C. In the 4th row, 25 is greater than 24 (a C), but less than 40, so goes onto the C peg. Similarly, disc 3 is on A, 2 is on A, and disc 1 is on A. Thus, discs 2 and 3 are on peg A, while 1, 4, and 5 are on peg C. - Gary W. Adamson, Jun 22 2012
Shares with arrays A253551 and A254053 the property that A001511(n) = k for all terms n on row k and when going downward in each column, terms grow by doubling. - Antti Karttunen, Feb 03 2015
Let P be the infinite palindromic word having initial word 0 and midword sequence (1,2,3,4,...) = A000027. Row n of the array A135764 gives the positions of n-1 in S. ("Infinite palindromic word" is defined at A260390.) - Clark Kimberling, Aug 13 2015
The probability distribution series 1 = 2/3 + 4/15 + 16/255 + 256/65535 + ... + A001146(n-1)/A051179(n) governs the proportions of terms in A001511 from row n of the array. In A001511(1..15) there are ((2/3) * 15) = ten terms from row one of the array, ((4/15) * 15) = four terms from row two, and ((16/255) * 15) = one (rounded), giving one term from row three (a 4). - Gary W. Adamson, Dec 16 2021
From Gary W. Adamson, Dec 30 2021: (Start)
Subarrays representing the number of divisors of an integer can be mapped on the table. For 60, write the odd divisors on the top row: 1, 3, 5, 15. Since 60 has 12 divisors, let the left column equal 1, 2, 4, where 4 is the highest power of 2 dividing 60. Multiplying top row terms by left column terms, we get the result:
1 3 5 15
2 6 10 30
4 12 20 60. (End)

Examples

			The table begins
   1,  3,   5,   7,   9,  11,  13,  15,  17,  19,  21,  23, ...
   2,  6,  10,  14,  18,  22,  26,  30,  34,  38,  42,  46, ...
   4, 12,  20,  28,  36,  44,  52,  60,  68,  76,  84,  92, ...
   8, 24,  40,  56,  72,  88, 104, 120, 136, 152, 168, 184, ...
  16, 48,  80, 112, 144, 176, 208, 240, 272, 304, 336, 368, ...
  32, 96, 160, 224, 288, 352, 416, 480, 544, 608, 672, 736, ...
etc.
For n = 6, we have [A002260(6), A004736(6)] = [3, 1] (i.e., 6 corresponds to location 3,1 (row,col) in above table) and A(3,1) = A000079(3-1) * A005408(1-1) = 2^2 * 1 = 4.
For n = 13, we have [A002260(13), A004736(13)] = [3, 3] (13 corresponds to location 3,3 (row,col) in above table) and A(3,3) = A000079(3-1) * A005408(3-1) = 2^2 * 5 = 20.
For n = 23, we have [A002260(23), A004736(23)] = [2, 6] (23 corresponds to location 2,6) and A(2,6) = A000079(2-1) * A005408(6-1) = 2^1 * 11 = 22.
		

Crossrefs

Transpose: A054582.
Inverse permutation: A249725.
Column 1: A000079.
Row 1: A005408.
Cf. A001511 (row index), A003602 (column index, both one-based).
Related arrays: A135765, A253551, A254053, A254055.
Cf. also permutations A246675, A246676, A249741, A249811, A249812.
Cf. A260390.

Programs

  • Maple
    seq(seq(2^(j-1)*(2*(i-j)+1),j=1..i),i=1..20); # Robert Israel, Feb 03 2015
  • Mathematica
    f[n_] := Block[{i, j}, {1}~Join~Flatten@ Last@ Reap@ For[j = 1, j <= n, For[i = j, i > 0, Sow[2^(j - i - 1)*(2 i + 1)], i--], j++]]; f@ 10 (* Michael De Vlieger, Feb 03 2015 *)
  • PARI
    a(n) = {s = ceil((1 + sqrt(1 + 8*n)) / 2); r = n - binomial(s-1, 2) - 1;k = s - r - 2; 2^r * (2 * k + 1) } \\ David A. Corneth, Feb 05 2015
  • Scheme
    (define (A135764 n) (A135764bi (A002260 n) (A004736 n)))
    (define (A135764bi row col) (* (A000079 (- row 1)) (+ -1 col col)))
    ;; Antti Karttunen, Feb 03 2015
    

Formula

From Antti Karttunen, Feb 03 2015: (Start)
A(row, col) = 2^(row-1) * ((2*col)-1) = A000079(row-1) * A005408(col-1).
A(row,col) = A064989(A135765(row,A249746(col))).
A(row,col) = A(row+1,col)/2 [discarding the topmost row and halving the rest of terms gives the array back].
A(row,col) = A(row,col+1) - A000079(row) [discarding the leftmost column and subtracting 2^{row number} from the rest of terms gives the array back].
(End)
G.f.: ((2*x+1)*Sum_{i>=0} 2^i*x^(i*(i+1)/2) + 2*(1-2*x)*Sum_{i>=0} i*x^(i*(i+1)/2) + (1-6*x)*Sum_{i>=0} x^(i*(i+1)/2) - 1 - 2*x)*x/(1-2*x)^2. These sums are related to Jacobi theta functions. - Robert Israel, Feb 03 2015

Extensions

More terms from Sean A. Irvine, Nov 23 2010
Name amended and the illustration of array in the example section transposed by Antti Karttunen, Feb 03 2015

A329050 Square array A(n,k) = prime(n+1)^(2^k), read by descending antidiagonals (0,0), (0,1), (1,0), (0,2), (1,1), (2,0), ...; Fermi-Dirac primes (A050376) in matrix form, sorted into rows by their prime divisor.

Original entry on oeis.org

2, 4, 3, 16, 9, 5, 256, 81, 25, 7, 65536, 6561, 625, 49, 11, 4294967296, 43046721, 390625, 2401, 121, 13, 18446744073709551616, 1853020188851841, 152587890625, 5764801, 14641, 169, 17, 340282366920938463463374607431768211456, 3433683820292512484657849089281, 23283064365386962890625, 33232930569601, 214358881, 28561, 289, 19
Offset: 0

Views

Author

Antti Karttunen and Peter Munn, Nov 02 2019

Keywords

Comments

This sequence is a permutation of A050376, so every positive integer is the product of a unique subset, S_factors, of its terms. If we restrict S_factors to be chosen from a subset, S_0, consisting of numbers from specified rows and/or columns of this array, there are notable sequences among those that may be generated. See the examples. Other notable sequences can be generated if we restrict the intersection of S_factors with specific rows/columns to have even cardinality. In any of the foregoing cases, the numbers in the resulting sequence form a group under the binary operation A059897(.,.).
Shares with array A246278 the property that columns grow downward by iterating A003961, and indeed, this array can be obtained from A246278 by selecting its columns 1, 2, 8, 128, ..., 2^((2^k)-1), for k >= 0.
A(n,k) is the image of the lattice point with coordinates X=n and Y=k under the inverse of the bijection f defined in the first comment of A306697. This geometric relationship can be used to construct an isomorphism from the polynomial ring GF(2)[x,y] to a ring over the positive integers, using methods similar to those for constructing A297845 and A306697. See A329329, the ring's multiplicative operator, for details.

Examples

			The top left 5 X 5 corner of the array:
  n\k |   0     1       2           3                   4
  ----+-------------------------------------------------------
   0  |   2,    4,     16,        256,              65536, ...
   1  |   3,    9,     81,       6561,           43046721, ...
   2  |   5,   25,    625,     390625,       152587890625, ...
   3  |   7,   49,   2401,    5764801,     33232930569601, ...
   4  |  11,  121,  14641,  214358881,  45949729863572161, ...
Column 0 continues as a list of primes, column 1 as a list of their squares, column 2 as a list of their 4th powers, and so on.
Every nonnegative power of 2 (A000079) is a product of a unique subset of numbers from row 0; every squarefree number (A005117) is a product of a unique subset of numbers from column 0. Likewise other rows and columns generate the sets of numbers from sequences:
Row 1:                 A000244 Powers of 3.
Column 1:              A062503 Squares of squarefree numbers.
Row 2:                 A000351 Powers of 5.
Column 2:              A113849 4th powers of squarefree numbers.
Union of rows 0 and 1:     A003586 3-smooth numbers.
Union of columns 0 and 1:  A046100 Biquadratefree numbers.
Union of row 0 / column 0: A122132 Oddly squarefree numbers.
Row 0 excluding column 0:  A000302 Powers of 4.
Column 0 excluding row 0:  A056911 Squarefree odd numbers.
All rows except 0:         A005408 Odd numbers.
All columns except 0:      A000290\{0} Positive squares.
All rows except 1:         A001651 Numbers not divisible by 3.
All columns except 1:      A252895 (have odd number of square divisors).
If, instead of restrictions on choosing individual factors of the product, we restrict the product to be of an even number of terms from each row of the array, we get A262675. The equivalent restriction applied to columns gives us A268390; applied only to column 0, we get A028260 (product of an even number of primes).
		

Crossrefs

Transpose: A329049.
Permutation of A050376.
Rows 1-4: A001146, A011764, A176594, A165425 (after the two initial terms).
Antidiagonal products: A191555.
Subtable of A182944, A242378, A246278, A329332.
A000290, A003961, A225546 are used to express relationship between terms of this sequence.
Related binary operations: A059897, A306697, A329329.
See also the table in the example section.

Programs

  • Mathematica
    Table[Prime[#]^(2^k) &[m - k + 1], {m, 0, 7}, {k, m, 0, -1}] // Flatten (* Michael De Vlieger, Dec 28 2019 *)
  • PARI
    up_to = 105;
    A329050sq(n,k) = (prime(1+n)^(2^k));
    A329050list(up_to) = { my(v = vector(up_to), i=0); for(a=0, oo, for(col=0, a, i++; if(i > up_to, return(v)); v[i] = A329050sq(col, a-col))); (v); };
    v329050 = A329050list(up_to);
    A329050(n) = v329050[1+n];
    for(n=0,up_to-1,print1(A329050(n),", ")); \\ Antti Karttunen, Nov 06 2019

Formula

A(0,k) = 2^(2^k), and for n > 0, A(n,k) = A003961(A(n-1,k)).
A(n,k) = A182944(n+1,2^k).
A(n,k) = A329332(2^n,2^k).
A(k,n) = A225546(A(n,k)).
A(n,k+1) = A000290(A(n,k)) = A(n,k)^2.

Extensions

Example annotated for clarity by Peter Munn, Feb 12 2020

A060803 a(n) = Sum_{k=0..n} 2^(2^k).

Original entry on oeis.org

2, 6, 22, 278, 65814, 4295033110, 18446744078004584726, 340282366920938463481821351509772796182, 115792089237316195423570985008687907853610267032561502502939405359422902436118
Offset: 0

Views

Author

Varol Akman (akman(AT)cs.bilkent.edu.tr), Apr 28 2001

Keywords

Comments

Partial sums of A001146.
Number of Boolean functions with up to n arguments. - Paul Tarau (paul.tarau(AT)gmail.com), Jun 06 2008

Examples

			a(3) = 278 because a(3) = 2^2^0 + 2^2^1 + 2^2^2 + 2^2^3 = 2 + 4 + 16 + 256.
		

Crossrefs

Programs

  • Haskell
    -- code generating the infinite sequence:
    scanl (+) 2 (map (\x->2^2^x) [1..]) - Paul Tarau (paul.tarau(AT)gmail.com), Jun 06 2008
    
  • Mathematica
    Accumulate[2^(2^Range[0,10])] (* Harvey P. Dale, Sep 25 2023 *)
  • PARI
    { for (n=0, 11, write("b060803.txt", n, " ", sum(k=0, n, 2^(2^k))); ) } \\ Harry J. Smith, Jul 12 2009

Formula

a(0) = 2 and a(n) - a(n-1) = 2^2^n, n > 0.

Extensions

More terms from Benoit Cloitre, May 13 2002
Edited by N. J. A. Sloane, Jun 07 2008

A077585 a(n) = 2^(2^n - 1) - 1.

Original entry on oeis.org

0, 1, 7, 127, 32767, 2147483647, 9223372036854775807, 170141183460469231731687303715884105727, 57896044618658097711785492504343953926634992332820282019728792003956564819967
Offset: 0

Views

Author

Henry Bottomley, Nov 07 2002

Keywords

Comments

a(2), a(3), a(5) and a(7) are prime; a(11) is not.
Let S be a set of n elements. First we perform a set partition on S. Let SU be the set of all nonempty subsets of S. As is well known, 2^n - 1 is the number of nonempty subsets of a set with n elements (see A000225). That is, 2^n - 1 is the number of elements |SU| of SU. In the second step, we select k elements from SU. We want to know how many different selections are possible. Let W be the resulting set of selections formed from SU. Then the number of elements |W| of W is |W| = Sum_{k=1..2^n-1} binomial(2^n-1,k) = 2^(2^n-1) - 1 = A077585. Example: |W(n)| = a(n=2) = 7, because W = {[[1, 2]], [[1]], [[1, 2], [1]], [[1, 2], [2], [1]], [[1, 2], [2]], [[2], [1]], [[2]]}. - Thomas Wieder, Nov 08 2007
These are so-called double Mersenne numbers, sometimes denoted MM(n) where M = A000225. MM(n) is prime iff M(n) = 2^n-1 is a Mersenne prime exponent (A000043), which isn't possible unless n itself is also in A000043. Primes of this form are called double Mersenne primes MM(p). For all Mersenne exponents between 7 and 61, factors of MM(p) are known. MM(61) is far too large for any currently known primality test, but a distributed search for factors of this and other MM(p) is ongoing, see the doublemersenne.org web site. - M. F. Hasler, Mar 05 2020
This is a strong divisibility sequence, that is, gcd(a(n),a(m)) = a(gcd(n,m)) for n, m >= 1. - Peter Bala, Dec 05 2022

Examples

			a(5) = 2^(2^5 - 1) - 1 = 2^31 - 1 = 2147483647.
		

Crossrefs

Programs

  • Maple
    a:= n-> 2^(2^n-1)-1:
    seq(a(n), n=0..8); # Thomas Wieder, Nov 08 2007
  • Mathematica
    2^(2^Range[0, 9] - 1) - 1 (* Vladimir Joseph Stephan Orlovsky, Jun 22 2011 *)
  • PARI
    a(n)=if(n<1,0,-1+2*(1+a(n-1))^2)
    
  • PARI
    apply( {A077585(n)=1<<(1<M. F. Hasler, Mar 05 2020
    
  • Python
    def A077585(n): return (1<<(1<Chai Wah Wu, Mar 14 2023

Formula

a(n) = A000225(A000225(n)).
a(n) = A058891(n+1) - 1. - corrected by Maurizio De Leo, Feb 25 2015
a(n) = (A001146(n) - 2)/2.
a(n) = A056220(1+a(n-1)).
a(n) = Sum_{k=1..2^n-1} binomial(2^n-1,k). - Thomas Wieder, Nov 08 2007
a(n) = 2*a(n-1)^2 + 4*a(n-1) + 1. - Roderick MacPhee, Oct 05 2012

Extensions

Corrected by Lekraj Beedassy, Jan 02 2007

A326941 Number of T_0 sets of subsets of {1..n}.

Original entry on oeis.org

2, 4, 14, 224, 64210, 4294322204, 18446744009291513774, 340282366920938463075992982725615419816, 115792089237316195423570985008687907843742078391854287068939455414919611614210
Offset: 0

Views

Author

Gus Wiseman, Aug 07 2019

Keywords

Comments

The dual of a multiset partition has, for each vertex, one block consisting of the indices (or positions) of the blocks containing that vertex, counted with multiplicity. For example, the dual of {{1,2},{2,3}} is {{1},{1,2},{2}}. The T_0 condition means that the dual is strict (no repeated edges).

Examples

			The a(0) = 2 through a(2) = 14 sets of subsets:
  {}    {}        {}
  {{}}  {{}}      {{}}
        {{1}}     {{1}}
        {{},{1}}  {{2}}
                  {{},{1}}
                  {{},{2}}
                  {{1},{2}}
                  {{1},{1,2}}
                  {{2},{1,2}}
                  {{},{1},{2}}
                  {{},{1},{1,2}}
                  {{},{2},{1,2}}
                  {{1},{2},{1,2}}
                  {{},{1},{2},{1,2}}
		

Crossrefs

The non-T_0 version is A001146.
The covering case is A326939.
The case without empty edges is A326940.
The unlabeled version is A326949.

Programs

  • Mathematica
    dual[eds_]:=Table[First/@Position[eds,x],{x,Union@@eds}];
    Table[Length[Select[Subsets[Subsets[Range[n]]],UnsameQ@@dual[#]&]],{n,0,3}]

Formula

a(n) = 2 * A326940(n).
Binomial transform of A326939.

Extensions

a(5)-a(8) from Andrew Howroyd, Aug 14 2019

A384247 The number of integers from 1 to n whose largest divisor that is an infinitary divisor of n is 1.

Original entry on oeis.org

1, 1, 2, 3, 4, 2, 6, 4, 8, 4, 10, 6, 12, 6, 8, 15, 16, 8, 18, 12, 12, 10, 22, 8, 24, 12, 18, 18, 28, 8, 30, 16, 20, 16, 24, 24, 36, 18, 24, 16, 40, 12, 42, 30, 32, 22, 46, 30, 48, 24, 32, 36, 52, 18, 40, 24, 36, 28, 58, 24, 60, 30, 48, 48, 48, 20, 66, 48, 44, 24
Offset: 1

Views

Author

Amiram Eldar, May 23 2025

Keywords

Comments

Analogous to A047994, as A064380 is analogous to A116550.

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := p^e*(1 - 1/p^(2^(IntegerExponent[e, 2]))); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n)); n * prod(i = 1, #f~, (1 - 1/f[i,1]^(1 << valuation(f[i,2], 2))));}

Formula

Multiplicative with a(p^e) = p^e * (1 - 1/p^A006519(e)).
a(n) >= A091732(n), with equality if and only if n is in A138302.
a(n) <= A047994(n), with equality if and only if n is in A138302.
a(n) >= A000010(n), with equality if and only if n is an exponentially odd number (A268335).
a(n) is odd if and only if n = 1 or 2^(2^k) for k >= 0 (A001146). a(2^(2^k)) = 2^(2^k)-1.
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = Product_{p prime} f(1/p) = 0.66718130416373472394..., and f(x) = 1 - (1-x)*Sum_{k>=1} x^(2^k)/(1-x^(2^k)).

A003181 Number of P-equivalence classes of nondegenerate Boolean functions of n variables.

Original entry on oeis.org

2, 2, 8, 68, 3904, 37329264, 25626412300941056, 67516342973185974302549277749387264, 2871827610052485009904013737758920847602293486924450772201235462734479360
Offset: 0

Views

Author

Keywords

Comments

Also the number of non-isomorphic sets of subsets of {1..n} with union {1..n}. - Gus Wiseman, Aug 05 2019

Examples

			From _Gus Wiseman_, Aug 05 2019: (Start)
Non-isomorphic representatives of the a(0) = 2 through a(2) = 8 sets of subsets:
  {}    {{1}}     {{1,2}}
  {{}}  {{},{1}}  {{1},{2}}
                  {{},{1,2}}
                  {{2},{1,2}}
                  {{},{1},{2}}
                  {{},{2},{1,2}}
                  {{1},{2},{1,2}}
                  {{},{1},{2},{1,2}}
(End)
		

References

  • S. Muroga, Threshold Logic and Its Applications. Wiley, NY, 1971, p. 38 and 214.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Maple
    b:= proc(n, i, l) `if`(n=0, 2^(w-> add(mul(2^igcd(t, l[h]),
          h=1..nops(l)), t=1..w)/w)(ilcm(l[])), `if`(i<1, 0,
          add(b(n-i*j, i-1, [l[], i$j])/j!/i^j, j=0..n/i)))
        end:
    a:= n-> `if`(n=0, 2, b(n$2, [])-b(n-1$2, [])):
    seq(a(n), n=0..8);  # Alois P. Heinz, Aug 14 2019
  • Mathematica
    b[n_, i_, l_] := If[n == 0, 2^Function[w, Sum[Product[2^GCD[t, l[[h]]], {h, 1, Length[l]}], {t, 1, w}]/w][If[l == {}, 1, LCM @@ l]], If[i < 1, 0, Sum[b[n - i*j, i - 1, Join[l, Table[i, {j}]]]/j!/i^j, {j, 0, n/i}]]];
    a[n_] := If[n == 0, 2, b[n, n, {}] - b[n - 1, n - 1, {}]];
    a /@ Range[0, 8] (* Jean-François Alcover, Apr 11 2020, after Alois P. Heinz *)

Formula

a(n) = A003180(n) - A003180(n-1), for n >= 1. - Christian Sievers, Jul 22 2016
a(n) = 2 * A055621(n). - Gus Wiseman, Aug 05 2019

Extensions

More terms from Christian Sievers, Jul 22 2016
Definition clarified by Ivo Timoteo, Mar 14 2017
Previous Showing 11-20 of 116 results. Next