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.

Previous Showing 51-60 of 60 results.

A176850 a(n,k) is the number of ways to choose integers i,j from {0,1,...,k} such that the inequalities |i-j| <= n <= i+j are satisfied.

Original entry on oeis.org

1, 2, 3, 1, 3, 6, 6, 3, 1, 4, 9, 11, 10, 6, 3, 1, 5, 12, 16, 17, 15, 10, 6, 3, 1, 6, 15, 21, 24, 24, 21, 15, 10, 6, 3, 1, 7, 18, 26, 31, 33, 32, 28, 21, 15, 10, 6, 3, 1, 8, 21, 31, 38, 42, 43, 41, 36, 28, 21, 15, 10, 6, 3, 1, 9, 24, 36, 45, 51, 54, 54, 51, 45, 36, 28, 21, 15, 10, 6, 3, 1, 10, 27, 41, 52, 60, 65, 67, 66, 62, 55, 45, 36, 28, 21, 15, 10, 6, 3
Offset: 0

Views

Author

Sean Murray, Apr 27 2010

Keywords

Comments

The rows are of length 1, 3, 5, 7, ...
a(n,k) is also the number of independent rank n tensor operators to appear in the tensor product of two spaces each spanned by k+1 tensor operators of ranks 0 to k,
{Y_{l,m},l=0,1,...,k, m:-l,-l+1,...,l} times {Y'_{l'm'}, l'=0,1,...,k, m':-l,-l+1,...,l}.
Basis elements of the tensor product space are given by
psi^{l,l'}{p,q} = Sum{m,m'} C^{ll'p}{mm'q} Y{l,m}Y'_{l'm'}
for all l,l' = 0,1,...,k and where p = |l-l'|, |l-l'|+1, ..., l+l' is the rank, q=-p, -p+1,...,p and where C^{ll'p}_{mm'q} are the Clebsch-Gordon coefficients.
Sum_{k=0..2*n+1} a(n,k)*(2*k+1) = (n+1)^4. - L. Edson Jeffery, Oct 29 2012
Sum_{k=0..2*n+1} (a(n,k) - a(n-1,k))*(2*k+1) = n^4 - (n-1)^4 = A005917(n+1), for n > 0. - L. Edson Jeffery, Nov 02 2012

Examples

			Triangle begins
   1;
   2,  3,  1;
   3,  6,  6,  3,  1;
   4,  9, 11, 10,  6,  3,  1;
   5, 12, 16, 17, 15, 10,  6,  3,  1;
   6, 15, 21, 24, 24, 21, 15, 10,  6,  3,  1;
   7, 18, 26, 31, 33, 32, 28, 21, 15, 10,  6,  3,  1;
   8, 21, 31, 38, 42, 43, 41, 36, 28, 21, 15, 10,  6,  3,  1;
   9, 24, 36, 45, 51, 54, 54, 51, 45, 36, 28, 21, 15, 10,  6,  3,  1;
  10, 27, 41, 52, 60, 65, 67, 66, 62, 55, 45, 36, 28, 21, 15, 10,  6,  3,  1;
		

Crossrefs

Cf. A005917.

Programs

  • Maple
    Seq:=[]: for k from 0 to 15 do for n from 0 to k do Seq:= [op(Seq), -(3/2)*n^2+2*k*n+(1/2)*n+k+1] end do; for n from k+1 to 2*k do Seq:= [op(Seq), (1/2)*(2*k-n+1)*(2*k-n+2)] end do; end do; Seq;
  • Mathematica
    Table[If[n <= k, -(3/2)*n^2 + 2*k*n + n/2 + k + 1, (2*k - n + 1)*(2*k - n + 2)/2], {k, 0, 8}, {n, 0, 2 k}] // Flatten (* Michael De Vlieger, Jul 10 2022 *)

Formula

a(n,k) = -(3/2)*n^2 + 2*k*n + n/2 + k + 1 for n=0,1,...,k, a(n) = (2*k-n+1)*(2*k-n+2)/2 for n = k+1,...,2*k.

Extensions

Edited by Sean Murray, Oct 05 2011

A194359 Triangle of divisors of 210^n, each number occurring once.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 14, 15, 21, 30, 35, 42, 70, 105, 210, 4, 9, 12, 18, 20, 25, 28, 36, 45, 49, 50, 60, 63, 75, 84, 90, 98, 100, 126, 140, 147, 150, 175, 180, 196, 225, 245, 252, 294, 300, 315, 350, 420, 441, 450, 490, 525, 588, 630, 700, 735, 882, 900
Offset: 0

Views

Author

T. D. Noe, Aug 26 2011

Keywords

Comments

The length of row k is A005917, the rhombic dodecahedral numbers, (k+1)^4 - k^4. The triangle has rows beginning with 2^k and ending with 210^k.

Crossrefs

Programs

  • Mathematica
    Join[{{1}}, Table[Complement[Divisors[210^n], Divisors[210^(n-1)]], {n, 9}]]
    Take[DeleteDuplicates[Flatten[Divisors/@(210^Range[5])]],100] (* Harvey P. Dale, Sep 03 2020 *)

A272297 a(n) = n^4 + 64.

Original entry on oeis.org

64, 65, 80, 145, 320, 689, 1360, 2465, 4160, 6625, 10064, 14705, 20800, 28625, 38480, 50689, 65600, 83585, 105040, 130385, 160064, 194545, 234320, 279905, 331840, 390689, 457040, 531505, 614720, 707345, 810064, 923585, 1048640, 1185985, 1336400, 1500689, 1679680, 1874225, 2085200
Offset: 0

Views

Author

Bruno Berselli, Apr 25 2016

Keywords

Comments

This is the case k=2 of Sophie Germain's Identity n^4+(2*k^2)^2 = ((n-k)^2+k^2)*((n+k)^2+k^2).

Crossrefs

Cf. A005917.
Subsequence of A227855.
Cf. A000583 (k=0), A057781 (k=1), A272298 (k=3).

Programs

  • Magma
    [n^4+64: n in [0..40]];
    
  • Mathematica
    Table[n^4 + 64, {n, 0, 40}]
  • Maxima
    makelist(n^4+64, n, 0, 40);
    
  • PARI
    vector(40, n, n--; n^4+64)
    
  • Python
    [n**4+64 for n in range(40)]
    
  • Python
    for n in range(0,10**5):print(n**4+64) # Soumil Mandal, Apr 30 2016
  • Sage
    [n^4+64 for n in (0..40)]
    

Formula

O.g.f.: (64 - 255*x + 395*x^2 - 245*x^3 + 65*x^4)/(1 - x)^5.
E.g.f.: (64 + x + 7*x^2 + 6*x^3 + x^4)*exp(x).
a(n) = (n^2 - 8)^2 + (4*n)^2.

A341050 Cube array read by upward antidiagonals ignoring zero and empty terms: T(n, k, r) is the number of n-ary strings of length k, containing r consecutive 0's.

Original entry on oeis.org

1, 1, 1, 3, 1, 1, 3, 1, 5, 8, 1, 1, 3, 1, 5, 8, 1, 7, 21, 19, 1, 1, 3, 1, 5, 8, 1, 7, 21, 20, 1, 9, 40, 81, 43, 1, 1, 3, 1, 5, 8, 1, 7, 21, 20, 1, 9, 40, 81, 47, 1, 11, 65, 208, 295, 94, 1, 1, 3, 1, 5, 8, 1, 7, 21, 20, 1, 9, 40, 81, 48, 1, 11, 65, 208, 297, 107, 1, 13, 96, 425, 1024, 1037, 201
Offset: 2

Views

Author

Robert P. P. McKone, Feb 04 2021

Keywords

Examples

			For n = 5, k = 6 and r = 4, there are 65 strings: {000000, 000001, 000002, 000003, 000004, 000010, 000011, 000012, 000013, 000014, 000020, 000021, 000022, 000023, 000024, 000030, 000031, 000032, 000033, 000034, 000040, 000041, 000042, 000043, 000044, 010000, 020000, 030000, 040000, 100000, 100001, 100002, 100003, 100004, 110000, 120000, 130000, 140000, 200000, 200001, 200002, 200003, 200004, 210000, 220000, 230000, 240000, 300000, 300001, 300002, 300003, 300004, 310000, 320000, 330000, 340000, 400000, 400001, 400002, 400003, 400004, 410000, 420000, 430000, 440000}
The first seven slices of the tetrahedron (or pyramid) are:
-----------------Slice 1-----------------
  1
-----------------Slice 2-----------------
    1
  1  3
-----------------Slice 3-----------------
      1
    1  3
  1  5  8
-----------------Slice 4-----------------
        1
      1  3
    1  5   8
  1  7  21  19
-----------------Slice 5-----------------
          1
        1  3
      1  5   8
    1  7  21  20
  1  9  40  81  43
-----------------Slice 6-----------------
              1
           1    3
        1    5     8
      1   7    21    20
    1   9   40    81    47
  1  11  65   208   295   94
-----------------Slice 7-----------------
                 1
              1     3
           1     5     8
         1    7     21    20
      1    9    40     81      48
    1   11   65    208     297     107
  1  13   96   425    1024    1037    201
		

Crossrefs

Cf. A340156 (r=2), A340242 (r=3).
Cf. A008466 (n=2, r=2), A186244 (n=3, r=2), A050231 (n=2, r=3), A231430 (n=3, r=3).
Cf. A000567 [(k=4, r=2),(k=5, r=3),(k=6, r=4),...,(k=x, r=x-2)].
Cf. A103532 [(k=6, r=3),(k=7, r=4),(k=8, r=5),...,(k=x, r=x-3)].

Programs

  • Mathematica
    m[r_, n_] := Normal[With[{p = 1/n}, SparseArray[{Band[{1, 2}] -> p, {i_, 1} /; i <= r -> 1 - p, {r + 1, r + 1} -> 1}]]]; T[n_, k_, r_] := MatrixPower[m[r, n], k][[1, r + 1]]*n^k; DeleteCases[Transpose[PadLeft[Reverse[Table[T[n, k, r], {k, 2, 8}, {r, 2, k}, {n, 2, r}], 2]], 2 <-> 3], 0, 3] // Flatten

A346263 Irregular triangle of numbers T(n,k) read by rows where each T(n,k) is the number of large or small squares that are used to tile elementary squares of type 2 whose length of side is A344332(n).

Original entry on oeis.org

9, 9, 9, 9, 25, 9, 9, 9, 9, 25, 9, 9, 9, 49, 9, 9, 25, 9, 9, 9, 9, 25, 9, 9, 9, 9, 25, 9, 9, 49, 9, 81, 9, 9, 25, 9, 9, 9, 9, 25, 9, 9, 9, 9, 25, 9, 49, 9, 9, 9, 9, 25, 9, 9, 9, 9, 25, 9, 121, 9, 9, 49, 9, 25, 9, 9, 81, 9, 9, 9, 25, 9, 9, 9, 9, 25, 9, 9, 49, 9, 9
Offset: 1

Views

Author

Bernard Schott, Jul 13 2021

Keywords

Comments

An elementary square of type 2 is the smallest square that can be tiled with squares of two different sides a < b satisfying a^2+b^2 = c^2 and so that the numbers of small and large squares are equal.
Every term is an odd square >= 9 and each odd square is present infinitely many times.
Notation: s_p (resp. s_e) = side of a primitive (resp. elementary) tiled square, a = side of small squares and b = side of large squares used to tile a primitive square, and z_p (z_e) = number of small squares = number of large squares used to tile a primitive (resp. elementary) square.
A primitive square with side s_p = a*c/(c-b) is tiled with z_p small and z_p large squares with sides a and b, and z_p = (a/(c-b))^2.
Each elementary square with a side s_e = k*s_p, k>0, is tiled with z_e small and z_e large squares with sides k*a and k*b, and z_e = z_p = (a/(c-b))^2.
When an elementary side A344332(n) is a multiple of m distinct primitive sides s_p, then there are m different values T(n,1), ..., T(n,m) in the row n (see example).

Examples

			The triangle T begins:
   n\k 1    2    3    4    5
   1:  9
   2:  9
   3:  9
   4:  9
   5: 25
   6:  9
   7:  9
   8:  9
   9:  9
  10: 25
  11:  9
  12:  9
  13:  9
  14: 49
  15:  9
  16:  9,   25
  17:  9
  ...
The first elementary side that is a multiple of two primitive sides (15 and 65) is A344332(16) = 195 = 13*15 = 3*65.
As 195 = 13*15, the number z of squares with side a = 13*3 = 39 and b = 13*4 = 52 to tile this elementary square is T(16,1) = (39/(65-52))^2 = 9.
As 195 = 3 * 65, the number z of squares with side a = 3*5 = 15 and b = 3*12 = 36 to tile this elementary square is T(16,2) = (15/(39-36))^2 = 25.
Hence, the elementary square with side A344332(16) = 195 has two different possible tilings: with T(16,1) = 9 squares of sides (a,b) = (39,52) or with T(16,2) = 25 squares of sides (a,b) = (15,36).
Elementary square 195 X 195 with a = 39, b = 52, s = 195, z = 9:
     ________ ________ ________ _____
    |        |        |        |     |
    |        |        |        |     |
    |        |        |        |_____|
    |________|________|________|     |
    |        |        |        |     |
    |        |        |        |_____|
    |        |        |        |     |
    |________|________|________|     |
    |        |        |        |_____|
    |        |        |        |     |
    |        |        |        |     |
    |_____ __|___ ____|_ ______|_____|
    |     |      |      |      |     |
    |     |      |      |      |     |
    |_____|______|______|______|_____|
		

Crossrefs

Cf. A345286 (similar for type 1).

A365206 Centered octachoral numbers.

Original entry on oeis.org

1, 49, 337, 1249, 3361, 7441, 14449, 25537, 42049, 65521, 97681, 140449, 195937, 266449, 354481, 462721, 594049, 751537, 938449, 1158241, 1414561, 1711249, 2052337, 2442049, 2884801, 3385201, 3948049, 4578337, 5281249, 6062161, 6926641, 7880449, 8929537
Offset: 1

Views

Author

Léo Cymrot Cymbalista, Aug 25 2023

Keywords

Comments

A octachoral number is a centered figurate number that represents a octachoron, which is a four-dimensional regular polytope composed of 8 cells (also known as tesseract).
One of the 6 centered regular polichoral (centered pentachoral, centered hexadecachoral, centered octachoral, centered icositetrachoral, centered hexacosichoral and centered hecatonicosachoral) numbers.

Crossrefs

Programs

  • Mathematica
    Table[8*n^4 - 16*n^3 + 16*n^2 - 8n + 1, {n, 1, 100}]

Formula

a(n) = 8*n^4 - 16*n^3 + 16*n^2 - 8n + 1.
G.f.: x*(1 + 44*x + 102*x^2 + 44*x^3 + x^4)/(1 - x)^5. - Stefano Spezia, Aug 26 2023

A260550 a(n) is the number of 2 X 2 matrices with entries in {1, ..., n} that are not the product of two 2 X 2 positive integer matrices.

Original entry on oeis.org

1, 15, 75, 237, 559, 1157, 2055, 3471, 5449, 8131, 11633, 16361, 22041, 29349, 38329, 48839, 61325, 76479, 93957, 114717, 138041, 164153, 194505, 229625, 268259, 311031, 359719, 413245, 472145, 537835, 608837, 688121, 774877, 867549, 971403, 1080637, 1198233, 1326059, 1467029, 1617451, 1777881, 1948219, 2132381, 2329081, 2539351
Offset: 1

Views

Author

Aldo González Lorenzo, Jul 29 2015

Keywords

Comments

a(n) <= A000583(n), which is the number of 2 X 2 matrices with entries in {1, ..., n}.
a(n) >= A005917(n), which is the number of 2 X 2 matrices with entries in {1, ..., n} that contain the element 1. All such matrices are not decomposable as a product of 2 X 2 positive integer matrices.
This definition is a generalization of the notion of prime numbers to the family of 2 X 2 positive integer matrices. Since the matrices do not contain 0, max(A*B) > max(A) and max(A*B) > max(B). Thus, for every matrix there is a finite number of possible decompositions to check.

Examples

			The matrix [2,2;3,3] is decomposable: [2,2;3,3] = [1,1;1,2] * [1,1;1,1]. However, the matrix [2,3;3;2] is not decomposable.
		

Crossrefs

Programs

  • Python
    # See Branicky link.

A272850 a(n) = (n^2 + (n+1)^2)*(n^2 + (n+1)^2 + 2*n*(n+1)).

Original entry on oeis.org

1, 45, 325, 1225, 3321, 7381, 14365, 25425, 41905, 65341, 97461, 140185, 195625, 266085, 354061, 462241, 593505, 750925, 937765, 1157481, 1413721, 1710325, 2051325, 2440945, 2883601, 3383901, 3946645, 4576825, 5279625
Offset: 0

Views

Author

Matthew Badley, May 07 2016

Keywords

Comments

Larger of pair of integers whose Pythagorean means are all integers.
The smaller of the pairs are: (A001844).
The arithmetic means are: (A007204)
The geometric means are: (A005917)
The harmonic means are: (A016754).
Subtracting terms in A016754 from A007204 gives complementary harmonics (A060300).

Crossrefs

Programs

Formula

a(n) = (2*n^2 + 2*n + 1)*(4*n^2 + 4*n + 1).
From Colin Barker, May 24 2016: (Start)
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5) for n>4.
G.f.: (1 + 40*x + 110*x^2 + 40*x^3 + x^4) / (1-x)^5. (End)

A346265 a(n) is the number of distinct possible tilings of type 1 (A344331) or of type 2 (A344332) for a square whose side is A344330(n).

Original entry on oeis.org

1, 1, 2, 5, 3, 2, 2, 9, 1, 1, 2, 2, 1, 4, 9, 4, 2, 2, 13, 5, 3, 2, 4, 10, 2, 5, 2, 2, 1, 16, 2, 4, 6, 2, 10, 4, 1, 4, 2, 2, 17, 3, 2, 9, 13, 3, 6, 2, 3, 19, 2, 3, 4, 6, 2, 10, 6, 2, 7, 23, 2, 2, 3, 4, 18, 8, 4, 4, 2, 18, 2, 2, 6, 2, 18, 2, 4, 2, 4, 2, 2, 21, 3, 4, 6, 11, 14, 6, 2, 23
Offset: 1

Views

Author

Bernard Schott, Aug 11 2021

Keywords

Comments

These squares with side = A344330(n) can be tiled with squares of two different sizes so that the numbers of large or small squares are equal.
Notation: s = side of the tiled squares, a = side of small squares, b = side of large squares, and z = number of small squares = number of large squares.

Examples

			-> A344330(1) = A344331(1) = 10 and there is no k_2 such that A344330(1) = A344332(k_2) = 10, then a(1) = A345287(1) = 1 (example below of type 1):
   Primitive square 10 X 10 corresponding to a(1) = 1 with
    a = 1, b = 2, s = 10, z = 20:
      ___ ___ _ ___ ___ _
     |   |   |_|   |   |_|
     |___|___|_|___|___|_|
     |   |   |_|   |   |_|
     |___|___|_|___|___|_|
     |   |   |_|   |   |_|
     |___|___|_|___|___|_|
     |   |   |_|   |   |_|
     |___|___|_|___|___|_|
     |   |   |_|   |   |_|
     |___|___|_|___|___|_|
-> A344330(2) = A344332(1) = 15 and there is no k_1 such that A344330(2) = A344331(k_1) = 15, then a(2) = A346264(1) = 1 (example below of type 2):
   Primitive square 15 X 15 corresponding to a(2) = 1 with
     a = 3, b = 4, c = 5, s = 15, z = 9:
        ________ ________ ________ ______
       |        |        |        |      |
       |        |        |        |      |
       |        |        |        |______|
       |_______ |________|________|      |
       |        |        |        |      |
       |        |        |        |______|
       |        |        |        |      |
       |________|________|________|      |
       |        |        |        |______|
       |        |        |        |      |
       |        |        |        |      |
       |_____ __|___ ____|_ ______|______|
       |     |      |      |      |      |
       |     |      |      |      |      |
       |_____|______|______|______|______|
-> A344330(4) = A344331(3) = A344332(2) = 30, then a(4) = A345287(3) + A346264(2) = 3+2 = 5 (see link with the corresponding 5 distinct tilings).
-> A344330(6) = A344332(3) = 45 and there is no k_1 such that A344330(6) = A344331(k_1) = 45, then a(6) = A346264(3) = 2.
		

References

  • Ivan Yashchenko, Invitation to a Mathematical Festival, pp. 10 and 102, MSRI, Mathematical Circles Library, 2013.

Crossrefs

Programs

  • PARI
    \\ isok1 from A344331 and isok2 from A344332
    isok3(s) = {if (!(s % 2) && ispower(s/2, 4), return (0)); my(m = sqrtnint(s, 3)); vecsearch(setbinop((x, y)->if (gcd(x, y)==1, x*y*(x^2+y^2), 0), [1..m]), s); } \\ A344333
    sd7(x) = sumdiv(x, d, if (isok3(d), numdiv(x/d))); \\ A345287
    isok7(k) = my(kk= sqrtnint(k\4, 3)+2); vecsearch(vector(kk, i, (i+1)^4 - i^4), k); \\ A005917
    sd4(x) = sumdiv(x, d, if (isok7(d), numdiv(x/d))); \\ A346264
    lista(nn) = {for (n=1, nn, my(b1 = isok1(n), b2 = isok2(n)); if (b1 || b2, my(x = 0); if (b1, x += sd7(n)); if (b2, x += sd4(n)); print1(x, ", ");););} \\ Michel Marcus, Dec 23 2021

Formula

If A344330(n) = A344331(k_1) and there is no k_2 such that A344330(n) = A344332(k_2) then a(n) = A345287(k_1).
If A344330(n) = A344332(k_2) and there is no k_1 such that A344330(n) = A344331(k_1) then a(n) = A346264(k_2).
If A344330(n) = A344331(k_1) = A344332(k_2) then a(n) = A345287(k_1) + A346264(k_2).

Extensions

a(19),a(59),a(86),a(87) corrected by Bernard Schott, Dec 23 2021

A353890 a(n) is the period of the binary sequence {b(m)} defined by b(m) = 1 if (m+1)^n - m^n and (m+2)^n - 2*(m+1)^n + m^n are coprime, 0 otherwise.

Original entry on oeis.org

1, 1, 5, 11, 91, 1247, 3485, 263017, 852841, 1241058127, 74966255, 243641132605417, 181556731572385303, 718802057694183783881, 6582662048285, 943422576750791493013356207217, 487331778345355477261, 607088607861933740557075591887834842297
Offset: 2

Views

Author

Samuel Harkness, May 09 2022

Keywords

Comments

For any n, consecutive n-th powers will never share a divisor > 1, so now consider the second differences. Specifically, each m > 0, define the binary sequence {b(m)} as follows: b(m) = 1 if the first difference (m+1)^n - m^n and the second difference (m+2)^n - 2*(m+1)^n + m^n are coprime, 0 otherwise. I conjecture that {b(m)} is periodic with period a(n).
If m^n mod p == (m+1)^n mod p == (m+2)^n mod p, then p is in the prime factorization of a(n).
All primes p >= 5 belong to a prime factorization for a(n). p will always belong to the prime factorization of n=p-1 due to Fermat's Little Theorem.
I conjecture that the greatest prime factor for any prime n >= 5 is phi(2^n+1)/2 + 1 = Jacobsthal(n). n*A069925 + 1 = A001045(n).
I conjecture that all prime factors "f" are f=n*k+1, unless n is composite, in which case additionally all prime factors for any divisor of n will also be included in the prime factorization for a(n).

Examples

			For n=2 and n=3, the first and second differences are coprime for all m. Each of their sequences {b(m)} consist only of 1's, which can be described trivially as [1] with a period of 1, so a(2) = a(3) = 1.
For n > 3, the first and second differences are coprime for some m values, but not for all. Each repeating periodic sequence {b(m)} begins at m=1, and can be used to predict what b(m) will be at any higher m value for that power n.
n=4 has the 5-term repeating sequence, beginning at m=1:
  [0 0 1 1 1], so a(4) = 5.
The sequence is repeating, so for example, f(41)..f(45) is also [0 0 1 1 1].
n=5 has the 11-term repeating sequence
  [1 1 0 1 1 0 1 1 1 1 1]
so a(5) = 11.
n=6 has the 91-term repeating sequence
  [0 0 0 0 0 0 1 0 0 0 0 1 1
   1 0 0 0 0 0 1 1 0 0 0 0 1
   1 1 0 0 0 0 1 1 1 0 0 0 0
   1 1 1 0 0 0 0 1 1 1 0 0 0
   0 1 1 1 0 0 0 0 1 1 1 0 0
   0 0 1 1 0 0 0 0 0 1 1 1 0
   0 0 0 1 0 0 0 0 0 0 1 1 1]
so a(6) = 91.
The period for higher n values has yet to be found. If they exist, it seems they would be quite large given the large expansion from 5, 11, to 91.
Example: the 233rd term in the sequence of values for n=6 is calculated by using m=233 and n=6. Define the first difference for the 233rd term as 234^6 - 233^6 = 4164782373647. The second difference for the 233rd term is 235^6 - 2*234^6 + 233^6 = 89948228762. The terms 4164782373647 and 89948228762 share a common factor, so the 233rd term of the sequence for 6th powered terms is denoted 0 (not coprime). Because the 6th powered terms repeat their tendency of being coprime or not every 91 terms, we could instead look at 233 mod 91 = 51, and from the table for n=6 above, the 51st term is 0.
		

Crossrefs

Programs

  • MATLAB
    See Links section.

Extensions

a(7)-a(19) from Jon E. Schoenfield, May 10 2022
Previous Showing 51-60 of 60 results.