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

A340747 Numbers in array A322744 that do not have a unique decomposition into numbers of A307002.

Original entry on oeis.org

24, 40, 60, 67, 72, 88, 96, 100, 120, 132, 136, 144, 147, 150, 160, 168, 180, 184, 200, 204, 216, 220, 227, 232, 240, 264, 267, 276, 280, 288, 300, 307, 312, 323, 328, 330, 340, 348, 352, 360, 367, 376, 384, 387, 396, 400, 408, 420, 424
Offset: 1

Views

Author

David Lovler, Jan 20 2021

Keywords

Comments

For i >= 2, A322744(i, a(n)) is in the sequence.
There are numbers in array A322744 that have three decompositions of the form A322744(4,p) = A322744(7,q) = A322744(10,r). In these cases, p = q + r. p, q and r need not be in A307002. There are two situations. (a) For n > 0, 60n = A322744(4,10n) = A322744(7,6n) = A322744(10,4n); (b) For n >= 0, 60n+40 = A322744(4,10n+7) = A322744(7,6n+4) = A322744(10,4n+3).
A proof of p = q + r. q must be even because A322744(7,q) = even. p and r must be both odd or both even, otherwise there is the contradiction that p gets equated with a fraction. When p and r are odd, (3*4*p - 4)/2 = (3*7*q - q)/2 = (3*10*r - 10)/2. Solving for p in terms of q, and p in terms of r gives p = (5/3)*q + 1/3 and p = (5/2)*r - 1/2. Multiplying the latter by 2/3 and adding the two equations gives (5/3)*p = (5/3)*q + (5/3)*r, thus p = q + r. When p and r are even, (3*4*p)/2 = (3*7*q - q)/2 = (3*10*r)/2, and the same follows.

Examples

			60 = A322744(4,10). Also 60 = A322744(6,7) and 60 = A322744(2,20). These decompositions are the same but different from A322744(4,10) as follows. 6 = A322744(2,2) and 20 = A322744(2,7), making 60 = A322744(A322744(2,2), 7) and 60 = A322744(2, A322744(2,7)). Thus 60 can be written as A322744(2,2,7), a well-defined composition because A322744(n,k) is associative. 2,4,7 and 10 are in A307002, thus A322744(4,10) and A322744(2,2,7) are different decompositions of 60, so 60 is in the sequence.
88 is in the sequence because 88 = A322744(3,22) = A322744(4,15) and 3,4,15 and 22 are in A307002.
Examples of A322744(4,p) = A322744(7,q) = A322744(10,r) with p = q + r:
60*1 + 40 = 100 = A322744(4,17) = A322744(7,10) = A322744(10,7) and 17 = 10 + 7, which works by commuting one of the decompositions. Note that 60 also works this way. 60 = A322744(4,10) = A322744(7,6) = A322744(10,4) and 10 = 6 + 4.
60*3 = 180 = A322744(4,30) = A322744(7,18) = A322744(10,12) and 30 = 18 + 12.
60*3 + 40 = 220 = A322744(4,37) = A322744(7,22) = A322744(10,15) and 37 = 22 + 15.
See A340746 for more examples.
		

Crossrefs

A340746 Numbers in array A322744 that do not have a unique decomposition into numbers of A307002 and are not equal to A322744(n,k), n > 1, k in the sequence.

Original entry on oeis.org

24, 40, 60, 67, 88, 100, 132, 136, 147, 150, 184, 204, 220, 227, 232, 276, 307, 323, 328, 330, 340, 348, 367, 376, 387, 424, 460, 472, 484, 492, 499, 510, 547, 550, 564, 567, 568, 580, 627, 636, 664, 675, 690, 707, 708, 712, 726, 748, 767
Offset: 1

Views

Author

David Lovler, Jan 18 2021

Keywords

Comments

While array A322744 has many properties of the multiplication table, one way the numbers that sieve out of the array fail to be prime numbers is that unique factorization does not hold. Some numbers have two or more decompositions.
The numbers in this sequence are primitive in the sense that they are not A322744 multiples of an earlier number in the sequence.

Examples

			24 = A322744(4,4). Also 24 = A322744(6,3) and 24 = A322744(2,8). These two decompositions are the same but they differ from A322744(4,4) as follows. 6 = A322744(2,2) and 8 = A322744(2,3), making 24 = A322744(A322744(2,2), 3) and 24 = A322744(2, A322744(2,3)). Thus 24 can be written as A322744(2,2,3), a well-defined composition because A322744(n,k) is associative. 2,3 and 4 are in A307002, thus A322744(4,4) and A322744(2,2,3) are different decompositions of 24, so 24 is in the sequence.
40 is in the sequence because 40 = A322744(3,10) = A322744(4,7) and 3,4,7 and 10 are in A307002.
67 is in the sequence because 67 = A322744(3,17) = A322744(7,7) and 3,7 and 17 are in A307002.
220 has three decompositions. 220 = A322744(4,37) = A322744(7,22) = A322744(10,15) and 4,7,10,15,22 and 37 are in A307002.
72 = A322744(2,2,2,3) = A322744(2,4,4) is not in the sequence because 72 = A322744(2,24) and 24 is in the sequence.
		

Crossrefs

Programs

  • Mathematica
    T[n_,k_]:=T[n,k]=(3*n*k-If[OddQ@n,If[OddQ@k,n+k-1,k],If[OddQ@k,n,0]])/2;F[d_]:=If[(q=Union[Sort/@(Position[Table[T[n,k],{n,2,Ceiling[d/3]},{k,2,Ceiling[d/3]}],d]+1)])=={},{{d}},q];FC[x_]:=FixedPoint[Union[Sort/@Flatten[Flatten/@Tuples[#]&/@((F/@#&/@#)&[#]),1]]&,F[x]];list={};Do[If[Length@FC@i>1&&ContainsNone[list,Flatten@F@i],AppendTo[list,i]],{i,500}];list (* Giorgos Kalogeropoulos, Apr 11 2021 *)

A327263 Array T(n,k) in which the i-th row consists of numbers > 1 not in array U(i;n,k) = (i*n*k - (i-2)*A319929(n,k))/2 where i >= 1, n >= 1 and k >= 1, read by antidiagonals.

Original entry on oeis.org

3, 5, 2, 9, 3, 2, 13, 5, 3, 2, 21, 7, 4, 3, 2, 25, 11, 5, 4, 3, 2, 33, 13, 7, 5, 4, 3, 2, 37, 17, 9, 6, 5, 4, 3, 2, 45, 19, 10, 7, 6, 5, 4, 3, 2, 57, 23, 13, 9, 7, 6, 5, 4, 3, 2, 61, 29, 15, 11, 8, 7, 6, 5, 4, 3, 2, 73, 31, 17, 12, 9, 8, 7, 6, 5, 4, 3, 2
Offset: 1

Views

Author

David Lovler, Oct 15 2019

Keywords

Comments

All the U(i;n,k) mimic the ordinary multiplication table in that they are commutative, associative, have identity element 1 and have 0. However (except when i=2) they are partially distributive, meaning that distributivity works except if an even number is partitioned into a sum of two odd numbers. Only when i=2, the odd-even-dependent A319929 term disappears and normal distributivity holds.
U(0;n,k) = A319929(n,k);
U(1;n,k) = A322630(n,k);
U(2;n,k) = n*k;
U(3;n,k) = A322744(n,k);
U(4;n,k) = A327259(n,k);
U(i;n,k) = 2i*floor(n/2)*floor(k/2) + A319929(n,k).
Row 1 is 2p-1 where p is a prime number (A076274 without 1).
Row 2 is the prime numbers.
Row 3 is A307002.
Row 4 is A327261.
The i-th row of T(n,k) consists of numbers that sieve out of the array U(i;n,k) = (i*n*k - (i-2)*A319929(n,k))/2, in numerical order.
From David Lovler, Sep 02 2020: (Start)
Row 1 has no even numbers. Row 2 has one even number. Generally, the even numbers of the i-th row start with i-1 consecutive even numbers (from 2). This is because U(i;2,2) = 2*i gives the first even number not in row i.
Row 3 seems to have even numbers that, after 2, coincide with A112774 which has an infinite number of terms. For i > 3, as i increases, row i has a denser presence of even numbers, thus each row has an infinite number of even terms.
Generalization of the twin prime conjecture: Since row 2 is the prime numbers, we can observe the twin prime conjecture that after the first three odd primes, the sprinkling of pairs of consecutive prime numbers never ends. Concerning just odd terms, a similar conjecture can be stated for rows i >= 3. Row 3 starts with four odd numbers then the sprinkling of three consecutive odd number never ends. Row 4 starts with five odd numbers then the sprinkling of four consecutive odd numbers never ends. The pattern continues as row i starts with i+1 odd numbers then the sprinkling of i consecutive odd numbers never ends. We can take this back to row 1 which starts with two odd numbers then continues with isolated odd numbers.
Studying the even terms, there is an analog to the above generalization of the twin prime conjecture. Row 3 starts with two even numbers then continues with isolated even numbers. Row 4 starts with three even numbers then the sprinkling of pairs of consecutive even numbers never ends. Row 5 starts with four even numbers then the sprinkling of three consecutive even numbers never ends. The pattern continues as row i starts with i-1 even numbers then the sprinkling of i-2 consecutive even numbers never ends.
(End)

Examples

			3  5  9  13  21  25  33  37  45  57  61  73  81  85  93 105 117 121 133 141 145 ...
2  3  5   7  11  13  17  19  23  29  31  37  41  43  47  53  59  61  67  71  73 ...
2  3  4   5   7   9  10  13  15  17  21  22  23  25  29  31  34  37  39  41  45 ...
2  3  4   5   6   7   9  11  12  14  15  17  19  21  22  25  27  28  29  31  35 ...
2  3  4   5   6   7   8   9  11  13  14  16  17  18  19  21  23  25  26  28  29 ...
2  3  4   5   6   7   8   9  10  11  13  15  16  18  19  20  21  22  23  25  27 ...
2  3  4   5   6   7   8   9  10  11  12  13  15  17  18  20  21  22  23  24  25 ...
2  3  4   5   6   7   8   9  10  11  12  13  14  15  17  19  20  22  23  24  25 ...
2  3  4   5   6   7   8   9  10  11  12  13  14  15  16  17  19  21  22  24  25 ...
2  3  4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19  21  23  24 ...
...
		

Crossrefs

Programs

  • Mathematica
    row=12;max=200;U[i_,n_,k_]:=(i*n*k-(i-2)If[OddQ@n,If[OddQ@k,n+k-1,k],If[OddQ@k,n,0]])/2;t=Table[c=Union@Flatten@Table[U[i,n,k],{n,2,max},{k,2,max}];Complement[Range[2,max],c][[;;row]],{i,row}];Flatten@Table[t[[m,k-m+1]],{k,row},{m,k}] (* Giorgos Kalogeropoulos, Jun 08 2021 *)

Formula

With one exception there are likely no formulas for the rows of T(n,k) since their creation is based on a sieving process like the familiar prime number sieve. The exception is T(1,k) = 2*T(2,k)-1.

A319929 Minimal arithmetic table similar to multiplication with different rules for odd and even products, read by antidiagonals.

Original entry on oeis.org

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

Views

Author

David Lovler, Dec 17 2018

Keywords

Comments

This table is akin to multiplication in that it is associative, 1 is the identity and 0 takes any number to 0. Associativity is proved by checking eight cases of three ordered odd and even numbers. Distributivity works except if an even number is partitioned into a sum of two odd numbers.

Examples

			T(3,5) = 3 + 5 - 1 = 7, T(4,7) = 4, T(8,8) = 0.
Array T(n,k) begins:
   1  2  3  4  5  6  7  8  9 10
   2  0  2  0  2  0  2  0  2  0
   3  2  5  4  7  6  9  8 11 10
   4  0  4  0  4  0  4  0  4  0
   5  2  7  4  9  6 11  8 13 10
   6  0  6  0  6  0  6  0  6  0
   7  2  9  4 11  6 13  8 15 10
   8  0  8  0  8  0  8  0  8  0
   9  2 11  4 13  6 15  8 17 10
  10  0 10  0 10  0 10  0 10  0
		

Crossrefs

Programs

  • Mathematica
    Table[Function[n, If[OddQ@ n, If[OddQ@ k, n + k - 1, k], If[OddQ@ k, n, 0]]][m - k + 1], {m, 12}, {k, m}] // Flatten (* Michael De Vlieger, Mar 24 2019 *)
  • PARI
    T(n,k) = if (n%2, if (k%2, n+k-1, k), if (k%2, n, 0));
    matrix(6, 6, n, k, T(n,k)) \\ Michel Marcus, Dec 22 2018

Formula

T(n,k) = n + k - 1 if n is odd and k is odd;
T(n,k) = n if n is even and k is odd;
T(n,k) = k if n is odd and k is even;
T(n,k) = 0 if n is even and k is even.

A327259 Array T(n,k) = 2*n*k - A319929(n,k), n >= 1, k >= 1, read by antidiagonals.

Original entry on oeis.org

1, 2, 2, 3, 8, 3, 4, 10, 10, 4, 5, 16, 13, 16, 5, 6, 18, 20, 20, 18, 6, 7, 24, 23, 32, 23, 24, 7, 8, 26, 30, 36, 36, 30, 26, 8, 9, 32, 33, 48, 41, 48, 33, 32, 9, 10, 34, 40, 52, 54, 54, 52, 40, 34, 10, 11, 40, 43, 64, 59, 72, 59, 64, 43, 40, 11, 12, 42, 50, 68, 72, 78, 78, 72, 68, 50, 42, 12
Offset: 1

Views

Author

David Lovler, Aug 27 2019

Keywords

Comments

Associative multiplication-like table whose values depend on whether n and k are odd or even.
Associativity is proved by checking the formula with eight cases of three odd and even arguments. T(n,k) is distributive as long as partitioning an even number into two odd numbers is not allowed.
T(n,k) has the same group structure as A319929, A322630 and A322744. For those arrays, position (3,3) is 5, 7 and 11 respectively. T(3,3) = 13. If we didn't have the formula for these arrays, their entries could be computed knowing one position and applying the arithmetic rules.

Examples

			Array T(n,k) begins:
   1   2   3   4   5   6   7   8   9  10
   2   8  10  16  18  24  26  32  34  40
   3  10  13  20  23  30  33  40  43  50
   4  16  20  32  36  48  52  64  68  80
   5  18  23  36  41  54  59  72  77  90
   6  24  30  48  54  72  78  96 102 120
   7  26  33  52  59  78  85 104 111 130
   8  32  40  64  72  96 104 128 136 160
   9  34  43  68  77 102 111 136 145 170
  10  40  50  80  90 120 130 160 170 200
		

Crossrefs

Equals A322744 + A322630 - A319929.
Equals 4*A322630 - 3*A319929.
0 and diagonal is A354595.

Programs

  • Mathematica
    T[n_,k_]:=2n*k-If[Mod[n,2]==1,If[Mod[k,2]==1,n+k-1,k],If[Mod[k,2]==1,n,0]]; MatrixForm[Table[T[n,k],{n,1,10},{k,1,10}]] (* Stefano Spezia, Sep 05 2019 *)
  • PARI
    T(n,k) = 2*n*k - if (n%2, if (k%2, n+k-1, k), if (k%2, n, 0));
    matrix(8, 8, n, k, T(n,k)) \\ Michel Marcus, Sep 04 2019

Formula

T(n,k) = 2*n*k - n - k + 1 if n is odd and k is odd;
T(n,k) = 2*n*k - n if n is even and k is odd;
T(n,k) = 2*n*k - k if n is odd and k is even;
T(n,k) = 2*n*k if n is even and k is even.
T(n,k) = 8*floor(n/2)*floor(k/2) + A319929(n,k).
T(n,n) = A354595(n). - David Lovler, Jul 09 2022
Writing T(n,k) as (4*n*k - 2*A319929(n,k))/2 shows that the array is U(4;n,k) of A327263. - David Lovler, Jan 15 2022

A307002 Numbers > 1 not of the form (3n*k - A319929(n,k))/2 where n and k > 1.

Original entry on oeis.org

2, 3, 4, 5, 7, 9, 10, 13, 15, 17, 21, 22, 23, 25, 29, 31, 34, 37, 39, 41, 45, 46, 49, 53, 55, 57, 58, 63, 65, 69, 71, 73, 77, 79, 81, 82, 85, 93, 94, 95, 97, 101, 105, 106, 109, 111, 118, 119, 121, 125, 129, 133, 135, 137, 141, 142, 143, 149, 151, 153, 157
Offset: 1

Views

Author

David Lovler, Mar 19 2019

Keywords

Comments

All even terms > 2 appear to be semiprimes of the form 6m+4 (A112774).
The subsequence of odd terms is A307001. - David Lovler, Jan 17 2022

Crossrefs

Third row of array A327263.

Programs

  • PARI
    T319929(n, k) = if (n%2, if (k%2, n+k-1, k), if (k%2, n, 0));
    T(n, k) = (3*n*k - T319929(n, k))/2; \\ A322744
    lista(nn) = {my(list = List()); for (n=2, nn, for (k=2, nn\n, listput(list, T(n, k)););); setminus([2..nn], Set(list));} \\ Michel Marcus, Jan 24 2021

Extensions

Name amended by David Lovler, Jan 25 2022

A340748 Numbers m > 3 such that m-1, m, m+1 belong to A307002.

Original entry on oeis.org

4, 22, 94, 142, 262, 334, 358, 694, 862, 934, 1174, 1678, 1822, 2182, 2854, 3022, 3862, 3958, 4054, 4702, 4894, 5062, 5398, 5854, 6022, 6238, 6382, 6694, 7534, 7558, 7822, 8038, 8422, 9502, 9934, 10078, 10342, 10558, 11062, 11758, 12574, 12622, 13942, 14038, 14254, 14374, 15094, 16438, 16462
Offset: 1

Views

Author

David Lovler, Jan 30 2021

Keywords

Comments

All terms m == 1 (mod 3). To prove this, look at gaps of 3 in row 2 of array A322744(n,k). The longest strings of consecutive numbers not in A322744(n,k) can occur only for these 3 numbers. The number following such a gap is A322744(2,e) = (3*2*e)/2 = 3e for some even e. The middle of the string, m = A322744(2,e) - 2 = 3e - 2. Thus m == 1 (mod 3). After the first term, all terms m == 2 (mod 4).- David Lovler, Nov 29 2021

Crossrefs

Programs

  • PARI
    T319929(n, k) = if (n%2, if (k%2, n+k-1, k), if (k%2, n, 0));
    T(n, k) = (3*n*k - T319929(n, k))/2; \\ A322744
    list(nn) = {my(list = List()); for (n=2, nn, for (k=2, nn\n, listput(list, T(n, k)); ); ); setminus([1..nn], Set(list)); } \\ A307002
    lista(nn) = {my(v=Vec(list(nn))); for (m=4, #v-1, my(x=v[m]); if (vecsearch(v, x-1) && vecsearch(v, x+1), print1(x, ", ")););} \\ Michel Marcus, Apr 02 2021

A354594 a(n) = n^2 + 2*floor(n/2)^2.

Original entry on oeis.org

0, 1, 6, 11, 24, 33, 54, 67, 96, 113, 150, 171, 216, 241, 294, 323, 384, 417, 486, 523, 600, 641, 726, 771, 864, 913, 1014, 1067, 1176, 1233, 1350, 1411, 1536, 1601, 1734, 1803, 1944, 2017, 2166, 2243, 2400, 2481, 2646, 2731, 2904
Offset: 0

Views

Author

David Lovler, Jun 01 2022

Keywords

Comments

The first bisection is A033581, the second bisection is A080859. - Bernard Schott, Jun 07 2022

Crossrefs

Programs

  • Mathematica
    a[n_] := n^2 + 2 Floor[n/2]^2
    Table[a[n], {n, 0, 90}]    (* A354594 *)
    LinearRecurrence[{1, 2, -2, -1, 1}, {0, 1, 6, 11, 24}, 60]
  • PARI
    a(n) = n^2 + 2*(n\2)^2;

Formula

a(n) = a(n-1) + 2*a(n-2) - 2*a(n-3) - a(n-4) + a(n-5), n >= 5.
a(n) = A000290(n) + 2*A008794(n).
G.f.: x*(1 + 5*x + 3*x^2 + 3*x^3)/((1 - x)^3*(1 + x)^2).
E.g.f.: (x*(1 + 3*x)*cosh(x) + (1 + 3*x + 3*x^2)*sinh(x))/2. - Stefano Spezia, Jun 07 2022

A354596 Array T(n,k) = k^2 + (2n-4)*floor(k/2)^2, n >= 0, k >= 0, read by descending antidiagonals.

Original entry on oeis.org

0, 1, 0, 0, 1, 0, 5, 2, 1, 0, 0, 7, 4, 1, 0, 9, 8, 9, 6, 1, 0, 0, 17, 16, 11, 8, 1, 0, 13, 18, 25, 24, 13, 10, 1, 0, 0, 31, 36, 33, 32, 15, 12, 1, 0, 17, 32, 49, 54, 41, 40, 17, 14, 1, 0, 0, 49, 64, 67, 72, 49, 48, 19, 16, 1, 0, 21, 50, 81, 96, 85, 90, 57, 56, 21, 18, 1, 0
Offset: 0

Views

Author

David Lovler, Jun 01 2022

Keywords

Comments

Column k is an arithmetic progression with difference 2*A008794(k).
Odd rows of A133728 triangle are contained in row 0.
For i = 0 through 4, row i is 0 and the diagonal of A319929, A322630 = A213037, A003991, A322744, and A327259, respectively. In general, row i is 0 and the diagonal of array U(i;n,k) described in A327263.

Examples

			T(n,k) begins:
  0,   1,   0,   5,   0,   9,   0,  13, ...
  0,   1,   2,   7,   8,  17,  18,  31, ...
  0,   1,   4,   9,  16,  25,  36,  49, ...
  0,   1,   6,  11,  24,  33,  54,  67, ...
  0,   1,   8,  13,  32,  41,  72,  85, ...
  0,   1,  10,  15,  40,  49,  90, 103, ...
  0,   1,  12,  17,  48,  57, 108, 121, ...
  ...
		

Crossrefs

Programs

  • Mathematica
    T[n_, k_] := k^2 + (2*n - 4)*Floor[k/2]^2; Table[T[n - k, k], {n, 0, 10}, {k, n, 0, -1}] // Flatten (* Amiram Eldar, Jun 20 2022 *)
  • PARI
    T(n,k) = k^2 + (2*n-4)*(k\2)^2;

Formula

T(n,k) = U(n;k,k) (see A327263).
For each row, T(n,k) = T(n,k-1) + 2*T(n,k-2) - 2*T(n,k-3) - T(n,k-4) + T(n,k-5), k >= 5.
G.f. for row n: x*(1 + (2*n-1)*x + 3*x^2 + (2*n-3)*x^3)/((1 - x)^3*(1 + x)^2). When n = 2, this reduces to x*(1 + x)/(1 - x)^3.
E.g.f. for row n: (((4-n)*x + n*x^2)*cosh(x) + (n-2 + n*x + n*x^2)*sinh(x))/2. When n = 2, this reduces to (x + x^2)*cosh(x) + (x + x^2)*sinh(x) = (x + x^2)*exp(x).

A307001 Odd numbers > 1 not of the form (3n*k - n - k + 1)/2 where n and k are odd numbers > 1.

Original entry on oeis.org

3, 5, 7, 9, 13, 15, 17, 21, 23, 25, 29, 31, 37, 39, 41, 45, 49, 53, 55, 57, 63, 65, 69, 71, 73, 77, 79, 81, 85, 93, 95, 97, 101, 105, 109, 111, 119, 121, 125, 129, 133, 135, 137, 141, 143, 149, 151, 153, 157, 161, 169, 175, 177, 181, 183, 185, 189, 193, 197
Offset: 1

Views

Author

David Lovler, Mar 19 2019

Keywords

Comments

Terms are the odd numbers > 1 not appearing in array A322744 with its first row and column omitted.
They are the odd numbers in A307002. - David Lovler, Jan 17 2022

Crossrefs

Programs

  • Mathematica
    Select[2 Range[100]-1, FindInstance[# == 1 + 2*n + k (2 + 6 n) && n>0 && k>0, {n,k}, Integers] === {} &] (* Giovanni Resta, May 06 2019 *)
  • PARI
    isok(n) = {my(kj, tij); if (n % 2, forstep (i=3, oo, 2, kj = 0; forstep (j=3, i, 2, tij = (3*i*j - i - j +1)/2; if (tij == n, return (0)); if (tij > n, kj = j; break);); if ((kj == 3) && (tij > n), break);); return (n>1));} \\ Michel Marcus, Apr 24 2019 and Jan 25 2022

Extensions

Definition amended by David Lovler, Jan 25 2022
Showing 1-10 of 12 results. Next