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 11 results. Next

A351415 Intersection of Beatty sequences for (1+sqrt(5))/2 and sqrt(5).

Original entry on oeis.org

4, 6, 8, 11, 17, 22, 24, 29, 33, 35, 38, 40, 42, 46, 51, 53, 55, 58, 64, 67, 69, 71, 76, 80, 82, 84, 87, 93, 98, 100, 105, 111, 114, 116, 118, 122, 127, 129, 131, 134, 140, 145, 147, 152, 156, 158, 160, 163, 165, 169, 174, 176, 181, 187, 190, 192, 194, 199
Offset: 1

Views

Author

Clark Kimberling, Feb 10 2022

Keywords

Comments

Conjecture: every term of the difference sequence is in {2,3,4,5,6}, and each occurs infinitely many times.
From Clark Kimberling, Jul 29 2022: (Start)
This is the first of four sequences that partition the positive integers. Starting with a general overview, suppose that u = (u(n)) and v = (v(n)) are increasing sequences of positive integers. Let u' and v' be their complements, and assume that the following four sequences are infinite:
(1) u ^ v = intersection of u and v (in increasing order);
(2) u ^ v';
(3) u' ^ v;
(4) u' ^ v'.
Every positive integer is in exactly one of the four sequences. For A351415, u, v, u', v', are the Beatty sequences given by u(n) = floor(n*(1+sqrt(5))/2) and v(n) = floor(n*sqrt(5)), so that r = (1+sqrt(5))/2, s = sqrt(5), r' = (3+sqrt(5))/2, s' = (5 + sqrt(5))/4.
(1) u ^ v = (4, 6, 8, 11, 17, 22, 24, 29, 33, 35, 38, 40, 42, ...) = A351415
(2) u ^ v' = (1, 3, 9, 12, 14, 16, 19, 21, 25, 27, 30, 32, ...) = A356101
(3) u' ^ v = (2, 13, 15, 20, 26, 31, 44, 49, 60, 62, 73, 78, ...) = A356102
(4) u' ^ v' = (5, 7, 10, 18, 23, 28, 34, 36, 39, 41, 47, 52, ...) = A356103
(End)

Examples

			The two Beatty sequences are (1,3,4,6,8,9,11,12,14,...) and (2,4,6,8,11,13,15,17,...), with common terms forming the sequence (4,6,8,11,...).
		

Crossrefs

Cf. A001950, A108598, A356101, A356102, A356103, A356104 (results of composition instead of intersections), A190509 (composites, reversed order).

Programs

  • Mathematica
    z = 200;
    r = (1 + Sqrt[5])/2; u = Table[Floor[n*r], {n, 1, z}]  (* A000201 *)
    u1 = Take[Complement[Range[1000], u], z]  (* A001950 *)
    r1 = Sqrt[5]; v = Table[Floor[n*r1], {n, 1, z}]  (* A022839 *)
    v1 = Take[Complement[Range[1000], v], z]  (* A108598 *)
    Intersection[u, v]    (* A351415 *)
    Intersection[u, v1]   (* A356101 *)
    Intersection[u1, v]   (* A356102 *)
    Intersection[u1, v1]  (* A356103 *)

A356217 a(n) = A022839(A000201(n)).

Original entry on oeis.org

2, 6, 8, 13, 17, 20, 24, 26, 31, 35, 38, 42, 46, 49, 53, 55, 60, 64, 67, 71, 73, 78, 82, 84, 89, 93, 96, 100, 102, 107, 111, 114, 118, 122, 125, 129, 131, 136, 140, 143, 147, 149, 154, 158, 160, 165, 169, 172, 176, 178, 183, 187, 190, 194, 196, 201, 205, 207
Offset: 1

Views

Author

Clark Kimberling, Oct 02 2022

Keywords

Comments

This is the first of four sequences that partition the positive integers. Suppose that u = (u(n)) and v = (v(n)) are increasing sequences of positive integers. Let u' and v' be their (increasing) complements, and consider these four sequences:
(1) v o u, defined by (v o u)(n) = v(u(n));
(2) v' o u;
(3) v o u';
(4) v' o u'.
Every positive integer is in exactly one of the four sequences. For the reverse composites, u o v, u o v', u' o v, u' o v', see A356104 to A356107.
Assume that if w is any of the sequences u, v, u', v', then lim_{n->oo} w(n)/n exists and defines the (limiting) density of w. For w = u,v,u',v', denote the densities by r,s,r',s'. Then the densities of sequences (1)-(4) exist, and
1/(r*r') + 1/(r*s') + 1/(s*s') + 1/(s*r') = 1.
For A356217 u, v, u', v', are the Beatty sequences given by u(n) = floor(n*(1+sqrt(5))/2) and v(n) = floor(n*sqrt(5)), so that r = (1+sqrt(5))/2, s = sqrt(5), r' = (3+sqrt(5))/2, s' = (5 + sqrt(5))/4.

Examples

			(1)  v o u = (2, 6, 8, 13, 17, 20, 24, 26, 31, 35, 38, 42, ...) = A356217
(2)  v' o u = (1, 5, 7, 10, 14, 16, 19, 21, 25, 28, 30, 34, ...) = A356218
(3)  v o u' = (4, 11, 15, 22, 29, 33, 40, 44, 51, 58, 62, 76, ...) = A190509
(4)  v' o u' = (3, 9, 12, 18, 23, 27, 32, 36, 41, 47, 50, 56, ...) = A356220
		

Crossrefs

Cf. u = A000201, u' = A001950, v = A022839, v' = A108598, A351415 (intersections), A356104 (reverse composites), A356218, A190509, A356220.

Programs

  • Mathematica
    z = 1000;
    u = Table[Floor[n*(1 + Sqrt[5])/2], {n, 1, z}];  (* A000201 *)
    u1 = Complement[Range[Max[u]], u];  (* A001950 *)
    v = Table[Floor[n*Sqrt[5]], {n, 1, z}];  (* A022839 *)
    v1 = Complement[Range[Max[v]], v];  (* A108598 *)
    Table[v[[u[[n]]]], {n, 1, z/4}]   (* A356217 *)
    Table[v1[[u[[n]]]], {n, 1, z/4}]  (* A356218 *)
    Table[v[[u1[[n]]]], {n, 1, z/4}]  (* A190509 *)
    Table[v1[[u1[[n]]]], {n, 1, z/4}] (* A356220 *)
  • Python
    from math import isqrt
    def A356217(n): return isqrt(5*(n+isqrt(5*n**2)>>1)**2) # Chai Wah Wu, Oct 14 2022

A356220 a(n) = A108598(A001950(n)).

Original entry on oeis.org

3, 9, 12, 18, 23, 27, 32, 36, 41, 47, 50, 56, 61, 65, 70, 74, 79, 85, 88, 94, 97, 103, 108, 112, 117, 123, 126, 132, 135, 141, 146, 150, 155, 161, 164, 170, 173, 179, 184, 188, 193, 197, 202, 208, 211, 217, 222, 226, 231, 235, 240, 246, 249, 255, 258, 264
Offset: 1

Views

Author

Clark Kimberling, Nov 13 2022

Keywords

Comments

This is the fourth of four sequences that partition the positive integers. See A356217.

Examples

			(1)  v o u = (2, 6, 8, 13, 17, 20, 24, 26, 31, 35, 38, 42, ...) = A356217
(2)  v' o u = (1, 5, 7, 10, 14, 16, 19, 21, 25, 28, 30, 34, ...) = A356218
(3)  v o u' = (4, 11, 15, 22, 29, 33, 40, 44, 51, 58, 62, 76, ...) = A190509
(4)  v' o u' = (3, 9, 12, 18, 23, 27, 32, 36, 41, 47, 50, 56, ...) = A356220
		

Crossrefs

Cf. A000201, A001950, A022839, A108598, A351415 (intersections), A356104 (reverse composites), A356217, A356218, A356219.

Programs

  • Mathematica
    z = 1000;
    u = Table[Floor[n*(1 + Sqrt[5])/2], {n, 1, z}];  (* A000201 *)
    u1 = Complement[Range[Max[u]], u];  (* A001950 *)
    v = Table[Floor[n*Sqrt[5]], {n, 1, z}];  (* A022839 *)
    v1 = Complement[Range[Max[v]], v];  (* A108598 *)
    zz = 120;
    Table[v[[u[[n]]]], {n, 1, z/4}]   (* A356217 *)
    Table[v1[[u[[n]]]], {n, 1, z/4}]  (* A356218 *)
    Table[v[[u1[[n]]]], {n, 1, z/4}]  (* A190509 *)
    Table[v1[[u1[[n]]]], {n, 1, z/4}] (* A356220 *)

A190508 a(n) = n+[n*s/r]+[n*t/r]+[n*u/r]; r=golden ratio, s=r^2, t=r^3, u=r^4.

Original entry on oeis.org

8, 18, 26, 36, 47, 55, 65, 73, 84, 94, 102, 112, 123, 131, 141, 149, 160, 170, 178, 188, 196, 207, 217, 225, 235, 246, 254, 264, 272, 283, 293, 301, 311, 322, 330, 340, 348, 358, 369, 377, 387, 395, 406, 416, 424, 434, 445, 453, 463, 471, 482, 492, 500, 510, 518, 529, 539, 547, 557, 568, 576, 586, 594, 605, 615, 623, 633, 644
Offset: 1

Views

Author

Clark Kimberling, May 11 2011

Keywords

Comments

This is one of four sequences that partition the positive integers. In general, suppose that r, s, t, u are positive real numbers for which the sets {i/r: i>=1}, {j/s: j>=1}, {k/t: k>=1}, {h/u: h>=1} are pairwise disjoint. Let a(n) be the rank of n/r when all the numbers in the four sets are jointly ranked. Define b(n), c(n), d(n) as the ranks of n/s, n/t, n/u, respectively.
It is easy to prove that
a(n)=n+[n*s/r]+[n*t/r]+[n*u/r],
b(n)=n+[n*r/s]+[n*t/s]+[n*u/s],
c(n)=n+[n*r/t]+[n*s/t]+[n*u/t],
d(n)=n+[n*r/u]+[n*s/u]+[n*t/u], where []=floor.
Taking r=golden ratio, s=r^2, t=r^3, u=r^4 gives a=A190508, b=A190509, c=A054770, d=A190511.

Crossrefs

Cf. A190509, A054770, A190511 (the other three sequences in the partition of N).

Programs

  • Mathematica
    r=GoldenRatio; s=r^2; t=r^3; u=r^4;
    a[n_] := n + Floor[n*s/r] + Floor[n*t/r]+Floor[n*u/r];
    b[n_] := n + Floor[n*r/s] + Floor[n*t/s]+Floor[n*u/s];
    c[n_] := n + Floor[n*r/t] + Floor[n*s/t]+Floor[n*u/t];
    d[n_] := n + Floor[n*r/u] + Floor[n*s/u]+Floor[n*t/u];
    Table[a[n], {n, 1, 120}]  (*A190508*)
    Table[b[n], {n, 1, 120}]  (*A190509*)
    Table[c[n], {n, 1, 120}]  (*A054770*)
    Table[d[n], {n, 1, 120}]  (*A190511*)

Formula

A190508: a(n)=n+[n*r]+[n*r^2]+[n*r^3]
A190509: b(n)=[n/r]+n+[n*r]+[n*r^2]
A054770: c(n)=[n/r^2]+[n/r]+n+[n*r]
A190511: d(n)=[n/r^3]+[n/r^2]+[n/r]+n

A190511 n+[nr/u]+[ns/u]+[nt/u]; r=golden ratio, s=r^2, t=r^3, u=r^4.

Original entry on oeis.org

1, 3, 5, 7, 10, 12, 14, 16, 19, 21, 23, 25, 28, 30, 32, 34, 37, 39, 41, 43, 45, 48, 50, 52, 54, 57, 59, 61, 63, 66, 68, 70, 72, 75, 77, 79, 81, 83, 86, 88, 90, 92, 95, 97, 99, 101, 104, 106, 108, 110, 113, 115, 117, 119, 121, 124, 126, 128, 130, 133, 135, 137, 139, 142, 144, 146, 148, 151, 153, 155, 157, 159, 162, 164, 166, 168, 171, 173
Offset: 1

Views

Author

Clark Kimberling, May 11 2011

Keywords

Comments

A190508: a(n)=n+[nr]+[nr^2]+[nr^3]
A190509: b(n)=[n/r]+n+[nr]+[nr^2]
A054770: c(n)=[n/r^2]+[n/r]+n+[nr]
A190511: d(n)=[n/r^3]+[n/r^2]+[n/r]+n
Are the differences between successive terms always 2 or 3? - Harvey P. Dale, Apr 03 2025

Crossrefs

Programs

  • Mathematica
    (See A190508.)
    Module[{g=GoldenRatio,g4},g4=g^4;Table[n+Floor[(n*g)/g4]+Floor[(n*g^2)/g4]+Floor[(n*g^3)/g4],{n,80}]] (* Harvey P. Dale, Apr 03 2025 *)

Formula

a(n) = A022839(n)-1. - Michel Dekking, May 04 2019

A356101 Intersection of A000201 and A022839.

Original entry on oeis.org

1, 3, 9, 12, 14, 16, 19, 21, 25, 27, 30, 32, 37, 43, 45, 48, 50, 56, 59, 61, 63, 66, 72, 74, 77, 79, 85, 88, 90, 92, 95, 97, 101, 103, 106, 108, 110, 113, 119, 121, 124, 126, 132, 135, 137, 139, 142, 144, 148, 150, 153, 155, 161, 166, 168, 171, 173, 177, 179
Offset: 1

Views

Author

Clark Kimberling, Sep 04 2022

Keywords

Comments

This is the second of four sequences that partition the positive integers. See A351415.

Examples

			Starting with a general overview, suppose that u = (u(n)) and v = (v(n)) are increasing sequences of positive integers.  Let u' and v' be their complements, and assume that the following four sequences are infinite:
(1)  u ^ v = intersection of u and v (in increasing order);
(2)  u ^ v';
(3)  u' ^ v;
(4)  u' ^ v'.
Every positive integer is in exactly one of the four sequences. For A351415, u, v, u', v', are the Beatty sequences given by u(n) = floor(n*(1+sqrt(5))/2) and v(n) = floor(n*sqrt(5)), so r = (1+sqrt(5))/2, s = sqrt(5), r' = (3+sqrt(5))/2, s' = (5 + sqrt(5))/4.
(1)  u ^ v = (4, 6, 8, 11, 17, 22, 24, 29, 33, 35, 38, 40, 42, ...) =  A351415
(2)  u ^ v' = (1, 3, 9, 12, 14, 16, 19, 21, 25, 27, 30, 32, 37, ...) =  A356101
(3)  u' ^ v = (2, 13, 15, 20, 26, 31, 44, 49, 60, 62, 73, 78, ...) = A356102
(4)  u' ^ v' = (5, 7, 10, 18, 23, 28, 34, 36, 39, 41, 47, 52, 54, ...) = A356103
		

Crossrefs

Cf. u = A000201, u' = A001950, v = A022839, v' = A108598, A351415, A356102, A356103, A356104 (results of composition instead of intersections), A190509 (composites, reversed order).

Programs

  • Mathematica
    z = 200;
    r = (1 + Sqrt[5])/2; u = Table[Floor[n*r], {n, 1, z}]  (* A000201 *)
    u1 = Take[Complement[Range[1000], u], z]  (* A001950 *)
    r1 = Sqrt[5]; v = Table[Floor[n*r1], {n, 1, z}]  (* A022839 *)
    v1 = Take[Complement[Range[1000], v], z]  (* A108598 *)
    Intersection[u, v]   (* A351415 *)
    Intersection[u, v1]  (* A356101 *)
    Intersection[u1, v]  (* A356102 *)
    Intersection[u1, v1] (* A356103 *)

A356102 Intersection of A001950 and A022839.

Original entry on oeis.org

2, 13, 15, 20, 26, 31, 44, 49, 60, 62, 73, 78, 89, 91, 96, 102, 107, 109, 120, 125, 136, 138, 143, 149, 154, 167, 172, 178, 183, 185, 196, 201, 212, 214, 219, 225, 230, 243, 248, 259, 261, 272, 277, 290, 295, 301, 306, 308, 319, 324, 326, 328, 330, 333, 335
Offset: 1

Views

Author

Clark Kimberling, Sep 04 2022

Keywords

Comments

This is the third of four sequences that partition the positive integers. See A351415.

Examples

			Starting with a general overview, suppose that u = (u(n)) and v = (v(n)) are increasing sequences of positive integers.  Let u' and v' be their complements, and assume that the following four sequences are infinite:
(1)  u ^ v = intersection of u and v (in increasing order);
(2)  u ^ v';
(3)  u' ^ v;
(4)  u' ^ v'.
Every positive integer is in exactly one of the four sequences. For A351415, u, v, u', v', are the Beatty sequences given by u(n) = floor(n*(1+sqrt(5))/2) and v(n) = floor(n*sqrt(5)), so that r = (1+sqrt(5))/2, s = sqrt(5), r' = (3+sqrt(5))/2, s' = (5 + sqrt(5))/4.
(1)  u ^ v = (4, 6, 8, 11, 17, 22, 24, 29, 33, 35, 38, 40, 42, ...) =  A351415
(2)  u ^ v' = (1, 3, 9, 12, 14, 16, 19, 21, 25, 27, 30, 32, 37, ...) =  A356101
(3)  u' ^ v = (2, 13, 15, 20, 26, 31, 44, 49, 60, 62, 73, 78, ...) = A356102
(4)  u' ^ v' = (5, 7, 10, 18, 23, 28, 34, 36, 39, 41, 47, 52, 54, ...) = A356103
		

Crossrefs

Cf. u = A000201, u' = A001950, v = A022839, v' = A108598, A351415, A356101, A356103, A356104 (results of composition instead of intersections), A190509 (composites, reversed order).

Programs

  • Mathematica
    z = 200;
    r = (1 + Sqrt[5])/2; u = Table[Floor[n*r], {n, 1, z}]  (* A000201 *)
    u1 = Take[Complement[Range[1000], u], z]  (* A001950 *)
    r1 = Sqrt[5]; v = Table[Floor[n*r1], {n, 1, z}]  (* A022839 *)
    v1 = Take[Complement[Range[1000], v], z]  (* A108598 *)
    Intersection[u, v]   (* A351415 *)
    Intersection[u, v1]  (* A356101 *)
    Intersection[u1, v]  (* A356102 *)
    Intersection[u1, v1] (* A356103 *)

A356103 Intersection of A001950 and A108598.

Original entry on oeis.org

5, 7, 10, 18, 23, 28, 34, 36, 39, 41, 47, 52, 54, 57, 65, 68, 70, 75, 81, 83, 86, 94, 99, 104, 112, 115, 117, 123, 128, 130, 133, 141, 146, 151, 157, 159, 162, 164, 170, 175, 180, 188, 191, 193, 198, 204, 206, 209, 217, 222, 227, 233, 235, 238, 240, 246, 251
Offset: 1

Views

Author

Clark Kimberling, Sep 04 2022

Keywords

Comments

This is the fourth of four sequences that partition the positive integers. See A351415.

Examples

			Starting with a general overview, suppose that u = (u(n)) and v = (v(n)) are increasing sequences of positive integers.  Let u' and v' be their complements, and assume that the following four sequences are infinite:
(1)  u ^ v = intersection of u and v (in increasing order);
(2)  u ^ v';
(3)  u' ^ v;
(4)  u' ^ v'.
Every positive integer is in exactly one of the four sequences. For A351415, u, v, u', v', are the Beatty sequences given by u(n) = floor(n*(1+sqrt(5))/2) and v(n) = floor(n*sqrt(5)), so that r = (1+sqrt(5))/2, s = sqrt(5), r' = (3+sqrt(5))/2, s' = (5 + sqrt(5))/4.
(1)  u ^ v = (4, 6, 8, 11, 17, 22, 24, 29, 33, 35, 38, 40, 42, ...) =  A351415
(2)  u ^ v' = (1, 3, 9, 12, 14, 16, 19, 21, 25, 27, 30, 32, 37, ...) =  A356101
(3)  u' ^ v = (2, 13, 15, 20, 26, 31, 44, 49, 60, 62, 73, 78, ...) = A356102
(4)  u' ^ v' = (5, 7, 10, 18, 23, 28, 34, 36, 39, 41, 47, 52, 54, ...) = A356103
		

Crossrefs

Cf. u = A000201, u' = A001950, v = A022839, v' = A108598, A351415, A356101, A356102, A356104 (results of composition instead of intersections), A190509 (composites, reversed order).

Programs

  • Mathematica
    z = 200;
    r = (1 + Sqrt[5])/2; u = Table[Floor[n*r], {n, 1, z}]  (* A000201 *)
    u1 = Take[Complement[Range[1000], u], z]  (* A001950 *)
    r1 = Sqrt[5]; v = Table[Floor[n*r1], {n, 1, z}]  (* A022839 *)
    v1 = Take[Complement[Range[1000], v], z]  (* A108598 *)
    Intersection[u, v]   (* A351415 *)
    Intersection[u, v1]  (* A356101 *)
    Intersection[u1, v]  (* A356102 *)
    Intersection[u1, v1] (* A356103 *)

A356218 a(n) = A108598(A000201(n)).

Original entry on oeis.org

1, 5, 7, 10, 14, 16, 19, 21, 25, 28, 30, 34, 37, 39, 43, 45, 48, 52, 54, 57, 59, 63, 66, 68, 72, 75, 77, 81, 83, 86, 90, 92, 95, 99, 101, 104, 106, 110, 113, 115, 119, 121, 124, 128, 130, 133, 137, 139, 142, 144, 148, 151, 153, 157, 159, 162, 166, 168, 171
Offset: 1

Views

Author

Clark Kimberling, Oct 02 2022

Keywords

Comments

This is the second of four sequences that partition the positive integers. See A356217.

Examples

			(1)  v o u = (2, 6, 8, 13, 17, 20, 24, 26, 31, 35, 38, 42, ...) = A356217
(2)  v' o u = (1, 5, 7, 10, 14, 16, 19, 21, 25, 28, 30, 34, ...) = A356218
(3)  v o u' = (4, 11, 15, 22, 29, 33, 40, 44, 51, 58, 62, 76, ...) = A190509
(4)  v' o u' = (3, 9, 12, 18, 23, 27, 32, 36, 41, 47, 50, 56, ...) = A356220
		

Crossrefs

Cf. u = A000201, u' = A001950, v = A022839, v' = A108598, A351415 (intersections), A356104 (reverse composites), A356217, A190509, A356220.

Programs

  • Mathematica
    z = 1000;
    u = Table[Floor[n*(1 + Sqrt[5])/2], {n, 1, z}];  (* A000201 *)
    u1 = Complement[Range[Max[u]], u];  (* A001950 *)
    v = Table[Floor[n*Sqrt[5]], {n, 1, z}];  (* A022839 *)
    v1 = Complement[Range[Max[v]], v];  (* A108598 *)
    zz = 120;
    Table[v[[u[[n]]]], {n, 1, z/4}]   (* A356217 *)
    Table[v1[[u[[n]]]], {n, 1, z/4}]  (* A356218 *)
    Table[v[[u1[[n]]]], {n, 1, z/4}]  (* A190509 *)
    Table[v1[[u1[[n]]]], {n, 1, z/4}] (* A356220 *)

A340429 Array T(n, k) is the number x such that frac(x*phi) + frac(n*phi)*frac(k*phi) = 1 where phi is the golden ratio A001622 and frac(y) is the fractional part of y, read by antidiagonals.

Original entry on oeis.org

1, 3, 3, 4, 8, 4, 6, 11, 11, 6, 8, 16, 15, 16, 8, 9, 21, 22, 22, 21, 9, 11, 24, 29, 32, 29, 24, 11, 12, 29, 33, 42, 42, 33, 29, 12, 14, 32, 40, 48, 55, 48, 40, 32, 14, 16, 37, 44, 58, 63, 63, 58, 44, 37, 16, 17, 42, 51, 64, 76, 72, 76, 64, 51, 42, 17
Offset: 1

Views

Author

Michel Marcus, Jan 07 2021

Keywords

Examples

			Array begins:
  1  3  4  6  8 ...
  3  8 11 16 21 ...
  4 11 15 22 29 ...
  6 16 22 32 42 ...
  8 21 29 42 55 ...
  ...
		

Crossrefs

Cf. A000201 (row 1), A003623 (row 2), A190509 (row 3), A371388 (main diagonal).

Programs

  • Maple
    h := n -> ceil(2*n / (sqrt(5) + 3)):
    T := (n, k) -> 3*n*k - n*h(k) - k*h(n):
    seq(lprint(seq(T(n, k), k = 1..9)), n = 1..7);  # Peter Luschny, Mar 21 2024
  • Mathematica
    A340429[n_, k_] := Floor[n * GoldenRatio] * k + Floor[k * GoldenRatio] * n - n * k;
    Table[A340429[n - k + 1, k], {n, 15}, {k, n}] (* Paolo Xausa, Mar 21 2024 *)
  • PARI
    f(n) = 2*floor(n*(1+sqrt(5))/2) - 3*n; \\ A339765
    T(n, k) = 2*n*k + f(n)*k/2 + f(k)*n/2;

Formula

T(n, k) = 2*n*k + A339765(n)*k/2 + A339765(k)*n/2.
T(n, k) = T(k, n), array is symmetric.
T(n, k) = 3*n*k - n*h(k) - k*h(n) where h(n) = ceiling(2*n / (sqrt(5) + 3)) = A189663(n + 1). - Peter Luschny, Mar 21 2024
Showing 1-10 of 11 results. Next