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-3 of 3 results.

A223701 Irregular triangle of numbers k such that prime(n) is the largest prime factor of k^2 - 1.

Original entry on oeis.org

3, 2, 5, 7, 17, 4, 9, 11, 19, 26, 31, 49, 161, 6, 8, 13, 15, 29, 41, 55, 71, 97, 99, 127, 244, 251, 449, 4801, 8749, 10, 21, 23, 34, 43, 65, 76, 89, 109, 111, 197, 199, 241, 351, 485, 769, 881, 1079, 6049, 19601, 12, 14, 25, 27, 51, 53, 64, 79, 129, 131, 155
Offset: 1

Views

Author

T. D. Noe, Apr 03 2013

Keywords

Comments

Note that the first number of each row forms the sequence 3, 2, 4, 6, 10, 12,..., which is A039915. The first 25 rows, except the first, are in A181447-A181470.

Examples

			Irregular triangle:
  {3},
  {2, 5, 7, 17},
  {4, 9, 11, 19, 26, 31, 49, 161},
  {6, 8, 13, 15, 29, 41, 55, 71, 97, 99, 127, 244, 251, 449, 4801, 8749}
		

Crossrefs

Row 26 is A181568.
Cf. A039915 (first terms), A175607 (last terms), A181471 (row lengths), A379344 (row sums).
Cf. A223702, A223703, A223704 (related tables).

Programs

  • Mathematica
    t = Table[FactorInteger[n^2 - 1][[-1,1]], {n, 2, 10^5}]; Table[1 + Flatten[Position[t, Prime[n]]], {n, 6}]

A181452 Numbers k such that 17 is the largest prime factor of k^2 - 1.

Original entry on oeis.org

16, 33, 35, 50, 67, 69, 101, 103, 118, 120, 169, 188, 239, 271, 307, 339, 441, 511, 545, 577, 749, 883, 1121, 1189, 1376, 1429, 1665, 1871, 2024, 2177, 2311, 2449, 2549, 3401, 4115, 4861, 4999, 5201, 9827, 11663, 24751, 28799, 57121, 62425, 74359, 388961, 672281
Offset: 1

Views

Author

Artur Jasinski, Oct 21 2010

Keywords

Comments

Numbers k such that A076605(k) = 17.
Sequence is finite, for proof see A175607.
Search for terms can be restricted to the range from 2 to A175607(7) = 672281; primepi(17) = 7.

Crossrefs

Programs

  • Magma
    [ n: n in [2..350000] | m eq 17 where m is D[#D] where D is PrimeDivisors(n^2-1) ]; // Klaus Brockhaus, Feb 19 2011
    
  • Magma
    p:=(97*89*83*79*73*71)^5 *(67*61*59*53*47*43*41)^6 *(37*31*29)^7 *(23*19*17)^8 *13^9 *11^10 *7^13 *5^15 *3^23 *2^36; [ n: n in [2..700000] | p mod (n^2-1) eq 0 and (D[#D] eq 17 where D is PrimeDivisors(n^2-1)) ]; // Klaus Brockhaus, Feb 24 2011
    
  • Mathematica
    jj = 2^36*3^23*5^15*7^13*11^10*13^9*17^8*19^8*23^8*29^7*31^7*37^7*41^6 *43^6*47^6*53^6*59^6*61^6*67^6*71^5*73^5*79^5*83^5*89^5*97^5; rr = {}; n = 2; While[n < 700000, If[GCD[jj, n^2 - 1] == n^2 - 1, k = FactorInteger[n^2 - 1]; kk = Last[k][[1]]; If[kk == 17, AppendTo[rr, n]]]; n++ ]; rr (* Artur Jasinski *)
    Select[Range[680000], FactorInteger[#^2-1][[-1, 1]]==17&]
  • PARI
    is(n)=n=n^2-1; forprime(p=2, 13, n/=p^valuation(n, p)); n>1 && 17^valuation(n, 17)==n \\ Charles R Greathouse IV, Jul 01 2013

A181454 Numbers k such that 23 is the largest prime factor of k^2 - 1.

Original entry on oeis.org

22, 24, 45, 47, 91, 116, 137, 139, 183, 208, 229, 254, 298, 321, 323, 344, 415, 461, 505, 551, 599, 645, 781, 783, 919, 967, 1013, 1057, 1126, 1151, 1310, 1471, 1519, 1749, 1793, 2186, 2209, 2276, 2393, 2575, 2874, 2991, 3704, 3725, 4047, 4049, 4369
Offset: 1

Views

Author

Artur Jasinski, Oct 21 2010

Keywords

Comments

Numbers k such that A076605(k) = 23.
Sequence is finite, for proof see A175607.
Search for terms can be restricted to the range from 2 to A175607(9) = 10285001; primepi(23) = 9.

Crossrefs

Programs

  • Magma
    [ n: n in [2..300000] | m eq 23 where m is D[#D] where D is PrimeDivisors(n^2-1) ]; // Klaus Brockhaus, Feb 18 2011
    
  • Magma
    p:=(97*89*83*79*73*71)^5 *(67*61*59*53*47*43*41)^6 *(37*31*29)^7 *(23*19*17)^8 *13^9 *11^10 *7^13 *5^15 *3^23 *2^36; [ n: n in [2..10300000] | p mod (n^2-1) eq 0 and (D[#D] eq 23 where D is PrimeDivisors(n^2-1)) ]; // Klaus Brockhaus, Feb 24 2011
    
  • Mathematica
    jj=2^36*3^23*5^15*7^13*11^10*13^9*17^8*19^8*23^8*29^7*31^7*37^7*41^6 *43^6*47^6*53^6*59^6*61^6*67^6*71^5*73^5*79^5*83^5*89^5*97^5; rr ={};n = 2; While[n < 14000000, If[GCD[jj, n^2 - 1] == n^2 - 1, k = FactorInteger[n^2 - 1]; kk = Last[k][[1]]; If[kk == 23, AppendTo[rr, n]]]; n++ ]; rr
    Select[Range[300000], FactorInteger[#^2-1][[-1, 1]]==23&]
  • PARI
    is(n)=n=n^2-1; forprime(p=2, 19, n/=p^valuation(n, p)); n>1 && 23^valuation(n, 23)==n \\ Charles R Greathouse IV, Jul 01 2013
Showing 1-3 of 3 results.