A028362 -id:A028362 - cp's OEIS frontend

A048675 If n = p_i^e_i * ... * p_k^e_k, p_i < ... < p_k primes (with p_i = prime(i)), then a(n) = (1/2) * (e_i * 2^i + ... + e_k * 2^k).

0, 1, 2, 2, 4, 3, 8, 3, 4, 5, 16, 4, 32, 9, 6, 4, 64, 5, 128, 6, 10, 17, 256, 5, 8, 33, 6, 10, 512, 7, 1024, 5, 18, 65, 12, 6, 2048, 129, 34, 7, 4096, 11, 8192, 18, 8, 257, 16384, 6, 16, 9, 66, 34, 32768, 7, 20, 11, 130, 513, 65536, 8, 131072, 1025, 12, 6, 36, 19

Offset: 1

Antti Karttunen, Jul 14 1999

nonn

The original motivation for this sequence was to encode the prime factorization of n in the binary representation of a(n), each such representation being unique as long as this map is restricted to A005117 (squarefree numbers, resulting a permutation of nonnegative integers A048672) or any of its subsequence, resulting an injective function like A048623 and A048639.
However, also the restriction to A260443 (not all terms of which are squarefree) results a permutation of nonnegative integers, namely A001477, the identity permutation.
When a polynomial with nonnegative integer coefficients is encoded with the prime factorization of n (e.g., as in A206296, A260443), then a(n) gives the evaluation of that polynomial at x=2.
The primitive completely additive integer sequence that satisfies a(n) = a(A225546(n)), n >= 1. By primitive, we mean that if b is another such sequence, then there is an integer k such that b(n) = k * a(n) for all n >= 1. - Peter Munn, Feb 03 2020
If the binary rank of an integer partition y is given by Sum_i 2^(y_i-1), and the Heinz number is Product_i prime(y_i), then a(n) is the binary rank of the integer partition with Heinz number n. Note the function taking a set s to Sum_i 2^(s_i-1) is the inverse of A048793 (binary indices), and the function taking a multiset m to Product_i prime(m_i) is the inverse of A112798 (prime indices). - Gus Wiseman, May 22 2024

			From _Gus Wiseman_, May 22 2024: (Start)
The A018819(7) = 6 cases of binary rank 7 are the following, together with their prime indices:
   30: {1,2,3}
   40: {1,1,1,3}
   54: {1,2,2,2}
   72: {1,1,1,2,2}
   96: {1,1,1,1,1,2}
  128: {1,1,1,1,1,1,1}
(End)

Antti Karttunen, Table of n, a(n) for n = 1..1024
Hans Havermann, Factorization of the first 10000 terms, in format [[primes], [exponents]]

Row 2 of A104244.
Similar logarithmic functions: A001414, A056239, A090880, A289506, A293447.
Left inverse of the following sequences: A000079, A019565, A038754, A068911, A134683, A260443, A332824.
A003961, A028234, A032742, A055396, A064989, A067029, A225546, A297845 are used to express relationship between terms of this sequence.
Cf. also A048623, A048676, A099884, A277896 and tables A277905, A285325.
Cf. A297108 (Möbius transform), A332813 and A332823 [= a(n) mod 3].
Pairs of sequences (f,g) that satisfy a(f(n)) = g(n), possibly with offset change: (A000203,A331750), (A005940,A087808), (A007913,A248663), (A007947,A087207), (A097248,A048675), (A206296,A000129), (A248692,A056239), (A283477,A005187), (A284003,A006068), (A285101,A028362), (A285102,A068052), (A293214,A001065), (A318834,A051953), (A319991,A293897), (A319992,A293898), (A320017,A318674), (A329352,A069359), (A332461,A156552), (A332462,A156552), (A332825,A000010) and apparently (A163511,A135529).
See comments/formulas in A277333, A331591, A331740 giving their relationship to this sequence.
The formula section details how the sequence maps the terms of A329050, A329332.
A277892, A322812, A322869, A324573, A324575 give properties of the n-th term of this sequence.
The term k appears A018819(k) times.
The inverse transformation is A019565 (Heinz number of binary indices).
The version for distinct prime indices is A087207.
Numbers k such that a(k) is prime are A277319, counts A372688.
Grouping by image gives A277905.
A014499 lists binary indices of prime numbers.
A061395 gives greatest prime index, least A055396.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.
Binary indices:
- listed A048793, sum A029931
- reversed A272020
- opposite A371572, sum A230877
- length A000120, complement A023416
- min A001511, opposite A000012
- max A070939, opposite A070940
- complement A368494, sum A359400
- opposite complement A371571, sum A359359
Cf. A000720, A035100, A071814, A096111, A225620, A243055, A304818, A355536, A358136, A372890.
Cf. A087207, A097248.

nthprime := proc(n) local i; if(isprime(n)) then for i from 1 to 1000000 do if(ithprime(i) = n) then RETURN(i); fi; od; else RETURN(0); fi; end; # nthprime(2) = 1, nthprime(3) = 2, nthprime(5) = 3, etc. - this is also A049084.
A048675 := proc(n) local s,d; s := 0; for d in ifactors(n)[ 2 ] do s := s + d[ 2 ]*(2^(nthprime(d[ 1 ])-1)); od; RETURN(s); end;
# simpler alternative
f:= n -> add(2^(numtheory:-pi(t[1])-1)*t[2], t=ifactors(n)[2]):
map(f, [$1..100]); # Robert Israel, Oct 10 2016

a[1] = 0; a[n_] := Total[ #[[2]]*2^(PrimePi[#[[1]]]-1)& /@ FactorInteger[n] ]; Array[a, 100] (* Jean-François Alcover, Mar 15 2016 *)

a(n) = my(f = factor(n)); sum(k=1, #f~, f[k,2]*2^primepi(f[k,1]))/2; \\ Michel Marcus, Oct 10 2016

\\ The following program reconstructs terms (e.g. for checking purposes) from the factorization file prepared by Hans Havermann:
v048675sigs = readvec("a048675.txt");
A048675(n) = if(n<=2,n-1,my(prsig=v048675sigs[n],ps=prsig[1],es=prsig[2]); prod(i=1,#ps,ps[i]^es[i])); \\ Antti Karttunen, Feb 02 2020

from sympy import factorint, primepi
def a(n):
    if n==1: return 0
    f=factorint(n)
    return sum([f[i]*2**(primepi(i) - 1) for i in f])
print([a(n) for n in range(1, 51)]) # Indranil Ghosh, Jun 19 2017

a(1) = 0, a(n) = 1/2 * (e1*2^i1 + e2*2^i2 + ... + ez*2^iz) if n = p_{i1}^e1*p_{i2}^e2*...*p_{iz}^ez, where p_i is the i-th prime. (e.g. p_1 = 2, p_2 = 3).
Totally additive with a(p^e) = e * 2^(PrimePi(p)-1), where PrimePi(n) = A000720(n). [Missing factor e added to the comment by Antti Karttunen, Jul 29 2015]
From Antti Karttunen, Jul 29 2015: (Start)
a(1) = 0; for n > 1, a(n) = 2^(A055396(n)-1) + a(A032742(n)). [Where A055396(n) gives the index of the smallest prime dividing n and A032742(n) gives the largest proper divisor of n.]
a(1) = 0; for n > 1, a(n) = (A067029(n) * (2^(A055396(n)-1))) + a(A028234(n)).
Other identities. For all n >= 0:
a(A019565(n)) = n.
a(A260443(n)) = n.
a(A206296(n)) = A000129(n).
a(A005940(n+1)) = A087808(n).
a(A007913(n)) = A248663(n).
a(A007947(n)) = A087207(n).
a(A283477(n)) = A005187(n).
a(A284003(n)) = A006068(n).
a(A285101(n)) = A028362(1+n).
a(A285102(n)) = A068052(n).
Also, it seems that a(A163511(n)) = A135529(n) for n >= 1. (End)
a(1) = 0, a(2n) = 1+a(n), a(2n+1) = 2*a(A064989(2n+1)). - Antti Karttunen, Oct 11 2016
From Peter Munn, Jan 31 2020: (Start)
a(n^2) = a(A003961(n)) = 2 * a(n).
a(A297845(n,k)) = a(n) * a(k).
a(n) = a(A225546(n)).
a(A329332(n,k)) = n * k.
a(A329050(n,k)) = 2^(n+k).
(End)
From Antti Karttunen, Feb 02-25 2020, Feb 01 2021: (Start)
a(n) = Sum_{d|n} A297108(d) = Sum_{d|A225546(n)} A297108(d).
a(n) = a(A097248(n)).
For n >= 2:
A001221(a(n)) = A322812(n), A001222(a(n)) = A277892(n).
A000203(a(n)) = A324573(n), A033879(a(n)) = A324575(n).
For n >= 1, A331750(n) = a(A000203(n)).
For n >= 1, the following chains hold:
A293447(n) >= a(n) >= A331740(n) >= A331591(n).
a(n) >= A087207(n) >= A248663(n).
(End)
a(n) = A087207(A097248(n)). - Flávio V. Fernandes, Jul 16 2025

Entry revised by Antti Karttunen, Jul 29 2015

More linking formulas added by Antti Karttunen, Apr 18 2017

A005329 a(n) = Product_{i=1..n} (2^i - 1). Also called 2-factorial numbers.

Original entry on oeis.org

1, 1, 3, 21, 315, 9765, 615195, 78129765, 19923090075, 10180699028325, 10414855105976475, 21319208401933844325, 87302158405919092510875, 715091979502883286756577125, 11715351900195736886933003038875, 383876935713713710574133710574817125

Offset: 0

N. J. A. Sloane

Conjecture: this sequence is the inverse binomial transform of A075272 or, equivalently, the inverse binomial transform of the BinomialMean transform of A075271. - John W. Layman, Sep 12 2002
To win a game, you must flip n+1 heads in a row, where n is the total number of tails flipped so far. Then the probability of winning for the first time after n tails is A005329 / A006125. The probability of having won before n+1 tails is A114604 / A006125. - Joshua Zucker, Dec 14 2005
Number of upper triangular n X n (0,1)-matrices with no zero rows. - Vladeta Jovovic, Mar 10 2008
Equals the q-Fibonacci series for q = (-2), and the series prefaced with a 1: (1, 1, 1, 3, 21, ...) dot (1, -2, 4, -8, ...) if n is even, and (-1, 2, -4, 8, ...) if n is odd. For example, a(3) = 21 = (1, 1, 1, 3) dot (-1, 2, -4, 8) = (-1, 2, -4, 24) and a(4) = 315 = (1, 1, 1, 3, 21) dot (1, -2, 4, -8 16) = (1, -2, 4, -24, 336). - Gary W. Adamson, Apr 17 2009
Number of chambers in an A_n(K) building where K=GF(2) is the field of two elements. This is also the number of maximal flags in an n-dimensional vector space over a field of two elements. - Marcos Spreafico, Mar 22 2012
Given probability p = 1/2^n that an outcome will occur at the n-th stage of an infinite process, then starting at n=1, A114604(n)/A006125(n+2) = 1-a(n)/A006125(n+1) is the probability that the outcome has occurred up to and including the n-th iteration. The limiting ratio is 1-A048651 ~ 0.7112119. These observations are a more formal and generalized statement of Joshua Zucker's Dec 14, 2005 comment. - Bob Selcoe, Mar 02 2016
Also the number of dominating sets in the n-triangular honeycomb rook graph. - Eric W. Weisstein, Jul 14 2017
Empirical: Letting Q denote the Hall-Littlewood Q basis of the symmetric functions over the field of fractions of the univariate polynomial ring in t over the field of rational numbers, and letting h denote the complete homogeneous basis, a(n) is equal to the absolute value of 2^A000292(n) times the coefficient of h_{1^(n*(n+1)/2)} in Q_{(n, n-1, ..., 1)} with t evaluated at 1/2. - John M. Campbell, Apr 30 2018
The series f(x) = Sum_{n>=0} x^(2^n-1)/a(n) satisfies f'(x) = f(x^2), f(0) = 1. - Lucas Larsen, Jan 05 2022

			G.f. = 1 + x + 3*x^2 + 21*x^3 + 315*x^4 + 9765*x^5 + 615195*x^6 + 78129765*x^7 + ...

Annie Cuyt, Vigdis Brevik Petersen, Brigitte Verdonk, Haakon Waadeland, and William B. Jones, Handbook of continued fractions for special functions, Springer, New York, 2008. (see 19.2.1)
Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, p. 358.
Mark Ronan, Lectures on Buildings (Perspectives in Mathematics; Vol. 7), Academic Press Inc., 1989.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 0..50
E. Andresen and K. Kjeldsen, On certain subgraphs of a complete transitively directed graph, Discrete Math., Vol. 14, No. 2 (1976), pp. 103-119.
Paul Barry, Centered polygon numbers, heptagons and nonagons, and the Robbins numbers, arXiv:2104.01644 [math.CO], 2021.
Geoffrey Critzer, Combinatorics of Vector Spaces over Finite Fields, Master's thesis, Emporia State University, 2018.
Hsien-Kuei Hwang, Emma Yu Jin, and Michael J. Schlosser, Asymptotics and statistics on Fishburn Matrices: dimension distribution and a conjecture of Stoimenow, arXiv:2012.13570 [math.CO], 2020.
Kent E. Morrison, Integer Sequences and Matrices Over Finite Fields, Journal of Integer Sequences, Vol. 9 (2006), Article 06.2.1.
Eric Weisstein's World of Mathematics, Dominating Set.
Eric Weisstein's World of Mathematics, q-Factorial.
Index entries for sequences related to factorial numbers.

Cf. A000225, A005321, A006125, A114604, A006088, A028362, A027871 (3-fac), A027872 (5-fac), A027873 (6-fac), A048651, A048652, A075271, A075272, A032085, A122746.
Cf. A048651, A079555, A152476 (inverse binomial transform).
Column q=2 of A069777.

List([0..15],n->Product([1..n],i->2^i-1)); # Muniru A Asiru, May 18 2018

[1] cat [&*[ 2^k-1: k in [1..n] ]: n in [1..16]]; // Vincenzo Librandi, Dec 24 2015

A005329 := proc(n) option remember; if n<=1 then 1 else (2^n-1)*procname(n-1); end if; end proc: seq(A005329(n), n=0..15);

a[0] = 1; a[n_] := a[n] = (2^n-1)*a[n-1]; a /@ Range[0,14] (* Jean-François Alcover, Apr 22 2011 *)
FoldList[Times, 1, 2^Range[15] - 1] (* Harvey P. Dale, Dec 21 2011 *)
Table[QFactorial[n, 2], {n, 0, 14}] (* Arkadiusz Wesolowski, Oct 30 2012 *)
QFactorial[Range[0, 10], 2] (* Eric W. Weisstein, Jul 14 2017 *)
a[ n_] := If[ n < 0, 0, (-1)^n QPochhammer[ 2, 2, n]]; (* Michael Somos, Jan 28 2018 *)

a(n)=polcoeff(sum(m=0,n,2^(m*(m+1)/2)*x^m/prod(k=0,m,1+2^k*x+x*O(x^n))),n) \\ Paul D. Hanna, Sep 17 2009

Dx(n,F)=local(D=F);for(i=1,n,D=deriv(D));D
a(n)=local(A=1+x+x*O(x^n));for(i=1,n,A=1+sum(k=1,n,x^k/k!*Dx(k,x*A+x*O(x^n) ))); polcoeff(A,n) \\ Paul D. Hanna, Apr 21 2012

{a(n) = if( n<0, 0, prod(k=1, n, 2^k - 1))}; /* Michael Somos, Jan 28 2018 */

{a(n) = if( n<0, 0, (-1)^n * sum(k=0, n+1, (-1)^k * 2^(k*(k+1)/2) * prod(j=1, k, (2^(n+1-j) - 1) / (2^j - 1))))}; /* Michael Somos, Jan 28 2018 */

a(n)/2^(n*(n+1)/2) -> c = 0.2887880950866024212788997219294585937270... (see A048651, A048652).
From Paul D. Hanna, Sep 17 2009: (Start)
G.f.: Sum_{n>=0} 2^(n*(n+1)/2) * x^n / (Product_{k=0..n} (1+2^k*x)).
Compare to: 1 = Sum_{n>=0} 2^(n*(n+1)/2) * x^n/(Product_{k=1..n+1} (1+2^k*x)). (End)
G.f. satisfies: A(x) = 1 + Sum_{n>=1} x^n/n! * d^n/dx^n x*A(x). - Paul D. Hanna, Apr 21 2012
a(n) = 2^(binomial(n+1,2))*(1/2; 1/2){n}, where (a;q){n} is the q-Pochhammer symbol. - G. C. Greubel, Dec 23 2015
a(n) = Product_{i=1..n} A000225(i). - Michel Marcus, Dec 27 2015
From Peter Bala, Nov 10 2017: (Start)
O.g.f. as a continued fraction of Stieltjes' type: A(x) = 1/(1 - x/(1 - 2*x/(1 - 6*x/(1 - 12*x/(1 - 28*x/(1 - 56*x/(1 - ... -(2^n - 2^floor(n/2))*x/(1 - ... )))))))) (follows from Heine's continued fraction for the ratio of two q-hypergeometric series at q = 2. See Cuyt et al. 19.2.1).
A(x) = 1/(1 + x - 2*x/(1 - (2 - 1)^2*x/(1 + x - 2^3*x/(1 - (2^2 - 1)^2*x/(1 + x - 2^5*x/(1 - (2^3 - 1)^2*x/(1 + x - 2^7*x/(1 - (2^4 - 1)^2*x/(1 + x - ... ))))))))). (End)
0 = a(n)*(a(n+1) - a(n+2)) + 2*a(n+1)^2 for all n>=0. - Michael Somos, Feb 23 2019
From Amiram Eldar, Feb 19 2022: (Start)
Sum_{n>=0} 1/a(n) = A079555.
Sum_{n>=0} (-1)^n/a(n) = A048651. (End)

Better definition from Leslie Ann Goldberg (leslie(AT)dcs.warwick.ac.uk), Dec 11 1999

A053632 Irregular triangle read by rows giving coefficients in expansion of Product_{k=1..n} (1 + x^k).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 1, 1, 1, 1, 1, 1, 2, 2, 3, 3, 3, 3, 3, 3, 2, 2, 1, 1, 1, 1, 1, 1, 2, 2, 3, 4, 4, 4, 5, 5, 5, 5, 4, 4, 4, 3, 2, 2, 1, 1, 1, 1, 1, 1, 2, 2, 3, 4, 5, 5, 6, 7, 7, 8, 8, 8, 8, 8, 7, 7, 6, 5, 5, 4, 3, 2, 2, 1, 1, 1, 1, 1, 1, 2, 2, 3, 4

Offset: 0

N. J. A. Sloane, Mar 22 2000

Or, triangle T(n,k) read by rows, giving number of subsets of {1,2,...,n} with sum k. - Roger CUCULIERE (cuculier(AT)imaginet.fr), Nov 19 2000
Row n consists of A000124(n) terms. These are also the successive vectors (their nonzero elements) when one starts with the infinite vector (of zeros) with 1 inserted somewhere and then shifts it one step (right or left) and adds to the original, then shifts the result two steps and adds, three steps and adds, etc. - Antti Karttunen, Feb 13 2002
T(n,k) = number of partitions of k into distinct parts <= n. Triangle of distribution of Wilcoxon's signed rank statistic. - Mitch Harris, Mar 23 2006
T(n,k) = number of binary words of length n in which the sum of the positions of the 0's is k. Example: T(4,5)=2 because we have 0110 (sum of the positions of the 0's is 1+4=5) and 1001 (sum of the positions of the 0's is 2+3=5). - Emeric Deutsch, Jul 23 2006
A fair coin is flipped n times. You receive i dollars for a "success" on the i-th flip, 1<=i<=n. T(n,k)/2^n is the probability that you will receive exactly k dollars. Your expectation is n(n+1)/4 dollars. - Geoffrey Critzer, May 16 2010
From Gus Wiseman, Jan 02 2023: (Start)
With offset 1, also the number of integer compositions of n whose partial sums add up to k for k = n..n(n+1)/2. For example, row n = 6 counts the following compositions:
  6 15 24 33  42  51  141  231  321  411  1311 2211  3111  12111 21111 111111
          114 123 132 222  312  1131 1221 2121 11121 11211
                  213 1113 1122 1212 2112 1111
(End)

			Triangle begins:
  1;
  1, 1;
  1, 1, 1, 1;
  1, 1, 1, 2, 1, 1, 1;
  1, 1, 1, 2, 2, 2, 2, 2, 1, 1, 1;
  1, 1, 1, 2, 2, 3, 3, 3, 3, 3, 3, 2, 2, 1, 1, 1;
  1, 1, 1, 2, 2, 3, 4, 4, 4, 5, 5, 5, 5, 4, 4, 4, 3, 2, 2, 1, 1, 1;
  ...
Row n = 4 counts the following binary words, where k = sum of positions of zeros:
  1111  0111  1011  0011  0101  0110  0001  0010  0100  1000  0000
                    1101  1110  1001  1010  1100
Row n = 5 counts the following strict partitions of k with all parts <= n (0 is the empty partition):
  0  1  2  3  4  5  42  43  53  54  532  542  543  5431 5432 54321
           21 31 32 51  52  431 432 541  5321 5421
                 41 321 421 521 531 4321

A. V. Yurkin, New binomial and new view on light theory, (book), 2013, 78 pages, no publisher listed.

Alois P. Heinz, Rows n = 0..40, flattened
S. R. Finch, Signum equations and extremal coefficients.
Steven R. Finch, Signum equations and extremal coefficients, February 7, 2009. [Cached copy, with permission of the author]
FindStat - Combinatorial Statistic Finder, The major index of an integer composition
Alexander Rosa and Štefan Znám, A combinatorial problem in the theory of congruences. (Russian), Mat.-Fys. Casopis Sloven. Akad. Vied 15 1965 49-59. [Annotated scanned copy] See Table 1.
F. Wilcoxon, Individual Comparisons by Ranking Methods, Biometrics Bulletin, v. 1, no. 6 (1945), pp. 80-83.
A. V. Yurkin, On similarity of systems of geometrical and arithmetic triangles, in Mathematics, Computing, Education Conference XIX, 2012.
A. V. Yurkin, New view on the diffraction discovered by Grimaldi and Gaussian beams, arXiv preprint arXiv:1302.6287 [physics.optics], 2013.
A. V. Yurkin, About the evident description of distribution of beams and "wavy geometrical trajectories" in long thin pipes, 2014 (original in Russian).
A. V. Yurkin, Symmetric triangle of Pascal and non-linear arithmetic parallelepiped, Book Manuscript, Research Gate 2015.

Cf. A053633, A068009.
Rows reduced modulo 2 and interpreted as binary numbers: A068052, A068053. Rows converge towards A000009.
Row sums give A000079.
Cf. A001788, A028362.
Cf. A285101 (multiplicative encoding of each row), A285103 (number of odd terms on row n), A285105 (number of even terms).
Row lengths are A000124.
A reciprocal version is (A033999, A219977, A291983, A291984, A291985, ...).
A negative version is A231599.
A version for partitions is A358194, reversed partitions A264034.
Cf. A025591, A029931, A063865, A152947, A318283, A359042.

with(gfun,seriestolist); map(op,[seq(seriestolist(series(mul(1+(z^i), i=1..n),z,binomial(n+1,2)+1)), n=0..10)]); # Antti Karttunen, Feb 13 2002
# second Maple program:
g:= proc(n) g(n):= `if`(n=0, 1, expand(g(n-1)*(1+x^n))) end:
T:= n-> seq(coeff(g(n), x, k), k=0..degree(g(n))):
seq(T(n), n=0..10);  # Alois P. Heinz, Nov 19 2012

Table[CoefficientList[ Series[Product[(1 + t^i), {i, 1, n}], {t, 0, 100}], t], {n, 0, 8}] // Grid (* Geoffrey Critzer, May 16 2010 *)

From Mitch Harris, Mar 23 2006: (Start)
T(n,k) = T(n-1, k) + T(n-1, k-n), T(0,0)=1, T(0,k) = 0, T(n,k) = 0 if k < 0 or k > (n+1 choose 2).
G.f.: (1+x)*(1+x^2)*...*(1+x^n). (End)
Sum_{k>=0} k * T(n,k) = A001788(n). - Alois P. Heinz, Feb 09 2017
max_{k>=0} T(n,k) = A025591(n). - Alois P. Heinz, Jan 20 2023

A374848 Obverse convolution A000045**A000045; see Comments.

Original entry on oeis.org

0, 1, 2, 16, 162, 3600, 147456, 12320100, 2058386904, 701841817600, 488286500625000, 696425232679321600, 2038348954317776486400, 12259459134020160144810000, 151596002479762016373851690400, 3855806813438155578522841251840000

Offset: 0

Clark Kimberling, Jul 31 2024

nonn

The obverse convolution of sequences
  s = (s(0), s(1), ...) and t = (t(0), t(1), ...)
is introduced here as the sequence s**t given by
s**t(n) = (s(0)+t(n)) * (s(1)+t(n-1)) * ... * (s(n)+t(0)).
Swapping * and + in the representation s(0)*t(n) + s(1)*t(n-1) + ... + s(n)*t(0)
of ordinary convolution yields s**t.
If x is an indeterminate or real (or complex) variable, then for every sequence t of real (or complex) numbers, s**t is a sequence of polynomials p(n) in x, and the zeros of p(n) are the numbers -t(0), -t(1), ..., -t(n).
Following are abbreviations in the guide below for triples (s, t, s**t):
F = (0,1,1,2,3,5,...) = A000045, Fibonacci numbers
L = (2,1,3,4,7,11,...) = A000032, Lucas numbers
P = (2,3,5,7,11,...) = A000040, primes
T = (1,3,6,10,15,...) = A000217, triangular numbers
C = (1,2,6,20,70, ...) = A000984, central binomial coefficients
LW = (1,3,4,6,8,9,...) = A000201, lower Wythoff sequence
UW = (2,5,7,10,13,...) = A001950, upper Wythoff sequence
[ ] = floor
In the guide below, sequences s**t are identified with index numbers Axxxxxx; in some cases, s**t and Axxxxxx differ in one or two initial terms.
Table 1. s = A000012 = (1,1,1,1...) = (1);
t = A000012; 1                 s**t = A000079; 2^(n+1)
t = A000027; n                 s**t = A000142; (n+1)!
t = A000040, P                 s**t = A054640
t = A000040, P           (1/3) s**t = A374852
t = A000079, 2^n               s**t = A028361
t = A000079, 2^n         (1/3) s**t = A028362
t = A000045, F                 s**t = A082480
t = A000032, L                 s**t = A374890
t = A000201, LW                s**t = A374860
t = A001950, UW                s**t = A374864
t = A005408, 2*n+1             s**t = A000165, 2^n*n!
t = A016777, 3*n+1             s**t = A008544
t = A016789, 3*n+2             s**t = A032031
t = A000142, n!                s**t = A217757
t = A000051, 2^n+1             s**t = A139486
t = A000225, 2^n-1             s**t = A006125
t = A032766, [3*n/2]           s**t = A111394
t = A034472, 3^n+1             s**t = A153280
t = A024023, 3^n-1             s**t = A047656
t = A000217, T                 s**t = A128814
t = A000984, C                 s**t = A374891
t = A279019, n^2-n             s**t = A130032
t = A004526, 1+[n/2]           s**t = A010551
t = A002264, 1+[n/3]           s**t = A264557
t = A002265, 1+[n/4]           s**t = A264635
Sequences (c)**L, for c=2..4: A374656 to A374661
Sequences (c)**F, for c=2..6: A374662, A374662, A374982 to A374855
The obverse convolutions listed in Table 1 are, trivially, divisibility sequences. Likewise, if s = (-1,-1,-1,...) instead of s = (1,1,1,...), then s**t is a divisibility sequence for every choice of t; e.g. if s = (-1,-1,-1,...) and t = A279019, then s**t = A130031.
Table 2. s = A000027 = (0,1,2,3,4,5,...) = (n);
t = A000027, n                 s**t = A007778, n^(n+1)
t = A000290, n^2               s**t = A374881
t = A000040, P                 s**t = A374853
t = A000045, F                 s**t = A374857
t = A000032, L                 s**t = A374858
t = A000079, 2^n               s**t = A374859
t = A000201, LW                s**t = A374861
t = A005408, 2*n+1             s**t = A000407, (2*n+1)! / n!
t = A016777, 3*n+1             s**t = A113551
t = A016789, 3*n+2             s**t = A374866
t = A000142, n!                s**t = A374871
t = A032766, [3*n/2]           s**t = A374879
t = A000217, T                 s**t = A374892
t = A000984, C                 s**t = A374893
t = A038608, n*(-1)^n          s**t = A374894
Table 3. s = A000290 = (0,1,4,9,16,...) = (n^2);
t = A000290, n^2               s**t = A323540
t = A002522, n^2+1             s**t = A374884
t = A000217, T                 s**t = A374885
t = A000578, n^3               s**t = A374886
t = A000079, 2^n               s**t = A374887
t = A000225, 2^n-1             s**t = A374888
t = A005408, 2*n+1             s**t = A374889
t = A000045, F                 s**t = A374890
Table 4. s = t;
s = t = A000012, 1             s**s = A000079; 2^(n+1)
s = t = A000027, n             s**s = A007778, n^(n+1)
s = t = A000290, n^2           s**s = A323540
s = t = A000045, F             s**s = this sequence
s = t = A000032, L             s**s = A374850
s = t = A000079, 2^n           s**s = A369673
s = t = A000244, 3^n           s**s = A369674
s = t = A000040, P             s**s = A374851
s = t = A000201, LW            s**s = A374862
s = t = A005408, 2*n+1         s**s = A062971
s = t = A016777, 3*n+1         s**s = A374877
s = t = A016789, 3*n+2         s**s = A374878
s = t = A032766, [3*n/2]       s**s = A374880
s = t = A000217, T             s**s = A375050
s = t = A005563, n^2-1         s**s = A375051
s = t = A279019, n^2-n         s**s = A375056
s = t = A002398, n^2+n         s**s = A375058
s = t = A002061, n^2+n+1       s**s = A375059
If n = 2*k+1, then s**s(n) is a square; specifically,
s**s(n) = ((s(0)+s(n))*(s(1)+s(n-1))*...*(s(k)+s(k+1)))^2.
If n = 2*k, then s**s(n) has the form 2*s(k)*m^2, where m is an integer.
Table 5. Others
s = A000201, LW      t = A001950, UW         s**t = A374863
s = A000045, F       t = A000032, L          s**t = A374865
s = A005843, 2*n     t = A005408, 2*n+1      s**t = A085528, (2*n+1)^(n+1)
s = A016777, 3*n+1   t = A016789, 3*n+2      s**t = A091482
s = A005408, 2*n+1   t = A000045, F          s**t = A374867
s = A005408, 2*n+1   t = A000032, L          s**t = A374868
s = A005408, 2*n+1   t = A000079, 2^n        s**t = A374869
s = A000027, n       t = A000142, n!         s**t = A374871
s = A005408, 2*n+1   t = A000142, n!         s**t = A374872
s = A000079, 2^n     t = A000142, n!         s**t = A374874
s = A000142, n!      t = A000045, F          s**t = A374875
s = A000142, n!      t = A000032, L          s**t = A374876
s = A005408, 2*n+1   t = A016777, 3*n+1      s**t = A352601
s = A005408, 2*n+1   t = A016789, 3*n+2      s**t = A064352
Table 6. Arrays of coefficients of s(x)**t(x), where s(x) and t(x) are polynomials
s(x)              t(x)            s(x)**t(x)
 n                 x              A132393
n^2                x              A269944
x+1               x+1             A038220
x+2               x+2             A038244
 x                x+3             A038220
nx                x+1             A094638
 1              x^2+x+1           A336996
n^2 x             x+1             A375041
n^2 x            2x+1             A375042
n^2 x             x+2             A375043
2^n x             x+1             A375044
2^n              2x+1             A375045
2^n              x+2              A375046
x+1              F(n)             A375047
x+1             x+F(n)            A375048
x+F(n)          x+F(n)            A375049

			a(0) = 0 + 0 = 0
a(1) = (0+1) * (1+0) = 1
a(2) = (0+1) * (1+1) * (1+0) = 2
a(3) = (0+2) * (1+1) * (1+1) * (2+0) = 16
As noted above, a(2*k+1) is a square for k>=0. The first 5 squares are 1, 16, 3600, 12320100, 701841817600, with corresponding square roots 1, 4, 60, 3510, 837760.
If n = 2*k, then s**s(n) has the form 2*F(k)*m^2, where m is an integer and F(k) is the k-th Fibonacci number; e.g., a(6) = 2*F(3)*(192)^2.

Cf. A000045, A374850-to-A374881, A374884-to-A374894, A375041-to-A375049.

a:= n-> (F-> mul(F(n-j)+F(j), j=0..n))(combinat[fibonacci]):
seq(a(n), n=0..15);  # Alois P. Heinz, Aug 02 2024

s[n_] := Fibonacci[n]; t[n_] := Fibonacci[n];
u[n_] := Product[s[k] + t[n - k], {k, 0, n}];
Table[u[n], {n, 0, 20}]

a(n)=prod(k=0, n, fibonacci(k) + fibonacci(n-k)) \\ Andrew Howroyd, Jul 31 2024

a(n) ~ c * phi^(3*n^2/4 + n) / 5^((n+1)/2), where c = QPochhammer(-1, 1/phi^2)^2/2 if n is even and c = phi^(1/4) * QPochhammer(-phi, 1/phi^2)^2 / (phi + 1)^2 if n is odd, and phi = A001622 is the golden ratio. - Vaclav Kotesovec, Aug 01 2024

A028361 Number of totally isotropic spaces of index n in orthogonal geometry of dimension 2n.

Original entry on oeis.org

1, 2, 6, 30, 270, 4590, 151470, 9845550, 1270075950, 326409519150, 167448083323950, 171634285407048750, 351678650799042888750, 1440827432323678715208750, 11804699153027899713705288750, 193419995622362136809061156168750, 6338179836549184861096125026493768750

Offset: 0

N. J. A. Sloane

These numbers appear in first column of A155103. - Mats Granvik, Jan 20 2009
Equals row sums of unsigned triangle A158474. - Gary W. Adamson, Mar 20 2009
a(n) = (n+2) terms in the sequence (1, 1, 2, 4, 8, 16, ...) dot (n+2) terms in the sequence (1, 1, 2, 6, 30, 270, ...). Example: a(4) = 4590 = (1, 2, 4, 8, 16) dot (1, 1, 2, 6, 30, 270) = (1 + 1 + 4 + 24 + 240 + 4320). - Gary W. Adamson, Aug 02 2010
a(n) is the right-hand side of the mass formula used to classify Type II Self Dual Binary Linear Codes of length 2(n+1). a(n) is the number of distinct Type II Self Dual Binary Linear codes of length 2(n+1) when 2(n+1) = 0 MOD 8.  It is important to note that Type II codes are only possible when the length is a multiple of 8. In short, this sequence only applies to Type II codes when 2(n+1) = 0 MOD 8, else the right hand side of the mass formula is zero. - Nathan J. Russell, Mar 04 2016
This is almost certainly the sequence of number of Carlyle circles needed for the construction of regular polygons using straightedge and compass mentioned on page 107 of DeTemple (1991). - N. J. A. Sloane, Aug 05 2021
a(n) is also the number of Sp(oo, F2)-orbits of V^n, where V is the countable-dimensional symplectic vector space over the two-element field. - Jingjie Yang, Jul 30 2025

W. Cary Huffman and Vera Pless, Fundamentals of Error Correcting Codes, Cambridge University Press, 2003, Page 366. - Nathan J. Russell, Mar 04 2016

T. D. Noe, Table of n, a(n) for n=0..50
C. Bachoc and P. Gaborit, On extremal additive F_4 codes of length 10 to 18, J. Théorie Nombres Bordeaux, 12 (2000), 255-271.
Julien Clément and Antoine Genitrini, Binary Decision Diagrams: from Tree Compaction to Sampling, arXiv:1907.06743 [cs.DS], 2019.
John P. D'Angelo, Counting Tournament Brackets, J. Int. Seq., Vol. 25 (2022), Article 22.6.8.
Duane W. DeTemple, Carlyle circles and the Lemoine simplicity of polygon constructions, The American Mathematical Monthly 98.2 (1991): 97-108.
Davide Gaiotto and Justin Kulp, Orbifold groupoids, arXiv:2008.05960 [hep-th], 2020, page 42.

Cf. A006125, A028362, A155103, A158474, A323716 (product of 3^i+1).

[1] cat [ (&*[1+2^j: j in [0..n-1]]): n in [1..20]]; // G. C. Greubel, Jun 06 2020

seq( mul((1+2^j), j=0..n-1), n = 0..20); # G. C. Greubel, Jun 06 2020

Table[QPochhammer[-1, 2, n], {n, 0, 15}] (* Arkadiusz Wesolowski, Oct 29 2012 *)
Table[Product[2^i + 1, {i, 0, n/2 - 2}], {n, 2, 32, 2}] (* Nathan J. Russell, Mar 04 2016 *)
Table[Product[2^i + 1, {i, 0, n - 1}], {n, 0, 15}] (* Nathan J. Russell, Mar 04 2016 *)
FoldList[Times,1,2^Range[0,20]+1] (* Harvey P. Dale, Apr 11 2016 *)

PARI
```
{a(n) = prod(k=0, n-1, 2^k + 1)};
```

{a(n)=polcoeff(sum(m=0,n,2^(m*(m-1)/2)*x^m/prod(k=0,m,1-2^k*x+x*O(x^n))),n)} /* Paul D. Hanna, May 02 2012 */

for n in range(2,50,2):
  product = 1
  for i in range(0,n//2-2 + 1):
    product *= (2**i+1)
  print(product)
# Nathan J. Russell, Mar 01 2016

[product( 1+2^j for j in (0..n-1)) for n in (0..20)] # G. C. Greubel, Jun 06 2020

a(n) = Product_{i=0..n-1} ( 2^i + 1 ).
Asymptotic to C*2^(n*(n-1)/2) where C = A081845 = 4.76846205806274344829979857... = Product_{k>=0} (1 + 1/2^k). - Benoit Cloitre, Apr 09 2003
It appears that a(n) = 2^((1/2)*(n - 1)*n) * Product_{k>=0} (1 + 1/(2^k)) / Product_{k>=0} (1 + 1/(2^(n + k))). - Peter Moxey (pmoxey(AT)live.com), Mar 21 2010
G.f.: Sum_{n>=0} 2^(n*(n-1)/2) * x^n / Product_{k=0..n} (1 - 2^k*x). - Paul D. Hanna, May 02 2012
a(n) = (a(n-2)^3 + a(n-1) * a(n-3) * (a(n-1) - 2 * a(n-2))) * a(n-1) / (a(n-2)^2 * (a(n-2) - a(n-3))) if n>2. - Michael Somos, Aug 21 2012
0 = a(n)*(+a(n+1) + a(n+2)) + a(n+1)*(-2*a(n+1)) for all n>=0. - Michael Somos, Oct 10 2014
Sum_{k=0..n} 2^k/a(k) = 3-2/a(n) and Sum_{k=0..n} 4^k/a(k) = 9-(4*(1+2^n))/a(n) for n >= 0. - Werner Schulte, Dec 25 2016
G.f. A(x) satisfies: A(x) = (1 + x * A(2*x)) / (1 - x). - Ilya Gutkovskiy, Jun 06 2020
a(n) = Sum_{k=0..n} q_binomial(n, k, q=2) * 2^(k*(k-1)/2). - Jingjie Yang, Jul 30 2025

A003179 Number of self-dual binary codes of length 2n (up to column permutation equivalence).

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 3, 4, 7, 9, 16, 25, 55, 103, 261, 731, 3295, 24147, 519492

Offset: 0

N. J. A. Sloane

The length 36 binary self dual codes have been classified. - Nathan J. Russell, Feb 14 2016

This is number of binary self-dual codes of length 2n up to column permutation equivalence. Sequence A028362 gives an actual count of all possible binary self-dual codes of length 2n. - Nathan J. Russell, Nov 25 2018

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

R. T. Bilous and G. H. J. van Rees, An enumeration of binary self-dual codes of length 32, preprint.
R. T. Bilous and G. H. J. van Rees, An enumeration of binary self-dual codes of length 32, Designs, Codes Crypt., 26 (2002), 61-86.
J. H. Conway and V. S. Pless, On the enumeration of self-dual codes, J. Comb. Theory, A28 (1980), 26-53. MR0558873
J. H. Conway, V. Pless and N. J. A. Sloane, The Binary Self-Dual Codes of Length Up to 32: A Revised Enumeration, J. Comb. Theory, A28 (1980), 26-53 (Abstract, pdf, ps, Table A, Table D).
Masaaki Harada and Akihiro Munemasa, Classification of Self-Dual Codes of Length 36, arXiv:1012.5464 [math.CO], 2010-2012.
W. C. Huffman, On the classification and enumeration of self-dual codes, Finite Fields Applic. 11 (2005), 451-490.
W. Cary Huffman and Vera Pless, Fundamentals of Error Correcting Codes, Cambridge University Press, 2003, Pages 7,252-282,338-393.
G. Nebe, E. M. Rains and N. J. A. Sloane, Self-Dual Codes and Invariant Theory, Springer, Berlin, 2006.
E. M. Rains and N. J. A. Sloane, Self-dual codes, pp. 177-294 of Handbook of Coding Theory, Elsevier, 1998 (Abstract, pdf, ps).

Cf. A003178, A028362, A028363, A105685. Equals A106163 + A106165.

a(18) from Nathan J. Russell, Feb 14 2016

Name clarified by Nathan J. Russell, Nov 26 2018

A079555 Decimal expansion of Product_{k>=1} (1 + 1/2^k) = 2.384231029031371...

Original entry on oeis.org

2, 3, 8, 4, 2, 3, 1, 0, 2, 9, 0, 3, 1, 3, 7, 1, 7, 2, 4, 1, 4, 9, 8, 9, 9, 2, 8, 8, 6, 7, 8, 3, 9, 7, 2, 3, 8, 7, 7, 1, 6, 1, 9, 5, 1, 6, 5, 0, 8, 4, 3, 3, 4, 5, 7, 6, 9, 2, 1, 0, 1, 5, 0, 7, 9, 8, 9, 1, 8, 1, 2, 9, 3, 0, 3, 6, 0, 3, 7, 2, 5, 5, 1, 8, 6, 5, 3, 5, 2, 1, 0, 3, 6, 5, 6, 8, 0, 5, 2, 0, 0, 0, 2, 6, 8

Offset: 1

Benoit Cloitre, Jan 25 2003

			2.38423102903137172414989928867839723877161951650843345769...

G. C. Greubel, Table of n, a(n) for n = 1..1200
Richard J. McIntosh, Some Asymptotic Formulae for q-Hypergeometric Series, Journal of the London Mathematical Society, Vol. 51, No. 1 (1995), pp. 120-136; alternative link.

Cf. A028362, A048651, A081845, A098844, A100220, A132019-A132026, A132034-A132038, A132265-A132268, A132323-A132326, A132269, A132270, A000593.

Cf. A141848. - Gleb Koloskov, Apr 04 2021

digits = 105; NProduct[(1 + 1/2^k), {k, 1, Infinity}, WorkingPrecision -> digits+10, NProductFactors -> 200] // RealDigits[#, 10, digits]& // First (* Jean-François Alcover, Feb 14 2013 *)
N[QPochhammer[-1/2,1/2]] (* G. C. Greubel, Dec 05 2015 *)
1/N[QPochhammer[1/2, 1/4]] (* Gleb Koloskov, Apr 04 2021 *)

prodinf(n=1,1+2.^-n) \\ Charles R Greathouse IV, May 27 2015

1/prodinf(n=0, 1-2^(-2*n-1)) \\ Gleb Koloskov, Apr 04 2021

(1/2)*lim sup Product_{k=0..floor(log_2(n)), (1 + 1/floor(n/2^k))} for n-->oo. - Hieronymus Fischer, Aug 20 2007
(1/2)*lim sup A132369(n)/A098844(n) for n-->oo. - Hieronymus Fischer, Aug 20 2007
(1/2)*lim sup A132269(n)/n^((1+log_2(n))/2) for n-->oo. - Hieronymus Fischer, Aug 20 2007
(1/2)*lim sup A132270(n)/n^((log_2(n)-1)/2) for n-->oo. - Hieronymus Fischer, Aug 20 2007
exp(sum{n>0, 2^(-n)*sum{k|n, -(-1)^k/k}})=exp(sum{n>0, A000593(n)/(n*2^n)}). - Hieronymus Fischer, Aug 20 2007
(1/2)*lim sup A132269(n+1)/A132269(n)=2.3842310290313717241498992886... for n-->oo. - Hieronymus Fischer, Aug 20 2007
Equals (-1/2; 1/2){infinity}, where (a;q){infinity} is the q-Pochhammer symbol. - G. C. Greubel, Dec 05 2015
2 + Sum_{k>1} 1/(Product_{i=2..k} (2^i-1)) = 2 + 1/3 + 1/(3*7) + 1/(3*7*15) + 1/(3*7*15*31) + 1/(3*7*15*31*63) + ... (conjecture). - Werner Schulte, Dec 22 2016
From Peter Bala, Dec 15 2020: (Start)
The above conjecture of Schulte follows by setting x = 1/2 and t = -1 in the identity Product_{k >= 1} (1 - t*x^k) = Sum_{n >= 0} (-1)^n*x^(n*(n+1)/2)*t^n/( Product_{k = 1..n} 1 - x^k ), due to Euler.
Constant C = 1 + Sum_{n >= 0} (1/2)^(n+1)*Product_{k = 1..n} (1 + 1/2^k).
C = 2 + Sum_{n >= 0} (1/4)^(n+1)*Product_{k = 1..n} (1 + 1/2^k).
3*C = 7 + Sum_{n >= 0} (1/8)^(n+1)*Product_{k = 1..n} (1 + 1/2^k).
3*7*C = 50 + Sum_{n >= 0} (1/16)^(n+1)*Product_{k = 1..n} (1 + 1/2^k).
3*7*15*C = 751 + Sum_{n >= 0} (1/32)^(n+1)*Product_{k = 1..n} (1 + 1/2^k).
(End)
Equals 1/(1-P), where P is the Pell constant from A141848. - Gleb Koloskov, Apr 04 2021
Equals Sum_{k>=0} A000009(k)/2^k. - Vaclav Kotesovec, Sep 15 2021
From Amiram Eldar, Feb 19 2022: (Start)
Equals (sqrt(2)/2) * exp(log(2)/24 + Pi^2/(12*log(2))) * Product_{k>=1} (1 - exp(-2*(2*k-1)*Pi^2/log(2))) (McIntosh, 1995).
Equals (1/2) * A081845.
Equals Sum_{n>=0} 1/A005329(n). (End)

A003178 Number of indecomposable self-dual binary codes of length 2n.

Original entry on oeis.org

1, 1, 0, 0, 1, 0, 1, 1, 2, 2, 6, 8, 26, 45, 148, 457, 2523, 20786

Offset: 0

N. J. A. Sloane

R. T. Bilous, Enumeration of binary self-dual codes of length 34, Preprint, 2005.
R. T. Bilous and G. H. J. van Rees, An enumeration of binary self-dual codes of length 32, Designs, Codes Crypt., 26 (2002), 61-86.
J. H. Conway and V. S. Pless, On the enumeration of self-dual codes, J. Comb. Theory, A28 (1980), 26-53.
V. S. Pless, The children of the (32,16) doubly even codes, IEEE Trans. Inform. Theory, 24 (1978), 738-746.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

G. Nebe, E. M. Rains and N. J. A. Sloane, Self-Dual Codes and Invariant Theory, Springer, Berlin, 2006.
J. H. Conway, V. Pless and N. J. A. Sloane, The Binary Self-Dual Codes of Length Up to 32: A Revised Enumeration, J. Comb. Theory, A28 (1980), 26-53 (Abstract, pdf, ps, Table A, Table D).
E. M. Rains and N. J. A. Sloane, Self-dual codes, pp. 177-294 of Handbook of Coding Theory, Elsevier, 1998 (Abstract, pdf, ps).

Cf. A003179, A028362, A028363. Equals A106162 + A106164.

a(16) corrected and a(17) added by N. J. A. Sloane, based on data in Bilous's paper, Sep 06 2005

A028363 Total number of doubly-even self-dual binary codes of length 8n.

Original entry on oeis.org

1, 30, 9845550, 171634285407048750, 193419995622362136809061156168750, 14272693289804307141953423466197932293533748208968750

Offset: 0

N. J. A. Sloane

F. J. MacWilliams and N. J. A. Sloane, The Theory of Error-Correcting Codes, Elsevier-North Holland, 1978, p. 631.

Cf. A003178, A003179, A028362.

Join[{1},Table[2*Product[2^i+1,{i,4n-2}],{n,6}]] (* Harvey P. Dale, May 08 2013 *)
Table[Product[2^i + 1, {i, 0, n/2 - 2}], {n, 8, 40, 8}] (* Nathan J. Russell, Mar 04 2016 *)

for n in range(8, 50, 8):
    product = 1
    for i in range(n//2 - 1):
        product *= 2**i + 1
    print(product, end=", ")
# Nathan J. Russell, Mar 01 2016

a(n) = 2*Product_{i=1..4n-2} (2^i + 1).

There is an error in Eq. (75) of F. J. MacWilliams and N. J. A. Sloane, the lower subscript should be 1 not 0.

Formula corrected by N. J. A. Sloane, May 07 2013 following a suggestion from Harvey P. Dale

A068052 Start from 1, shift one left and sum mod 2 (bitwise-XOR) to get 3 (11 in binary), then shift two steps left and XOR to get 15 (1111 in binary), then three steps and XOR to get 119 (1110111 in binary), then four steps and so on.

Original entry on oeis.org

1, 3, 15, 119, 1799, 59367, 3743271, 481693095, 123123509927, 62989418816679, 64491023022979239, 132015402419352060071, 540829047855347718631591, 4430403202865824763042320551, 72583450474242118015031400337575, 2378466805556971511916001231449723047

Offset: 0

Antti Karttunen, Feb 13 2002

a(n) = each row of A053632 reduced mod 2 and interpreted as a binary number.

Antti Karttunen, Table of n, a(n) for n = 0..64

Same sequence shown in binary: A068053.

Cf. A000120, A003987, A028362 (using + instead of XOR), A048675, A053632, A080791, A248663, A285101, A285102, A285103, A285105.

with(gfun,seriestolist); [seq(foo(map(`mod`,seriestolist(series(mul(1+(z^i),i=1..n),z,binomial(n+1,2)+1)),2)), n=0..20)];
foo := proc(a) local i; add(a[i]*2^(i-1),i=1..nops(a)); end;
# second Maple program:
a:= proc(n) option remember; `if`(n=0, 1,
      (t-> Bits[Xor](2^n*t, t))(a(n-1)))
    end:
seq(a(n), n=0..16);  # Alois P. Heinz, Mar 07 2024

FoldList[BitXor[#, #*#2]&, 1, 2^Range[20]] (* Paolo Xausa, Mar 07 2024 *)

a(n) = if(n<1, 1, bitxor(a(n - 1), 2^n*a(n - 1))); \\ Indranil Ghosh, Apr 15 2017, after formula by Antti Karttunen

a(0) = 1; for n > 0, a(n) = a(n-1) XOR (2^n)*a(n-1), where XOR is bitwise-XOR (A003987).
a(n) = A248663(A285101(n)) = A048675(A285102(n)).
A000120(a(n)) = A285103(n). [Number of ones in binary representation.]
A080791(a(n)) = A285105(n). [Number of nonleading zeros.]

Formulas added by Antti Karttunen, Apr 15 2017

cp's OEIS Frontend

A048675 If n = p_i^e_i * ... * p_k^e_k, p_i < ... < p_k primes (with p_i = prime(i)), then a(n) = (1/2) * (e_i * 2^i + ... + e_k * 2^k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A005329 a(n) = Product_{i=1..n} (2^i - 1). Also called 2-factorial numbers.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

Formula

Extensions

A053632 Irregular triangle read by rows giving coefficients in expansion of Product_{k=1..n} (1 + x^k).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A374848 Obverse convolution A000045**A000045; see Comments.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A028361 Number of totally isotropic spaces of index n in orthogonal geometry of dimension 2n.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

Python

Sage

Formula

A003179 Number of self-dual binary codes of length 2n (up to column permutation equivalence).