A047202 -id:A047202 - cp's OEIS frontend

A191732 Dispersion of A047202, (numbers >1 and congruent to 2 or 3 or 4 mod 5), by antidiagonals.

1, 2, 5, 3, 8, 6, 4, 13, 9, 10, 7, 22, 14, 17, 11, 12, 37, 23, 28, 18, 15, 19, 62, 38, 47, 29, 24, 16, 32, 103, 63, 78, 48, 39, 27, 20, 53, 172, 104, 129, 79, 64, 44, 33, 21, 88, 287, 173, 214, 132, 107, 73, 54, 34, 25, 147, 478, 288, 357, 219, 178, 122, 89

Offset: 1

Clark Kimberling, Jun 14 2011

For a background discussion of dispersions and their fractal sequences, see A191426.  For dispersions of congruence sequences mod 3, mod 4, or mod 5, see A191655, A191663, A191667, A191702.
...
Suppose that {2,3,4,5,6} is partitioned as {x1, x2} and {x3,x4,x5}.  Let S be the increasing sequence of numbers >1 and congruent to x1 or x2 mod 5, and let T be the increasing sequence of numbers >1 and congruent to x3 or x4 or x5 mod 5.  There are 10 sequences in S, each matched by a (nearly) complementary sequence in T.  Each of the 20 sequences generates a dispersion, as listed here:
...
A191722=dispersion of A008851 (0, 1 mod 5 and >1)
A191723=dispersion of A047215 (0, 2 mod 5 and >1)
A191724=dispersion of A047218 (0, 3 mod 5 and >1)
A191725=dispersion of A047208 (0, 4 mod 5 and >1)
A191726=dispersion of A047216 (1, 2 mod 5 and >1)
A191727=dispersion of A047219 (1, 3 mod 5 and >1)
A191728=dispersion of A047209 (1, 4 mod 5 and >1)
A191729=dispersion of A047221 (2, 3 mod 5 and >1)
A191730=dispersion of A047211 (2, 4 mod 5 and >1)
A191731=dispersion of A047204 (3, 4 mod 5 and >1)
...
A191732=dispersion of A047202 (2,3,4 mod 5 and >1)
A191733=dispersion of A047206 (1,3,4 mod 5 and >1)
A191734=dispersion of A032793 (1,2,4 mod 5 and >1)
A191735=dispersion of A047223 (1,2,3 mod 5 and >1)
A191736=dispersion of A047205 (0,3,4 mod 5 and >1)
A191737=dispersion of A047212 (0,2,4 mod 5 and >1)
A191738=dispersion of A047222 (0,2,3 mod 5 and >1)
A191739=dispersion of A008854 (0,1,4 mod 5 and >1)
A191740=dispersion of A047220 (0,1,3 mod 5 and >1)
A191741=dispersion of A047217 (0,1,2 mod 5 and >1)
...
For further information about these 20 dispersions, see A191722.
...
Regarding the dispersions A191722-A191741, there are general formulas for sequences of the type "(a or b mod m)" and "(a or b or c mod m)" used in the relevant Mathematica programs.

			Northwest corner:
1....2....3....4....7
5....8....13...22...37
6....9....14...23...38
10...17...28...47...78
11...18...29...48...79
15...24...39...64...107

Ivan Neretin, Table of n, a(n) for n = 1..5050

Cf. A008851, A047202, A191722, A191426.

(* Program generates the dispersion array t of the increasing sequence f[n] *)
r = 40; r1 = 12;  c = 40; c1 = 12;
a=2; b=3; c2=4; m[n_]:=If[Mod[n,3]==0,1,0];
f[n_]:=a*m[n+2]+b*m[n+1]+c2*m[n]+5*Floor[(n-1)/3]
Table[f[n], {n, 1, 30}]  (* A047202 *)
mex[list_] := NestWhile[#1 + 1 &, 1, Union[list][[#1]] <= #1 &, 1, Length[Union[list]]]
rows = {NestList[f, 1, c]};
Do[rows = Append[rows, NestList[f, mex[Flatten[rows]], r]], {r}];
t[i_, j_] := rows[[i, j]];
TableForm[Table[t[i, j], {i, 1, 10}, {j, 1, 10}]] (* A191732 *)
Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191732  *)

A191708 Dispersion of A047202, (numbers >1 and congruent to 0, 2, 3, or 4 mod 5), by antidiagonals.

Original entry on oeis.org

1, 2, 6, 3, 8, 11, 4, 10, 14, 16, 5, 13, 18, 20, 21, 7, 17, 23, 25, 27, 26, 9, 22, 29, 32, 34, 33, 31, 12, 28, 37, 40, 43, 42, 39, 36, 15, 35, 47, 50, 54, 53, 49, 45, 41, 19, 44, 59, 63, 68, 67, 62, 57, 52, 46, 24, 55, 74, 79, 85, 84, 78, 72, 65, 58, 51, 30

Offset: 1

Clark Kimberling, Jun 12 2011

For a background discussion of dispersions and their fractal sequences, see A191426.  For dispersions of congruence sequences mod 3 or mod 4, see A191655, A191663, A191667.
...
Each of the sequences (5n, n>1), (5n+1, n>1), (5n+2, n>=0), (5n+3, n>=0), (5n+4, n>=0), generates a dispersion.  Each complement (beginning with its first term >1) also generates a dispersion.  The ten sequences and dispersions are listed here:
...
A191702=dispersion of A008587 (5k, k>=1)
A191703=dispersion of A016861 (5k+1, k>=1)
A191704=dispersion of A016873 (5k+2, k>=0)
A191705=dispersion of A016885 (5k+3, k>=0)
A191706=dispersion of A016897 (5k+4, k>=0)
A191707=dispersion of A047201 (1, 2, 3, 4 mod 5 and >1)
A191708=dispersion of A047202 (0, 2, 3, 4 mod 5 and >1)
A191709=dispersion of A047207 (0, 1, 3, 4 mod 5 and >1)
A191710=dispersion of A032763 (0, 1, 2, 4 mod 5 and >1)
A191711=dispersion of A001068 (0, 1, 2, 3 mod 5 and >1)
...
EXCEPT for at most 2 initial terms (so that column 1 always starts with 1):
A191702 has 1st col A047201, all else A008587
A191703 has 1st col A047202, all else A016861
A191704 has 1st col A047207, all else A016873
A191705 has 1st col A032763, all else A016885
A191706 has 1st col A001068, all else A016897
A191707 has 1st col A008587, all else A047201
A191708 has 1st col A042968, all else A047203
A191709 has 1st col A042968, all else A047207
A191710 has 1st col A042968, all else A032763
A191711 has 1st col A042968, all else A001068
...
Regarding the dispersions A191670-A191673, there is a formula for sequences of the type "(a or b or c or d mod m)", (as in the relevant Mathematica programs):
...
If f(n)=(n mod 3), then (a,b,c,d,a,b,c,d,a,b,c,d,...) is given by a*f(n+3)+b*f(n+2)+c*f(n+1)+d*f(n); so that for n>=1, "(a, b, c, d mod m)" is given by
a*f(n+3)+b*f(n+2)+c*f(n+1)+d*f(n)+m*floor((n-1)/4)).

			Northwest corner:
1....2....3....4....5
6....8....10...13...17
11...14...18...23...29
16...20...25...32...40
21...27...34...43...54
26...33...42...53...67

Ivan Neretin, Table of n, a(n) for n = 1..5050 (first 100 antidiagonals, flattened)

Cf. A047202, A016861, A191703, A191426.

(* Program generates the dispersion array T of the increasing sequence f[n] *)
r = 40; r1 = 12;  c = 40; c1 = 12;
a=2; b=3; c2=4; d=5; m[n_]:=If[Mod[n,4]==0,1,0];
f[n_]:=a*m[n+3]+b*m[n+2]+c2*m[n+1]+d*m[n]+5*Floor[(n-1)/4]
Table[f[n], {n, 1, 30}]  (* A047202 *)
mex[list_] := NestWhile[#1 + 1 &, 1, Union[list][[#1]] <= #1 &, 1, Length[Union[list]]]
rows = {NestList[f, 1, c]};
Do[rows = Append[rows, NestList[f, mex[Flatten[rows]], r]], {r}];
t[i_, j_] := rows[[i, j]];
TableForm[Table[t[i, j], {i, 1, 10}, {j, 1, 10}]] (* A191708 *)
Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191708  *)

A191709 Dispersion of A047202, (numbers >1 and congruent to 0, 1, 3, or 4 mod 5), by antidiagonals.

Original entry on oeis.org

1, 3, 2, 5, 4, 7, 8, 6, 10, 12, 11, 9, 14, 16, 17, 15, 13, 19, 21, 23, 22, 20, 18, 25, 28, 30, 29, 27, 26, 24, 33, 36, 39, 38, 35, 32, 34, 31, 43, 46, 50, 49, 45, 41, 37, 44, 40, 55, 59, 64, 63, 58, 53, 48, 42, 56, 51, 70, 75, 81, 80, 74, 68, 61, 54, 47, 71

Offset: 1

Clark Kimberling, Jun 12 2011

For a background discussion of dispersions and their fractal sequences, see A191426.  For dispersions of congruence sequences mod 3 or mod 4, see A191655, A191663, A191667.
...
Each of the sequences (5n, n>1), (5n+1, n>1), (5n+2, n>=0), (5n+3, n>=0), (5n+4, n>=0), generates a dispersion.  Each complement (beginning with its first term >1) also generates a dispersion.  The ten sequences and dispersions are listed here:
...
A191702=dispersion of A008587 (5k, k>=1)
A191703=dispersion of A016861 (5k+1, k>=1)
A191704=dispersion of A016873 (5k+2, k>=0)
A191705=dispersion of A016885 (5k+3, k>=0)
A191706=dispersion of A016897 (5k+4, k>=0)
A191707=dispersion of A047201 (1, 2, 3, 4 mod 5 and >1)
A191708=dispersion of A047202 (0, 2, 3, 4 mod 5 and >1)
A191709=dispersion of A047207 (0, 1, 3, 4 mod 5 and >1)
A191710=dispersion of A032763 (0, 1, 2, 4 mod 5 and >1)
A191711=dispersion of A001068 (0, 1, 2, 3 mod 5 and >1)
...
EXCEPT for at most 2 initial terms (so that column 1 always starts with 1):
A191702 has 1st col A047201, all else A008587
A191703 has 1st col A047202, all else A016861
A191704 has 1st col A047207, all else A016873
A191705 has 1st col A032763, all else A016885
A191706 has 1st col A001068, all else A016897
A191707 has 1st col A008587, all else A047201
A191708 has 1st col A042968, all else A047203
A191709 has 1st col A042968, all else A047207
A191710 has 1st col A042968, all else A032763
A191711 has 1st col A042968, all else A001068
...
Regarding the dispersions A191670-A191673, there is a formula for sequences of the type "(a or b or c or d mod m)", (as in the relevant Mathematica programs):
...
If f(n)=(n mod 3), then (a,b,c,d,a,b,c,d,a,b,c,d,...) is given by a*f(n+3)+b*f(n+2)+c*f(n+1)+d*f(n); so that for n>=1, "(a, b, c, d mod m)" is given by
a*f(n+3)+b*f(n+2)+c*f(n+1)+d*f(n)+m*floor((n-1)/4)).

			Northwest corner:
1....3....5....8....11
2....4....6....9....13
7....10...14...19...25
12...16...21...28...36
17...23...30...39...50
22...29...38...49...63

Ivan Neretin, Table of n, a(n) for n = 1..5050 (first 100 antidiagonals, flattened)

Cf. A047207, A016873, A191704, A191426.

(* Program generates the dispersion array T of the increasing sequence f[n] *)
r = 40; r1 = 12;  c = 40; c1 = 12;
a=3; b=4; c2=5; d=6; m[n_]:=If[Mod[n,4]==0,1,0];
f[n_]:=a*m[n+3]+b*m[n+2]+c2*m[n+1]+d*m[n]+5*Floor[(n-1)/4]
Table[f[n], {n, 1, 30}]  (* A047207 *)
mex[list_] := NestWhile[#1 + 1 &, 1, Union[list][[#1]] <= #1 &, 1, Length[Union[list]]]
rows = {NestList[f, 1, c]};
Do[rows = Append[rows, NestList[f, mex[Flatten[rows]], r]], {r}];
t[i_, j_] := rows[[i, j]];
TableForm[Table[t[i, j], {i, 1, 10}, {j, 1, 10}]] (* A191709 *)
Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191709  *)

A016813 a(n) = 4*n + 1.

Original entry on oeis.org

1, 5, 9, 13, 17, 21, 25, 29, 33, 37, 41, 45, 49, 53, 57, 61, 65, 69, 73, 77, 81, 85, 89, 93, 97, 101, 105, 109, 113, 117, 121, 125, 129, 133, 137, 141, 145, 149, 153, 157, 161, 165, 169, 173, 177, 181, 185, 189, 193, 197, 201, 205, 209, 213, 217, 221, 225, 229, 233, 237

Offset: 0

N. J. A. Sloane

Apart from initial term(s), dimension of the space of weight 2n cusp forms for Gamma_0( 23 ).
Apart from initial term(s), dimension of the space of weight 2n cuspidal newforms for Gamma_0( 64 ).
Numbers k such that k and (k+1) have the same binary digital sum. - Benoit Cloitre, Jun 05 2002
Numbers k such that (1 + sqrt(k))/2 is an algebraic integer. - Alonso del Arte, Jun 04 2012
Numbers k such that 2 is the only prime p that satisfies the relationship p XOR k = p + k. - Brad Clardy, Jul 22 2012
This may also be interpreted as the array T(n,k) = A001844(n+k) + A008586(k) read by antidiagonals:
    1,   9,  21,  37,  57,  81, ...
    5,  17,  33,  53,  77, 105, ...
   13,  29,  49,  73, 101, 133, ...
   25,  45,  69,  97, 129, 165, ...
   41,  65,  93, 125, 161, 201, ...
   61,  89, 121, 157, 197, 241, ...
   ...
- R. J. Mathar, Jul 10 2013
With leading term 2 instead of 1, 1/a(n) is the largest tolerance of form 1/k, where k is a positive integer, so that the nearest integer to (n - 1/k)^2 and to (n + 1/k)^2 is n^2. In other words, if interval arithmetic is used to square [n - 1/k, n + 1/k], every value in the resulting interval of length 4n/k rounds to n^2 if and only if k >= a(n). - Rick L. Shepherd, Jan 20 2014
Odd numbers for which the number of prime factors congruent to 3 (mod 4) is even. - Daniel Forgues, Sep 20 2014
For the Collatz conjecture, we identify two types of odd numbers. This sequence contains all the descenders: where (3*a(n) + 1) / 2 is even and requires additional divisions by 2. See A004767 for the ascenders. - Fred Daniel Kline, Nov 29 2014 [corrected by Jaroslav Krizek, Jul 29 2016]
a(n-1), n >= 1, is also the complex dimension of the manifold M(S), the set of all conjugacy classes of irreducible representations of the fundamental group pi_1(X,x_0) of rank 2, where S = {a_1, ..., a_{n}, a_{n+1} = oo}, a subset of P^1 = C U {oo}, X = X(S) = P^1 \ S, and x_0 a base point in X. See the Iwasaki et al. reference, Proposition 2.1.4. p. 150. - Wolfdieter Lang, Apr 22 2016
For n > 3, also the number of (not necessarily maximal) cliques in the n-sunlet graph. - Eric W. Weisstein, Nov 29 2017
For integers k with absolute value in A047202, also exponents of the powers of k having the same unit digit of k in base 10. - Stefano Spezia, Feb 23 2021
Starting with a(1) = 5, numbers ending with 01 in base 2. - John Keith, May 09 2022

			From _Leo Tavares_, Jul 02 2021: (Start)
Illustration of initial terms:
                                        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
(End)

K. Iwasaki, H. Kimura, S. Shimomura and M. Yoshida, From Gauss to Painlevé, Vieweg, 1991. p. 150.

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Colin Defant and Noah Kravitz, Loops and Regions in Hitomezashi Patterns, arXiv:2201.03461 [math.CO], 2022. Theorem 1.3.
Gennady Eremin, Partitioning the set of natural numbers into Mersenne trees and into arithmetic progressions; Natural Matrix and Linnik's constant, arXiv:2405.16143 [math.CO], 2024. See pp. 5, 8, 14.
Gennady Eremin, Infinite matrix of odd natural numbers. A bit about Sophie Germain prime numbers, arXiv:2501.17090 [math.GM], 2025. See pp. 3, 11.
L. B. W. Jolley, Summation of Series, Dover Publications, 1961, p. 16.
Tanya Khovanova, Recursive Sequences
Konrad Knopp, Theorie und Anwendung der unendlichen Reihen, Berlin, J. Springer, 1922. (Original German edition of "Theory and Application of Infinite Series")
Konrad Knopp, Theory and Application of Infinite Series, Dover, p. 269.
Luis Manuel Rivera, Integer sequences and k-commuting permutations, arXiv preprint arXiv:1406.3081 [math.CO], 2014-2015.
William A. Stein's The Modular Forms Database, PARI-readable dimension tables for Gamma_0(N)
Eric Weisstein's World of Mathematics, Clique
Eric Weisstein's World of Mathematics, Hilbert Number
Eric Weisstein's World of Mathematics, Sunlet Graph
Wikipedia, Interval arithmetic
Index entries for linear recurrences with constant coefficients, signature (2,-1).

Subsequence of A042963 and of A079523.
a(n) = A093561(n+1, 1), (4, 1)-Pascal column.
Cf. A016921, A017281, A017533, A047202, A158057, A161705, A161709, A161714, A128470. - Reinhard Zumkeller, Jun 17 2009
Cf. A004772 (complement).
Cf. A017557.

List([0..70],n->4*n+1); # Muniru A Asiru, Aug 08 2018

a016813 = (+ 1) . (* 4)
a016813_list = [1, 5 ..]  -- Reinhard Zumkeller, Feb 14 2012

Magma
```
[n: n in [1..250 by 4]];
```

seq(4*k+1, k=0..100); # Wesley Ivan Hurt, Sep 28 2013

Range[1, 237, 4] (* Vladimir Joseph Stephan Orlovsky, May 26 2011 *)
Table[4 n + 1, {n, 0, 20}] (* Eric W. Weisstein, Nov 29 2017 *)
4 Range[0, 20] + 1 (* Eric W. Weisstein, Nov 29 2017 *)
LinearRecurrence[{2, -1}, {5, 9}, {0, 20}] (* Eric W. Weisstein, Nov 29 2017 *)
CoefficientList[Series[(1 + 3 x)/(-1 + x)^2, {x, 0, 20}], x] (* Eric W. Weisstein, Nov 29 2017 *)

a(n)=4*n+1 \\ Charles R Greathouse IV, Mar 22 2013

x='x+O('x^100); Vec((1+3*x)/(1-x)^2) \\ Altug Alkan, Oct 22 2015

(0 to 59).map(4 *  + 1) // _Alonso del Arte, Aug 08 2018

a(n) = A005408(2*n).
Sum_{n>=0} (-1)^n/a(n) = (1/(4*sqrt(2)))*(Pi+2*log(sqrt(2)+1)) = A181048 [Jolley]. - Benoit Cloitre, Apr 05 2002 [corrected by Amiram Eldar, Jul 30 2023]
G.f.: (1+3*x)/(1-x)^2. - Paul Barry, Feb 27 2003 [corrected for offset 0 by Wolfdieter Lang, Oct 03 2014]
(1 + 5*x + 9*x^2 + 13*x^3 + ...) = (1 + 2*x + 3*x^2 + ...) / (1 - 3*x + 9*x^2 - 27*x^3 + ...). - Gary W. Adamson, Jul 03 2003
a(n) = A001969(n) + A000069(n). - Philippe Deléham, Feb 04 2004
a(n) = A004766(n-1). - R. J. Mathar, Oct 26 2008
a(n) = 2*a(n-1) - a(n-2); a(0)=1, a(1)=5. a(n) = 4 + a(n-1). - Philippe Deléham, Nov 03 2008
A056753(a(n)) = 3. - Reinhard Zumkeller, Aug 23 2009
A179821(a(n)) = a(A179821(n)). - Reinhard Zumkeller, Jul 31 2010
a(n) = 8*n - 2 - a(n-1) for n > 0, a(0) = 1. - Vincenzo Librandi, Nov 20 2010
The identity (4*n+1)^2 - (4*n^2+2*n)*(2)^2 = 1 can be written as a(n)^2 - A002943(n)*2^2 = 1. - Vincenzo Librandi, Mar 11 2009 - Nov 25 2012
A089911(6*a(n)) = 8. - Reinhard Zumkeller, Jul 05 2013
a(n) = A004767(n) - 2. - Jean-Bernard François, Sep 27 2013
a(n) = A058281(3n+1). - Eli Jaffe, Jun 07 2016
From Ilya Gutkovskiy, Jul 29 2016: (Start)
E.g.f.: (1 + 4*x)*exp(x).
a(n) = Sum_{k = 0..n} A123932(k).
a(A005098(k)) = x^2 + y^2.
Inverse binomial transform of A014480. (End)
Dirichlet g.f.: 4*Zeta(-1 + s) + Zeta(s). - Stefano Spezia, Nov 02 2018

A191702 Dispersion of A008587 (5,10,15,20,25,30,...), by antidiagonals.

Original entry on oeis.org

1, 5, 2, 25, 10, 3, 125, 50, 15, 4, 625, 250, 75, 20, 6, 3125, 1250, 375, 100, 30, 7, 15625, 6250, 1875, 500, 150, 35, 8, 78125, 31250, 9375, 2500, 750, 175, 40, 9, 390625, 156250, 46875, 12500, 3750, 875, 200, 45, 11, 1953125, 781250, 234375, 62500, 18750

Offset: 1

Clark Kimberling, Jun 12 2011

For a background discussion of dispersions and their fractal sequences, see A191426. For dispersions of congruence sequences mod 3 or mod 4, see A191655, A191663, A191667.
...
Each of the sequences (5n, n>1), (5n+1, n>1), (5n+2, n>=0), (5n+3, n>=0), (5n+4, n>=0), generates a dispersion. Each complement (beginning with its first term >1) also generates a dispersion. The ten sequences and dispersions are listed here:
...
A191702=dispersion of A008587 (5k, k>=1)
A191703=dispersion of A016861 (5k+1, k>=1)
A191704=dispersion of A016873 (5k+2, k>=0)
A191705=dispersion of A016885 (5k+3, k>=0)
A191706=dispersion of A016897 (5k+4, k>=0)
A191707=dispersion of A047201 (1, 2, 3, 4 mod 5 and >1)
A191708=dispersion of A047202 (0, 2, 3, 4 mod 5 and >1)
A191709=dispersion of A047207 (0, 1, 3, 4 mod 5 and >1)
A191710=dispersion of A032763 (0, 1, 2, 4 mod 5 and >1)
A191711=dispersion of A001068 (0, 1, 2, 3 mod 5 and >1)
...
EXCEPT for at most 2 initial terms (so that column 1 always starts with 1):
A191702 has 1st col A047201, all else A008587
A191703 has 1st col A047202, all else A016861
A191704 has 1st col A047207, all else A016873
A191705 has 1st col A032763, all else A016885
A191706 has 1st col A001068, all else A016897
A191707 has 1st col A008587, all else A047201
A191708 has 1st col A042968, all else A047203
A191709 has 1st col A042968, all else A047207
A191710 has 1st col A042968, all else A032763
A191711 has 1st col A042968, all else A001068
...
Regarding the dispersions A191670-A191673, there is a formula for sequences of the type
"(a or b or c or d mod m)", used in Mathematica programs for A191707-A191711):  if f(n)=(n mod 3), then
  (a,b,c,d,a,b,c,d,a,b,c,d,...) is given by
   a*f(n+3)+b*f(n+2)+c*f(n+1)+d*f(n), so that for n>=1,
"(a, b, c, d mod m)" is given by
   a*f(n+3)+b*f(n+2)+c*f(n+1)+d*f(n)+m*floor((n-1)/4)).

			Northwest corner:
  1...5....25....125...625
  2...10...50....250...1250
  3...15...75....375...1875
  4...20...100...500...2500
  6...30...150...750...3750

Ivan Neretin, Table of n, a(n) for n = 1..5050 (first 100 antidiagonals, flattened)

Cf. A047201, A008587, A191707, A191426.

(* Program generates the dispersion array T of the increasing sequence f[n] *)
r = 40; r1 = 12;  c = 40; c1 = 12;
f[n_] := 5n
Table[f[n], {n, 1, 30}]  (* A008587 *)
mex[list_] := NestWhile[#1 + 1 &, 1, Union[list][[#1]] <= #1 &, 1, Length[Union[list]]]
rows = {NestList[f, 1, c]};
Do[rows = Append[rows, NestList[f, mex[Flatten[rows]], r]], {r}];
t[i_, j_] := rows[[i, j]];
TableForm[Table[t[i, j], {i, 1, 10}, {j, 1, 10}]]
(* A191702 *)
Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191702  *)

T(i,j) = T(i,1)*T(1,j) = (i-1+floor((i+3)/4))*5^(j-1), i>=1, j>=1.

A010685 Period 2: repeat (1,4).

Original entry on oeis.org

1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1, 4, 1

Offset: 0

N. J. A. Sloane

Continued fraction of (1 + sqrt(2))/2. - R. J. Mathar, Nov 21 2011

This sequence can be generated by an infinite number of formulas all having the form a^(b*n) mod c subject to the following conditions. The number a can be congruent to either 2,3, or 4 mod 5 (A047202). If a is congruent to 2 or 3 mod 5, then b can be any number of the form 4k+2 and c = 5 or 15. If a is congruent to 4 mod 5, then b can be any number of the form 2k+1 and c = 5. For example: a(n) = 29^(13*n) mod 5, a(n) = 24^(11*n) mod 5, and a(n) = 22^(10*n) mod 15. - Gary Detlefs, May 19 2014

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

Cf. sequences listed in Comments section of A283393.

Cf. A047202.

[Modexp(4,n,5): n in [0..100]]; // G. C. Greubel, Nov 22 2021

A010685 := proc(n)
    if type(n,'even') then
        1 ;
    else
        4;
    end if;
end proc: # R. J. Mathar, Aug 03 2015

Table[(5-3(-1)^n)/2, {n, 0, 100}] (* Wesley Ivan Hurt, Mar 26 2014 *)
PadRight[{},120,{1,4}] (* Harvey P. Dale, Aug 08 2022 *)

values(m)=my(v=[]);for(i=1,m,v=concat([1,4],v));v; /* Anders Hellström, Aug 03 2015 */

[power_mod(4,n,5)for n in range(0,81)] # Zerinvary Lajos, Nov 26 2009

a(2n) = 1, a(2n+1) = 4.
From Paul Barry, Jun 03 2003: (Start)
G.f.: (1+4*x)/((1-x)*(1+x)).
E.g.f.: (5*exp(x) - 3*exp(-x))/2.
a(n) = (5 - 3*(-1)^n)/2.
a(n) = 4^((1-(-1)^n)/2) = 2^(1-(-1)^n) = 2/(2^((-1)^n)).
a(n) = 4^(ceiling(n/2) - floor(n/2)). (End)
a(n) = gcd((n-1)^2, (n+1)^2). - Paul Barry, Sep 16 2004
a(n) = A160700(A000302(n)). - Reinhard Zumkeller, Jun 10 2009
a(n) = 4^n mod 5. - Zerinvary Lajos, Nov 26 2009
a(n) = 4^(n mod 2). - Wesley Ivan Hurt, Mar 29 2014

A191722 Dispersion of A008851, (numbers >1 and congruent to 0 or 1 mod 5), by antidiagonals.

Original entry on oeis.org

1, 5, 2, 15, 6, 3, 40, 16, 10, 4, 101, 41, 26, 11, 7, 255, 105, 66, 30, 20, 8, 640, 265, 166, 76, 51, 21, 9, 1601, 665, 416, 191, 130, 55, 25, 12, 4005, 1665, 1041, 480, 326, 140, 65, 31, 13, 10015, 4165, 2605, 1201, 816, 351, 165, 80, 35, 14, 25040, 10415

Offset: 1

Clark Kimberling, Jun 13 2011

For a background discussion of dispersions and their fractal sequences, see A191426.  For dispersions of congruence sequences mod 3, mod 4, or mod 5, see A191655, A191663, A191667, A191702.
...
Suppose that {2,3,4,5,6} is partitioned as {x1, x2} and {x3,x4,x5}.  Let S be the increasing sequence of numbers >1 and congruent to x1 or x2 mod 5, and let T be the increasing sequence of numbers >1 and congruent to x3 or x4 or x5 mod 5.  There are 10 sequences in S, each matched by a (nearly) complementary sequence in T.  Each of the 20 sequences generates a dispersion, as listed here:
...
A191722=dispersion of A008851 (0, 1 mod 5 and >1)
A191723=dispersion of A047215 (0, 2 mod 5 and >1)
A191724=dispersion of A047218 (0, 3 mod 5 and >1)
A191725=dispersion of A047208 (0, 4 mod 5 and >1)
A191726=dispersion of A047216 (1, 2 mod 5 and >1)
A191727=dispersion of A047219 (1, 3 mod 5 and >1)
A191728=dispersion of A047209 (1, 4 mod 5 and >1)
A191729=dispersion of A047221 (2, 3 mod 5 and >1)
A191730=dispersion of A047211 (2, 4 mod 5 and >1)
A191731=dispersion of A047204 (3, 4 mod 5 and >1)
...
A191732=dispersion of A047202 (2,3,4 mod 5 and >1)
A191733=dispersion of A047206 (1,3,4 mod 5 and >1)
A191734=dispersion of A032793 (1,2,4 mod 5 and >1)
A191735=dispersion of A047223 (1,2,3 mod 5 and >1)
A191736=dispersion of A047205 (0,3,4 mod 5 and >1)
A191737=dispersion of A047212 (0,2,4 mod 5 and >1)
A191738=dispersion of A047222 (0,2,3 mod 5 and >1)
A191739=dispersion of A008854 (0,1,4 mod 5 and >1)
A191740=dispersion of A047220 (0,1,3 mod 5 and >1)
A191741=dispersion of A047217 (0,1,2 mod 5 and >1)
...
EXCEPT for at most 2 initial terms (so that column 1 always starts with 1):
A191722 has 1st col A047202, all else A008851
A191723 has 1st col A047206, all else A047215
A191724 has 1st col A032793, all else A047218
A191725 has 1st col A047223, all else A047208
A191726 has 1st col A047205, all else A047216
A191727 has 1st col A047212, all else A047219
A191728 has 1st col A047222, all else A047209
A191729 has 1st col A008854, all else A047221
A191730 has 1st col A047220, all else A047211
A191731 has 1st col A047217, all else A047204
...
A191732 has 1st col A000851, all else A047202
A191733 has 1st col A047215, all else A047206
A191734 has 1st col A047218, all else A032793
A191735 has 1st col A047208, all else A047223
A191736 has 1st col A047216, all else A047205
A191737 has 1st col A047219, all else A047212
A191738 has 1st col A047209, all else A047222
A191739 has 1st col A047221, all else A008854
A191740 has 1st col A047211, all else A047220
A191741 has 1st col A047204, all else A047217
...
Regarding the dispersions A191722-A191741, there are general formulas for sequences of the type "(a or b mod m)" and "(a or b or c mod m)" used in the relevant Mathematica programs.

			Northwest corner:
1....5....15...40...101
2....6....16...41...105
3....10...26...66...166
4....11...30...76...191
7....20...51...130..326
8....21...55...140..351

Ivan Neretin, Table of n, a(n) for n = 1..5050

Cf. A047202, A008851, A191732, A191702, A191426.

(* Program generates the dispersion array t of the increasing sequence f[n] *)
r = 40; r1 = 12;  c = 40; c1 = 12;
a=5; b=6; m[n_]:=If[Mod[n,2]==0,1,0];
f[n_]:=a*m[n+1]+b*m[n]+5*Floor[(n-1)/2]
Table[f[n], {n, 1, 30}]  (* A008851 *)
mex[list_] := NestWhile[#1 + 1 &, 1, Union[list][[#1]] <= #1 &, 1, Length[Union[list]]]
rows = {NestList[f, 1, c]};
Do[rows = Append[rows, NestList[f, mex[Flatten[rows]], r]], {r}];
t[i_, j_] := rows[[i, j]];
TableForm[Table[t[i, j], {i, 1, 10}, {j, 1, 10}]]
(* A191722 *)
Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191722  *)

A191723 Dispersion of A047215, (numbers >1 and congruent to 0 or 2 mod 5), by antidiagonals.

Original entry on oeis.org

1, 2, 3, 5, 7, 4, 12, 17, 10, 6, 30, 42, 25, 15, 8, 75, 105, 62, 37, 20, 9, 187, 262, 155, 92, 50, 22, 11, 467, 655, 387, 230, 125, 55, 27, 13, 1167, 1637, 967, 575, 312, 137, 67, 32, 14, 2917, 4092, 2417, 1437, 780, 342, 167, 80, 35, 16, 7292, 10230, 6042

Offset: 1

Clark Kimberling, Jun 13 2011

For a background discussion of dispersions and their fractal sequences, see A191426.  For dispersions of congruence sequences mod 3, mod 4, or mod 5, see A191655, A191663, A191667, A191702.
...
Suppose that {2,3,4,5,6} is partitioned as {x1, x2} and {x3,x4,x5}.  Let S be the increasing sequence of numbers >1 and congruent to x1 or x2 mod 5, and let T be the increasing sequence of numbers >1 and congruent to x3 or x4 or x5 mod 5.  There are 10 sequences in S, each matched by a (nearly) complementary sequence in T.  Each of the 20 sequences generates a dispersion, as listed here:
...
A191722=dispersion of A008851 (0, 1 mod 5 and >1)
A191723=dispersion of A047215 (0, 2 mod 5 and >1)
A191724=dispersion of A047218 (0, 3 mod 5 and >1)
A191725=dispersion of A047208 (0, 4 mod 5 and >1)
A191726=dispersion of A047216 (1, 2 mod 5 and >1)
A191727=dispersion of A047219 (1, 3 mod 5 and >1)
A191728=dispersion of A047209 (1, 4 mod 5 and >1)
A191729=dispersion of A047221 (2, 3 mod 5 and >1)
A191730=dispersion of A047211 (2, 4 mod 5 and >1)
A191731=dispersion of A047204 (3, 4 mod 5 and >1)
...
A191732=dispersion of A047202 (2,3,4 mod 5 and >1)
A191733=dispersion of A047206 (1,3,4 mod 5 and >1)
A191734=dispersion of A032793 (1,2,4 mod 5 and >1)
A191735=dispersion of A047223 (1,2,3 mod 5 and >1)
A191736=dispersion of A047205 (0,3,4 mod 5 and >1)
A191737=dispersion of A047212 (0,2,4 mod 5 and >1)
A191738=dispersion of A047222 (0,2,3 mod 5 and >1)
A191739=dispersion of A008854 (0,1,4 mod 5 and >1)
A191740=dispersion of A047220 (0,1,3 mod 5 and >1)
A191741=dispersion of A047217 (0,1,2 mod 5 and >1)
...
For further information about these 20 dispersions, see A191722.
...
Regarding the dispersions A191722-A191741, there are general formulas for sequences of the type "(a or b mod m)" and "(a or b or c mod m)" used in the relevant Mathematica programs.

			Northwest corner:
1....2....5....12....30
3....7....17...42....105
4....10...25...62....155
6....15...37...92....230
8....20...50...125...312
9....22...55...137...342

Ivan Neretin, Table of n, a(n) for n = 1..5050

Cf. A047206, A047215, A191733, A191722, A191426.

(* Program generates the dispersion array t of the increasing sequence f[n] *)
r = 40; r1 = 12;  c = 40; c1 = 12;
a=2; b=5; m[n_]:=If[Mod[n,2]==0,1,0];
f[n_]:=a*m[n+1]+b*m[n]+5*Floor[(n-1)/2]
Table[f[n], {n, 1, 30}]  (* A047215 *)
mex[list_] := NestWhile[#1 + 1 &, 1, Union[list][[#1]] <= #1 &, 1, Length[Union[list]]]
rows = {NestList[f, 1, c]};
Do[rows = Append[rows, NestList[f, mex[Flatten[rows]], r]], {r}];
t[i_, j_] := rows[[i, j]];
TableForm[Table[t[i, j], {i, 1, 10}, {j, 1, 10}]] (* A191722 *)
Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191722  *)

A191724 Dispersion of A047218, (numbers >1 and congruent to 0 or 3 mod 5), by antidiagonals.

Original entry on oeis.org

1, 3, 2, 8, 5, 4, 20, 13, 10, 6, 50, 33, 25, 15, 7, 125, 83, 63, 38, 18, 9, 313, 208, 158, 95, 45, 23, 11, 783, 520, 395, 238, 113, 58, 28, 12, 1958, 1300, 988, 595, 283, 145, 70, 30, 14, 4895, 3250, 2470, 1488, 708, 363, 175, 75, 35, 16, 12238, 8125, 6175

Offset: 1

Clark Kimberling, Jun 13 2011

For a background discussion of dispersions and their fractal sequences, see A191426.  For dispersions of congruence sequences mod 3, mod 4, or mod 5, see A191655, A191663, A191667, A191702.
...
Suppose that {2,3,4,5,6} is partitioned as {x1, x2} and {x3,x4,x5}.  Let S be the increasing sequence of numbers >1 and congruent to x1 or x2 mod 5, and let T be the increasing sequence of numbers >1 and congruent to x3 or x4 or x5 mod 5.  There are 10 sequences in S, each matched by a (nearly) complementary sequence in T.  Each of the 20 sequences generates a dispersion, as listed here:
...
A191722=dispersion of A008851 (0, 1 mod 5 and >1)
A191723=dispersion of A047215 (0, 2 mod 5 and >1)
A191724=dispersion of A047218 (0, 3 mod 5 and >1)
A191725=dispersion of A047208 (0, 4 mod 5 and >1)
A191726=dispersion of A047216 (1, 2 mod 5 and >1)
A191727=dispersion of A047219 (1, 3 mod 5 and >1)
A191728=dispersion of A047209 (1, 4 mod 5 and >1)
A191729=dispersion of A047221 (2, 3 mod 5 and >1)
A191730=dispersion of A047211 (2, 4 mod 5 and >1)
A191731=dispersion of A047204 (3, 4 mod 5 and >1)
...
A191732=dispersion of A047202 (2,3,4 mod 5 and >1)
A191733=dispersion of A047206 (1,3,4 mod 5 and >1)
A191734=dispersion of A032793 (1,2,4 mod 5 and >1)
A191735=dispersion of A047223 (1,2,3 mod 5 and >1)
A191736=dispersion of A047205 (0,3,4 mod 5 and >1)
A191737=dispersion of A047212 (0,2,4 mod 5 and >1)
A191738=dispersion of A047222 (0,2,3 mod 5 and >1)
A191739=dispersion of A008854 (0,1,4 mod 5 and >1)
A191740=dispersion of A047220 (0,1,3 mod 5 and >1)
A191741=dispersion of A047217 (0,1,2 mod 5 and >1)
...
For further information about these 20 dispersions, see A191722.
...
Regarding the dispersions A191722-A191741, there are general formulas for sequences of the type "(a or b mod m)" and "(a or b or c mod m)" used in the relevant Mathematica programs.

			Northwest corner:
1....3....8....20....50
2....5....13...33....83
4....10...25...63....158
6....15...38...95....238
7....18...45...113...283
9....23...58...145...363

Ivan Neretin, Table of n, a(n) for n = 1..5050

Cf. A032793, A047218, A191734, A191722, A191426.

(* Program generates the dispersion array t of the increasing sequence f[n] *)
r = 40; r1 = 12;  c = 40; c1 = 12;
a=3; b=5; m[n_]:=If[Mod[n,2]==0,1,0];
f[n_]:=a*m[n+1]+b*m[n]+5*Floor[(n-1)/2]
Table[f[n], {n, 1, 30}]  (* A047218 *)
mex[list_] := NestWhile[#1 + 1 &, 1, Union[list][[#1]] <= #1 &, 1, Length[Union[list]]]
rows = {NestList[f, 1, c]};
Do[rows = Append[rows, NestList[f, mex[Flatten[rows]], r]], {r}];
t[i_, j_] := rows[[i, j]];
TableForm[Table[t[i, j], {i, 1, 10}, {j, 1, 10}]] (* A191724 *)
Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191724  *)

A191725 Dispersion of A047208, (numbers >1 and congruent to 0 or 4 mod 5), by antidiagonals.

Original entry on oeis.org

1, 4, 2, 10, 5, 3, 25, 14, 9, 6, 64, 35, 24, 15, 7, 160, 89, 60, 39, 19, 8, 400, 224, 150, 99, 49, 20, 11, 1000, 560, 375, 249, 124, 50, 29, 12, 2500, 1400, 939, 624, 310, 125, 74, 30, 13, 6250, 3500, 2349, 1560, 775, 314, 185, 75, 34, 16, 15625, 8750, 5874

Offset: 1

Clark Kimberling, Jun 13 2011

For a background discussion of dispersions and their fractal sequences, see A191426.  For dispersions of congruence sequences mod 3, mod 4, or mod 5, see A191655, A191663, A191667, A191702.
...
Suppose that {2,3,4,5,6} is partitioned as {x1, x2} and {x3,x4,x5}.  Let S be the increasing sequence of numbers >1 and congruent to x1 or x2 mod 5, and let T be the increasing sequence of numbers >1 and congruent to x3 or x4 or x5 mod 5.  There are 10 sequences in S, each matched by a (nearly) complementary sequence in T.  Each of the 20 sequences generates a dispersion, as listed here:
...
A191722=dispersion of A008851 (0, 1 mod 5 and >1)
A191723=dispersion of A047215 (0, 2 mod 5 and >1)
A191724=dispersion of A047218 (0, 3 mod 5 and >1)
A191725=dispersion of A047208 (0, 4 mod 5 and >1)
A191726=dispersion of A047216 (1, 2 mod 5 and >1)
A191727=dispersion of A047219 (1, 3 mod 5 and >1)
A191728=dispersion of A047209 (1, 4 mod 5 and >1)
A191729=dispersion of A047221 (2, 3 mod 5 and >1)
A191730=dispersion of A047211 (2, 4 mod 5 and >1)
A191731=dispersion of A047204 (3, 4 mod 5 and >1)
...
A191732=dispersion of A047202 (2,3,4 mod 5 and >1)
A191733=dispersion of A047206 (1,3,4 mod 5 and >1)
A191734=dispersion of A032793 (1,2,4 mod 5 and >1)
A191735=dispersion of A047223 (1,2,3 mod 5 and >1)
A191736=dispersion of A047205 (0,3,4 mod 5 and >1)
A191737=dispersion of A047212 (0,2,4 mod 5 and >1)
A191738=dispersion of A047222 (0,2,3 mod 5 and >1)
A191739=dispersion of A008854 (0,1,4 mod 5 and >1)
A191740=dispersion of A047220 (0,1,3 mod 5 and >1)
A191741=dispersion of A047217 (0,1,2 mod 5 and >1)
...
For further information about these 20 dispersions, see A191722.
...
Regarding the dispersions A191722-A191741, there are general formulas for sequences of the type "(a or b mod m)" and "(a or b or c mod m)" used in the relevant Mathematica programs.

			Northwest corner:
1....4....10....25....64
2....5....14....35...89
3....9....24...60...150
6....15...39...99...249
7....19...49...124..310
8....20...50...125...314

Ivan Neretin, Table of n, a(n) for n = 1..5050

Cf. A047223, A047208, A191735, A191722, A191426.

(* Program generates the dispersion array t of the increasing sequence f[n] *)
r = 40; r1 = 12;  c = 40; c1 = 12;
a=4; b=5; m[n_]:=If[Mod[n,2]==0,1,0];
f[n_]:=a*m[n+1]+b*m[n]+5*Floor[(n-1)/2]
Table[f[n], {n, 1, 30}]  (* A047208 *)
mex[list_] := NestWhile[#1 + 1 &, 1, Union[list][[#1]] <= #1 &, 1, Length[Union[list]]]
rows = {NestList[f, 1, c]};
Do[rows = Append[rows, NestList[f, mex[Flatten[rows]], r]], {r}];
t[i_, j_] := rows[[i, j]];
TableForm[Table[t[i, j], {i, 1, 10}, {j, 1, 10}]] (* A191725 *)
Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191725  *)

cp's OEIS Frontend

A191732 Dispersion of A047202, (numbers >1 and congruent to 2 or 3 or 4 mod 5), by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A191708 Dispersion of A047202, (numbers >1 and congruent to 0, 2, 3, or 4 mod 5), by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A191709 Dispersion of A047202, (numbers >1 and congruent to 0, 1, 3, or 4 mod 5), by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A016813 a(n) = 4*n + 1.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

PARI

Scala

Formula

A191702 Dispersion of A008587 (5,10,15,20,25,30,...), by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A010685 Period 2: repeat (1,4).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

A191722 Dispersion of A008851, (numbers >1 and congruent to 0 or 1 mod 5), by antidiagonals.

Views

Author

Keywords

Comments