A015552 -id:A015552 - cp's OEIS frontend

A206369 a(p^k) = p^k - p^(k-1) + p^(k-2) - ... +- 1, and then extend by multiplicativity.

1, 1, 2, 3, 4, 2, 6, 5, 7, 4, 10, 6, 12, 6, 8, 11, 16, 7, 18, 12, 12, 10, 22, 10, 21, 12, 20, 18, 28, 8, 30, 21, 20, 16, 24, 21, 36, 18, 24, 20, 40, 12, 42, 30, 28, 22, 46, 22, 43, 21, 32, 36, 52, 20, 40, 30, 36, 28, 58, 24, 60, 30, 42, 43, 48, 20, 66, 48, 44, 24, 70, 35

Offset: 1

N. J. A. Sloane, Feb 06 2012

For more information see the Comments in A061020.
a(n) is the number of integers j such that 1 <= j <= n and gcd(n,j) is a perfect square. For example, a(12) = 6 because |{1,4,5,7,8,11}|=6 and the respective GCDs with 12 are 1,4,1,1,4,1, which are squares. - Geoffrey Critzer, Feb 16 2015
If m is squarefree (A005117), then a(m) = A000010(m) where A000010 is the Euler totient function. - Michel Marcus, Nov 08 2017
Also it appears that the primorials (A002110) is the sequence of indices of minimum records for a(n)/n, and these records are A038110(n)/A060753(n). - Michel Marcus, Nov 09 2017
Also called rho(n). When rho(n) | n, then n is called k-imperfect, with k = n/rho(n), cf. A127724. - M. F. Hasler, Feb 13 2020

P. J. McCarthy, Introduction to Arithmetical Functions, Springer Verlag, 1986, page 25.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
László Tóth, A survey of the alternating sum-of-divisors function, arXiv:1111.4842 [math.NT], 2011-2014.

Cf. A061020, A206368.
Cf. A027748 row, A124010, A206475 (first differences).
Cf. A078429.
Cf. A127724 (k-imperfect), A127725 (2-imperfect), A127726 (3-imperfect).

a206369 n = product $
   zipWith h (a027748_row n) (map toInteger $ a124010_row n) where
           h p e = sum $ take (fromInteger e + 1) $
                         iterate ((* p) . negate) (1 - 2 * (e `mod` 2))
-- Reinhard Zumkeller, Feb 08 2012

a:= n-> mul(add(i[1]^(i[2]-j)*(-1)^j, j=0..i[2]), i=ifactors(n)[2]):
seq(a(n), n=1..100);  # Alois P. Heinz, Nov 03 2017

Table[Length[Select[Range[n], IntegerQ[GCD[n, #]^(1/2)] &]], {n, 72}] (* Geoffrey Critzer, Feb 16 2015 *)
a[n_] := n*DivisorSum[n, LiouvilleLambda[#]/#&]; Array[a, 72] (* Jean-François Alcover, Dec 04 2017, after Enrique Pérez Herrero *)
f[p_,e_] := Sum[(-1)^(e-k)*p^k, {k,0,e}]; a[1] = 1; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Jan 01 2020 *)

a(n) = sum(k=1, n, issquare(gcd(n, k)));

ak(p,e)=my(s=1); for(i=1,e, s=s*p + (-1)^i); s
a(n)=my(f=factor(n)); prod(i=1,#f~, ak(f[i,1],f[i,2])) \\ Charles R Greathouse IV, Dec 27 2016

a(n) = sumdiv(n, d, eulerphi(n/d) * issquare(d)); \\ Daniel Suteu, Jun 27 2018

apply( {A206369(n)=vecprod([f[1]^(f[2]+1)\/(f[1]+1)|f<-factor(n)~])}, [1..99]) \\ M. F. Hasler, Feb 13 2020

from math import prod
from sympy import factorint
def A206369(n): return prod((lambda x:x[0]+int((x[1]<<1)>=p+1))(divmod(p**(e+1),p+1)) for p, e in factorint(n).items()) # Chai Wah Wu, Mar 05 2024

a(n) = abs(A061020(n)).
a(n) = n*Sum_{d|n} lambda(d)/d, where lambda(n) is A008836(n). - Enrique Pérez Herrero, Sep 23 2012
Dirichlet g.f.: zeta(s - 1)*zeta(2*s)/zeta(s). - Geoffrey Critzer, Feb 25 2015
From Michel Marcus, Nov 05 2017: (Start)
a(2^n) = A001045(n+1);
a(3^n) = A015518(n+1);
a(5^n) = A015531(n+1);
a(7^n) = A015552(n+1);
a(11^n) = A015592(n+1). (End)
a(p^k) = p^k - a(p^(k - 1)) for k > 0 and prime p. - David A. Corneth, Nov 09 2017
a(n) = Sum_{d|n, d is a perfect square} phi(n/d), where phi(k) is the Euler totient function. - Daniel Suteu, Jun 27 2018
a(p^k) = A071324(p^k), for k >= 0 and prime p. - Michel Marcus, Aug 11 2018
Sum_{k=1..n} a(k) ~ Pi^2 * n^2 / 30. - Vaclav Kotesovec, Feb 07 2019
G.f.: Sum_{k>=1} lambda(k)*x^k/(1 - x^k)^2. - Ilya Gutkovskiy, May 23 2019
a(n) = Sum_{i=1..n} A010052(gcd(n,i)). - Ridouane Oudra, Nov 24 2019
a(p^k) = round(p^(k+1)/(p+1)). - M. F. Hasler, Feb 13 2020

A015565 a(n) = 7a(n-1) + 8a(n-2), a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 7, 57, 455, 3641, 29127, 233017, 1864135, 14913081, 119304647, 954437177, 7635497415, 61083979321, 488671834567, 3909374676537, 31274997412295, 250199979298361, 2001599834386887, 16012798675095097, 128102389400760775, 1024819115206086201, 8198552921648689607

Offset: 0

Olivier Gérard

A linear 2nd order recurrence. A Jacobsthal number sequence.
Binomial transform of A053573 (preceded by zero). - Paul Barry, Apr 09 2003
Second binomial transform of A080424. Binomial transform of A053573, with leading zero. Binomial transform is 0,1,9,81,729,....(9^n - 0^n)/9. Second binomial transform is 0,1,11,111,1111,... (A002275: repunits). - Paul Barry, Mar 14 2004
Number of walks of length n between any two distinct nodes of the complete graph K_9. Example: a(2)=7 because the walks of length 2 between the nodes A and B of the complete graph ABCDEFGHI are: ACB, ADB, AEB, AFB, AGB, AHB and AIB. - Emeric Deutsch, Apr 01 2004
Unsigned version of A014990. - Philippe Deléham, Feb 13 2007
General form: k=8^n-k. Also: A001045, A078008, A097073, A115341, A015518, A054878, A015521, A109499, A015531, A109500, A109501, A015552. - Vladimir Joseph Stephan Orlovsky, Dec 11 2008
The ratio a(n+1)/a(n) converges to 8 as n approaches infinity. - Felix P. Muga II, Mar 09 2014

			G.f. = x + 7*x^2 + 57*x^3 + 455*x^4 + 3641*x^5 + 29127*x^6 + 233017*x^7 + ...

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Jean-Paul Allouche, Jeffrey Shallit, Zhixiong Wen, Wen Wu, and Jiemeng Zhang, Sum-free sets generated by the period-k-folding sequences and some Sturmian sequences, arXiv:1911.01687 [math.CO], 2019.
M. Dukes and C. D. White, Web Matrices: Structural Properties and Generating Combinatorial Identities, arXiv:1603.01589 [math.CO], 2016.
Dale Gerdemann, Fractal generated from (7,8) recursion, YouTube Video, Dec 5, 2014.
Index entries for linear recurrences with constant coefficients, signature (7,8).

Cf. A002275, A014990, A053573, A080424, A082311, A082365.
Cf. A062160, A099322, A106566.
Cf. A001045, A078008, A097073, A115341, A015518, A054878, A015521, A109499, A015531, A109500, A109501, A015552. - Vladimir Joseph Stephan Orlovsky, Dec 11 2008

[Round(8^n/9): n in [0..30]]; // Vincenzo Librandi, Jun 24 2011

seq(round(8^n/9),n=0..25); # Mircea Merca, Dec 28 2010

k=0;lst={k};Do[k=8^n-k;AppendTo[lst, k], {n, 0, 5!}];lst (* Vladimir Joseph Stephan Orlovsky, Dec 11 2008 *)
LinearRecurrence[{7,8},{0,1},30] (* Harvey P. Dale, Mar 04 2016 *)

x='x+O('x^30); concat([0], Vec(x/(1-7*x-8*x^2))) \\ G. C. Greubel, Dec 30 2017

[lucas_number1(n,7,-8) for n in range(0, 21)] # Zerinvary Lajos, Apr 24 2009

From Paul Barry, Apr 09 2003: (Start)
a(n) = (8^n - (-1)^n)/9.
a(n) = J(3*n)/3 = A001045(3*n)/3. (End)
From Emeric Deutsch, Apr 01 2004: (Start)
a(n) = 8^(n-1) - a(n-1).
G.f.: x/(1-7*x-8*x^2). (End)
a(n) = Sum_{k = 0..n} A106566(n,k)*A099322(k). - Philippe Deléham, Oct 30 2008
a(n) = round(8^n/9). - Mircea Merca, Dec 28 2010
From Peter Bala, May 31 2024: (Start)
G.f: A(x) = x/(1 - x^2) o x/(1 - x^2), where o denotes the black diamond product of power series as defined by Dukes and White. Cf. A054878.
The black diamond product A(x) o A(x) is the g.f. for the number of walks of length n between any two distinct nodes of the complete graph K_81.
Row 8 of A062160. (End)
E.g.f.: exp(-x)*(exp(9*x) - 1)/9. - Elmo R. Oliveira, Aug 17 2024

A015585 a(n) = 9a(n-1) + 10a(n-2).

Original entry on oeis.org

0, 1, 9, 91, 909, 9091, 90909, 909091, 9090909, 90909091, 909090909, 9090909091, 90909090909, 909090909091, 9090909090909, 90909090909091, 909090909090909, 9090909090909091, 90909090909090909, 909090909090909091, 9090909090909090909, 90909090909090909091

Offset: 0

Olivier Gérard

Number of walks of length n between any two distinct nodes of the complete graph K_11. Example: a(2)=9 because the walks of length 2 between the nodes A and B of the complete graph ABCDEFGHIJK are: ACB, ADB, AEB, AFB, AGB, AHB, AIB, AJB and AKB. - Emeric Deutsch, Apr 01 2004
Beginning with n=1 and a(1)=1, these are the positive integers whose balanced base-10 representations (A097150) are the first n digits of 1,-1,1,-1,.... Also, a(n) = (-1)^(n-1)*A014992(n) = |A014992(n)| for n >= 1. - Rick L. Shepherd, Jul 30 2004
General form: k=10^n-k. Also: A001045, A078008, A097073, A115341, A015518, A054878, A015521, A109499, A015531, A109500, A109501, A015552, A093134, A015565, A015577. - Vladimir Joseph Stephan Orlovsky, Dec 11 2008

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Jean-Paul Allouche, Jeffrey Shallit, Zhixiong Wen, Wen Wu, Jiemeng Zhang, Sum-free sets generated by the period-k-folding sequences and some Sturmian sequences, arXiv:1911.01687 [math.CO], 2019.
Index entries for linear recurrences with constant coefficients, signature (9,10).

Cf. A014992 (q-integers for q=-10), A097150.
Cf. A001045, A078008, A097073, A115341, A015518, A054878, A015521, A109499, A015531, A109500, A109501, A015552, A093134, A015565, A015577. - Vladimir Joseph Stephan Orlovsky, Dec 11 2008
Cf. A094028, A098610, A098611.

[Round(10^n/11): n in [0..30]]; // Vincenzo Librandi, Jun 24 2011

k=0;lst={k};Do[k=10^n-k;AppendTo[lst, k], {n, 0, 5!}];lst (* Vladimir Joseph Stephan Orlovsky, Dec 11 2008 *)
LinearRecurrence[{9,10},{0,1},30] (* Harvey P. Dale, Aug 08 2014 *)

a(n)=10^n\/11 \\ Charles R Greathouse IV, Jun 24 2011

[lucas_number1(n,9,-10) for n in range(0, 19)] # Zerinvary Lajos, Apr 26 2009

[abs(gaussian_binomial(n,1,-10)) for n in range(0,19)] # Zerinvary Lajos, May 28 2009

a(n) = 9*a(n-1) + 10*a(n-2).
From Emeric Deutsch, Apr 01 2004: (Start)
a(n) = 10^(n-1) - a(n-1).
G.f.: x/(1 - 9x - 10x^2). (End)
From Henry Bottomley, Sep 17 2004: (Start)
a(n) = round(10^n/11).
a(n) = (10^n - (-1)^n)/11.
a(n) = A098611(n)/11 = 9*A094028(n+1)/A098610(n). (End)
E.g.f.: exp(-x)*(exp(11*x) - 1)/11. - Elmo R. Oliveira, Aug 17 2024

Extended by T. D. Noe, May 23 2011

A090018 a(n) = 6a(n-1) + 3a(n-2) for n > 2, a(0)=1, a(1)=6.

Original entry on oeis.org

1, 6, 39, 252, 1629, 10530, 68067, 439992, 2844153, 18384894, 118841823, 768205620, 4965759189, 32099171994, 207492309531, 1341251373168, 8669985167601, 56043665125110, 362271946253463, 2341762672896108, 15137391876137037, 97849639275510546, 632510011281474387

Offset: 0

Paul Barry, Nov 19 2003

From Johannes W. Meijer, Aug 09 2010: (Start)

a(n) represents the number of n-move routes of a fairy chess piece starting in a given corner or side square on a 3 X 3 chessboard. This fairy chess piece behaves like a white queen on the eight side and corner squares but on the central square the queen explodes with fury and turns into a red queen, see A180032. The central square leads to A180028. (End)

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Index entries for sequences related to Chebyshev polynomials.
Index entries for linear recurrences with constant coefficients, signature (6,3).

Cf. A141041, A180028, A180032.

Sequences with g.f. of the form 1/(1 - 6*x - k*x^2): A106392 (k=-10), A027471 (k=-9), A006516 (k=-8), A081179 (k=-7), A030192 (k=-6), A003463 (k=-5), A084326 (k=-4), A138395 (k=-3), A154244 (k=-2), A001109 (k=-1), A000400 (k=0), A005668 (k=1), A135030 (k=2), this sequence (k=3), A135032 (k=4), A015551 (k=5), A057089 (k=6), A015552 (k=7), A189800 (k=8), A189801 (k=9), A190005 (k=10), A015553 (k=11).

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

a:= n-> (<<0|1>, <3|6>>^n. <<1,6>>)[1,1]:
seq(a(n), n=0..30);  # Alois P. Heinz, Jan 17 2011

Join[{a=1,b=6},Table[c=6*b+3*a;a=b;b=c,{n,100}]] (* Vladimir Joseph Stephan Orlovsky, Jan 16 2011 *)
LinearRecurrence[{6,3}, {1,6}, 41] (* G. C. Greubel, Oct 10 2022 *)

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

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

a(n) = (3+2*sqrt(3))^n*(sqrt(3)/4+1/2) + (1/2-sqrt(3)/4)*(3-2*sqrt(3))^n.
a(n) = (-i*sqrt(3))^n * ChebyshevU(n, isqrt(3)), i^2=-1.
From Johannes W. Meijer, Aug 09 2010: (Start)
G.f.: 1/(1 - 6*x - 3*x^2).
Limit_{k->oo} a(n+k)/a(k) = A141041(n) + A090018(n-1)*sqrt(12) for n >= 1.
Limit_{n->oo} A141041(n)/A090018(n-1) = sqrt(12). (End)
a(n) = Sum_{k=0..n} A099089(n,k)*3^k. - Philippe Deléham, Nov 21 2011
E.g.f.: exp(3*x)*(2*cosh(2*sqrt(3)*x) + sqrt(3)*sinh(2*sqrt(3)*x))/2. - Stefano Spezia, Apr 23 2025

Typo in Mathematica program corrected by Vincenzo Librandi, Nov 15 2011

A015577 a(n+1) = 8a(n) + 9a(n-1), a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 8, 73, 656, 5905, 53144, 478297, 4304672, 38742049, 348678440, 3138105961, 28242953648, 254186582833, 2287679245496, 20589113209465, 185302018885184, 1667718169966657, 15009463529699912, 135085171767299209, 1215766545905692880, 10941898913151235921

Offset: 0

Olivier Gérard

Binomial transform is A011557, with a leading zero. - Paul Barry, Jul 09 2003
Number of walks of length n between any two distinct nodes of the complete graph K_10. Example: a(2) = 8 because the walks of length 2 between the nodes A and B of the complete graph ABCDEFGHIJ are: ACB, ADB, AEB, AFB, AGB, AHB, AIB and AJB. - Emeric Deutsch, Apr 01 2004
The ratio a(n+1)/a(n) converges to 9 as n approaches infinity. - Felix P. Muga II, Mar 09 2014

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Jean-Paul Allouche, Jeffrey Shallit, Zhixiong Wen, Wen Wu, Jiemeng Zhang, Sum-free sets generated by the period-k-folding sequences and some Sturmian sequences, arXiv:1911.01687 [math.CO], 2019.
Martin Burtscher, Igor Szczyrba, Rafał Szczyrba, Analytic Representations of the n-anacci Constants and Generalizations Thereof, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.5.
Index entries for linear recurrences with constant coefficients, signature (8,9).

Cf. A001045, A078008, A097073, A115341, A015518, A054878, A015521, A109499, A015531, A109500, A109501, A015552, A093134, A015565. - Vladimir Joseph Stephan Orlovsky, Dec 11 2008

[Round(9^n/10): n in [0..30]]; // Vincenzo Librandi, Jun 24 2011

seq(round(9^n/10),n=0..25); # Mircea Merca, Dec 28 2010

k=0;lst={k};Do[k=9^n-k;AppendTo[lst, k], {n, 0, 5!}];lst (* Vladimir Joseph Stephan Orlovsky, Dec 11 2008 *)
Table[(9^n - (-1)^n)/10, {n,0,30}] (* or *) LinearRecurrence[{8,9}, {0,1}, 30] (* G. C. Greubel, Jan 06 2018 *)

a[0]:0$
a[n]:=9^(n-1)-a[n-1]$
A015577(n):=a[n]$
makelist(A015577(n),n,0,30); /* Martin Ettl, Nov 05 2012 */

A015577_vec(N=20)=Vec(O(x^N)+1/(1-8*x-9*x^2), -N-1) \\ M. F. Hasler, Jun 14 2008, edited Oct 25 2019

for(n=0,30, print1((9^n - (-1)^n)/10, ", ")) \\ G. C. Greubel, Jan 06 2018

apply( {A015577(n)=9^n\/10}, [0..25]) \\ M. F. Hasler, Oct 25 2019

[lucas_number1(n,8,-9) for n in range(0, 19)] # Zerinvary Lajos, Apr 25 2009

From Paul Barry, Jul 09 2003: (Start)
G.f.: x/((1+x)*(1-9*x)).
E.g.f. exp(4*x)*sinh(5*x)/5.
a(n) = (9^n - (-1)^n)/10. (End)
a(n) = 9^(n-1)-a(n-1). - Emeric Deutsch, Apr 01 2004
a(n) = round(9^n/10). - Mircea Merca, Dec 28 2010

Extended by T. D. Noe, May 23 2011

A062160 Square array T(n,k) = (n^k - (-1)^k)/(n+1), n >= 0, k >= 0, read by falling antidiagonals.

Original entry on oeis.org

0, 1, 0, -1, 1, 0, 1, 0, 1, 0, -1, 1, 1, 1, 0, 1, 0, 3, 2, 1, 0, -1, 1, 5, 7, 3, 1, 0, 1, 0, 11, 20, 13, 4, 1, 0, -1, 1, 21, 61, 51, 21, 5, 1, 0, 1, 0, 43, 182, 205, 104, 31, 6, 1, 0, -1, 1, 85, 547, 819, 521, 185, 43, 7, 1, 0, 1, 0, 171, 1640, 3277, 2604, 1111, 300, 57, 8, 1, 0, -1, 1, 341, 4921, 13107, 13021, 6665, 2101, 455, 73, 9, 1, 0

Offset: 0

Henry Bottomley, Jun 08 2001

For n >= 1, T(n, k) equals the number of walks of length k between any two distinct vertices of the complete graph K_(n+1). - Peter Bala, May 30 2024

			From _Seiichi Manyama_, Apr 12 2019: (Start)
Square array begins:
   0, 1, -1,  1,  -1,    1,    -1,      1, ...
   0, 1,  0,  1,   0,    1,     0,      1, ...
   0, 1,  1,  3,   5,   11,    21,     43, ...
   0, 1,  2,  7,  20,   61,   182,    547, ...
   0, 1,  3, 13,  51,  205,   819,   3277, ...
   0, 1,  4, 21, 104,  521,  2604,  13021, ...
   0, 1,  5, 31, 185, 1111,  6665,  39991, ...
   0, 1,  6, 43, 300, 2101, 14706, 102943, ... (End)

Seiichi Manyama, Antidiagonals n = 0..139, flattened
M. Dukes and C. D. White, Web Matrices: Structural Properties and Generating Combinatorial Identities, arXiv:1603.01589 [math.CO], 2016.

Rows include A062157, A000035, A001045, A015518, A015521, A015531, A015540, A015552, A015565, A015577, A015585, A015592, A015609.
Columns include A000004, A000012, A023443, A002061, A062158, A060884, A062159, A060888.
Related to repunits in negative bases (cf. A055129 for positive bases).
Main diagonal gives A081216.
Cf. A109502.

seq(print(seq((n^k - (-1)^k)/(n+1), k = 0..10)), n = 0..10); # Peter Bala, May 31 2024

T[n_,k_]:=(n^k - (-1)^k)/(n+1); Join[{0},Table[Reverse[Table[T[n-k,k],{k,0,n}]],{n,12}]]//Flatten (* Stefano Spezia, Feb 20 2024 *)

T(n, k) = n^(k-1) - n^(k-2) + n^(k-3) - ... + (-1)^(k-1) = n^(k-1) - T(n, k-1) = n*T(n, k-1) - (-1)^k = (n - 1)*T(n, k-1) + n*T(n, k-2) = round[n^k/(n+1)] for n > 1.
T(n, k) = (-1)^(k+1) * resultant( n*x + 1, (x^k-1)/(x-1) ). - Max Alekseyev, Sep 28 2021
G.f. of row n: x/((1+x) * (1-n*x)). - Seiichi Manyama, Apr 12 2019
E.g.f. of row n: (exp(n*x) - exp(-x))/(n+1). - Stefano Spezia, Feb 20 2024
From Peter Bala, May 31 2024: (Start)
Binomial transform of the m-th row: Sum_{k = 0..n} binomial(n, k)*T(m, k) = (m + 1)^(n-1) for n >= 1.
Let R(m, x) denote the g.f. of the m-th row of the square array. Then R(m_1, x) o R(m_2, x) = R(m_1 + m_2 + m_1*m_2, x),  where o denotes the black diamond product of power series as defined by Dukes and White. Cf. A109502.
T(m_1 + m_2 + m_1*m_2, k) = Sum_{i = 0..k} Sum_{j = i..k} binomial(k, i)* binomial(k-i, j-i)*T(m_1, j)*T(m_2, k-i).  (End)

A093134 A Jacobsthal trisection.

Original entry on oeis.org

1, 0, 8, 56, 456, 3640, 29128, 233016, 1864136, 14913080, 119304648, 954437176, 7635497416, 61083979320, 488671834568, 3909374676536, 31274997412296, 250199979298360, 2001599834386888, 16012798675095096, 128102389400760776, 1024819115206086200, 8198552921648689608

Offset: 0

Paul Barry, Mar 23 2004

Counts closed walks at a vertex of the complete graph on 9 nodes K_9.

Second binomial transform is A047855.

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (7,8).

Other sequences with a(n+1) = 8^n - a(n) are A001045, A078008, A097073, A115341, A015518, A054878, A015521, A109499, A015531, A109500, A109501, A015552, A015565. - Vladimir Joseph Stephan Orlovsky, Dec 11 2008

Cf. A047855.

[(8^n/9+8*(-1)^n/9): n in [0..20]]; // Vincenzo Librandi, Oct 11 2011

k=0;lst={1, k};Do[k=8^n-k;AppendTo[lst, k], {n, 1, 5!}];lst (* Vladimir Joseph Stephan Orlovsky, Dec 11 2008 *)
Table[(8^n + 8*(-1)^n)/9, {n,0,30}] (* or *) LinearRecurrence[{7,8}, {1,0}, 30] (* G. C. Greubel, Jan 06 2018 *)

for(n=0,30, print1((8^n + 8*(-1)^n)/9, ", ")) \\ G. C. Greubel, Jan 06 2018

G.f.: (1-7*x)/(1 - 7*x - 8*x^2).
a(n) = (8^n + 8*(-1)^n)/9.
a(n) = 8*A001045(3*n-3)/3.
From Elmo R. Oliveira, Aug 17 2024: (Start)
E.g.f.: exp(-x)*(exp(9*x) + 8)/9.
a(n) = 7*a(n-1) + 8*a(n-2) for n > 1. (End)

A033117 Base 7 digits are, in order, the first n terms of the periodic sequence with initial period 1,0.

Original entry on oeis.org

1, 7, 50, 350, 2451, 17157, 120100, 840700, 5884901, 41194307, 288360150, 2018521050, 14129647351, 98907531457, 692352720200, 4846469041400, 33925283289801, 237476983028607, 1662338881200250, 11636372168401750, 81454605178812251, 570182236251685757, 3991275653761800300

Offset: 1

Clark Kimberling

Partial sums of round(7^n/8), A015552. - Mircea Merca, Dec 28 2010

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (7,1,-7).

Cf. A015552.

[Floor((7*7^n-1)/48): n in [1..30]]; // Vincenzo Librandi, Jun 25 2011

I:=[1, 7, 50]; [n le 3 select I[n] else 7*Self(n-1)+Self(n-2)-7*Self(n-3): n in [1..30]]; // Vincenzo Librandi, Mar 26 2014

A033117 := proc(n) add( round(7^i/8),i=0..n) ; end proc:

Join[{a=1,b=7},Table[c=6*b+7*a+1;a=b;b=c,{n,60}]] (* Vladimir Joseph Stephan Orlovsky, Feb 06 2011 *)
Module[{nn=30,c},c=PadRight[{},nn,{1,0}];Table[FromDigits[Take[c,n],7],{n,nn}]] (* or *) LinearRecurrence[{7,1,-7},{1,7,50},30] (* Harvey P. Dale, Feb 13 2014 *)
CoefficientList[Series[1/((1 - x) (1 - 7 x) (1 + x)), {x, 0, 50}], x] (* Vincenzo Librandi, Mar 26 2014 *)

G.f.: x / ((1-x)*(1-7*x)*(1+x)).
a(n) = 7*a(n-1) + a(n-2) - 7*a(n-3).
a(n) = (7*7^n - 4 - 3*(-1)^n)/48. - Bruno Berselli, Jan 19 2011
a(n) = (1/6)*floor(7^(n+1)/8) = floor((7*7^n-1)/48) = ceiling((7*7^n-7)/48) = round((7*7^n-7)/48) = round((7*7^n-4)/48); a(n) = a(n-2) + 7^(n-1), n > 2. - Mircea Merca, Dec 28 2010

A015592 a(n) = 10a(n-1) + 11a(n-2).

Original entry on oeis.org

0, 1, 10, 111, 1220, 13421, 147630, 1623931, 17863240, 196495641, 2161452050, 23775972551, 261535698060, 2876892678661, 31645819465270, 348104014117971, 3829144155297680, 42120585708274481, 463326442791019290, 5096590870701212191, 56062499577713334100

Offset: 0

Olivier Gérard

Number of walks of length n between any two distinct nodes of the complete graph K_12. Example: a(2)=10 because the walks of length 2 between the nodes A and B of the complete graph ABCDEFGHIJKL are ACB, ADB, AEB, AFB, AGB, AHB, AIB, AJB, AKB and ALB. - Emeric Deutsch, Apr 01 2004

General form: k=11^n-k. Also: A001045, A078008, A097073, A115341, A015518, A054878, A015521, A109499, A015531, A109500, A109501, A015552, A093134, A015565, A015577, A015585. - Vladimir Joseph Stephan Orlovsky, Dec 11 2008

Vincenzo Librandi, Table of n, a(n) for n = 0..900
Index entries for linear recurrences with constant coefficients, signature (10,11).

Cf. A001045, A078008, A097073, A115341, A015518, A054878, A015521, A109499, A015531, A109500, A109501, A015552, A093134, A015565, A015577, A015585. - Vladimir Joseph Stephan Orlovsky, Dec 11 2008

[-(1/12)*(-1)^n+(1/12)*11^n: n in [0..20]]; // Vincenzo Librandi, Oct 11 2011

k=0;lst={k};Do[k=11^n-k;AppendTo[lst, k], {n, 0, 5!}];lst (* Vladimir Joseph Stephan Orlovsky, Dec 11 2008 *)

[lucas_number1(n,10,-11) for n in range(0, 18)] # Zerinvary Lajos, Apr 26 2009

a(n) = 11^(n-1) - a(n-1). G.f.: x/(1 - 10x - 11x^2). - Emeric Deutsch, Apr 01 2004
From Elmo R. Oliveira, Aug 17 2024: (Start)
E.g.f.: exp(5*x)*sinh(6*x)/6.
a(n) = (11^n - (-1)^n)/12. (End)

A189800 a(n) = 6a(n-1) + 8a(n-2), with a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 6, 44, 312, 2224, 15840, 112832, 803712, 5724928, 40779264, 290475008, 2069084160, 14738305024, 104982503424, 747801460736, 5326668791808, 37942424436736, 270267896954880, 1925146777223168, 13713023838978048, 97679317251653632, 695780094221746176

Offset: 0

Vladimir Joseph Stephan Orlovsky, May 24 2011

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (6,8).

Sequences of the form a(n) = c*a(n-1) + d*a(n-2), with a(0)=0, a(1)=1:
c/d...1.......2.......3.......4.......5.......6.......7.......8.......9......10
.A000045,A001045,A006130,A006131,A015440,A015441,A015442,A015443,A015445,A015446
.A000129,A002605,A015518,A063727,A002532,A083099,A015519,A003683,A002534,A083102
.A006190,A007482,A030195,A015521,A015523,A083858,A015524,A015525,A099012,A015528
.A001076,A090017,A015530,A057087,A015531,A085939,A015532,A180222,A015533,A180226
.A052918,A015535,A015536,A015537,A057088,A015540,A015541,A015544,A015545,A180250
.A005668,A135030,A090018,A135032,A015551,A057089,A015552,A189800,A189801,A190005
.A054413,A015555,A015559,A015561,A015562,A015564,A057090,A015565,A015566,A015568
.A041025,A190331,A015574,A190510,A015575,A190560,A015576,A057091,A015577,A190953
.A099371,A015579,A181353,A015580,A015581,A153191,A015583,A015584,A057092,A015585
A041041,A191014,A015588,A190954,A190955,A190956,A015589,A190957,A015591,A057093

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

LinearRecurrence[{6, 8}, {0, 1}, 50]
CoefficientList[Series[-(x/(-1+6 x+8 x^2)),{x,0,50}],x] (* Harvey P. Dale, Jul 26 2011 *)

a(n)=([0,1; 8,6]^n*[0;1])[1,1] \\ Charles R Greathouse IV, Oct 03 2016

G.f.: x/(1 - 2*x*(3+4*x)). - Harvey P. Dale, Jul 26 2011

cp's OEIS Frontend

A206369 a(p^k) = p^k - p^(k-1) + p^(k-2) - ... +- 1, and then extend by multiplicativity.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

PARI

Python

Formula

A015565 a(n) = 7*a(n-1) + 8*a(n-2), a(0) = 0, a(1) = 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

A015585 a(n) = 9*a(n-1) + 10*a(n-2).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Sage

Sage

Formula

Extensions

A090018 a(n) = 6*a(n-1) + 3*a(n-2) for n > 2, a(0)=1, a(1)=6.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

Extensions

A015577 a(n+1) = 8*a(n) + 9*a(n-1), a(0) = 0, a(1) = 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Maxima

PARI

PARI

PARI

Sage

A015565 a(n) = 7a(n-1) + 8a(n-2), a(0) = 0, a(1) = 1.

A015585 a(n) = 9a(n-1) + 10a(n-2).

A090018 a(n) = 6a(n-1) + 3a(n-2) for n > 2, a(0)=1, a(1)=6.

A015577 a(n+1) = 8a(n) + 9a(n-1), a(0) = 0, a(1) = 1.

A015592 a(n) = 10a(n-1) + 11a(n-2).

A189800 a(n) = 6a(n-1) + 8a(n-2), with a(0)=0, a(1)=1.