A024023 -id:A024023 - cp's OEIS frontend

A062318 Numbers of the form 3^m - 1 or 2*3^m - 1; i.e., the union of sequences A048473 and A024023.

0, 1, 2, 5, 8, 17, 26, 53, 80, 161, 242, 485, 728, 1457, 2186, 4373, 6560, 13121, 19682, 39365, 59048, 118097, 177146, 354293, 531440, 1062881, 1594322, 3188645, 4782968, 9565937, 14348906, 28697813, 43046720, 86093441, 129140162

Offset: 1

Ahmed Fares (ahmedfares(AT)my-deja.com), Jul 05 2001

WARNING: The offset of this sequence has been changed from 0 to 1 without correcting the formulas and programs, many of them correspond to the original indexing a(0)=0, a(1)=1, ... - M. F. Hasler, Oct 06 2014
Numbers n such that no entry in n-th row of Pascal's triangle is divisible by 3, i.e., such that A062296(n) = 0.
The base 3 representation of these numbers is 222...222 or 122...222.
a(n+1) is the smallest number with ternary digit sum = n: A053735(a(n+1)) = n and A053735(m) <> n for m < a(n+1). - Reinhard Zumkeller, Sep 15 2006
A138002(a(n)) = 0. - Reinhard Zumkeller, Feb 26 2008
Also, number of terms in S(n), where S(n) is defined in A114482. - N. J. A. Sloane, Nov 13 2014
a(n+1) is also the Moore lower bound on the order of a (4,g)-cage. - Jason Kimberley, Oct 30 2011

			The first rows in Pascal's triangle with no multiples of 3 are:
row 0: 1;
row 1: 1, 1;
row 2: 1, 2,  1;
row 5: 1, 5, 10, 10,  5,  1;
row 8: 1, 8, 28, 56, 70, 56, 28, 8, 1;

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Daniel Birmajer, Juan B. Gil, Jordan O. Tirrell, and Michael D. Weiner, Pattern-avoiding stabilized-interval-free permutations, arXiv:2306.03155 [math.CO], 2023.
Sayan Dutta, Lorenz Halbeisen, and Norbert Hungerbühler, Properties of Hesse derivatives of cubic curves, arXiv:2309.05048 [math.AG], 2023. See p. 9.
Gyula Tasi and Fujio Mizukami, Quantum algebraic-combinatoric study of the conformational properties of n-alkanes, J. Math. Chemistry, 25, 1999, 55-64 (see p. 60).
Index entries for linear recurrences with constant coefficients, signature (1,3,-3).

Cf. A062296, A024023, A048473, A114482. Pairwise sums of A052993.
Moore lower bound on the order of a (k,g) cage: A198300 (square); rows: A000027 (k=2), A027383 (k=3), this sequence (k=4), A061547 (k=5), A198306 (k=6), A198307 (k=7), A198308 (k=8), A198309 (k=9), A198310 (k=10), A094626 (k=11); columns: A020725 (g=3), A005843 (g=4), A002522 (g=5), A051890 (g=6), A188377 (g=7). - Jason Kimberley, Oct 30 2011
Cf. A037233 (actual order of a (4,g)-cage).
Smallest number whose base b sum of digits is n: A000225 (b=2), this sequence (b=3), A180516 (b=4), A181287 (b=5), A181288 (b=6), A181303 (b=7), A165804 (b=8), A140576 (b=9), A051885 (b=10).

I:=[0,1,2]; [n le 3 select I[n] else Self(n-1)+3*Self(n-2) -3*Self(n-3): n in [1..40]]; // Vincenzo Librandi, Apr 20 2012

A062318 :=proc(n)
    if n mod 2 = 1 then
        3^((n-1)/2)-1
    else
        2*3^(n/2-1)-1
    fi
end proc:
seq(A062318(n), n=1..37); # Emeric Deutsch, Feb 03 2005, offset updated

CoefficientList[Series[x^2*(1+x)/((1-x)*(1-3*x^2)),{x,0,40}],x] (* Vincenzo Librandi, Apr 20 2012 *)
A062318[n_]:= (1/3)*(Boole[n==0] -3 +3^(n/2)*(2*Mod[n+1,2] +Sqrt[3] *Mod[n, 2]));
Table[A062318[n], {n, 50}] (* G. C. Greubel, Apr 17 2023 *)

PARI
```
a(n)=3^(n\2)<M. F. Hasler, Oct 06 2014
```

def A062318(n): return (1/3)*(int(n==0) - 3 + 2*((n+1)%2)*3^(n/2) + (n%2)*3^((n+1)/2))
[A062318(n) for n in range(1,41)] # G. C. Greubel, Apr 17 2023

a(n) = 2*3^(n/2-1)-1 if n is even; a(n) = 3^(n/2-1/2)-1 if n is odd. - Emeric Deutsch, Feb 03 2005, offset updated.
From Paul Curtz, Feb 21 2008: (Start)
a(n) = a(n-1) + 3*a(n-2) - 3*a(n-3).
Partial sums of A108411. (End)
G.f.: x^2*(1+x)/((1-x)*(1-3*x^2)). - Colin Barker, Apr 02 2012
a(2n+1) = 3*a(2n-1) + 2;  a(2n) = ( a(2n-1) + a(2n+1) )/2. See A060647 for case where a(1)= 1. - Richard R. Forberg, Nov 30 2013
a(n) = 2^((1+(-1)^n)/2) * 3^((2*n-3-(-1)^n)/4) - 1. - Luce ETIENNE, Aug 29 2014
a(n) = A052993(n-1) + A052993(n-2). - R. J. Mathar, Sep 10 2021
E.g.f.: (1 - 3*cosh(x) + 2*cosh(sqrt(3)*x) - 3*sinh(x) + sqrt(3)*sinh(sqrt(3)*x))/3. - Stefano Spezia, Apr 06 2022
a(n) = (1/3)*([n=0] - 3 + (1+(-1)^n)*3^(n/2) + ((1-(-1)^n)/2)*3^((n+1)/2)). - G. C. Greubel, Apr 17 2023

More terms from Emeric Deutsch, Feb 03 2005

Entry revised by N. J. A. Sloane, Jul 29 2011

A064079 Zsigmondy numbers for a = 3, b = 1: Zs(n, 3, 1) is the greatest divisor of 3^n - 1^n (A024023) that is relatively prime to 3^m - 1^m for all positive integers m < n.

Original entry on oeis.org

2, 1, 13, 5, 121, 7, 1093, 41, 757, 61, 88573, 73, 797161, 547, 4561, 3281, 64570081, 703, 581130733, 1181, 368089, 44287, 47071589413, 6481, 3501192601, 398581, 387440173, 478297, 34315188682441, 8401, 308836698141973, 21523361

Offset: 1

Jens Voß, Sep 04 2001

nonn

By Zsigmondy's theorem, the n-th Zsigmondy number for bases a and b is not 1 except in the three cases (1) a = 2, b = 1, n = 1, (2) a = 2, b = 1, n = 6, (3) n = 2 and a+b is a power of 2.

K. Zsigmondy, Zur Theorie der Potenzreste, Monatsh. f. Math. 3 (1892) 265-284.

Cf. A024023, A064078, A064080, A064081, A064082, A064083.

More terms from Vladeta Jovovic, Sep 06 2001

Definition corrected by Jerry Metzger, Nov 04 2009

A034472 a(n) = 3^n + 1.

Original entry on oeis.org

2, 4, 10, 28, 82, 244, 730, 2188, 6562, 19684, 59050, 177148, 531442, 1594324, 4782970, 14348908, 43046722, 129140164, 387420490, 1162261468, 3486784402, 10460353204, 31381059610, 94143178828, 282429536482, 847288609444, 2541865828330, 7625597484988

Offset: 0

N. J. A. Sloane

Companion numbers to A003462.
a(n) = A024101(n)/A024023(n). - Reinhard Zumkeller, Feb 14 2009
Mahler exhibits this sequence with n>=2 as a proof that there exists an infinite number of x coprime to 3, such that x belongs to A005836 and x^2 belong to A125293. - Michel Marcus, Nov 12 2012
a(n-1) is the number of n-digit base 3 numbers that have an even number of digits 0. - Yifan Xie, Jul 13 2024

			a(3)=28 because 4*a(2)-3*a(1)=4*10-3*4=28 (28 is also 3^3 + 1).
G.f. = 2 + 4*x + 10*x^2 + 28*x^3 + 82*x^4 + 244*x^5 + 730*x^5 + ...

Knuth, Donald E., Satisfiability, Fascicle 6, volume 4 of The Art of Computer Programming. Addison-Wesley, 2015, pages 148 and 220, Problem 191.
P. Ribenboim, The Little Book of Big Primes, Springer-Verlag, NY, 1991, pp. 35-36, 53.

T. D. Noe, Table of n, a(n) for n=0..200
T. A. Gulliver, Divisibility of sums of powers of odd integers, Int. Math. For. 5 (2010) 3059-3066, eq 5.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 454
Kurt Mahler, The representation of squares to the base 3, Acta Arith. Vol. 53, Issue 1 (1989), p. 99-106.
Burkard Polster, Special numbers in 3-coloring of Pascal's triangle, Mathologer video (2019).
Amelia Carolina Sparavigna, Some Groupoids and their Representations by Means of Integer Sequences, International Journal of Sciences (2019) Vol. 8, No. 10.
D. Suprijanto and I. W. Suwarno, Observation on Sums of Powers of Integers Divisible by 3k-1, Applied Mathematical Sciences, Vol. 8, 2014, no. 45, 2211 - 2217.
Eric Weisstein's World of Mathematics, Lucas Sequence
Index entries for linear recurrences with constant coefficients, signature (4,-3).

Cf. A003462, A000204, A007051, A000051, A052539, A034474, A062394, A034491, A062395, A062396, A062397, A007689, A063376, A063481, A074600 - A074624, A034524, A178248, A228081, A279396.

[3^n+1: n in [0..30]]; // Vincenzo Librandi, Jan 11 2017

ZL:= [S, {S=Union(Sequence(Z), Sequence(Union(Z, Z, Z)))}, unlabeled]: seq(combstruct[count](ZL, size=n), n=0..25); # Zerinvary Lajos, Jun 19 2008
g:=1/(1-3*z): gser:=series(g, z=0, 43): seq(coeff(gser, z, n)+1, n=0..31); # Zerinvary Lajos, Jan 09 2009

Mathematica
```
Table[3^n + 1, {n, 0, 24}]
```
PARI
```
a(n) = 3^n + 1
```

Vec(2*(1-2*x)/((1-x)*(1-3*x)) + O(x^50)) \\ Altug Alkan, Nov 15 2015

[lucas_number2(n,4,3) for n in range(27)] # Zerinvary Lajos, Jul 08 2008

[sigma(3,n) for n in range(27)] # Zerinvary Lajos, Jun 04 2009

[3^n+1 for n in range(30)] # Bruno Berselli, Jan 11 2017

a(n) = 3*a(n-1) - 2 = 4*a(n-1) - 3*a(n-2). (Lucas sequence, with A003462, associated to the pair (4, 3).)
G.f.: 2*(1-2*x)/((1-x)*(1-3*x)). Inverse binomial transforms yields 2,2,4,8,16,... i.e., A000079 with the first entry changed to 2. Binomial transform yields A063376 without A063376(-1). - R. J. Mathar, Sep 05 2008
E.g.f.: exp(x) + exp(3*x). - Mohammad K. Azarian, Jan 02 2009
a(n) = A279396(n+3,3). - Wolfdieter Lang, Jan 10 2017
a(n) = 2*A007051(n). - R. J. Mathar, Apr 07 2022

Additional comments from Rick L. Shepherd, Feb 13 2002

A052548 a(n) = 2^n + 2.

Original entry on oeis.org

3, 4, 6, 10, 18, 34, 66, 130, 258, 514, 1026, 2050, 4098, 8194, 16386, 32770, 65538, 131074, 262146, 524290, 1048578, 2097154, 4194306, 8388610, 16777218, 33554434, 67108866, 134217730, 268435458, 536870914, 1073741826, 2147483650

Offset: 0

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

The most "compact" sequence that satisfies Bertrand's Postulate. Begin with a(1) = 3 = n, then 2n - 2 = 4 = n_1, 2n_1 - 2 = 6 = n_2, 2n_2 - 2 = 10, etc. = a(n), hence there is guaranteed to be at least one prime between successive members of the sequence. - Andrew S. Plewe, Dec 11 2007

Number of 2-sided prudent polygons of area n, for n>0, see Beaton, p. 5. - Jonathan Vos Post, Nov 30 2010

Vincenzo Librandi, Table of n, a(n) for n = 0..240
Nicholas R. Beaton, Philippe Flajolet, and Anthony J. Guttmann, The Enumeration of Prudent Polygons by Area and its Unusual Asymptotics, arXiv:1011.6195 [math.CO], Nov 29, 2010.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 485
Popular Computing (Calabasas, CA), Sieves: Problem 43, Vol. 2 (No. 13, Apr 1974), pp. 6-7. This is Sieve #6 with K=2. [Annotated and scanned copy]
Eric Weisstein's World of Mathematics, Bertrand's Postulate
Index entries for sequences generated by sieves
Index entries for linear recurrences with constant coefficients, signature (3,-2).

Apart from initial term, same as A056469.
Cf. A003462, A007051, A034472, A024023, A067771, A029858, A134931, A115099, A100774, A079004, A058481, A100585, A100586, A058896, A000918, A173786.
Cf. also A000079, A000051, A100314.

a052548 = (+ 2) . a000079
a052548_list = iterate ((subtract 2) . (* 2)) 3
-- Reinhard Zumkeller, Sep 05 2015

[2^n + 2: n in [0..35]]; // Vincenzo Librandi, Apr 29 2011

spec := [S,{S=Union(Sequence(Union(Z,Z)),Sequence(Z),Sequence(Z))},unlabeled]: seq(combstruct[count](spec,size=n), n=0..20);

2^Range[0,40]+2 (* Harvey P. Dale, Jun 26 2012 *)

a(n)=1<Charles R Greathouse IV, Nov 20 2011

G.f.: (3-5*x)/((1-2*x)*(1-x)) = (3-5*x)/(1 - 3*x + 2*x^2) = 2/(1-x) + 1/(1-2*x).
a(0)=3, a(1)=4, a(n) = 3*a(n-1) - 2*a(n-2).
a(n) = A058896(n)/A000918(n), for n>0. - Reinhard Zumkeller, Feb 14 2009
a(n) = A173786(n,1), for n>0. - Reinhard Zumkeller, Feb 28 2010
a(n)*A000918(n) = A028399(2*n), for n>0. - Reinhard Zumkeller, Feb 28 2010
a(0)=3, a(n) = 2*a(n-1) - 2. - Vincenzo Librandi, Aug 06 2010
E.g.f.: (2 + exp(x))*exp(x). - Ilya Gutkovskiy, Aug 16 2016

More terms from James Sellers, Jun 06 2000

A013609 Triangle of coefficients in expansion of (1+2*x)^n.

Original entry on oeis.org

1, 1, 2, 1, 4, 4, 1, 6, 12, 8, 1, 8, 24, 32, 16, 1, 10, 40, 80, 80, 32, 1, 12, 60, 160, 240, 192, 64, 1, 14, 84, 280, 560, 672, 448, 128, 1, 16, 112, 448, 1120, 1792, 1792, 1024, 256, 1, 18, 144, 672, 2016, 4032, 5376, 4608, 2304, 512, 1, 20, 180, 960, 3360, 8064, 13440, 15360, 11520, 5120, 1024

Offset: 0

N. J. A. Sloane

T(n,k) is the number of lattice paths from (0,0) to (n,k) with steps (1,0) and two kinds of steps (1,1). The number of paths with steps (1,0) and s kinds of steps (1,1) corresponds to the expansion of (1+s*x)^n. - Joerg Arndt, Jul 01 2011
Also sum of rows in A046816. - Lior Manor, Apr 24 2004
Also square array of unsigned coefficients of Chebyshev polynomials of second kind. - Philippe Deléham, Aug 12 2005
The rows give the number of k-simplices in the n-cube. For example, 1, 6, 12, 8 shows that the 3-cube has 1 volume, 6 faces, 12 edges and 8 vertices. - Joshua Zucker, Jun 05 2006
Triangle whose (i, j)-th entry is binomial(i, j)*2^j.
With offset [1,1] the triangle with doubled numbers, 2*a(n,m), enumerates sequences of length m with nonzero integer entries n_i satisfying sum(|n_i|) <= n. Example n=4, m=2: [1,3], [3,1], [2,2] each in 2^2=4 signed versions: 2*a(4,2) = 2*6 = 12. The Sum over m (row sums of 2*a(n,m)) gives 2*3^(n-1), n >= 1. See the W. Lang comment and a K. A. Meissner reference under A024023. - Wolfdieter Lang, Jan 21 2008
n-th row of the triangle = leftmost column of nonzero terms of X^n, where X = an infinite bidiagonal matrix with (1,1,1,...) in the main diagonal and (2,2,2,...) in the subdiagonal. - Gary W. Adamson, Jul 19 2008
Numerators of a matrix square-root of Pascal's triangle A007318, where the denominators for the n-th row are set to 2^n. - Gerald McGarvey, Aug 20 2009
From Johannes W. Meijer, Sep 22 2010: (Start)
The triangle sums (see A180662 for their definitions) link the Pell-Jacobsthal triangle, whose mirror image is A038207, with twenty-four different sequences; see the crossrefs.
This triangle may very well be called the Pell-Jacobsthal triangle in view of the fact that A000129 (Kn21) are the Pell numbers and A001045 (Kn11) the Jacobsthal numbers.
(End)
T(n,k) equals the number of n-length words on {0,1,2} having n-k zeros. - Milan Janjic, Jul 24 2015
T(n-1,k-1) is the number of 2-compositions of n with zeros having k positive parts; see Hopkins & Ouvry reference. - Brian Hopkins, Aug 16 2020
T(n,k) is the number of chains 0=x_0Geoffrey Critzer, Oct 01 2022
Excluding the initial 1, T(n,k) is the number of k-faces of a regular n-cross polytope. See A038207 for n-cube and A135278 for n-simplex. - Mohammed Yaseen, Jan 14 2023

			Triangle begins:
  1;
  1,  2;
  1,  4,   4;
  1,  6,  12,    8;
  1,  8,  24,   32,   16;
  1, 10,  40,   80,   80,    32;
  1, 12,  60,  160,  240,   192,    64;
  1, 14,  84,  280,  560,   672,   448,    128;
  1, 16, 112,  448, 1120,  1792,  1792,   1024,    256;
  1, 18, 144,  672, 2016,  4032,  5376,   4608,   2304,    512;
  1, 20, 180,  960, 3360,  8064, 13440,  15360,  11520,   5120,  1024;
  1, 22, 220, 1320, 5280, 14784, 29568,  42240,  42240,  28160, 11264,  2048;
  1, 24, 264, 1760, 7920, 25344, 59136, 101376, 126720, 112640, 67584, 24576, 4096;
From _Peter Bala_, Apr 20 2012: (Start)
The triangle can be written as the matrix product A038207*(signed version of A013609).
  |.1................||.1..................|
  |.2...1............||-1...2..............|
  |.4...4...1........||.1..-4...4..........|
  |.8..12...6...1....||-1...6...-12...8....|
  |16..32..24...8...1||.1..-8....24.-32..16|
  |..................||....................|
(End)

B. N. Cyvin et al., Isomer enumeration of unbranched catacondensed polygonal systems with pentagons and heptagons, Match, No. 34 (Oct 1996), pp. 109-121.
G. Hotz, Zur Reduktion von Schaltkreispolynomen im Hinblick auf eine Verwendung in Rechenautomaten, El. Datenverarbeitung, Folge 5 (1960), pp. 21-27.

T. D. Noe, Rows n=0..50 of triangle, flattened
Feryal Alayont and Evan Henning, Edge Covers of Caterpillars, Cycles with Pendants, and Spider Graphs, J. Int. Seq. (2023) Vol. 26, Art. 23.9.4.
H. J. Brothers, Pascal's Prism: Supplementary Material, PDF version.
John Cartan, Cartan's triangle shows the relationship to the n-cube.
R. Ehrenborg and M. Readdy, Characterization of the factorial functions of Eulerian binomial and Sheffer posets
J. Goldman and J. Haglund, Generalized rook polynomials, J. Combin. Theory A91 (2000), 509-530, 1-rook coefficients of k rooks on the 2xn board, all heights 2.
W. G. Harter, Representations of multidimensional symmetries in networks, J. Math. Phys., 15 (1974), 2016-2021.
Brian Hopkins and Stéphane Ouvry, Combinatorics of Multicompositions, arXiv:2008.04937 [math.CO], 2020.
D. A. Zaitsev, A generalized neighborhood for cellular automata, Theoretical Computer Science, 666 (2017), 21-35.

Cf. A007318, A013610, etc.
Appears in A167580 and A167591. - Johannes W. Meijer, Nov 23 2009
From Johannes W. Meijer, Sep 22 2010: (Start)
Triangle sums (see the comments): A000244 (Row1); A000012 (Row2); A001045 (Kn11); A026644 (Kn12); 4*A011377 (Kn13); A000129 (Kn21); A094706 (Kn22); A099625 (Kn23); A001653 (Kn3); A007583 (Kn4); A046717 (Fi1); A007051 (Fi2); A077949 (Ca1); A008998 (Ca2); A180675 (Ca3); A092467 (Ca4); A052942 (Gi1); A008999 (Gi2); A180676 (Gi3); A180677 (Gi4); A140413 (Ze1); A180678 (Ze2); A097117 (Ze3); A055588 (Ze4).
(End)
Cf. A024023, A038207, A046816, A105728, A115068, A135278, A171977, A180662.
T(2n,n) gives A059304.

a013609 n = a013609_list !! n
a013609_list = concat $ iterate ([1,2] *) [1]
instance Num a => Num [a] where
   fromInteger k = [fromInteger k]
   (p:ps) + (q:qs) = p + q : ps + qs
   ps + qs         = ps ++ qs
   (p:ps) * qs'@(q:qs) = p * q : ps * qs' + [p] * qs
    *                = []
-- Reinhard Zumkeller, Apr 02 2011

a013609 n k = a013609_tabl !! n !! k
a013609_row n = a013609_tabl !! n
a013609_tabl = iterate (\row -> zipWith (+) ([0] ++ row) $
                                zipWith (+) ([0] ++ row) (row ++ [0])) [1]
-- Reinhard Zumkeller, Jul 22 2013, Feb 27 2013

[2^k*Binomial(n,k): k in [0..n], n in [0..15]]; // G. C. Greubel, Sep 17 2021

bin2:=proc(n,k) option remember; if k<0 or k>n then 0 elif k=0 then 1 else 2*bin2(n-1,k-1)+bin2(n-1,k); fi; end; # N. J. A. Sloane, Jun 01 2009

Flatten[Table[CoefficientList[(1 + 2*x)^n, x], {n, 0, 10}]][[1 ;; 59]] (* Jean-François Alcover, May 17 2011 *)
BinomialROW[n_, k_, t_] := Sum[Binomial[n, k]*Binomial[k, j]*(-1)^(k - j)*t^j, {j, 0, k}]; Column[Table[BinomialROW[n, k, 3], {n, 0, 10}, {k, 0, n}], Center] (* Kolosov Petro, Jan 28 2019 *)

a(n,k):=coeff(expand((1+2*x)^n),x^k);
create_list(a(n,k),n,0,6,k,0,n); /* Emanuele Munarini, Nov 21 2012 */

/* same as in A092566 but use */
steps=[[1,0], [1,1], [1,1]]; /* note double [1,1] */
/* Joerg Arndt, Jul 01 2011 */

flatten([[2^k*binomial(n,k) for k in (0..n)] for n in (0..15)]) # G. C. Greubel, Sep 17 2021

G.f.: 1 / (1 - x*(1+2*y)).
T(n,k) = 2^k*binomial(n,k).
T(n,k) = 2*T(n-1,k-1) + T(n-1,k). - Jon Perry, Nov 22 2005
Row sums are 3^n = A000244(n). - Joerg Arndt, Jul 01 2011
T(n,k) = Sum_{i=n-k..n} C(i,n-k)*C(n,i). - Mircea Merca, Apr 28 2012
E.g.f.: exp(2*y*x + x). - Geoffrey Critzer, Nov 12 2012
Riordan array (x/(1 - x), 2*x/(1 - x)). Exp(2*x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(2*x)*(1 + 6*x + 12*x^2/2! + 8*x^3/3!) = 1 + 8*x + 40*x^2/2! + 160*x^3/3! + 560*x^4/4! + .... The same property holds more generally for Riordan arrays of the form (f(x), 2*x/(1 - x)). - Peter Bala, Dec 21 2014
T(n,k) = Sum_{j=0..k} (-1)^(k-j) * binomial(n,k) * binomial(k,j) * 3^j. - Kolosov Petro, Jan 28 2019
T(n,k) = 2*(n+1-k)*T(n,k-1)/k, T(n,0) = 1. - Alexander R. Povolotsky, Oct 08 2023
For n >= 1, GCD(T(n,1), ..., T(n,n)) = GCD(T(n,1),T(n,n)) = GCD(2*n,2^n) = A171977(n). - Pontus von Brömssen, Nov 01 2024

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

A029858 a(n) = (3^n - 3)/2.

Original entry on oeis.org

0, 3, 12, 39, 120, 363, 1092, 3279, 9840, 29523, 88572, 265719, 797160, 2391483, 7174452, 21523359, 64570080, 193710243, 581130732, 1743392199, 5230176600, 15690529803, 47071589412, 141214768239

Offset: 1

Christian G. Bower

Also the number of 2-block covers of a labeled n-set. a(n) = A055154(n,2). Generally, number of k-block covers of a labeled n-set is T(n,k) = (1/k!)*Sum_{i = 1..k + 1} Stirling1(k + 1,i)*(2^(i - 1) - 1)^n. In particular, T(n,2) = (1/2!)*(3^n - 3), T(n,3) = (1/3!)*(7^n - 6*3^n + 11), T(n,4) = (1/4)!*(15^n - 10*7^n + 35*3^n - 50), ... - Vladeta Jovovic, Jan 19 2001
Conjectured to be the number of integers from 0 to 10^(n-1) - 1 that lack 0, 1, 2, 3, 4, 5 and 6 as a digit. - Alexandre Wajnberg, Apr 25 2005. This is easily verified to be true. - Renzo Benedetti, Sep 25 2008
Number of monic irreducible polynomials of degree 1 in GF(3)[x1,...,xn]. - Max Alekseyev, Jan 23 2006
Also, the greatest number of identical weights among which an odd one can be identified and it can be decided if the odd one is heavier or lighter, using n weighings with a comparing balance. If the odd one only needs to be identified, the sequence starts 4, 13, 40 and is A003462 (3^n - 1)/2, n > 1. - Tanya Khovanova, Dec 11 2006; corrected by Samuel E. Rhoads, Apr 18 2016
Binomial transform yields A134057. Inverse binomial transform yields A062510 with one additional 0 in front. - R. J. Mathar, Jun 18 2008
Numbers k where the recurrence s(0)=0, if s(k-1) >= k then s(k) = s(k-1) - k otherwise s(k) = s(k-1) + k produces s(k) = 0. - Hugo Pfoertner, Jan 05 2012
For n > 1: A008344(a(n)) = a(n). - Reinhard Zumkeller, May 09 2012
Also the number of edges in the (n-1)-Hanoi graph. - Eric W. Weisstein, Jun 18 2017
A level 1 Sierpiński triangle graph is a triangle. Level n+1 is formed from three copies of level n by identifying pairs of corner vertices of each pair of triangles.  a(n) is the number of degree 4 vertices in the level n Sierpinski triangle graph. - Allan Bickle, Jul 30 2020
Also the number of minimum vertex cuts in the n-Apollonian network. - Eric W. Weisstein, Dec 20 2020
Also the minimum number of turns in n-dimensional Euclidean space needed to visit all 3^n points of the grid {0, 1, 2}^n, moving in straight lines between turns (repeated visits and direction changes at non-grid points are allowed). - Marco Ripà, Aug 06 2025

			For the Sierpiński triangle, Level 1 is a triangle, so a(1) = 0.
Level 2 has three corners (degree 2) and three degree 4 vertices, so a(2) = 3.
The level 2 Hanoi graph has 3 triangles joined by 3 edges, so a(2+1) = 12.

Vincenzo Librandi, Table of n, a(n) for n = 1..300
Allan Bickle, Properties of Sierpinski Triangle Graphs, Springer PROMS 448 (2021) 295-303.
Natalia Agudelo Muñetón, Agustín Moreno Cañadas, Pedro Fernando Fernández Espinosa, and Isaías David Marín Gaviria, Brauer Configuration Algebras and Their Applications in Graph Energy Theory, Mathematics (2021) Vol. 9, 3042.
Alex Born, Cor A. J. Hukrnes, and Gerhard J. Woeginger, How to detect a counterfeit coin: adaptive versus non-adaptive solutions, Inf. Proc. Lett. 86 (2003) 137-141.
Gary Darby, The Counterfeit Coin
Madeleine Goertz and Aaron Williams, The Quaternary Gray Code and How It Can Be Used to Solve Ziggurat and Other Ziggu Puzzles, arXiv:2411.19291 [math.CO], 2024. See p. 17.
Lorenz Halbeisen and Norbert Hungerbuhler, The general counterfeit coin problem, Discr. Math 147 (1-3) (1995) 139-150, Theorem 1 with b=1.
Andreas Hinz, Sandi Klavzar, and Sara Sabrina Zemljic, A survey and classification of Sierpinski-type graphs, Discrete Applied Mathematics 217 3 (2017), 565-600.
Bennet Manvel, Counterfeit coin problems, Math. Mag. 50 (2) (1977) 90-92, theorem 2.
Marco Ripà, Solving the 106 years old 3^k points problem with the clockwise-algorithm, Journal of Fundamental Mathematics and Applications, 2020, 3(2), 84-97.
Allen Stenger and Jack Wert, The Twelve Coins (or Twelve bags of Gold)
Eric Weisstein's World of Mathematics, Apollonian Network
Eric Weisstein's World of Mathematics, Edge Count
Eric Weisstein's World of Mathematics, Hanoi Graph
Eric Weisstein's World of Mathematics, Minimum Vertex Cut
Index entries for linear recurrences with constant coefficients, signature (4,-3).

Cf. A055154, A003462, A007051, A034472, A024023.
Cf. A000203, A008776.
Cf. A007283, A029858, A067771, A233774, A233775, A246959 (Sierpiński triangle graphs).
Cf. A000225, A029858, A058809, A375256 (Hanoi graphs).

a029858 = (`div` 2) . (subtract 3) . (3 ^)
a029858_list = iterate ((+ 3) . (* 3)) 0
-- Reinhard Zumkeller, May 09 2012

[(3^n-3)/2: n in [1..30]]; // Vincenzo Librandi, Jun 05 2011

a:=n->sum(3^j,j=1..n): seq(a(n),n=0..23); # Zerinvary Lajos, Jun 27 2007

Table[(3^n - 3)/2, {n, 24}] (* Alonso del Arte, Dec 29 2014 *)

a(n)=(3^n-3)\2 \\ Charles R Greathouse IV, Apr 17 2012

a(n) = 3*a(n-1) + 3. - Alexandre Wajnberg, Apr 25 2005
O.g.f: 3*x^2/((1-x)*(1-3*x)). - R. J. Mathar, Jun 18 2008
a(n) = 3^(n-1) + a(n-1) (with a(1)=0). - Vincenzo Librandi, Nov 18 2010
a(n) = 3*A003462(n-1). - R. J. Mathar, Sep 10 2015
E.g.f.: 3*(-1 + exp(2*x))*exp(x)/2. - Ilya Gutkovskiy, Apr 19 2016
a(n) = A067771(n-1) - 3. - Allan Bickle, Jul 30 2020
a(n) = sigma(A008776(n-2)) for n>=2. - Flávio V. Fernandes, Apr 20 2021

Corrected by T. D. Noe, Nov 07 2006

A058481 a(n) = 3^n - 2.

Original entry on oeis.org

1, 7, 25, 79, 241, 727, 2185, 6559, 19681, 59047, 177145, 531439, 1594321, 4782967, 14348905, 43046719, 129140161, 387420487, 1162261465, 3486784399, 10460353201, 31381059607, 94143178825, 282429536479, 847288609441

Offset: 1

Vladeta Jovovic, Nov 26 2000

a(n) = number of 2 X n binary matrices with no zero rows or columns.
a(n)^2 + 2*a(n+1) + 1 is a square number, i.e., a(n)^2 + 2*a(n+1) + 1 = (a(n)+3)^2: for n=2, a(2)^2 + 2*a(3) + 1 = 7^2 + 2*25 + 1 = 100 = (7+3)^2; for n=3, a(3)^2 + 2*a(4) + 1 = 25^2 + 2*79 + 1 = 784 = (25+3)^2. - Bruno Berselli, Apr 23 2010
Sum of n-th row of triangle of powers of 3: 1; 3 1 3; 9 3 1 3 9; 27 9 3 1 3 9 27; ... . - Philippe Deléham, Feb 24 2014
a(n) = least k such that k*3^n + 1 is a square. Thus, the square is given by (3^n-1)^2. - Derek Orr, Mar 23 2014
Binomial transform of A058481: (1, 6, 12, 24, 48, 96, ...) and second binomial transform of (1, 5, 1, 5, 1, 5, ...). - Gary W. Adamson, Aug 24 2016
Number of ordered pairs of nonempty sets whose union is [n]. a(2) = 7: ({1,2},{1,2}), ({1,2},{1}), ({1,2},{2}), ({1},{1,2}), ({1},{2}), ({2},{1,2}), ({2},{1}). If "nonempty" is omitted we get A000244. - Manfred Boergens, Mar 29 2023

			G.f. = x + 7*x^2 + 25*x^3 + 79*x^4 + 241*x^5 + 727*x^6 + 2185*x^7 + 6559*x^8 + ...
a(1) = 1;
a(2) = 3 + 1 + 3 = 7;
a(3) = 9 + 3 + 1 + 3 + 9 = 25;
a(4) = 27 + 9 + 3 + 1 + 3 + 9 + 27 = 79; etc. - _Philippe Deléham_, Feb 24 2014

G. C. Greubel, Table of n, a(n) for n = 1..1000
John Elias, Illustration: Union of Triple-Sierpinski Triangles
Index entries for linear recurrences with constant coefficients, signature (4,-3).

Cf. A055602, A024206 (unlabeled case), A055609, A058482, A000244.

Cf. A003462, A007051, A034472, A024023, A067771, A029858, A134931, A115099, A100774, A079004.

A058481:=n->3^n-2; seq(A058481(n), n=1..30); # Wesley Ivan Hurt, Mar 24 2014

a=1;lst={a};Do[a=a*3+4;AppendTo[lst,a],{n,0,5!}];lst (* Vladimir Joseph Stephan Orlovsky, Dec 25 2008 *)
3^Range[30]-2  (* Harvey P. Dale, Mar 28 2011 *)
LinearRecurrence[{4, -3}, {1, 7}, 25] (* G. C. Greubel, Aug 25 2016 *)

a(n)=3^n-2 \\ Charles R Greathouse IV, Feb 06 2017

{a(n) = if( n<1, 0, 3^n - 2)}; /* Michael Somos, Feb 17 2017 */

Number of m X n binary matrices with no zero rows or columns is Sum_{j=0..m} (-1)^j*C(m, j)*(2^(m-j)-1)^n.
From Mohammad K. Azarian, Jan 14 2009: (Start)
G.f.: 1/(1-3*x)-2/(1-x)+1.
E.g.f.: e^(3*x)-2*(e^x)+1. (End)
a(n) = 3*a(n-1) + 4 (with a(1)=1). - Vincenzo Librandi, Aug 07 2010
a(n) = 4*a(n-1) - 3*a(n-2). - G. C. Greubel, Aug 25 2016

More terms from Larry Reeves (larryr(AT)acm.org), Dec 04 2000

cp's OEIS Frontend

A062318 Numbers of the form 3^m - 1 or 2*3^m - 1; i.e., the union of sequences A048473 and A024023.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

SageMath

Formula

Extensions

A064079 Zsigmondy numbers for a = 3, b = 1: Zs(n, 3, 1) is the greatest divisor of 3^n - 1^n (A024023) that is relatively prime to 3^m - 1^m for all positive integers m < n.

Views

Author

Keywords

Comments

Links

Crossrefs

Extensions

A034472 a(n) = 3^n + 1.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

Sage

Sage

Sage

Formula

Extensions

A052548 a(n) = 2^n + 2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Formula

Extensions

A013609 Triangle of coefficients in expansion of (1+2*x)^n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Haskell

Magma

Maple

Mathematica

Maxima

PARI

Sage

Formula

A374848 Obverse convolution A000045**A000045; see Comments.