A006190 -id:A006190 - cp's OEIS frontend

A054456 Convolution triangle of A000129(n) (Pell numbers).

1, 2, 1, 5, 4, 1, 12, 14, 6, 1, 29, 44, 27, 8, 1, 70, 131, 104, 44, 10, 1, 169, 376, 366, 200, 65, 12, 1, 408, 1052, 1212, 810, 340, 90, 14, 1, 985, 2888, 3842, 3032, 1555, 532, 119, 16, 1, 2378, 7813, 11784, 10716, 6482, 2709, 784, 152, 18, 1, 5741, 20892, 35223

Offset: 0

Wolfdieter Lang, Apr 27 2000 and May 08 2000

In the language of the Shapiro et al. reference (given in A053121) such a lower triangular (ordinary) convolution array, considered as a matrix, belongs to the Bell-subgroup of the Riordan-group.
The G.f. for the row polynomials p(n,x) (increasing powers of x) is Pell(z)/(1-x*z*Pell(z)) with Pell(x)=1/(1-2*x-x^2) = g.f. for A000129(n+1) (Pell numbers without 0).
Column sequences are A000129(n+1), A006645(n+1), A054457(n) for m=0..2.
Riordan array (1/(1-2x-x^2),x/(1-2x-x^2)). - Paul Barry, Mar 15 2005
As a Riordan array, this factors as (1/(1-x^2),x/(1-x^2))*(1/(1-2x),x/(1-2x)), [abs(A049310) times square of A007318, or A038207]. - Paul Barry, Jul 28 2005
Coefficients of polynomials defined by P(x, 0) = 1; P(x, 1) = 2 - x; P(x, n) = (2 - x)*P(x, n - 1) + P(x, n - 2). - Roger L. Bagula, Mar 24 2008
Subtriangle (obtained by dropping the first column) of the triangle given by (0, 2, 1/2, -1/2, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Feb 19 2013
T(n,k) is the number of words of length n over {0,1,2,3} having k letters 3 and avoiding runs of odd length of the letter 0. - Milan Janjic, Jan 14 2017

			Fourth row polynomial (n=3): p(3,x)= 12+14*x+6*x^2+x^3
Triangle begins:
{1},
{2, 1},
{5, 4, 1},
{12, 14, 6, 1},
{29, 44, 27, 8, 1},
{70, 131,104, 44, 10, 1},
{169, 376, 366, 200, 65, 12, 1},
{408, 1052, 1212, 810, 340, 90, 14, 1},
{985, 2888, 3842, 3032, 1555, 532, 119, 16, 1},
{2378, 7813, 11784, 10716, 6482, 2709, 784, 152, 18, 1},
{5741, 20892, 35223, 36248, 25235, 12432, 4396, 1104, 189, 20, 1},
The triangle (0, 2, 1/2, -1/2, 0, 0, 0, ...) DELTA (1, 0, 0, 0, 0, 0, ...) begins:
1
0, 1
0, 2, 1
0, 5, 4, 1
0, 12, 14, 6, 1
0, 29, 44, 27, 8, 1 - _Philippe Deléham_, Feb 19 2013

Milan Janjić, Words and Linear Recurrences, J. Int. Seq. 21 (2018), #18.1.4.

Cf. A000129. Row sums: A006190(n+1).

Cf. A129844.

G := 1/(1-(x+2)*z-z^2): Gser := simplify(series(G, z = 0, 18)): for n from 0 to 15 do P[n] := sort(coeff(Gser, z, n)) end do: for n from 0 to 15 do seq(coeff(P[n], x, j), j = 0 .. n) end do; # yields sequence in triangular form - Emeric Deutsch, Aug 30 2015
T := (n,k) -> `if`(n=0,1,2^(n-k)*binomial(n,k)*hypergeom([(k-n)/2,(k-n+1)/2],[-n], -1)): seq(seq(simplify(T(n,k)),k=0..n),n=0..10); # Peter Luschny, Apr 25 2016
# Uses function PMatrix from A357368. Adds column 1,0,0,0,... to the left.
PMatrix(10, A000129); # Peter Luschny, Oct 19 2022

P[x_, 0] := 1; P[x_, 1] := 2 - x; P[x_, n_] := P[x, n] = (2 - x) P[x, n - 1] + P[x, n - 2]; Table[Abs@ CoefficientList[P[x, n], x], {n, 0, 10}] // Flatten (* Roger L. Bagula, Mar 24 2008, edited by Michael De Vlieger, Apr 25 2018 *)

a(n, m) := ((n-m+1)*a(n, m-1) + (n+m)*a(n-1, m-1))/(4*m), n >= m >= 1, a(n, 0)= P(n+1)= A000129(n+1) (Pell numbers without P(0)), a(n, m) := 0 if n

G.f. for column m: Pell(x)*(x*Pell(x))^m, m >= 0, with Pell(x) G.f. for A000129(n+1).

Number triangle T(n, k) with T(n, 0)=A000129(n), T(1, 1)=1, T(n, k)=0 if k>n, T(n, k)=T(n-1, k-1)+T(n-2, k)+2T(n-1, k) otherwise; T(n, k)=if(k<=n, sum{j=0..floor((n-k)/2), C(n-j, k)C(n-k-j, j)2^(n-2j-k)}; - Paul Barry, Mar 15 2005

Bivariate g.f. G(x,z) = 1/[1 - (2 + x)z - z^2]. G.f. for column k = z^k/(1 - 2z - z^2)^{k+1} (k>=0). - Emeric Deutsch, Aug 30 2015

T(n,k) = 2^(n-k)*C(n,k)*hypergeom([(k-n)/2,(k-n+1)/2],[-n],-1)) for n>=1. - Peter Luschny, Apr 25 2016

A052924 Expansion of g.f.: (1-x)/(1 - 3*x - x^2).

Original entry on oeis.org

1, 2, 7, 23, 76, 251, 829, 2738, 9043, 29867, 98644, 325799, 1076041, 3553922, 11737807, 38767343, 128039836, 422886851, 1396700389, 4612988018, 15235664443, 50319981347, 166195608484, 548906806799, 1812916028881

Offset: 0

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

Euler encountered this sequence when finding the largest root of z^2 - 3z - 1 = 0. - V. Frederick Rickey (fred-rickey(AT)usma.edu), Aug 20 2003
Let M = a triangle with the Pell series A000129 (1, 2, 5, 12, ...) in each column, with the leftmost column shifted upwards one row. A052924 starting (1, 2, 7, 23, ...) = lim_{n->infinity} M^n, the left-shifted vector considered as a sequence. - Gary W. Adamson, Jul 31 2010
a(n) is the number of compositions of n when there are 2 types of 1 and 3 types of other natural numbers. - Milan Janjic, Aug 13 2010
Equals partial sums of A108300 prefaced with a 1: (1, 1, 5, 16, 53, 175, 578, ...). - Gary W. Adamson, Feb 15 2012

L. Euler, Introductio in analysin infinitorum, 1748, section 338. English translation by John D. Blanton, Introduction to Analysis of the Infinite, 1988, Springer, p. 286.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Sergio Falcón, The k-Fibonacci difference sequences, Chaos, Solitons & Fractals, Volume 87, June 2016, Pages 153-157.
Sergio Falcón and Ángel Plaza, On the Fibonacci k-numbers, Chaos, Solitons & Fractals 2007; 32(5): 1615-24.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 909
Tanya Khovanova, Recursive Sequences
Index entries for linear recurrences with constant coefficients, signature (3,1).

A108300 (first differences), A006190 (partial sums), A355981 (primes).

Cf. A000032, A001077, A078343, A164581, A179237, A180148, A329723.

a:=[1,2];; for n in [3..30] do a[n]:=3*a[n-1]+a[n-2]; od; a; # G. C. Greubel, Jun 09 2019

R:=PowerSeriesRing(Integers(), 30); Coefficients(R!( (1-x)/(1-3*x-x^2) )); // G. C. Greubel, Jun 09 2019

spec:= [S,{S=Sequence(Prod(Sequence(Z),Union(Z,Z,Prod(Z,Z))))}, unlabeled]: seq(combstruct[count](spec,size=n), n=0..30);
seq(coeff(series((1-x)/(1-3*x-x^2), x, n+1), x, n), n = 0..30); # G. C. Greubel, Oct 16 2019

CoefficientList[Series[(1-x)/(1-3*x-x^2), {x,0,30}], x] (* G. C. Greubel, Jun 09 2019 *)

Vec((1-x)/(1-3*x-x^2)+O(x^30)) \\ Charles R Greathouse IV, Nov 20 2011

((1-x)/(1-3*x-x^2)).series(x, 30).coefficients(x, sparse=False) # G. C. Greubel, Jun 09 2019

a(n) = 3*a(n-1) + a(n-2).
a(n) = Sum_{alpha=RootOf(-1+3*x+x^2)} (1/13)*(1+5*alpha)*alpha^(-1-n).
With offset 1: a(1)=1; for n > 1, a(n) = Sum_{i=1..3*n-4} a(ceiling(i/3)). - Benoit Cloitre, Jan 04 2004
Binomial transform of A006130. a(n) = (1/2 - sqrt(13)/26)*(3/2 - sqrt(13)/2)^n + (1/2 + sqrt(13)/26)*(3/2 + sqrt(13)/2)^n. - Paul Barry, Jul 20 2004
From Creighton Dement, Nov 04 2004: (Start)
a(n) = A006190(n+1) - A006190(n);
4*a(n) = 9*A006190(n+1) - A006497(n+1) - 2*A003688(n+1). (End)
Numerators in continued fraction [1, 2, 3, 3, 3, ...], where the latter = 0.69722436226...; the length of an inradius of a right triangle with legs 2 and 3. - Gary W. Adamson, Dec 19 2007
If p[1]=2, p[i]=3, (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-1) = det A. - Milan Janjic, Apr 29 2010
a(n) = A006190(n) + A003688(n). - R. J. Mathar, Jul 06 2012
a(n) = Sum_{k=0..n-2} A168561(n-2,k)*3^k + 2 * Sum_{k=0..n-1} A168561(n-1,k)*3^k, n>0. - R. J. Mathar, Feb 14 2024
From Peter Bala, Jul 08 2025: (Start)
The following series telescope:
Sum_{n >= 1} 1/(a(n) + 3*(-1)^(n+1)/a(n)) = 1/2, since 1/(a(n) + 3*(-1)^(n+1)/a(n)) = b(n) - b(n+1), where b(n) = (1/3) * (a(n) + a(n-1)) / (a(n)*a(n-1)).
Sum_{n >= 1} (-1)^(n+1)/(a(n) + 3*(-1)^(n+1)/a(n)) = 1/6, since 1/(a(n) + 3*(-1)^(n+1)/a(n)) = c(n) + c(n+1), where c(n) = (1/3) * (a(n) - a(n-1)) / (a(n)*a(n-1)). (End)

More terms from James Sellers, Jun 06 2000

A367300 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 3 + 2x, p(n,x) = up(n-1,x) + vp(n-2,x) for n >= 3, where u = p(2,x), v = 1 - 2x - x^2.

Original entry on oeis.org

1, 3, 2, 10, 10, 3, 33, 46, 22, 4, 109, 194, 131, 40, 5, 360, 780, 678, 296, 65, 6, 1189, 3036, 3228, 1828, 581, 98, 7, 3927, 11546, 14514, 10100, 4194, 1036, 140, 8, 12970, 43150, 62601, 51664, 26479, 8604, 1722, 192, 9, 42837, 159082, 261598, 249720, 152245, 61318, 16248, 2712, 255, 10

Offset: 1

Clark Kimberling, Dec 23 2023

Because (p(n,x)) is a strong divisibility sequence, for each integer k, the sequence (p(n,k)) is a strong divisibility sequence of integers.

			First eight rows:
     1
     3      2
    10     10      3
    33     46     22      4
   109    194    131     40     5
   360    780    678    296    65     6
  1189   3036   3228   1828   581    98    7
  3927  11546  14514  10100  4194  1036  140  8
Row 4 represents the polynomial p(4,x) = 33 + 46*x + 22*x^2 + 4*x^3, so (T(4,k)) = (33,46,22,4), k=0..3.

Rigoberto Flórez, Robinson Higuita, and Antara Mukherjee, Characterization of the strong divisibility property for generalized Fibonacci polynomials, Integers, 18 (2018), Paper No. A14.

Cf. A006190 (column 1); A000027 (p(n,n-1)); A107839 (row sums, p(n,1)); A001045 (alternating row sums, p(n,-1)); A030240 (p(n,2)); A039834 (signed Fibonacci numbers, p(n,-2)); A016130 (p(n,3)); A225883 (p(n,-3)); A099450 (p(n,-4)); A094440, A367208, A367209, A367210, A367211, A367297, A367298, A367299.

p[1, x_] := 1; p[2, x_] := 3 + 2 x; u[x_] := p[2, x]; v[x_] := 1 - 2 x - x^2;
p[n_, x_] := Expand[u[x]*p[n - 1, x] + v[x]*p[n - 2, x]]
Grid[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]
Flatten[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]

p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where p(1,x) = 1, p(2,x) = 3 + 2*x, u = p(2,x), and v = 1 - 2*x - x^2.

p(n,x) = k*(b^n - c^n), where k = -(1/sqrt(13 + 4*x)), b = (1/2) (2*x + 3 + 1/k), c = (1/2) (2*x + 3 - 1/k).

A006138 a(n) = a(n-1) + 3*a(n-2).

Original entry on oeis.org

1, 2, 5, 11, 26, 59, 137, 314, 725, 1667, 3842, 8843, 20369, 46898, 108005, 248699, 572714, 1318811, 3036953, 6993386, 16104245, 37084403, 85397138, 196650347, 452841761, 1042792802, 2401318085, 5529696491, 12733650746, 29322740219, 67523692457, 155491913114

Offset: 0

N. J. A. Sloane

The binomial transform of a(n) is b(n) = A006190(n+1), which satisfies b(n) = 3*b(n-1) + b(n-2). - Paul Barry, May 21 2006
Partial sums of A105476. - Paul Barry, Feb 02 2007
An elephant sequence, see A175654. For the corner squares four A[5] vectors, with decimal values 69, 261, 321 and 324, lead to this sequence. For the central square these vectors lead to the companion sequence A105476 (without the first leading 1). - Johannes W. Meijer, Aug 15 2010
Equals the INVERTi transform of A063782: (1, 3, 10, 32, 104, ...). - Gary W. Adamson, Aug 14 2010
Pisano period lengths: 1, 3, 1, 6, 24, 3, 24, 6, 3, 24, 120, 6, 156, 24, 24, 12, 16, 3, 90, 24, ... - R. J. Mathar, Aug 10 2012
The sequence is the INVERT transform of A016116: (1, 1, 2, 2, 4, 4, 8, 8, ...). - Gary W. Adamson, Jul 17 2015

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

Vincenzo Librandi, Table of n, a(n) for n = 0..200
N. T. Gridgeman, A new look at Fibonacci generalization, Fib,. Quart., 11 (1973), 40-55.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Index entries for linear recurrences with constant coefficients, signature (1,3).

Cf. A063782. - Gary W. Adamson, Aug 14 2010

Cf. A006130, A006190, A016116, A105476, A122950, A175654, A274977.

a:=[1,2];; for n in [3..40] do a[n]:=a[n-1]+3*a[n-2]; od; a; # G. C. Greubel, Nov 19 2019

[n le 2 select n else Self(n-1)+3*Self(n-2): n in [1..40]]; // Vincenzo Librandi, Sep 15 2016

A006138:=-(1+z)/(-1+z+3*z**2); # Simon Plouffe in his 1992 dissertation

CoefficientList[Series[(1+z)/(1-z-3*z^2), {z,0,40}], z] (* Vladimir Joseph Stephan Orlovsky, Jun 11 2011 *)
Table[(I*Sqrt[3])^(n-1)*(I*Sqrt[3]*ChebyshevU[n, 1/(2*I*Sqrt[3])] + ChebyshevU[n-1, 1/(2*I*Sqrt[3])]), {n, 0, 40}]//Simplify (* G. C. Greubel, Nov 19 2019 *)
LinearRecurrence[{1,3},{1,2},40] (* Harvey P. Dale, May 29 2025 *)

main(size)={my(v=vector(size),i);v[1]=1;v[2]=2;for(i=3,size,v[i]=v[i-1]+3*v[i-2]);return(v);} /* Anders Hellström, Jul 17 2015 */

def A006138_list(prec):
    P. = PowerSeriesRing(ZZ, prec)
    return P((1+x)/(1-x-3*x^2)).list()
A006138_list(40) # G. C. Greubel, Nov 19 2019

a(n) = Sum_{k=0..n+1} A122950(n+1,k)*2^(n+1-k). - Philippe Deléham, Jan 04 2008
G.f.: (1+x)/(1-x-3*x^2). - Paul Barry, May 21 2006
a(n) = Sum_{k=0..n} C(floor((2n-k)/2),n-k)*3^floor(k/2). - Paul Barry, Feb 02 2007
a(n) = A006130(n) + A006130(n-1). - R. J. Mathar, Jun 22 2011
a(n) = (i*sqrt(3))^(n-1)*(i*sqrt(3)*ChebyshevU(n, 1/(2*i*sqrt(3))) + ChebyshevU(n-1, 1/(2*i*sqrt(3)))), where i=sqrt(-1). - G. C. Greubel, Nov 19 2019

Typo in formula corrected by Johannes W. Meijer, Aug 15 2010

A057089 Scaled Chebyshev U-polynomials evaluated at i*sqrt(6)/2. Generalized Fibonacci sequence.

Original entry on oeis.org

1, 6, 42, 288, 1980, 13608, 93528, 642816, 4418064, 30365280, 208700064, 1434392064, 9858552768, 67757668992, 465697330560, 3200729997312, 21998563967232, 151195763787264, 1039165966526976, 7142170381885440

Offset: 0

Wolfdieter Lang, Aug 11 2000

a(n) gives the length of the word obtained after n steps with the substitution rule 0->1^6, 1->(1^6)0, starting from 0. The number of 1's and 0's of this word is 6*a(n-1) and 6*a(n-2), resp.

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Martin Burtscher, Igor Szczyrba, and Rafał Szczyrba, Analytic Representations of the n-anacci Constants and Generalizations Thereof, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.5.
A. F. Horadam, Special properties of the sequence W_n(a,b; p,q), Fib. Quart., 5.5 (1967), 424-434. Case n->n+1, a=0,b=1; p=6, q=6.
Tanya Khovanova, Recursive Sequences
Wolfdieter Lang, On polynomials related to powers of the generating function of Catalan's numbers, Fib. Quart. 38 (2000) 408-419. Eqs.(39) and (45),rhs, m=6.
Index entries for sequences related to Chebyshev polynomials.
Index entries for linear recurrences with constant coefficients, signature (6,6).

Cf. A001076, A006190, A007482, A015520, A015521, A015523, A015524, A015525, A015528, A015529, A015530, A015531, A015532, A015533, A015534, A015535, A015536, A015537, A015440, A015441, A015443, A015444, A015445, A015447, A015548, A030195, A053404, A057087, A057088, A083858, A085939, A090017, A091914, A099012, A135030, A135032, A180222, A180226, A180250.

I:=[1,6]; [n le 2 select I[n] else 6*Self(n-1)+6*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 14 2011

Join[{a=0,b=1},Table[c=6*b+6*a;a=b;b=c,{n,100}]] (* Vladimir Joseph Stephan Orlovsky, Jan 16 2011 *)
LinearRecurrence[{6,6},{1,6},40] (* Harvey P. Dale, Nov 05 2011 *)

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

[lucas_number1(n,6,-6) for n in range(1, 21)] # Zerinvary Lajos, Apr 24 2009

a(n) = 6*a(n-1) + 6*a(n-2); a(0)=1, a(1)=6.
a(n) = S(n, i*sqrt(6))*(-i*sqrt(6))^n with S(n, x) := U(n, x/2), Chebyshev's polynomials of the 2nd kind, A049310.
G.f.: 1/(1-6*x-6*x^2).
a(n) = Sum_{k=0..n} 5^k*A063967(n,k). - Philippe Deléham, Nov 03 2006

A084844 Denominators of the continued fraction n + 1/(n + 1/...) [n times].

Original entry on oeis.org

1, 2, 10, 72, 701, 8658, 129949, 2298912, 46866034, 1082120050, 27916772489, 795910114440, 24851643870041, 843458630403298, 30918112619119426, 1217359297034666112, 51240457936070359069, 2296067756927144738850, 109127748348241605689981

Offset: 1

Hollie L. Buchanan II, Jun 08 2003

The (n-1)-th term of the Lucas sequence U(n,-1). The numerator is the n-th term. Adjacent terms of the sequence U(n,-1) are relatively prime. - T. D. Noe, Aug 19 2004
From Flávio V. Fernandes, Mar 05 2021: (Start)
Also, the n-th term of the n-th metallic sequence (the diagonal through the array A073133, and its equivalents, which is rows formed by sequences beginning with A000045, A000129, A006190, A001076, A052918) as shown below (for n>=1):
   0    1    0    1     0      1  ... A000035
   0   [1]   1    2     3      5  ... A000045
   0    1   [2]   5    12     29  ... A000129
   0    1    3  [10]   33    109  ... A006190
   0    1    4   17   [72]   305  ... A001076
   0    1    5   26   135   [701] ... A052918. (End)

			a(4) = 72 since 4 + 1/(4 + 1/(4 + 1/4)) = 305/72.

Alois P. Heinz, Table of n, a(n) for n = 1..387
Eric Weisstein's World of Mathematics, Lucas Sequence

Cf. A084845 (numerators).

Cf. A000045, A097690, A097691, A117715, A290864 (primes in this sequence).

A084844 :=proc(n) combinat[fibonacci](n, n) end:
seq(A084844(n), n=1..30); # Zerinvary Lajos, Jan 03 2007

myList[n_] := Module[{ex = {n}}, Do[ex = {ex, n}, {n - 1}]; Flatten[ex]] Table[Denominator[FromContinuedFraction[myList[n]]], {n, 1, 20}]
Table[s=n; Do[s=n+1/s, {n-1}]; Denominator[s], {n, 20}] (* T. D. Noe, Aug 19 2004 *)
Table[Fibonacci[n, n], {n, 1, 20}] (* Vladimir Reshetnikov, May 07 2016 *)
Table[DifferenceRoot[Function[{y,m},{y[2+m]==n*y[1+m]+y[m],y[0]==0,y[1]==1}]][n],{n,1,20}] (* Benedict W. J. Irwin, Nov 03 2016 *)

from sympy import fibonacci
def a(n):
    return fibonacci(n, n)
print([a(n) for n in range(1, 31)]) # Indranil Ghosh, Aug 12 2017

a(n) = (s^n - (-s)^(-n))/(2*s - n), where s = (n + sqrt(n^2 + 4))/2. - Vladimir Reshetnikov, May 07 2016
a(n) = y(n,n), where y(m+2,n) = n*y(m+1,n) + y(m,n), with y(0,n)=0, y(1,n)=1 for all n. - Benedict W. J. Irwin, Nov 03 2016
a(n) ~ n^(n-1). - Vaclav Kotesovec, Jun 03 2017
a(n) = A117715(n,n). - Bobby Jacobs, Aug 12 2017
a(n) = [x^n] x/(1 - n*x - x^2). - Ilya Gutkovskiy, Oct 10 2017
a(n) == 0 (mod n) for even n and 1 (mod n) for odd n. - Flávio V. Fernandes, Dec 08 2020
a(n) == 0 (mod n) for even n and 1 (mod n^2) for odd n; see A065599. - Flávio V. Fernandes, Dec 25 2020
a(n) == 0 (mod 2*(n/2)^2) for even n and 1 (mod n^2) for odd n; see A129194. - Flávio V. Fernandes, Feb 06 2021

A090880 Suppose n=(p1^e1)(p2^e2)... where p1,p2,... are the prime numbers and e1,e2,... are nonnegative integers. Then a(n) = e1 + (e2)3 + (e3)9 + (e4)27 + ... + (ek)(3^(k-1)) + ...

Original entry on oeis.org

0, 1, 3, 2, 9, 4, 27, 3, 6, 10, 81, 5, 243, 28, 12, 4, 729, 7, 2187, 11, 30, 82, 6561, 6, 18, 244, 9, 29, 19683, 13, 59049, 5, 84, 730, 36, 8, 177147, 2188, 246, 12, 531441, 31, 1594323, 83, 15, 6562, 4782969, 7, 54, 19, 732, 245, 14348907, 10, 90, 30, 2190

Offset: 1

Sam Alexander, Dec 12 2003

Replace "3" with "x" and extend the definition of a to positive rationals and a becomes an isomorphism between positive rationals under multiplication and polynomials over Z under addition. This remark generalizes A001222, A048675 and A054841: evaluate said polynomial at x=1, x=2 and x=10, respectively.

For examples of such evaluations at x=3, see "Other identities" in the Formula section. - Antti Karttunen, Jul 31 2015

Joseph J. Rotman, The Theory of Groups: An Introduction, 2nd ed. Boston: Allyn and Bacon, Inc. 1973. Page 9, problem 1.26.

Antti Karttunen, Table of n, a(n) for n = 1..512
Sam Alexander, Post to sci.math.

Row 3 of A104244.

Cf. A001222, A006190, A048675, A054841, A090881, A090882, A090883, A090884, A178590, A206296, A260443.

(define (A090880 n) (if (= 1 n) (- n 1) (+ (A000244 (- (A055396 n) 1)) (A090880 (A032742 n))))) ;; Antti Karttunen, Jul 29 2015

a(1) = 0; for n > 1, a(n) = 3^(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.] - Antti Karttunen, Jul 29 2015
Other identities. For all n >= 0:
a(A206296(n)) = A006190(n).
a(A260443(n)) = A178590(n).

More terms from Ray Chandler, Dec 20 2003

Previous Showing 61-70 of 146 results. Next

cp's OEIS Frontend

A309584 Numbers k with 2 zeros in a fundamental period of A000129 mod k.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A309585 Numbers k with 4 zeros in a fundamental period of A000129 mod k.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A367209 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 1 + 4*x, p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where u = p(2,x), v = 1 - x - x^2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A054456 Convolution triangle of A000129(n) (Pell numbers).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A052924 Expansion of g.f.: (1-x)/(1 - 3*x - x^2).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

Extensions

A367300 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 3 + 2*x, p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where u = p(2,x), v = 1 - 2*x - x^2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A006138 a(n) = a(n-1) + 3*a(n-2).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

A367209 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 1 + 4x, p(n,x) = up(n-1,x) + v*p(n-2,x) for n >= 3, where u = p(2,x), v = 1 - x - x^2.

A367300 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 3 + 2x, p(n,x) = up(n-1,x) + vp(n-2,x) for n >= 3, where u = p(2,x), v = 1 - 2x - x^2.

A090880 Suppose n=(p1^e1)(p2^e2)... where p1,p2,... are the prime numbers and e1,e2,... are nonnegative integers. Then a(n) = e1 + (e2)3 + (e3)9 + (e4)27 + ... + (ek)(3^(k-1)) + ...