cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 20 results. Next

A033581 a(n) = 6*n^2.

Original entry on oeis.org

0, 6, 24, 54, 96, 150, 216, 294, 384, 486, 600, 726, 864, 1014, 1176, 1350, 1536, 1734, 1944, 2166, 2400, 2646, 2904, 3174, 3456, 3750, 4056, 4374, 4704, 5046, 5400, 5766, 6144, 6534, 6936, 7350, 7776, 8214, 8664, 9126, 9600, 10086, 10584, 11094, 11616
Offset: 0

Views

Author

Keywords

Comments

Number of edges of a complete 4-partite graph of order 4n, K_n,n,n,n. - Roberto E. Martinez II, Oct 18 2001
Number of edges of the complete bipartite graph of order 7n, K_n, 6n. - Roberto E. Martinez II, Jan 07 2002
Number of edges in the line graph of the product of two cycle graphs, each of order n, L(C_n x C_n). - Roberto E. Martinez II, Jan 07 2002
Total surface area of a cube of edge length n. See A000578 for cube volume. See A070169 and A071399 for surface area and volume of a regular tetrahedron and links for the other Platonic solids. - Rick L. Shepherd, Apr 24 2002
a(n) can represented as n concentric hexagons (see example). - Omar E. Pol, Aug 21 2011
Sequence found by reading the line from 0, in the direction 0, 6, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. Opposite numbers to the members of A003154 in the same spiral. - Omar E. Pol, Sep 08 2011
Together with 1, numbers m such that floor(2*m/3) and floor(3*m/2) are both squares. Example: floor(2*150/3) = 100 and floor(3*150/2) = 225 are both squares, so 150 is in the sequence. - Bruno Berselli, Sep 15 2014
a(n+1) gives the number of vertices in a hexagon-like honeycomb built from A003215(n) congruent regular hexagons (see link). Example: a hexagon-like honeycomb consisting of 7 congruent regular hexagons has 1 core hexagon inside a perimeter of six hexagons. The perimeter has 18 vertices. The core hexagon has 6 vertices. a(2) = 18 + 6 = 24 is the total number of vertices. - Ivan N. Ianakiev, Mar 11 2015
a(n) is the area of the Pythagorean triangle whose sides are (3n, 4n, 5n). - Sergey Pavlov, Mar 31 2017
More generally, if k >= 5 we have that the sequence whose formula is a(n) = (2*k - 4)*n^2 is also the sequence found by reading the line from 0, in the direction 0, (2*k - 4), ..., in the square spiral whose vertices are the generalized k-gonal numbers. In this case k = 5. - Omar E. Pol, May 13 2018
The sequence also gives the number of size=1 triangles within a match-made hexagon of size n. - John King, Mar 31 2019
For hexagons, the number of matches required is A045945; thus number of size=1 triangles is A033581; number of larger triangles is A307253 and total number of triangles is A045949. See A045943 for analogs for Triangles; see A045946 for analogs for Stars. - John King, Apr 04 2019

Examples

			From _Omar E. Pol_, Aug 21 2011: (Start)
Illustration of initial terms as concentric hexagons:
.
.                                 o o o o o o
.                                o           o
.              o o o o          o   o o o o   o
.             o       o        o   o       o   o
.   o o      o   o o   o      o   o   o o   o   o
.  o   o    o   o   o   o    o   o   o   o   o   o
.   o o      o   o o   o      o   o   o o   o   o
.             o       o        o   o       o   o
.              o o o o          o   o o o o   o
.                                o           o
.                                 o o o o o o
.
.    6            24                   54
.
(End)
		

Crossrefs

Bisection of A032528. Central column of triangle A001283.
Cf. A017593 (first differences).

Programs

Formula

a(n) = A000290(n)*6. - Omar E. Pol, Dec 11 2008
a(n) = A001105(n)*3 = A033428(n)*2. - Omar E. Pol, Dec 13 2008
a(n) = 12*n + a(n-1) - 6, with a(0)=0. - Vincenzo Librandi, Aug 05 2010
G.f.: 6*x*(1+x)/(1-x)^3. - Colin Barker, Feb 14 2012
For n > 0: a(n) = A005897(n) - 2. - Reinhard Zumkeller, Apr 27 2014
a(n) = 3*floor(1/(1-cos(1/n))) = floor(1/(1-n*sin(1/n))) for n > 0. - Clark Kimberling, Oct 08 2014
a(n) = t(4*n) - 4*t(n), where t(i) = i*(i+k)/2 for any k. Special case (k=1): a(n) = A000217(4*n) - 4*A000217(n). - Bruno Berselli, Aug 31 2017
From Amiram Eldar, Feb 03 2021: (Start)
Sum_{n>=1} 1/a(n) = Pi^2/36.
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/72 (A086729).
Product_{n>=1} (1 + 1/a(n)) = sqrt(6)*sinh(Pi/sqrt(6))/Pi.
Product_{n>=1} (1 - 1/a(n)) = sqrt(6)*sin(Pi/sqrt(6))/Pi. (End)
E.g.f.: 6*exp(x)*x*(1 + x). - Stefano Spezia, Aug 19 2022

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Nov 08 2001

A055874 a(n) = largest m such that 1, 2, ..., m divide n.

Original entry on oeis.org

1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2
Offset: 1

Views

Author

Leroy Quet, Jul 16 2000

Keywords

Comments

From Antti Karttunen, Nov 20 2013 & Jan 26 2014: (Start)
Differs from A232098 for the first time at n=840, where a(840)=8, while A232098(840)=7. A232099 gives all the differing positions. See also the comments at A055926 and A232099.
The positions where a(n) is an odd prime is given by A017593 up to A017593(34)=414 (so far all 3's), after which comes the first 7 at a(420). (A017593 gives the positions of 3's.)
(Continued on Jan 26 2014):
Only terms of A181062 occur as values.
A235921 gives such n where a(n^2) (= A235918(n)) differs from A071222(n-1) (= A053669(n)-1). (End)
a(n) is the largest m such that A003418(m) divides n. - David W. Wilson, Nov 20 2014
a(n) is the largest number of consecutive integers dividing n. - David W. Wilson, Nov 20 2014
A051451 gives indices where record values occur. - Gionata Neri, Oct 17 2015
Yuri Matiyasevich calls this the maximum inheritable divisor of n. - N. J. A. Sloane, Dec 14 2023

Examples

			a(12) = 4 because 1, 2, 3, 4 divide 12, but 5 does not.
		

Crossrefs

Programs

  • Haskell
    a055874 n = length $ takeWhile ((== 0) . (mod n)) [1..]
    -- Reinhard Zumkeller, Feb 21 2012, Dec 09 2010
    
  • Maple
    N:= 1000: # to get a(1) to a(N)
    A:= Vector(N,1);
    for m from 2 do
      Lm:= ilcm($1..m);
      if Lm > N then break fi;
      if Lm mod (m+1) = 0 then next fi;
      for k from 1 to floor(N/Lm) do
        A[k*Lm]:=m
      od
    od:
    convert(A,list); # Robert Israel, Nov 28 2014
  • Mathematica
    a[n_] := Module[{m = 1}, While[Divisible[n, m++]]; m - 2]; Array[a, 100] (* Jean-François Alcover, Mar 07 2016 *)
  • PARI
    a(n) = my(m = 1); while ((n % m) == 0, m++); m - 1; \\ Michel Marcus, Jan 17 2014
    
  • Python
    from itertools import count
    def A055874(n):
        for m in count(1):
            if n % m:
                return m-1 # Chai Wah Wu, Jan 02 2022
  • Scheme
    (define (A055874 n) (let loop ((m 1)) (if (not (zero? (modulo n m))) (- m 1) (loop (+ 1 m))))) ;; Antti Karttunen, Nov 18 2013
    

Formula

a(n) = A007978(n) - 1. - Antti Karttunen, Jan 26 2014
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = A064859 (Farhi, 2009). - Amiram Eldar, Jul 25 2022

A141419 Triangle read by rows: T(n, k) = A000217(n) - A000217(n - k) with 1 <= k <= n.

Original entry on oeis.org

1, 2, 3, 3, 5, 6, 4, 7, 9, 10, 5, 9, 12, 14, 15, 6, 11, 15, 18, 20, 21, 7, 13, 18, 22, 25, 27, 28, 8, 15, 21, 26, 30, 33, 35, 36, 9, 17, 24, 30, 35, 39, 42, 44, 45, 10, 19, 27, 34, 40, 45, 49, 52, 54, 55
Offset: 1

Views

Author

Roger L. Bagula, Aug 05 2008

Keywords

Comments

As a rectangle, the accumulation array of A051340.
From Clark Kimberling, Feb 05 2011: (Start)
Here all the weights are divided by two where they aren't in Cahn.
As a rectangle, A141419 is in the accumulation chain
... < A051340 < A141419 < A185874 < A185875 < A185876 < ...
(See A144112 for the definition of accumulation array.)
row 1: A000027
col 1: A000217
diag (1,5,...): A000326 (pentagonal numbers)
diag (2,7,...): A005449 (second pentagonal numbers)
diag (3,9,...): A045943 (triangular matchstick numbers)
diag (4,11,...): A115067
diag (5,13,...): A140090
diag (6,15,...): A140091
diag (7,17,...): A059845
diag (8,19,...): A140672
(End)
Let N=2*n+1 and k=1,2,...,n. Let A_{N,n-1} = [0,...,0,1; 0,...,0,1,1; ...; 0,1,...,1; 1,...,1], an n X n unit-primitive matrix (see [Jeffery]). Let M_n=[A_{N,n-1}]^4. Then t(n,k)=[M_n](1,k), that is, the n-th row of the triangle is given by the first row of M_n. - _L. Edson Jeffery, Nov 20 2011
Conjecture. Let N=2*n+1 and k=1,...,n. Let A_{N,0}, A_{N,1}, ..., A_{N,n-1} be the n X n unit-primitive matrices (again see [Jeffery]) associated with N, and define the Chebyshev polynomials of the second kind by the recurrence U_0(x) = 1, U_1(x) = 2*x and U_r(x) = 2*x*U_(r-1)(x) - U_(r-2)(x) (r>1). Define the column vectors V_(k-1) = (U_(k-1)(cos(Pi/N)), U_(k-1)(cos(3*Pi/N)), ..., U_(k-1)(cos((2*n-1)*Pi/N)))^T, where T denotes matrix transpose. Let S_N = [V_0, V_1, ..., V_(n-1)] be the n X n matrix formed by taking V_(k-1) as column k-1. Let X_N = [S_N]^T*S_N, and let [X_N](i,j) denote the entry in row i and column j of X_N, i,j in {0,...,n-1}. Then t(n,k) = [X_N](k-1,k-1), and row n of the triangle is given by the main diagonal entries of X_N. Remarks: Hence t(n,k) is the sum of squares t(n,k) = sum[m=1,...,n (U_(k-1)(cos((2*m-1)*Pi/N)))^2]. Finally, this sequence is related to A057059, since X_N = [sum_{m=1,...,n} A057059(n,m)*A_{N,m-1}] is also an integral linear combination of unit-primitive matrices from the N-th set. - L. Edson Jeffery, Jan 20 2012
Row sums: n*(n+1)*(2*n+1)/6. - L. Edson Jeffery, Jan 25 2013
n-th row = partial sums of n-th row of A004736. - Reinhard Zumkeller, Aug 04 2014
T(n,k) is the number of distinct sums made by at most k elements in {1, 2, ... n}, for 1 <= k <= n, e.g., T(6,2) = the number of distinct sums made by at most 2 elements in {1,2,3,4,5,6}. The sums range from 1, to 5+6=11. So there are 11 distinct sums. - Derek Orr, Nov 26 2014
A number n occurs in this sequence A001227(n) times, the number of odd divisors of n, see A209260. - Hartmut F. W. Hoft, Apr 14 2016
Conjecture: 2*n + 1 is composite if and only if gcd(t(n,m),m) != 1, for some m. - L. Edson Jeffery, Jan 30 2018
From Peter Munn, Aug 21 2019 in respect of the sequence read as a triangle: (Start)
A number m can be found in column k if and only if A286013(m, k) is nonzero, in which case m occurs in column k on row A286013(m, k).
The first occurrence of m is in row A212652(m) column A109814(m), which is the rightmost column in which m occurs. This occurrence determines where m appears in A209260. The last occurrence of m is in row m column 1.
Viewed as a sequence of rows, consider the subsequences (of rows) that contain every positive integer. The lexicographically latest of these subsequences consists of the rows with row numbers in A270877; this is the only one that contains its own row numbers only once.
(End)

Examples

			As a triangle:
   1,
   2,  3,
   3,  5,  6,
   4,  7,  9, 10,
   5,  9, 12, 14, 15,
   6, 11, 15, 18, 20, 21,
   7, 13, 18, 22, 25, 27, 28,
   8, 15, 21, 26, 30, 33, 35, 36,
   9, 17, 24, 30, 35, 39, 42, 44, 45,
  10, 19, 27, 34, 40, 45, 49, 52, 54, 55;
As a rectangle:
   1   2   3   4   5   6   7   8   9  10
   3   5   7   9  11  13  15  17  19  21
   6   9  12  15  18  21  24  27  30  33
  10  14  18  22  26  30  34  38  42  46
  15  20  25  30  35  40  45  50  55  60
  21  27  33  39  45  51  57  63  69  75
  28  35  42  49  56  63  70  77  84  91
  36  44  52  60  68  76  84  92 100 108
  45  54  63  72  81  90  99 108 117 126
  55  65  75  85  95 105 115 125 135 145
Since the odd divisors of 15 are 1, 3, 5 and 15, number 15 appears four times in the triangle at t(3+(5-1)/2, 5) in column 5 since 5+1 <= 2*3, t(5+(3-1)/2, 3), t(1+(15-1)/2, 2*1) in column 2 since 15+1 > 2*1, and t(15+(1-1)/2, 1). - _Hartmut F. W. Hoft_, Apr 14 2016
		

References

  • R. N. Cahn, Semi-Simple Lie Algebras and Their Representations, Dover, NY, 2006, ISBN 0-486-44999-8, p. 139.

Crossrefs

Cf. A000330 (row sums), A004736, A057059, A070543.
A144112, A051340, A141419, A185874, A185875, A185876 are accumulation chain related.
A141418 is a variant.
Cf. A001227, A209260. - Hartmut F. W. Hoft, Apr 14 2016
A109814, A212652, A270877, A286013 relate to where each natural number appears in this sequence.
A000027, A000217, A000326, A005449, A045943, A059845, A115067, A140090, A140091, A140672 are rows, columns or diagonals - refer to comments.

Programs

  • Haskell
    a141419 n k =  k * (2 * n - k + 1) `div` 2
    a141419_row n = a141419_tabl !! (n-1)
    a141419_tabl = map (scanl1 (+)) a004736_tabl
    -- Reinhard Zumkeller, Aug 04 2014
  • Maple
    a:=(n,k)->k*n-binomial(k,2): seq(seq(a(n,k),k=1..n),n=1..12); # Muniru A Asiru, Oct 14 2018
  • Mathematica
    T[n_, m_] = m*(2*n - m + 1)/2; a = Table[Table[T[n, m], {m, 1, n}], {n, 1, 10}]; Flatten[a]

Formula

t(n,m) = m*(2*n - m + 1)/2.
t(n,m) = A000217(n) - A000217(n-m). - L. Edson Jeffery, Jan 16 2013
Let v = d*h with h odd be an integer factorization, then v = t(d+(h-1)/2, h) if h+1 <= 2*d, and v = t(d+(h-1)/2, 2*d) if h+1 > 2*d; see A209260. - Hartmut F. W. Hoft, Apr 14 2016
G.f.: y*(-x + y)/((-1 + x)^2*(-1 + y)^3). - Stefano Spezia, Oct 14 2018
T(n, 2) = A060747(n) for n > 1. T(n, 3) = A008585(n - 1) for n > 2. T(n, 4) = A016825(n - 2) for n > 3. T(n, 5) = A008587(n - 2) for n > 4. T(n, 6) = A016945(n - 3) for n > 5. T(n, 7) = A008589(n - 3) for n > 6. T(n, 8) = A017113(n - 4) for n > 7.r n > 5. T(n, 7) = A008589(n - 3) for n > 6. T(n, 8) = A017113(n - 4) for n > 7. T(n, 9) = A008591(n - 4) for n > 8. T(n, 10) = A017329(n - 5) for n > 9. T(n, 11) = A008593(n - 5) for n > 10. T(n, 12) = A017593(n - 6) for n > 11. T(n, 13) = A008595(n - 6) for n > 12. T(n, 14) = A147587(n - 7) for n > 13. T(n, 15) = A008597(n - 7) for n > 14. T(n, 16) = A051062(n - 8) for n > 15. T(n, 17) = A008599(n - 8) for n > 16. - Stefano Spezia, Oct 14 2018
T(2*n-k, k) = A070543(n, k). - Peter Munn, Aug 21 2019

Extensions

Simpler name by Stefano Spezia, Oct 14 2018

A017641 a(n) = 12*n + 10.

Original entry on oeis.org

10, 22, 34, 46, 58, 70, 82, 94, 106, 118, 130, 142, 154, 166, 178, 190, 202, 214, 226, 238, 250, 262, 274, 286, 298, 310, 322, 334, 346, 358, 370, 382, 394, 406, 418, 430, 442, 454, 466, 478, 490, 502, 514, 526, 538, 550, 562, 574, 586, 598, 610, 622, 634
Offset: 0

Views

Author

Keywords

Comments

Exponents e such that x^e + x^2 - 1 is reducible.
If Y is a 4-subset of an (2n+1)-set X then, for n>=3, a(n-2) is the number of 3-subsets of X having at least two elements in common with Y. - Milan Janjic, Dec 16 2007

Crossrefs

Programs

Formula

A030132(a(n)) = 9. - Reinhard Zumkeller, Jul 04 2007
G.f.: 2*(5 + x)/(1 - x)^2. - Stefano Spezia, May 09 2021
Sum_{n>=0} (-1)^n/a(n) = Pi/12 - sqrt(3)*log(2 + sqrt(3))/12. - Amiram Eldar, Dec 12 2021
From Elmo R. Oliveira, Apr 04 2025: (Start)
E.g.f.: 2*exp(x)*(5 + 6*x).
a(n) = 2*A016969(n).
a(n) = 2*a(n-1) - a(n-2) for n >= 2. (End)

A140811 a(n) = 6*n^2 - 1.

Original entry on oeis.org

-1, 5, 23, 53, 95, 149, 215, 293, 383, 485, 599, 725, 863, 1013, 1175, 1349, 1535, 1733, 1943, 2165, 2399, 2645, 2903, 3173, 3455, 3749, 4055, 4373, 4703, 5045, 5399, 5765, 6143, 6533, 6935, 7349, 7775, 8213, 8663, 9125, 9599, 10085, 10583, 11093, 11615
Offset: 0

Views

Author

Paul Curtz, Jul 16 2008

Keywords

Comments

Also: The numerators in the j=2 column of the array a(i,j) defined in A140825, where the columns j=0 and j=1 are represented by A000012 and A005408. This could be extended to column j=3: 1, -1, 9, 55, 161, ... The common feature of these sequences derived from a(i,j) is that their j-th differences are constant sequences defined by A091137(j).
a(n) is the set of all k such that 6*k + 6 is a perfect square. - Gary Detlefs, Mar 04 2010
The identity (6*n^2 - 1)^2 - (9*n^2 - 3)*(2*n)^2 = 1 can be written as a(n+1)^2 - A157872(n)*A005843(n+1)^2 = 1. - Vincenzo Librandi, Feb 05 2012
Apart from first term, sequence found by reading the line from 5, in the direction 5, 23, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. - Omar E. Pol, Jul 18 2012
From Paul Curtz, Sep 17 2018: (Start)
Terms from center to right in the following spiral:
.
65--63--61--59
/ \
67 31--29--27 57
/ / \ \
69 33 9---7 25 55
/ / / \ \ \
71 35 11 -1===5==23==53==>
/ / / / / /
37 13 1---3 21 51
\ \ / /
39 15--17--19 49
\ /
41--43--45--47 (End)

References

  • P. Curtz, Intégration numérique des systèmes différentiels à conditions initiales, Note 12, Centre de Calcul Scientifique de l'Armement, Arcueil, 1969, 132 pages, pp. 28-36. CCSA, then CELAR. Now DGA Maitrise de l'Information 35131 Bruz.

Crossrefs

Programs

Formula

a(n) = 2*a(n-1) - a(n-2) + 12.
First differences: a(n+1) - a(n) = A017593(n).
Second differences: A071593(n+1) - A071593(n) = 12.
G.f.: (1-8*x-5*x^2)/(x-1)^3. - Jaume Oliver Lafont, Aug 30 2009
From Vincenzo Librandi, Feb 05 2012: (Start)
a(n) = a(n-1) + 12*n - 6.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). (End)
a(n) = A033581(n) - 1. - Omar E. Pol, Jul 18 2012
a(n) = A032528(2*n) - 1. - Adriano Caroli, Jul 21 2013
For n > 0, a(n) = floor(3/(cosh(1/n) - 1)) = floor(1/(n*sinh(1/n) - 1)); for similar formulas for cosine and sine, see A033581. - Clark Kimberling, Oct 19 2014, corrected by M. F. Hasler, Oct 21 2014
a(-n) = a(n). - Paul Curtz, Sep 17 2018
From Amiram Eldar, Feb 04 2021: (Start)
Sum_{n>=1} 1/a(n) = (1 - (Pi/sqrt(6))*cot(Pi/sqrt(6)))/2.
Sum_{n>=1} (-1)^(n+1)/a(n) = ((Pi/sqrt(6))*csc(Pi/sqrt(6)) - 1)/2.
Product_{n>=1} (1 + 1/a(n)) = (Pi/sqrt(6))*csc(Pi/sqrt(6)).
Product_{n>=1} (1 - 1/a(n)) = csc(Pi/sqrt(6))*sin(Pi/sqrt(3))/sqrt(2). (End)
a(n) = A003154(n+1) - 2*A016777(n). - Leo Tavares, May 13 2022
E.g.f.: exp(x)*(6*x^2 + 6*x - 1). - Elmo R. Oliveira, Jan 16 2025

Extensions

Edited and extended by R. J. Mathar, Aug 06 2008
Better description Ray Chandler, Feb 03 2009

A055926 Numbers k such that {largest m such that 1, 2, ..., m divide k} is different from {largest m such that m! divides k}; numbers k which are either odd multiples of 12 or the largest m such that (m-1)! divides k is a composite number > 5.

Original entry on oeis.org

12, 36, 60, 84, 108, 120, 132, 156, 180, 204, 228, 240, 252, 276, 300, 324, 348, 360, 372, 396, 420, 444, 468, 480, 492, 516, 540, 564, 588, 600, 612, 636, 660, 684, 708, 732, 756, 780, 804, 828, 840, 852, 876, 900, 924, 948, 960, 972, 996, 1020, 1044, 1068
Offset: 1

Views

Author

Leroy Quet, Jul 16 2000

Keywords

Comments

From Antti Karttunen, Nov 20 - Dec 06 2013: (Start)
This sequence has several interpretations:
Numbers k such that A055874(k) differs from A055881(k). [Leroy Quet's original definition of the sequence. Note that A055874(k) >= A055881(k) for all k.]
Numbers k such that {largest m such that m! divides k^2} is different from {largest m such that m! divides k}, i.e., numbers k for which A232098(k) > A055881(k).
Numbers k which are either 12 times an odd number (A073762) or the largest m such that (m-1)! divides k is a composite number > 5 (A232743).
Please see my attached notes for the proof of the equivalence of these interpretations.
Additional implications based on that proof:
A232099 is a subset of this sequence.
A055881(a(n))+1 is always composite. In the range n = 1..17712, only values 4, 6, 8, 9 and 10 occur.
The new definition can be also rephrased by saying that the sequence contains all the positive integers k whose factorial base representation of (A007623(k)) either ends as '...200' (in which case k is an odd multiple of 12, 12 = '200', 36 = '1200', 60 = '2200', ...) or the number of trailing zeros + 2 in that representation is a composite number greater than or equal to 6, e.g. 120 = '10000' (in other words, A055881(k) is one of the terms of A072668 after the initial 3). Together these conditions also imply that all the terms are divisible by 12.
(End)

Examples

			12 is included because 3! is the largest factorial to divide 12, but 1, 2, 3 and 4 all divide 12. Equally, 12 is included because it is one of the terms of A073762, or equally, because its factorial base representation ends with digits '...200': A007623(12) = 200.
840 (= 3*5*7*8) is included because the largest factorial which divides 840 is 5! (840 = 7*120), but all positive integers up to 8 divide 840. Equally, 840 is included because it is one of the terms of A232743 as 5+1 = 6 is a composite number larger than 5. Note that A007623(840) = 110000.
		

Crossrefs

Union of A073762 and A232743. Equivalently, setwise difference of A232742 and A017593. Subset: A232099.

Extensions

More terms from Antti Karttunen, Dec 01 2013

A017581 a(n) = 12*n + 5.

Original entry on oeis.org

5, 17, 29, 41, 53, 65, 77, 89, 101, 113, 125, 137, 149, 161, 173, 185, 197, 209, 221, 233, 245, 257, 269, 281, 293, 305, 317, 329, 341, 353, 365, 377, 389, 401, 413, 425, 437, 449, 461, 473, 485, 497, 509, 521, 533, 545, 557, 569, 581, 593, 605, 617, 629
Offset: 0

Views

Author

Keywords

Comments

Apart from initial term(s), dimension of the space of weight 2n cusp forms for Gamma_0(71).
A089911(2*a(n)) = 7. - Reinhard Zumkeller, Jul 05 2013
Equivalently, intersection of A016813 and A016789. - Bruno Berselli, Jan 24 2018

Crossrefs

Programs

Formula

a(n) = 2*a(n-1) - a(n-2) for n>1, a(0)=5, a(1)=17. - Vincenzo Librandi, Jun 08 2011
G.f.: x*(5 + 7*x)/(1 - x)^2. - Wolfdieter Lang, Jul 04 2023
E.g.f.: exp(x)*(5 + 12*x). - Stefano Spezia, Feb 21 2024
a(n) = A016969(2*n) = A016789(4*n+1). - Elmo R. Oliveira, Apr 10 2025

A141310 The odd numbers interlaced with the constant-2 sequence.

Original entry on oeis.org

1, 2, 3, 2, 5, 2, 7, 2, 9, 2, 11, 2, 13, 2, 15, 2, 17, 2, 19, 2, 21, 2, 23, 2, 25, 2, 27, 2, 29, 2, 31, 2, 33, 2, 35, 2, 37, 2, 39, 2, 41, 2, 43, 2, 45, 2, 47, 2, 49, 2, 51, 2, 53, 2, 55, 2, 57, 2, 59, 2, 61, 2, 63, 2, 65, 2, 67, 2, 69, 2, 71, 2, 73, 2, 75, 2, 77, 2, 79, 2, 81, 2, 83, 2, 85, 2, 87, 2, 89, 2, 91, 2, 93, 2, 95, 2, 97
Offset: 0

Views

Author

Paul Curtz, Aug 02 2008

Keywords

Comments

Similarly, the principle of interlacing a sequence and its first differences leads from A000012 and its differences A000004 to A059841, or from A140811 and its first differences A017593 to a sequence -1, 6, 5, 18, ...
If n is even then a(n) = n + 1 ; otherwise a(n) = 2. - Wesley Ivan Hurt, Jun 05 2013
Denominators of floor((n+1)/2) / (n+1), n > 0. - Wesley Ivan Hurt, Jun 14 2013
a(n) is also the number of minimum total dominating sets in the (n+1)-gear graph for n>1. - Eric W. Weisstein, Apr 11 2018
a(n) is also the number of minimum total dominating sets in the (n+1)-sun graph for n>1. - Eric W. Weisstein, Sep 09 2021
Denominators of Cesàro means sequence of A114112, corresponding numerators are in A354008. - Bernard Schott, May 14 2022
Also, denominators of Cesàro means sequence of A237420, corresponding numerators are in A354280. - Bernard Schott, May 22 2022

Crossrefs

Programs

  • Maple
    a:= n-> n+1-(n-1)*(n mod 2): seq(a(n), n=0..96); # Wesley Ivan Hurt, Jun 05 2013
  • Mathematica
    Flatten[Table[{2 n - 1, 2}, {n, 40}]] (* Alonso del Arte, Jun 15 2013 *)
    Riffle[Range[1, 79, 2], 2] (* Alonso del Arte, Jun 14 2013 *)
    Table[((-1)^n (n - 1) + n + 3)/2, {n, 0, 20}] (* Eric W. Weisstein, Apr 11 2018 *)
    Table[Floor[(n + 1)/2]/(n + 1), {n, 0, 20}] // Denominator (* Eric W. Weisstein, Apr 11 2018 *)
    LinearRecurrence[{0, 2, 0, -1}, {2, 3, 2, 5}, {0, 20}] (* Eric W. Weisstein, Apr 11 2018 *)
    CoefficientList[Series[(1 + 2 x + x^2 - 2 x^3)/(-1 + x^2)^2, {x, 0, 20}], x] (* Eric W. Weisstein, Apr 11 2018 *)
  • PARI
    A141310(n) = if(n%2,2,1+n); \\ (for offset=0 version) - Antti Karttunen, Oct 02 2018
    
  • PARI
    A141310off1(n) = if(n%2,n,2); \\ (for offset=1 version) - Antti Karttunen, Oct 02 2018
    
  • Python
    def A141310(n): return 2 if n % 2 else n + 1 # Chai Wah Wu, May 24 2022

Formula

a(2n) = A005408(n). a(2n+1) = 2.
First differences: a(n+1) - a(n) = (-1)^(n+1)*A109613(n-1), n > 0.
b(2n) = -A008586(n), and b(2n+1) = A060747(n), where b(n) = a(n+1) - 2*a(n).
a(n) = 2*a(n-2) - a(n-4). - R. J. Mathar, Feb 23 2009
G.f.: (1+2*x+x^2-2*x^3)/((x-1)^2*(1+x)^2). - R. J. Mathar, Feb 23 2009
From Wesley Ivan Hurt, Jun 05 2013: (Start)
a(n) = n + 1 - (n - 1)*(n mod 2).
a(n) = (n + 1) * (n - floor((n+1)/2))! / floor((n+1)/2)!.
a(n) = A000142(n+1) / A211374(n+1). (End)

Extensions

Edited by R. J. Mathar, Feb 23 2009
Term a(45) corrected, and more terms added by Antti Karttunen, Oct 02 2018

A030133 a(n+1) is the sum of digits of (a(n) + a(n-1)).

Original entry on oeis.org

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

Views

Author

Keywords

Comments

a(n) = A010888(A000032(n)). - Reinhard Zumkeller, Aug 20 2011
Similar to the digital roots of several Fibonacci sequences, this digital root sequence for Lucas numbers (A000032) has period 24 with digits summing to 117.
Decimal expansion of 23719213606865169775282 / 111111111111111111111111 = 0.[213472922461786527977538] (periodic). - Daniel Forgues, Feb 27 2017

Crossrefs

Programs

  • Haskell
    a030133 n = a030133_list !! n
    a030133_list =
       2 : 1 : map a007953 (zipWith (+) a030133_list $ tail a030133_list)
    -- Reinhard Zumkeller, Aug 20 2011
    
  • Mathematica
    Transpose[NestList[{Last[#],Total[IntegerDigits[Total[#]]]}&, {2,1}, 100]] [[1]] (* Harvey P. Dale, Jul 25 2011 *)
  • PARI
    V=[2,1];for(n=1,100,V=concat(V,sumdigits(V[n]+V[n+1])));V \\ Derek Orr, Feb 27 2017
    
  • PARI
    Vec((2 + x + 3*x^2 + 4*x^3 + 7*x^4 + 2*x^5 + 9*x^6 + 2*x^7 + 2*x^8 + 4*x^9 + 6*x^10 + x^11 + 7*x^12 + 8*x^13 + 6*x^14 + 5*x^15 + 2*x^16 + 7*x^17 + 9*x^18 + 7*x^19 + 7*x^20 + 5*x^21 + 3*x^22 + 8*x^23) / (1 - x^24)  + O(x^80)) \\ Colin Barker, Sep 25 2019

Formula

a(n+24) = a(n); a(A017593(n)) = 9. - Reinhard Zumkeller, Jul 04 2007
G.f.: (2 + x + 3*x^2 + 4*x^3 + 7*x^4 + 2*x^5 + 9*x^6 + 2*x^7 + 2*x^8 + 4*x^9 + 6*x^10 + x^11 + 7*x^12 + 8*x^13 + 6*x^14 + 5*x^15 + 2*x^16 + 7*x^17 + 9*x^18 + 7*x^19 + 7*x^20 + 5*x^21 + 3*x^22 + 8*x^23) / (1 - x^24). - Colin Barker, Sep 25 2019

A332514 Numbers k such that phi(k) == 6 (mod 12), where phi is the Euler totient function (A000010).

Original entry on oeis.org

7, 9, 14, 18, 19, 27, 31, 38, 43, 49, 54, 62, 67, 79, 81, 86, 98, 103, 127, 134, 139, 151, 158, 162, 163, 199, 206, 211, 223, 243, 254, 271, 278, 283, 302, 307, 326, 331, 343, 361, 367, 379, 398, 422, 439, 446, 463, 486, 487, 499, 523, 542, 547, 566, 571, 607, 614
Offset: 1

Views

Author

Amiram Eldar, Feb 14 2020

Keywords

Comments

Dence and Pomerance showed that the asymptotic number of the terms below x is ~ (3/8) * x/log(x).

Examples

			19 is a term since phi(19) = 18 == 6 (mod 12).
		

Crossrefs

Programs

  • Magma
    [k:k in [1..650]| EulerPhi(k) mod 12 eq 6]; // Marius A. Burtea, Feb 14 2020
  • Mathematica
    Select[Range[600], Mod[EulerPhi[#], 12] == 6 &]
Showing 1-10 of 20 results. Next