A000420 -id:A000420 - cp's OEIS frontend

A083098 a(n) = 2a(n-1) + 6a(n-2).

1, 1, 8, 22, 92, 316, 1184, 4264, 15632, 56848, 207488, 756064, 2757056, 10050496, 36643328, 133589632, 487039232, 1775616256, 6473467904, 23600633344, 86042074112, 313687948288, 1143628341248, 4169384372224, 15200538791936

Offset: 0

Mario Catalani (mario.catalani(AT)unito.it), Apr 22 2003

a(n+1) = a(n) + 7*A083099(n-1); a(n+1)/A083099(n) converges to sqrt(7).
Binomial transform of expansion of cosh(sqrt(7)x) (A000420 with interpolated zeros: 1, 0, 7, 0, 49, 0, 343, 0, ...).
The same sequence may be obtained by the following process. Starting a priori with the fraction 1/1, the numerators of fractions built according to the rule: add top and bottom to get the new bottom, add top and 7 times the bottom to get the new top. The limit of the sequence of fractions is sqrt(7). - Cino Hilliard, Sep 25 2005
a(n) is the number of compositions of n when there are 1 type of 1 and 7 types of other natural numbers. - Milan Janjic, Aug 13 2010

John Derbyshire, Prime Obsession, Joseph Henry Press, April 2004, see p. 16.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Project Euler, Problem 752, sequence alpha(n).
Index entries for linear recurrences with constant coefficients, signature (2,6).

The following sequences (and others) belong to the same family: A001333, A000129, A026150, A002605, A046717, A015518, A084057, A063727, A002533, A002532, A083098, A083099, A083100, A015519.

I:=[1,1]; [n le 2 select I[n] else 2*Self(n-1) + 6*Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 08 2018

CoefficientList[Series[(1+6x)/(1-2x-6x^2), {x, 0, 25}], x]
LinearRecurrence[{2, 6}, {1, 1}, 25] (* Sture Sjöstedt, Dec 06 2011 *)
a[n_] := Simplify[((1 + Sqrt[7])^n + (1 - Sqrt[7])^n)/2]; Array[a, 25, 0] (* Robert G. Wilson v, Sep 18 2013 *)

x='x+O('x^30); Vec((1-x)/(1-2*x-6*x^2)) \\ G. C. Greubel, Jan 08 2018

[lucas_number2(n,2,-6)/2 for n in range(0, 25)] # Zerinvary Lajos, Apr 30 2009

G.f.: (1-x)/(1-2*x-6*x^2).
a(n) = (1+sqrt(7))^n/2 + (1-sqrt(7))^n/2.
E.g.f.: exp(x)*cosh(sqrt(7)x).
a(n) = Sum_{k=0..n} A098158(n,k)*7^(n-k). - Philippe Deléham, Dec 26 2007
If p[1]=1, and p[i]=7, (i>1), and if A is Hessenberg matrix of order n defined by: A[i,j]=p[j-i+1], (i<=j), A[i,j]=-1, (i=j+1), and A[i,j]=0 otherwise. Then, for n>=1, a(n) = det A. - Milan Janjic, Apr 29 2010
G.f.: G(0)/2, where G(k)= 1 + 1/(1 - x*(7*k-1)/(x*(7*k+6) - 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 26 2013

A329332 Table of powers of squarefree numbers, powers of A019565(n) in increasing order in row n. Square array A(n,k) n >= 0, k >= 0 read by descending antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 8, 9, 6, 1, 1, 16, 27, 36, 5, 1, 1, 32, 81, 216, 25, 10, 1, 1, 64, 243, 1296, 125, 100, 15, 1, 1, 128, 729, 7776, 625, 1000, 225, 30, 1, 1, 256, 2187, 46656, 3125, 10000, 3375, 900, 7, 1, 1, 512, 6561, 279936, 15625, 100000, 50625, 27000, 49, 14

Offset: 0

Peter Munn, Nov 10 2019

The A019565 row order gives the table neat relationships with A003961, A003987, A059897, A225546, A319075 and A329050. See the formula section.

Transposition of this table, that is reflection about its main diagonal, has subtle symmetries. For example, consider the unique factorization of a number into powers of distinct primes. This can be restated as factorization into numbers from rows 2^n (n >= 0) with no more than one from each row. Reflecting about the main diagonal, this factorization becomes factorization (of a related number) into numbers from columns 2^k (k >= 0) with no more than one from each column. This is also unique and is factorization into powers of squarefree numbers with distinct exponents that are powers of two. See the example section.

			Square array A(n,k) begins:
n\k |  0   1     2      3        4          5           6             7
----+------------------------------------------------------------------
   0|  1   1     1      1        1          1           1             1
   1|  1   2     4      8       16         32          64           128
   2|  1   3     9     27       81        243         729          2187
   3|  1   6    36    216     1296       7776       46656        279936
   4|  1   5    25    125      625       3125       15625         78125
   5|  1  10   100   1000    10000     100000     1000000      10000000
   6|  1  15   225   3375    50625     759375    11390625     170859375
   7|  1  30   900  27000   810000   24300000   729000000   21870000000
   8|  1   7    49    343     2401      16807      117649        823543
   9|  1  14   196   2744    38416     537824     7529536     105413504
  10|  1  21   441   9261   194481    4084101    85766121    1801088541
  11|  1  42  1764  74088  3111696  130691232  5489031744  230539333248
  12|  1  35  1225  42875  1500625   52521875  1838265625   64339296875
Reflection of factorization about the main diagonal: (Start)
The canonical (prime power) factorization of 864 is 2^5 * 3^3 = 32 * 27. Reflecting the factors about the main diagonal of the table gives us 10 * 36 = 10^1 * 6^2 = 360. This is the unique factorization of 360 into powers of squarefree numbers with distinct exponents that are powers of two.
Reflection about the main diagonal is given by the self-inverse function A225546(.). Clearly, all positive integers are in the domain of A225546, whether or not they appear in the table. It is valid to start from 360, observe that A225546(360) = 864, then use 864 to derive 360's factorization into appropriate powers of squarefree numbers as above.
(End)

The range of values is A072774.
Rows (abbreviated list): A000079(1), A000244(2), A000400(3), A000351(4), A011557(5), A001024(6), A009974(7), A000420(8), A001023(9), A009965(10), A001020(16), A001022(32), A001026(64).
A019565 is column 1, A334110 is column 2, and columns that are sorted in increasing order (some without the 1) are: A005117(1), A062503(2), A062838(3), A113849(4), A113850(5), A113851(6), A113852(7).
Other subtables: A182944, A319075, A329050.
Re-ordered subtable of A297845, A306697, A329329.
A000290, A003961, A003987, A059897 and A225546 are used to express relationships between terms of this sequence.
Cf. A285322.

A(n,k) = A019565(n)^k.
A(k,n) = A225546(A(n,k)).
A(n,2k) = A000290(A(n,k)) = A(n,k)^2.
A(2n,k) = A003961(A(n,k)).
A(n,2k+1) = A(n,2k) * A(n,1).
A(2n+1,k) = A(2n,k) * A(1,k).
A(A003987(n,m), k) = A059897(A(n,k), A(m,k)).
A(n, A003987(m,k)) = A059897(A(n,m), A(n,k)).
A(2^n,k) = A319075(k,n+1).
A(2^n, 2^k) = A329050(n,k).
A(n,k) = A297845(A(n,1), A(1,k)) = A306697(A(n,1), A(1,k)), = A329329(A(n,1), A(1,k)).
Sum_{n>=0} 1/A(n,k) = zeta(k)/zeta(2*k), for k >= 2. - Amiram Eldar, Dec 03 2022

A076512 Denominator of cototient(n)/totient(n).

Original entry on oeis.org

1, 1, 2, 1, 4, 1, 6, 1, 2, 2, 10, 1, 12, 3, 8, 1, 16, 1, 18, 2, 4, 5, 22, 1, 4, 6, 2, 3, 28, 4, 30, 1, 20, 8, 24, 1, 36, 9, 8, 2, 40, 2, 42, 5, 8, 11, 46, 1, 6, 2, 32, 6, 52, 1, 8, 3, 12, 14, 58, 4, 60, 15, 4, 1, 48, 10, 66, 8, 44, 12, 70, 1, 72, 18, 8, 9, 60, 4, 78, 2, 2, 20, 82, 2, 64, 21

Offset: 1

Reinhard Zumkeller, Oct 15 2002

a(n)=1 iff n=A007694(k) for some k.
Numerator of phi(n)/n=Prod_{p|n} (1-1/p). - Franz Vrabec, Aug 26 2005
From Wolfdieter Lang, May 12 2011: (Start)
For n>=2, a(n)/A109395(n) = sum(((-1)^r)*sigma_r,r=0..M(n)) with the elementary symmetric functions (polynomials) sigma_r of the indeterminates {1/p_1,...,1/p_M(n)} if n = prod((p_j)^e(j),j=1..M(n)) where M(n)=A001221(n) and sigma_0=1.
This follows by expanding the above given product for phi(n)/n.
The n-th member of this rational sequence 1/2, 2/3, 1/2, 4/5, 1/3, 6/7, 1/2, 2/3, 2/5,... is also (2/n^2)*sum(k,with 1<=k=2.
Therefore, this scaled sum depends only on the distinct prime factors of n.
See also A023896. Proof via PIE (principle of inclusion and exclusion). (End)
In the sequence of rationals r(n)=eulerphi(n)/n: 1, 1/2, 2/3, 1/2, 4/5, 1/3, 6/7, 1/2, 2/3, 2/5, 10/11, 1/3, ... one can observe that new values are obtained for squarefree indices (A005117); while for a nonsquarefree number n (A013929), r(n) = r(A007947(n)), where A007947(n) is the squarefree kernel of n. - Michel Marcus, Jul 04 2015

T. D. Noe, Table of n, a(n) for n=1..1000

Cf. A076511 (numerator of cototient(n)/totient(n)), A051953.

Phi(m)/m = k: A000079 \ {1} (k=1/2), A033845 (k=1/3), A000244 \ {1} (k=2/3), A033846 (k=2/5), A000351 \ {1} (k=4/5), A033847 (k=3/7), A033850 (k=4/7), A000420 \ {1} (k=6/7), A033848 (k=5/11), A001020 \ {1} (k=10/11), A288162 (k=6/13), A001022 \ {1} (12/13), A143207 (k=4/15), A033849 (k=8/15), A033851 (k=24/35).

[Numerator(EulerPhi(n)/n): n in [1..100]]; // Vincenzo Librandi, Jul 04 2015

Table[Denominator[(n - EulerPhi[n])/EulerPhi[n]], {n, 80}] (* Alonso del Arte, May 12 2011 *)

vector(80, n, numerator(eulerphi(n)/n)) \\ Michel Marcus, Jul 04 2015

a(n) = A000010(n)/A009195(n).

A064628 a(n) = floor((4/3)^n).

Original entry on oeis.org

1, 1, 1, 2, 3, 4, 5, 7, 9, 13, 17, 23, 31, 42, 56, 74, 99, 133, 177, 236, 315, 420, 560, 747, 996, 1328, 1771, 2362, 3149, 4199, 5599, 7466, 9954, 13273, 17697, 23596, 31462, 41950, 55933, 74577, 99437, 132583, 176777, 235703, 314271, 419028, 558704

Offset: 0

Labos Elemer, Oct 01 2001

a(n) is the perimeter of a hexaflake (rounded down) after n iterations. The total number of holes = A000420(n) - 1. The total number of irregular polygon holes = A000420(n-1) - 1. The total number of triangle holes = 6*A000420(n-1). - Kival Ngaokrajang, Apr 18 2014

a(n) is composite infinitely often (Forman and Shapiro). More exactly, a(n) is divisible by at least one of 2, 3, 5 infinitely often (Dubickas and Novikas). - Tomohiro Yamada, Apr 15 2017

R. K. Guy, Unsolved Problems in Number Theory, E19.

Harry J. Smith, Table of n, a(n) for n=0,...,400
Arturas Dubickas and Aivaras Novikas, Integer parts of powers of rational numbers, Math. Z. 251 (2005), 635--648, available from the first author's page.
W. Forman and H. N. Shapiro, An arithmetic property of certain rational powers, Comm. Pure. Appl. Math. 20 (1967), 561-573.
Kival Ngaokrajang, Illustration of hexaflake for n = 0..3.
Eric Weisstein's World of Mathematics, Power Floors.
Wikipedia, n-flake.

Cf. A002379, A002380, A060692.
Cf. A094969 - A094500.
Cf. A046038, A070761, A070762, A067905 (Composites and Primes).

A064628:=n->floor(4^n/3^n); seq(A064628(n), n=0..30); # Wesley Ivan Hurt, Apr 19 2014

Table[Floor[(4/3)^n], {n, 0, 30}] (* Robert G. Wilson v *)

{ f=t=1; for (n=0, 400, write("b064628.txt", n, " ", f\t); f*=4; t*=3 ) } \\ Harry J. Smith, Sep 20 2009

def A064628(n): return floor((4/3)^n)
[A064628(n) for n in range(0, 46)] # Stefano Spezia, Oct 13 2024

More terms from Robert G. Wilson v, May 26 2004

OFFSET changed from 1 to 0 by Harry J. Smith, Sep 20 2009

A011754 Number of ones in the binary expansion of 3^n.

Original entry on oeis.org

1, 2, 2, 4, 3, 6, 6, 5, 6, 8, 9, 13, 10, 11, 14, 15, 11, 14, 14, 17, 17, 20, 19, 22, 16, 18, 24, 30, 25, 25, 25, 26, 26, 34, 29, 32, 27, 34, 36, 32, 28, 39, 38, 39, 34, 34, 45, 38, 41, 33, 41, 46, 42, 35, 39, 42, 39, 40, 42, 48, 56, 56, 49, 57, 56, 51, 45, 47, 55, 55, 64, 68, 58

Offset: 0

Allan C. Wechsler, Dec 11 1999

Conjecture: a(n)/n tends to log(3)/(2*log(2)) = 0.792481250... (A094148). - Ed Pegg Jr, Dec 05 2002
Senge & Straus prove that for every m, there is some N such that for all n > N, a(n) > m. Dimitrov & Howe make this effective, proving that for n > 25, a(n) > 22. - Charles R Greathouse IV, Aug 23 2021
Ed Pegg's conjecture means that about half of the bits of 3^n are nonzero. It appears that the same is true for 5^n (A000351, cf. A118738) and 7^n (A000420). - M. F. Hasler, Apr 17 2024

S. Wolfram, "A new kind of science", p. 903.

Hugo Pfoertner, Table of n, a(n) for n = 0..10000 (terms 0..1000 from T. D. Noe)
Vassil S. Dimitrov and Everett W. Howe, Powers of 3 with few nonzero bits and a conjecture of Erdős, arXiv:2105.06440 [math.NT], 2021.
Taylor Dupuy and David E. Weirich, Bits of 3^n in binary, Wieferich primes and a conjecture of Erdős, Journal of Number Theory, Volume 158, January 2016, Pages 268-280.
Hugo Pfoertner, Plot of a(n) - 0.79248*n, +-Pi*sqrt(n), n up to 10^6.
H. G. Senge and E. G. Straus, PV-numbers and sets of multiplicity, Periodica Mathematica Hungarica 3 (1973), pp. 93-100.

Cf. A007088, A000120 (Hamming weight), A000244 (3^n), A004656, A261009, A094148.

Cf. A118738 (same for 5^n).

a011754 = a000120 . a000244  -- Reinhard Zumkeller, Aug 14 2015

[&+Intseq(3^n, 2): n in [0..79]]; // Vincenzo Librandi, Nov 28 2018

f:= n -> convert(convert(3^n,base,2),`+`):
map(f, [$0..100]); # Robert Israel, Apr 17 2024

Table[DigitCount[3^n, 2][[1]], {n, 0, 100}] (* Stefan Steinerberger, Apr 03 2006 *)
DigitCount[3^Range[0,100],2,1] (* Harvey P. Dale, Apr 06 2012 *)

a(n)=hammingweight(3^n) \\ Charles R Greathouse IV, Feb 09 2017

A011754 = lambda n: (3**n).bit_count() # M. F. Hasler, Apr 17 2024

a(n) = A000120(3^n). - Benoit Cloitre, Dec 06 2002

a(n) = A000120(A000244(n)). - Reinhard Zumkeller, Aug 14 2015

More terms from Stefan Steinerberger, Apr 03 2006

A066003 Sum of digits of 7^n.

Original entry on oeis.org

1, 7, 13, 10, 7, 22, 28, 25, 31, 28, 43, 49, 37, 52, 58, 64, 52, 58, 73, 79, 76, 82, 97, 85, 73, 97, 112, 91, 133, 121, 118, 115, 103, 127, 142, 157, 136, 115, 130, 136, 142, 148, 136, 169, 175, 163, 187, 175, 136, 178, 184, 217, 196, 220, 217

Offset: 0

N. J. A. Sloane, Dec 11 2001

Harry J. Smith, Table of n, a(n) for n=0,...,1000

Cf. A000420 (7^n), A007953 (sum of digits).

Cf. sum of digits of k^n: A001370 (k=2), A004166 (k=3), A065713 (k=4), A066001 (k=5), A066002 (k=6), this sequence (k=7), A066004 (k=8), A065999 (k=9), A066005 (k=11), A066006 (k=12), A175527 (k=13).

Magma
```
[ &+Intseq(7^n): n in [0..60] ];
```

Table[Total[IntegerDigits[7^n]],{n,55}] (* Harvey P. Dale, Nov 22 2010 *)

a(n) = sumdigits(7^n); \\ Michel Marcus, Nov 01 2013

a(n) = A007953(A000420(n)). - Michel Marcus, Nov 01 2013

A027474 a(n) = 7^(n-2) * C(n,2).

Original entry on oeis.org

1, 21, 294, 3430, 36015, 352947, 3294172, 29647548, 259416045, 2219448385, 18643366434, 154231485954, 1259557135291, 10173346092735, 81386768741880, 645668365352248, 5084638377148953, 39779817891812397, 309398583602985310

Offset: 2

Olivier Gérard

7th binomial transform of (0,0,1,0,0,0,........). Starting at 1, the three-fold convolution of A000420 (powers of 7). - Paul Barry, Mar 08 2003

Vincenzo Librandi, Table of n, a(n) for n = 2..400
Index entries for linear recurrences with constant coefficients, signature (21,-147,343).

Third column of A027466.

Sequences similar to the form q^(n-2)*binomial(n, 2): A000217 (q=1), A001788 (q=2), A027472 (q=3), A038845 (q=4), A081135 (q=5), A081136 (q=6), this sequence (q=7), A081138 (q=8), A081139 (q=9), A081140 (q=10), A081141 (q=11), A081142 (q=12), A027476 (q=15).

[7^(n-2)* Binomial(n, 2): n in [2..20]]; /* Vincenzo Librandi, Oct 12 2011 */

seq(binomial(n, 2)*7^(n-2), n=2..30); # Zerinvary Lajos, Jun 12 2008

Table[7^(n-2) Binomial[n,2], {n,2,20}] (* Harvey P. Dale, Sep 25 2011 *)

a(n)=7^(n-2)*n*(n-1)/2 \\ Charles R Greathouse IV, Oct 07 2015

[7^(n-2)*binomial(n,2) for n in range(2, 21)] # Zerinvary Lajos, Mar 13 2009

From Paul Barry, Mar 08 2003: (Start)
G.f.: x^2 / (1-7*x)^3.
a(n) = 21*a(n-1) - 147*a(n-2) + 343*a(n-3), a(0) = a(1) = 0, a(2) = 1. (End)
Numerators of sequence a[3,n] in (a[i,j])^3 where a[i,j] = binomial(i-1, j-1)/2^(i-1) if j<=i, 0 if j>i.
E.g.f.: (x^2/2)*exp(7*x). - G. C. Greubel, May 13 2021
From Amiram Eldar, Jan 06 2022: (Start)
Sum_{n>=2} 1/a(n) = 14 - 84*log(7/6).
Sum_{n>=2} (-1)^n/a(n) = 112*log(8/7) - 14. (End)

Edited by Ralf Stephan, Dec 30 2004

A016177 a(n) = 8^n - 7^n.

Original entry on oeis.org

0, 1, 15, 169, 1695, 15961, 144495, 1273609, 11012415, 93864121, 791266575, 6612607849, 54878189535, 452866803481, 3719823438255, 30436810578889, 248242046141055, 2019169299698041, 16385984911571535, 132716292890482729, 1073129238309234975, 8664826172771491801

Offset: 0

N. J. A. Sloane

Number of ways to assign truth values to n ternary conjunctions connected by disjunctions such that the proposition is true. For example, a(2) = 15, since for the proposition '(a & b & c) v (d & e & f)' there are 15 assignments that make the proposition true. - Ori Milstein, Dec 22 2022
Equivalently, the number of length-n words over the alphabet {0,1,...,7} with at least one letter = 7. - Joerg Arndt, Jan 01 2023
a(n) is also the number of n-digit numbers whose smallest decimal digit is 2. - Stefano Spezia, Nov 15 2023

G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (15,-56).

Cf. A000420, A001018, A016140.

[8^n -7^n: n in [0..40]]; // G. C. Greubel, Nov 29 2024

Table[8^n - 7^n, {n, 0, 20}] (* Harvey P. Dale, Jan 31 2011 *)

a(n)=8^n-7^n \\ Charles R Greathouse IV, Sep 28 2015

def A016177(n): return pow(8,n) - pow(7,n)
print([A016177(n) for n in range(41)]) # G. C. Greubel, Nov 29 2024

G.f.: x/((1-7x)*(1-8x)).
a(n) = numerator(f(n-1)) where f(n) = Integral_{x=0..1/4} (1-x/2)^n dx. And denominator(f(n)) = 4*(n+1)*8^n. - Al Hakanson (hawkuu(AT)excite.com), Feb 22 2004 [corrected by Michel Marcus, Dec 23 2022]
a(n) = 15*a(n-1) - 56*a(n-2), n > 1. - Philippe Deléham, Jan 01 2009
E.g.f.: e^(8*x) - e^(7*x). - Mohammad K. Azarian, Jan 14 2009
a(n) = 8*a(n-1) + 7^(n-1), a(0)=0. - Vincenzo Librandi, Feb 09 2011

A024075 a(n) = 7^n-1.

Original entry on oeis.org

0, 6, 48, 342, 2400, 16806, 117648, 823542, 5764800, 40353606, 282475248, 1977326742, 13841287200, 96889010406, 678223072848, 4747561509942, 33232930569600, 232630513987206, 1628413597910448, 11398895185373142

Offset: 0

N. J. A. Sloane

In base 7 these are 0, 6, 66, 666, ... - David Rabahy, Dec 12 2016

Vincenzo Librandi, Table of n, a(n) for n = 0..500
Amelia Carolina Sparavigna, The groupoids of Mersenne, Fermat, Cullen, Woodall and other Numbers and their representations by means of integer sequences, Politecnico di Torino, Italy (2019), [math.NT].
Index entries for linear recurrences with constant coefficients, signature (8,-7).

Cf. A000420, A248724.

[7^n-1: n in [0..25]]; // Vincenzo Librandi, Jul 03 2011

CoefficientList[Series[1/(1 - 7 x) - 1/(1 - x), {x, 0, 19}], x] (* Michael De Vlieger, Jan 04 2020 *)

makelist(7^n-1, n, 0, 30); /* Martin Ettl, Nov 12 2012 */

a(n)=7^n-1 \\ Charles R Greathouse IV, Dec 12 2016

G.f.: 1/(1-7*x)-1/(1-x). - Mohammad K. Azarian, Jan 14 2009
E.g.f.: exp(7*x)-exp(x). - Mohammad K. Azarian, Jan 14 2009
a(n+1) = 7*a(n) + 6. - Reinhard Zumkeller, Nov 22 2009
a(n) = Sum_{i=1..n} 6^i*binomial(n,n-i) for n>0, a(0)=0. - Bruno Berselli, Nov 11 2015
a(n) = A000420(n) - 1. - Sean A. Irvine, Jun 19 2019
Sum_{n>=1} 1/a(n) = A248724. - Amiram Eldar, Nov 13 2020

A063265 Septinomial (also called heptanomial) coefficient array.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 4, 5, 6, 7, 6, 5, 4, 3, 2, 1, 1, 3, 6, 10, 15, 21, 28, 33, 36, 37, 36, 33, 28, 21, 15, 10, 6, 3, 1, 1, 4, 10, 20, 35, 56, 84, 116, 149, 180, 206, 224, 231, 224, 206, 180, 149, 116, 84, 56, 35

Offset: 0

Wolfdieter Lang, Jul 24 2001

The sequence of step width of this staircase array is [1,6,6,...], hence the degree sequence for the row polynomials is [0,6,12,18,...]= A008588.

The column sequences (without leading zeros) are for k=0..6 those of the lower triangular array A007318 (Pascal) and for k=7..9: A063267, A063417, A063418. Row sums give A000420 (powers of 7). Central coefficients give A025012.

			Triangle begins:
  {1};
  {1, 1, 1, 1, 1, 1, 1};
  {1, 2, 3, 4, 5, 6, 7, 6, 5, 4, 3, 2, 1};
  ...
N7(k,x)= 1 for k=0..6, N7(7,x)= 6-15*x+20*x^2-15*x^3+6*x^4-x^5 (from A063266).

L. Comtet, Advanced Combinatorics, Reidel, 1974, pp. 77,78.

T. D. Noe, Rows n = 0..25, flattened
Steven R. Finch, P. Sebah and Z.-Q. Bai, Odd Entries in Pascal's Trinomial Triangle, arXiv:0802.2654 [math.NT], 2008.

The q-nomial arrays are for q=2..8: A007318 (Pascal), A027907, A008287, A035343, A063260, A063265, A171890.

#Define the r-nomial coefficients for r = 1, 2, 3, ...
rnomial := (r,n,k) -> add((-1)^i*binomial(n,i)*binomial(n+k-1-r*i,n-1), i = 0..floor(k/r)):
#Display the 7-nomials as a table
r := 7:  rows := 10:
for n from 0 to rows do
seq(rnomial(r,n,k), k = 0..(r-1)*n)
end do;
# Peter Bala, Sep 07 2013

Flatten[Table[CoefficientList[(1 + x + x^2 + x^3 + x^4 + x^5 + x^6)^n, x], {n, 0, 25}]] (* T. D. Noe, Apr 04 2011 *)

a(n, k)=0 if n=-1 or k<0 or k >= 6*n; a(0, 0)=1; a(n, k)= sum(a(n-1, k-j), j=0..6) else.
G.f. for row n: (sum(x^j, j=0..6))^n.
G.f. for column k: (x^(ceiling(k/6)))*N7(k, x)/(1-x)^(k+1) with the row polynomials of the staircase array A063266(k, m).
T(n,k) = Sum_{i = 0..floor(k/7)} (-1)^i*binomial(n,i)*binomial(n+k-1-7*i,n-1) for n >= 0 and 0 <= k <= 6*n. - Peter Bala, Sep 07 2013

cp's OEIS Frontend

A083098 a(n) = 2*a(n-1) + 6*a(n-2).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Sage

Formula

A329332 Table of powers of squarefree numbers, powers of A019565(n) in increasing order in row n. Square array A(n,k) n >= 0, k >= 0 read by descending antidiagonals.

Views

Author

Keywords

Comments

Examples

Crossrefs

Formula

A076512 Denominator of cototient(n)/totient(n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A064628 a(n) = floor((4/3)^n).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Sage

Extensions

A011754 Number of ones in the binary expansion of 3^n.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Formula

Extensions

A066003 Sum of digits of 7^n.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A027474 a(n) = 7^(n-2) * C(n,2).

Views

Author

A083098 a(n) = 2a(n-1) + 6a(n-2).