A010876 -id:A010876 - cp's OEIS frontend

A130485 a(n) = Sum_{k=0..n} (k mod 7) (Partial sums of A010876).

0, 1, 3, 6, 10, 15, 21, 21, 22, 24, 27, 31, 36, 42, 42, 43, 45, 48, 52, 57, 63, 63, 64, 66, 69, 73, 78, 84, 84, 85, 87, 90, 94, 99, 105, 105, 106, 108, 111, 115, 120, 126, 126, 127, 129, 132, 136, 141, 147, 147, 148, 150, 153, 157, 162, 168, 168, 169, 171, 174, 178, 183

Offset: 0

Hieronymus Fischer, May 31 2007

Let A be the Hessenberg n X n matrix defined by A[1,j] = j mod 7, A[i,i]:=1, A[i,i-1]=-1. Then, for n >= 1, a(n)=det(A). - Milan Janjic, Jan 24 2010

Shawn A. Broyles, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (1,0,0,0,0,0,1,-1).

Cf. A010872, A010873, A010874, A010875, A010877, A130481, A130482, A130483, A130484.

a:=[0,1,3,6,10,15,21,21];; for n in [9..71] do a[n]:=a[n-1]+a[n-7]-a[n-8]; od; a; # G. C. Greubel, Aug 31 2019

I:=[0,1,3,6,10,15,21,21]; [n le 8 select I[n] else Self(n-1) + Self(n-7) - Self(n-8): n in [1..71]]; // G. C. Greubel, Aug 31 2019

a:=n->add(chrem( [n,j], [1,7] ),j=1..n):seq(a(n), n=1..70); # Zerinvary Lajos, Apr 07 2009

LinearRecurrence[{1,0,0,0,0,0,1,-1},{0,1,3,6,10,15,21,21},70] (* Harvey P. Dale, Jul 30 2017 *)

concat(0,Vec((1-7*x^6+6*x^7)/(1-x^7)/(1-x)^3+O(x^70))) \\ Charles R Greathouse IV, Dec 22 2011

def A130485_list(prec):
    P. = PowerSeriesRing(ZZ, prec)
    return P(x*(1-7*x^6+6*x^7)/((1-x^7)*(1-x)^3)).list()
A130485_list(70) # G. C. Greubel, Aug 31 2019

a(n) = 21*floor(n/7) + A010876(n)*(A010876(n) + 1)/2.
G.f.: (Sum_{k=1..6} k*x^k)/((1-x^7)*(1-x)).
G.f.: x*(1 - 7*x^6 + 6*x^7)/((1-x^7)*(1-x)^3).

A010888 Digital root of n (repeatedly add the digits of n until a single digit is reached).

Original entry on oeis.org

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

Offset: 0

N. J. A. Sloane

This is sometimes also called the additive digital root of n.
n mod 9 (A010878) is a very similar sequence.
Partial sums are given by A130487(n-1) + n (for n > 0). - Hieronymus Fischer, Jun 08 2007
Decimal expansion of 13717421/111111111 is 0.123456789123456789123456789... with period 9. - Eric Desbiaux, May 19 2008
Decimal expansion of 13717421 / 1111111110 = 0.0[123456789] (periodic) - Daniel Forgues, Feb 27 2017
a(A005117(n)) < 9. - Reinhard Zumkeller, Mar 30 2010
My friend Jahangeer Kholdi has found that 19 is the smallest prime p such that for each number n, a(p*n) = a(n). In fact we have: a(m*n) = a(a(m)*a(n)) so all numbers with digital root 1 (numbers of the form 9k + 1) have this property. See comment lines of A017173. Also we have a(m+n) = a(a(m) + a(n)). - Farideh Firoozbakht, Jul 23 2010

			The digits of 37 are 3 and 7, and 3 + 7 = 10. And the digits of 10 are 1 and 0, and 1 + 0 = 1, so a(37) = 1.

Martin Gardner, Mathematics, Magic and Mystery, 1956.

N. J. A. Sloane, Table of n, a(n) for n = 0..10000
Eric Weisstein's World of Mathematics, Digitaddition
Eric Weisstein's World of Mathematics, Digital Root
Wikipedia, Vedic square
Index entries for Colombian or self numbers and related sequences
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,0,1).

Cf. A007953, A007954, A031347, A113217, A113218, A010878 (n mod 9), A010872, A010873, A010874, A010875, A010876, A010877, A010879, A004526, A002264, A002265, A002266, A017173, A031286 (additive persistence of n), (multiplicative digital root of n), A031346 (multiplicative persistence of n).

a010888 = until (< 10) a007953
-- Reinhard Zumkeller, Oct 17 2011, May 12 2011

[n eq 0 select 0 else 1+(n-1) mod 9: n in [0..110]]; // Bruno Berselli, Mar 18 2016

A010888 := n->if n=0 then 0 else ((n-1) mod 9) + 1; fi; # N. J. A. Sloane, Feb 20 2013

Join[{0}, Array[Mod[ # - 1, 9] + 1 &, 104]] (* Robert G. Wilson v, Jan 04 2006 *)
Join[Range[0, 1], Table[n - 9 Floor[(n - 1) / 9], {n, 2, 100}]] (* José de Jesús Camacho Medina, Nov 10 2014 *) (* Corrected by Vincenzo Librandi, Nov 11 2014 *)
Join[{0},LinearRecurrence[{0, 0, 0, 0, 0, 0, 0, 0, 1},{1, 2, 3, 4, 5, 6, 7, 8, 9},104]] (* Ray Chandler, Aug 26 2015 *)
Table[FixedPoint[Total[IntegerDigits[#, 10]] &, n], {n, 0, 104}] (* IWABUCHI Yu(u)ki, Jun 03 2016 *)

A010888(n)=if(n,(n-1)%9+1) \\ M. F. Hasler, Jan 04 2011

def A010888(n):
    return 1 + (n - 1) % 9 if n else 0 # Chai Wah Wu, Aug 23 2014, Apr 23 2023

0 :: List.fill(10)(1 to 9).flatten // Alonso del Arte, Feb 01 2020

If n = 0 then a(n) = 0; otherwise a(n) = (n reduced mod 9), but if the answer is 0 change it to 9.
Equivalently, if n = 0 then a(n) = 0, otherwise a(n) = (n - 1 reduced mod 9) + 1.
If the initial 0 term is ignored, the sequence is periodic with period 9.
From Hieronymus Fischer, Jun 08 2007: (Start)
a(n) = A010878(n-1) + 1 (for n > 0).
G.f.: g(x) = x*(Sum_{k = 0..8}(k+1)*x^k)/(1 - x^9). Also: g(x) = x(9x^10 - 10x^9 + 1)/((1 - x^9)(1 - x)^2). (End)
a(n) = n - 9*floor((n-1)/9), for n > 0. - José de Jesús Camacho Medina, Nov 10 2014

A130481 a(n) = Sum_{k=0..n} (k mod 3) (i.e., partial sums of A010872).

Original entry on oeis.org

0, 1, 3, 3, 4, 6, 6, 7, 9, 9, 10, 12, 12, 13, 15, 15, 16, 18, 18, 19, 21, 21, 22, 24, 24, 25, 27, 27, 28, 30, 30, 31, 33, 33, 34, 36, 36, 37, 39, 39, 40, 42, 42, 43, 45, 45, 46, 48, 48, 49, 51, 51, 52, 54, 54, 55, 57, 57, 58, 60, 60, 61, 63, 63, 64, 66, 66, 67, 69, 69, 70, 72, 72

Offset: 0

Hieronymus Fischer, May 29 2007

Essentially the same as A092200. - R. J. Mathar, Jun 13 2008
Let A be the Hessenberg n X n matrix defined by: A[1,j]=j mod 3, A[i,i]:=1, A[i,i-1]=-1. Then, for n>=1, a(n)=det(A). - Milan Janjic, Jan 24 2010
2-adic valuation of A104537(n+1). - Gerry Martens, Jul 14 2015
Conjecture: a(n) is the exponent of the largest power of 2 that divides all the entries of the matrix {{3,1},{1,-1}}^n. - Greg Dresden, Sep 09 2018

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

Cf. A010872, A010873, A010874, A010875, A010876, A010877, A130482, A130483, A130484.

Cf. A092200, A104537.

List([0..80], n-> Int((n+1)/3) + Int(2*(n+1)/3)); # G. C. Greubel, Aug 31 2019

[Floor((n+1)/3) + Floor(2*(n+1)/3): n in [0..80]]; // G. C. Greubel, Aug 31 2019

seq(coeff(series(x*(1+2*x)/((1-x^3)*(1-x)), x, n+1), x, n), n = 0..80); # G. C. Greubel, Aug 31 2019

a[n_]:= Floor[(n+1)/3] + Floor[2(n+1)/3]; Table[a[n], {n, 0, 80}] (* Clark Kimberling, May 28 2012 *)
a[n_]:= IntegerExponent[A104537[n + 1], 2];
Table[a[n], {n, 0, 80}]  (* Gerry Martens, Jul 14 2015 *)
CoefficientList[Series[x(1+2x)/((1-x^3)(1-x)), {x, 0, 80}], x] (* Stefano Spezia, Sep 09 2018 *)
LinearRecurrence[{1,0,1,-1},{0,1,3,3},100] (* Harvey P. Dale, Jun 14 2021 *)

main(size)=my(n,k);vector(size,n,sum(k=0,n,k%3)) \\ Anders Hellström, Jul 14 2015

first(n)=my(s); concat(0, vector(n,k,s+=k%3)) \\ Charles R Greathouse IV, Jul 14 2015

a(n)=n\3*3+[0,1,3][n%3+1] \\ Charles R Greathouse IV, Jul 14 2015

def A130481_list(prec):
    P. = PowerSeriesRing(ZZ, prec)
    return P(x*(1+2*x)/((1-x^3)*(1-x))).list()
A130481_list(80) # G. C. Greubel, Aug 31 2019

a(n) = 3*floor(n/3) + A010872(n)*(A010872(n) + 1)/2.
G.f.: x*(1 + 2*x)/((1-x^3)*(1-x)).
a(n) = n + 1 - (Fibonacci(n+1) mod 2). - Gary Detlefs, Mar 13 2011
a(n) = floor((n+1)/3) + floor(2*(n+1)/3). - Clark Kimberling, May 28 2010
a(n) = n when n+1 is not a multiple of 3, and a(n) = n+1 when n+1 is a multiple of 3. - Dennis P. Walsh, Aug 06 2012
a(n) = n + 1 - sign((n+1) mod 3). - Wesley Ivan Hurt, Sep 25 2017
a(n) = n + (1-cos(2*(n+2)*Pi/3))/3 + sin(2*(n+2)*Pi/3)/sqrt(3). - Wesley Ivan Hurt, Sep 27 2017
a(n) = n + 1 - (n+1)^2 mod 3. - Ammar Khatab, Aug 14 2020
E.g.f.: ((1 + 3*x)*cosh(x) - (cos(sqrt(3)*x/2) + sqrt(3)*sin(sqrt(3)*x/2))*(cosh(x/2) - sinh(x/2)) + (1 + 3*x)*sinh(x))/3. - Stefano Spezia, May 28 2021
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(3*sqrt(3)) + log(2)/3. - Amiram Eldar, Sep 17 2022

A130482 a(n) = Sum_{k=0..n} (k mod 4) (Partial sums of A010873).

Original entry on oeis.org

0, 1, 3, 6, 6, 7, 9, 12, 12, 13, 15, 18, 18, 19, 21, 24, 24, 25, 27, 30, 30, 31, 33, 36, 36, 37, 39, 42, 42, 43, 45, 48, 48, 49, 51, 54, 54, 55, 57, 60, 60, 61, 63, 66, 66, 67, 69, 72, 72, 73, 75, 78, 78, 79, 81, 84, 84, 85, 87, 90, 90, 91, 93, 96, 96, 97, 99, 102, 102, 103, 105

Offset: 0

Hieronymus Fischer, May 29 2007

Let A be the Hessenberg n X n matrix defined by: A[1,j]=j mod 4, A[i,i]:=1, A[i,i-1]=-1. Then, for n>=1, a(n)=det(A). - Milan Janjic, Jan 24 2010

Colin Barker, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (1,0,0,1,-1).

Cf. A010872, A010874, A010875, A010876, A010877. A130481, A130483, A130484, A130485.

a:=[0,1,3,6,6];; for n in [6..71] do a[n]:=a[n-1]+a[n-4]-a[n-5]; od; a; # G. C. Greubel, Aug 31 2019

I:=[0,1,3,6,6]; [n le 5 select I[n] else Self(n-1) + Self(n-4) - Self(n-5): n in [1..71]]; // G. C. Greubel, Aug 31 2019

a:=n->add(chrem( [n,j], [1,4] ),j=1..n):seq(a(n), n=0..70); # Zerinvary Lajos, Apr 07 2009

Table[(6*n +(1-(-1)^n)*(1+2*I^(n+1)))/4, {n,0,70}] (* G. C. Greubel, Aug 31 2019 *)
LinearRecurrence[{1,0,0,1,-1},{0,1,3,6,6},80] (* Harvey P. Dale, Feb 16 2024 *)

a(n) = (1 - (-1)^n - (2*I)*(-I)^n + (2*I)*I^n + 6*n) / 4 \\ Colin Barker, Oct 15 2015

def A130482_list(prec):
    P. = PowerSeriesRing(ZZ, prec)
    return P(x*(1+2*x+3*x^2)/((1-x^4)*(1-x))).list()
A130482_list(70) # G. C. Greubel, Aug 31 2019

a(n) = 6*floor(n/4) + A010873(n)*(A010873(n)+1)/2.
G.f.: x*(1 + 2*x + 3*x^2)/((1-x^4)*(1-x)).
a(n) = (1 - (-1)^n - (2*i)*(-i)^n + (2*i)*i^n + 6*n) / 4 where i = sqrt(-1). - Colin Barker, Oct 15 2015
a(n) = 3*n/2 + (n mod 2)* ( (n-1) mod 4 ) - (n mod 2)/2. - Ammar Khatab, Aug 27 2020
E.g.f.: (3*x*exp(x) - 2*sin(x) + sinh(x))/2. - Stefano Spezia, Apr 22 2021
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(4*sqrt(3)) + log(3)/4. - Amiram Eldar, Sep 17 2022

A130483 a(n) = Sum_{k=0..n} (k mod 5) (Partial sums of A010874).

Original entry on oeis.org

0, 1, 3, 6, 10, 10, 11, 13, 16, 20, 20, 21, 23, 26, 30, 30, 31, 33, 36, 40, 40, 41, 43, 46, 50, 50, 51, 53, 56, 60, 60, 61, 63, 66, 70, 70, 71, 73, 76, 80, 80, 81, 83, 86, 90, 90, 91, 93, 96, 100, 100, 101, 103, 106, 110, 110, 111, 113, 116, 120, 120, 121, 123, 126, 130, 130

Offset: 0

Hieronymus Fischer, May 29 2007

Let A be the Hessenberg n X n matrix defined by: A[1,j]=j mod 5, A[i,i]=1, A[i,i-1]=-1. Then, for n>=1, a(n)=det(A). - Milan Janjic, Jan 24 2010

Shawn A. Broyles, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (1,0,0,0,1,-1).

Cf. A010872, A010873, A010875, A010876, A010877, A130481, A130482, A130484, A130485.

a:=[0,1,3,6,10,10];; for n in [7..71] do a[n]:=a[n-1]+a[n-5]-a[n-6]; od; a; # G. C. Greubel, Aug 31 2019

I:=[0,1,3,6,10,10]; [n le 6 select I[n] else Self(n-1) + Self(n-5) - Self(n-6): n in [1..71]]; // G. C. Greubel, Aug 31 2019

seq(coeff(series(x*(1+2*x+3*x^2+4*x^3)/((1-x^5)*(1-x)), x, n+1), x, n), n = 0..70); # G. C. Greubel, Aug 31 2019

Accumulate[Mod[Range[0,70],5]] (* or *) Accumulate[PadRight[{},70,{0,1,2,3,4}]] (* Harvey P. Dale, Nov 11 2016 *)

a(n) = sum(k=0, n, k % 5); \\ Michel Marcus, Apr 28 2018

def A130483_list(prec):
    P. = PowerSeriesRing(ZZ, prec)
    return P(x*(1+2*x+3*x^2+4*x^3)/((1-x^5)*(1-x))).list()
A130483_list(70) # G. C. Greubel, Aug 31 2019

a(n) = 10*floor(n/5) + A010874(n)*(A010874(n)+1)/2.
G.f.: x*(1 + 2*x + 3*x^2 + 4*x^3)/((1-x^5)*(1-x)).
From Wesley Ivan Hurt, Jul 23 2016: (Start)
a(n) = a(n-5) - a(n-6) for n>5; a(n) = a(n-5) + 10 for n>4.
a(n) = 10 + Sum_{k=1..4} k*floor((n-k)/5). (End)
a(n) = ((n mod 5)^2 - 3*(n mod 5) + 4*n)/2. - Ammar Khatab, Aug 13 2020

A010882 Period 3: repeat [1, 2, 3].

Original entry on oeis.org

1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3

Offset: 0

N. J. A. Sloane

Partial sums are given by A130481(n)+n+1. - Hieronymus Fischer, Jun 08 2007
41/333 = 0.123123123... - Eric Desbiaux, Nov 03 2008
Terms of the simple continued fraction for 3/(sqrt(37)-4). - Paolo P. Lava, Feb 16 2009
This is the lexicographically earliest sequence with no substring of more than 1 term being a palindrome. - Franklin T. Adams-Watters, Nov 24 2013

Index entries for linear recurrences with constant coefficients, signature (0,0,1).
Index entries for sequences that are fixed points of mappings

Cf. A010872, A010873, A010874, A010875, A010876, A004526, A002264, A002265, A002266, A177036 (decimal expansion of (4+sqrt(37))/7), A214090.

Cf. A130481, A180593.

a010882 = (+ 1) . (`mod` 3)
a010882_list = cycle [1,2,3]
-- Reinhard Zumkeller, Mar 20 2013

&cat[[1..3]^^30]; // Vincenzo Librandi, Feb 04 2016

seq(op([1, 2, 3]), n=0..50); # Wesley Ivan Hurt, Jul 05 2016

Nest[ Flatten[ # /. {1 -> {1, 2}, 2 -> {3, 1}, 3 -> {2, 3}}] &, {1}, 7] (* Robert G. Wilson v, Mar 08 2005 *)
PadRight[{},120,{1,2,3}] (* Harvey P. Dale, Apr 09 2018 *)

a(n) = 1 + n%3; \\ Michel Marcus, Feb 04 2016

G.f.: (1+2x+3x^2)/(1-x^3). - Paul Barry, May 25 2003
a(n) = 1 + (n mod 3). - Paolo P. Lava, Nov 21 2006
a(n) = A010872(n) + 1. - Hieronymus Fischer, Jun 08 2007
a(n) = 6 - a(n-1) - a(n-2) for n > 1. - Reinhard Zumkeller, Apr 13 2008
a(n) = n+1-3*floor(n/3) = floor(41*10^(n+1)/333)-floor(41*10^n/333)*10; a(n)-a(n-3)=0 with n>2. - Bruno Berselli, Jun 28 2010
a(n) = A180593(n+1)/3. - Reinhard Zumkeller, Oct 25 2010
a(n) = floor((4*n+3)/3) mod 4. - Gary Detlefs, May 15 2011
a(n) = -cos(2/3*Pi*n)-1/3*3^(1/2)*sin(2/3*Pi*n)+2. - Leonid Bedratyuk, May 13 2012
E.g.f.: 2*(3*exp(3*x/2) - sqrt(3)*cos(Pi/6-sqrt(3)*x/2))*exp(-x/2)/3. - Ilya Gutkovskiy, Jul 05 2016

A130484 a(n) = Sum_{k=0..n} (k mod 6) (Partial sums of A010875).

Original entry on oeis.org

0, 1, 3, 6, 10, 15, 15, 16, 18, 21, 25, 30, 30, 31, 33, 36, 40, 45, 45, 46, 48, 51, 55, 60, 60, 61, 63, 66, 70, 75, 75, 76, 78, 81, 85, 90, 90, 91, 93, 96, 100, 105, 105, 106, 108, 111, 115, 120, 120, 121, 123, 126, 130, 135, 135, 136, 138, 141, 145, 150, 150, 151, 153

Offset: 0

Hieronymus Fischer, May 31 2007

Let A be the Hessenberg n X n matrix defined by A[1,j] = j mod 6, A[i,i]=1, A[i,i-1]=-1. Then, for n >= 1, a(n)=det(A). - Milan Janjic, Jan 24 2010

Shawn A. Broyles, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (1,0,0,0,0,1,-1).

Cf. A010872, A010873, A010874, A010876, A010877, A130481, A130482, A130483, A130485.

a:=[0,1,3,6,10,15,15];; for n in [8..71] do a[n]:=a[n-1]+a[n-6]-a[n-7]; od; a; # G. C. Greubel, Aug 31 2019

I:=[0,1,3,6,10,15,15]; [n le 7 select I[n] else Self(n-1) + Self(n-6) - Self(n-7): n in [1..71]]; // G. C. Greubel, Aug 31 2019

seq(coeff(series(x*(1-6*x^5+5*x^6)/((1-x^6)*(1-x)^3), x, n+1), x, n), n = 0 .. 70); # G. C. Greubel, Aug 31 2019

Accumulate[Mod[Range[0,70],6]] (* or *) Accumulate[PadRight[ {},70, Range[0,5]]] (* Harvey P. Dale, Jul 12 2016 *)

a(n) = sum(k=0, n, k % 6); \\ Michel Marcus, Apr 28 2018

a(n)=n\6*15 + binomial(n%6+1,2) \\ Charles R Greathouse IV, Jan 24 2022

def A130484_list(prec):
    P. = PowerSeriesRing(ZZ, prec)
    return P(x*(1-6*x^5+5*x^6)/((1-x^6)*(1-x)^3)).list()
A130484_list(70) # G. C. Greubel, Aug 31 2019

a(n) = 15*floor(n/6) + A010875(n)*(A010875(n) + 1)/2.

G.f.: (Sum_{k=1..5} k*x^k)/((1-x^6)*(1-x)) = x*(1 - 6*x^5 + 5*x^6)/((1-x^6)*(1-x)^3).

A109720 Periodic sequence {0,1,1,1,1,1,1} or n^6 mod 7.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1

Offset: 0

Bruce Corrigan (scentman(AT)myfamily.com), Aug 09 2005

This sequence also represents n^12 mod 7; n^18 mod 7; (exponents are = 0 mod 6).
Characteristic sequence for numbers n>=1 to be relatively prime to 7. - Wolfdieter Lang, Oct 29 2008
a(n+4), n>=0, (periodic 1,1,1,0,1,1,1) is also the characteristic sequence for mod m reduced positive odd numbers (i.e., gcd(2*n+1,m)=1, n>=0) for each modulus m from 7*A003591 = [7,14,28,49,56,98,112,196,...]. [Wolfdieter Lang, Feb 04 2012]

Index entries for characteristic functions - _Reinhard Zumkeller_, Nov 30 2009
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,1)

Cf. A010876 = n mod 7; A053879 = n^2 mod 7; A070472 = m^3 mod 7; A070512 = n^4 mod 7; A070593 = n^5 mod 7.

Cf. A168185, A145568, A168184, A168182, A168181, A097325, A011558, A166486, A011655, A000035.

PadRight[{},120,{0,1,1,1,1,1,1}] (* Harvey P. Dale, Jul 09 2018 *)

a(n)=n^6%7 \\ Charles R Greathouse IV, Sep 24 2015

[power_mod(n,6,7)for n in range(0, 105)] # Zerinvary Lajos, Nov 06 2009

a(n) = 0 if n=0 mod 7; a(n)= 1 else.
G.f. = (x+x^2+x^3+x^4+x^5+x^6)/(1-x^7)= -x*(1+x)*(1+x+x^2)*(x^2-x+1) / ( (x-1)*(1+x+x^2+x^3+x^4+x^5+x^6) ).
a(n)=1-A082784(n); a(A047304(n))=1; a(A008589(n))=0; A033439(n) = SUM(a(k)*(n-k): 0<=k<=n). - Reinhard Zumkeller, Nov 30 2009
Multiplicative with a(p) = (if p=7 then 0 else 1), p prime. - Reinhard Zumkeller, Nov 30 2009
Dirichlet g.f. (1-7^(-s))*zeta(s). - R. J. Mathar, Mar 06 2011
For the general case: the characteristic function of numbers that are not multiples of m is a(n)=floor((n-1)/m)-floor(n/m)+1, m,n > 0. - Boris Putievskiy, May 08 2013

A010074 a(n) = sum of base-7 digits of a(n-1) + sum of base-7 digits of a(n-2).

Original entry on oeis.org

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

Offset: 0

N. J. A. Sloane, Leonid Broukhis

The digital sum analog (in base 7) of the Fibonacci recurrence. - Hieronymus Fischer, Jun 27 2007
a(n) and Fib(n)=A000045(n) are congruent modulo 6 which implies that (a(n) mod 6) is equal to (Fib(n) mod 6) = A082117(n-1) (for n>0). Thus (a(n) mod 6) is periodic with the Pisano period A001175(6)=24. - Hieronymus Fischer, Jun 27 2007
For general bases p>2, the inequality 2<=a(n)<=2p-3 holds (for n>2). Actually, a(n)<=11=A131319(7) for the base p=7. - Hieronymus Fischer, Jun 27 2007

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for Colombian or self numbers and related sequences

Cf. A000045, A010073, A010075, A010076, A010077, A131294, A131295, A131296, A131297, A131318, A131319, A131320.

nxt[{a_,b_}]:={b,Total[IntegerDigits[a,7]]+Total[IntegerDigits[b,7]]}; Transpose[NestList[nxt,{0,1},80]][[1]] (* Harvey P. Dale, Oct 12 2013 *)

Periodic from n=3 with period 24. - Franklin T. Adams-Watters, Mar 13 2006
From Hieronymus Fischer, Jun 27 2007: (Start)
a(n) = a(n-1)+a(n-2)-6*(floor(a(n-1)/7)+floor(a(n-2)/7)).
a(n) = floor(a(n-1)/7)+floor(a(n-2)/7)+(a(n-1)mod 7)+(a(n-2)mod 7).
a(n) = (a(n-1)+a(n-2)+6*(A010876(a(n-1))+A010876(a(n-2))))/7.
a(n) = Fib(n)-6*sum{1A000045(n). (End)

Incorrect comment removed by Michel Marcus, Apr 29 2018

A130486 a(n) = Sum_{k=0..n} (k mod 8) (Partial sums of A010877).

Original entry on oeis.org

0, 1, 3, 6, 10, 15, 21, 28, 28, 29, 31, 34, 38, 43, 49, 56, 56, 57, 59, 62, 66, 71, 77, 84, 84, 85, 87, 90, 94, 99, 105, 112, 112, 113, 115, 118, 122, 127, 133, 140, 140, 141, 143, 146, 150, 155, 161, 168, 168, 169, 171, 174, 178, 183, 189, 196, 196, 197, 199, 202, 206

Offset: 0

Hieronymus Fischer, May 31 2007

Let A be the Hessenberg n X n matrix defined by A[1,j] = j mod 8, A[i,i]:=1, A[i,i-1]=-1. Then, for n >= 1, a(n)=det(A). - Milan Janjic, Jan 24 2010

Shawn A. Broyles, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (1,0,0,0,0,0,0,1,-1).

Cf. A010872, A010873, A010874, A010875, A010876, A010878. A130481, A130482, A130483, A130484, A130485, A130487.

a:=[0,1,3,6,10,15,21,28,28];; for n in [10..71] do a[n]:=a[n-1]+a[n-8]-a[n-9]; od; a; # G. C. Greubel, Aug 31 2019

I:=[0,1,3,6,10,15,21,28,28]; [n le 9 select I[n] else Self(n-1) + Self(n-8) - Self(n-9): n in [1..71]]; // G. C. Greubel, Aug 31 2019

seq(coeff(series(x*(1-8*x^7+7*x^8)/((1-x^8)*(1-x)^3), x, n+1), x, n), n = 0 .. 40); # G. C. Greubel, Aug 31 2019

Array[28 Floor[#1/8] + #2 (#2 + 1)/2 & @@ {#, Mod[#, 8]} &, 61, 0] (* Michael De Vlieger, Apr 28 2018 *)
Accumulate[PadRight[{},100,Range[0,7]]] (* Harvey P. Dale, Dec 21 2018 *)

a(n) = sum(k=0, n, k % 8); \\ Michel Marcus, Apr 28 2018

def A130486_list(prec):
    P. = PowerSeriesRing(ZZ, prec)
    return P(x*(1-8*x^7+7*x^8)/((1-x^8)*(1-x)^3)).list()
A130486_list(70) # G. C. Greubel, Aug 31 2019

a(n) = 28*floor(n/8) + A010877(n)*(A010877(n) + 1)/2.
G.f.: (Sum_{k=1..7} k*x^k)/((1-x^8)*(1-x)).
G.f.: x*(1 - 8*x^7 + 7*x^8)/((1-x^8)*(1-x)^3).

cp's OEIS Frontend

A130485 a(n) = Sum_{k=0..n} (k mod 7) (Partial sums of A010876).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

A010888 Digital root of n (repeatedly add the digits of n until a single digit is reached).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Scala

Formula

A130481 a(n) = Sum_{k=0..n} (k mod 3) (i.e., partial sums of A010872).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

PARI

PARI

Sage

Formula

A130482 a(n) = Sum_{k=0..n} (k mod 4) (Partial sums of A010873).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

A130483 a(n) = Sum_{k=0..n} (k mod 5) (Partial sums of A010874).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage