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

A085152 All prime factors of n and n+1 are <= 5. (Related to the abc conjecture.)

Original entry on oeis.org

1, 2, 3, 4, 5, 8, 9, 15, 24, 80
Offset: 1

Views

Author

Benoit Cloitre, Jun 21 2003

Keywords

Comments

Equivalently: Numbers n such that n(n+1) is 5-smooth.
The ABC conjecture would imply that if the prime factors of A, B, C are prescribed in advance, then there is only a finite number of solutions to the equation A + B = C with gcd(A,B,C)=1 (indeed it would bound C to be no more than "roughly" the product of those primes). So in particular there ought to be only finitely many pairs of adjacent integers whose prime factors are limited to {2, 3, 5} (D. Rusin).
This sequence is complete by a theorem of Stormer. See A002071. - T. D. Noe, Mar 03 2008
This is the 3rd row of the table A138180. It has 10 = A002071(3) = A145604(1)+A145604(2)+A145604(3) terms and ends with A002072(3) = 80. It is the union of all terms in rows 1 through 3 of the table A145605. It is a subsequence of A252494 and A085153. - M. F. Hasler, Jan 16 2015

Crossrefs

Programs

  • Mathematica
    Select[Range[10000], FactorInteger[ # (# + 1)][[ -1,1]] <= 5 &] (* T. D. Noe, Mar 03 2008 *)
  • PARI
    for(n=1,99,vecmax(factor(n++)[,1])<6 && vecmax(factor(n--+(n<2))[,1])<6 && print1(n", ")) \\ This skips 2 if n+1 is not 5-smooth: twice as fast as the naive version. - M. F. Hasler, Jan 16 2015

Extensions

Edited by Dean Hickerson, Jun 30 2003

A138180 Irregular triangle read by rows: row n consists of all numbers x such that x and x+1 have no prime factor larger than prime(n).

Original entry on oeis.org

1, 1, 2, 3, 8, 1, 2, 3, 4, 5, 8, 9, 15, 24, 80, 1, 2, 3, 4, 5, 6, 7, 8, 9, 14, 15, 20, 24, 27, 35, 48, 49, 63, 80, 125, 224, 2400, 4374, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 14, 15, 20, 21, 24, 27, 32, 35, 44, 48, 49, 54, 55, 63, 80, 98, 99, 120, 125, 175, 224, 242, 384, 440, 539
Offset: 1

Views

Author

T. D. Noe, Mar 04 2008

Keywords

Comments

A number x is p-smooth if all prime factors of x are <= p. The length of row n is A002071(n). Row n begins with 1 and ends with A002072(n). Each term of row n-1 is in row n.
The n-th row is the union of the rows 1 to n of A145605. - M. F. Hasler, Jan 18 2015

Examples

			The table reads:
1,
1, 2, 3, 8,
1, 2, 3, 4, 5, 8, 9, 15, 24, 80,  (= A085152)
1, 2, 3, 4, 5, 6, 7, 8, 9, 14, 15, 20, 24, 27, 35, 48, 49, 63, 80, 125, 224, 2400, 4374, (= A085153)
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 14, 15, 20, 21, 24, 27, 32, 35, 44, 48, 49, 54, 55, 63, 80, 98, 99, 120, 125, 175, 224, 242, 384, 440, 539, 2400, 3024, 4374, 9800 (= A252494),
...
		

References

Crossrefs

Programs

  • Mathematica
    (* This program needs x maxima taken from A002072. *) xMaxima = A002072; smoothNumbers[p_, max_] := Module[{a, aa, k, pp, iter}, k = PrimePi[p]; aa = Array[a, k]; pp = Prime[Range[k]]; iter = Table[{a[j], 0, PowerExpand @ Log[pp[[j]], max/Times @@ (Take[pp, j-1]^Take[aa, j-1])]}, {j, 1, k}]; Table[Times @@ (pp^aa), Sequence @@ iter // Evaluate] // Flatten // Sort]; row[n_] := Module[{sn}, sn = smoothNumbers[Prime[n], xMaxima[[n]]+1]; Reap[Do[If[sn[[i]]+1 == sn[[i+1]], Sow[sn[[i]]]], {i, 1, Length[sn]-1}]][[2, 1]]]; Table[Print[n]; row[n], {n, 1, 10}] // Flatten (* Jean-François Alcover, Jan 16 2015, updated Nov 10 2016 *)
  • PARI
    A138180_row=[]; A138180(n,k)={if(k, A138180(n)[k], #A138180_rowA138180_row=concat(A138180_row,vector(n)); if(#A138180_row[n], A138180_row[n], k=0; p=prime(n); A138180_row[n]=vector(A002071(n),i, until( vecmax(factor(k++)[, 1])<=p && vecmax(factor(k--+(k<2))[, 1])<=p,k++); k)))} \\ A138180(n) (w/o 2nd arg. k) returns the whole row. - M. F. Hasler, Jan 16 2015

A074399 a(n) is the largest prime divisor of n(n+1).

Original entry on oeis.org

2, 3, 3, 5, 5, 7, 7, 3, 5, 11, 11, 13, 13, 7, 5, 17, 17, 19, 19, 7, 11, 23, 23, 5, 13, 13, 7, 29, 29, 31, 31, 11, 17, 17, 7, 37, 37, 19, 13, 41, 41, 43, 43, 11, 23, 47, 47, 7, 7, 17, 17, 53, 53, 11, 11, 19, 29, 59, 59, 61, 61, 31, 7, 13, 13, 67, 67, 23, 23, 71, 71, 73, 73, 37, 19
Offset: 1

Views

Author

N. J. A. Sloane, Nov 29 2002

Keywords

Comments

Størmer shows that a(n) tends to infinity with n. Pólya generalized this result to other polynomials.
Kotov shows that a(n) >> log log n. - Charles R Greathouse IV, Mar 26 2012
Keates and Schinzel give effective constants for the above; in particular the latter shows that lim inf a(n)/log log n >= 2/7. - Charles R Greathouse IV, Nov 12 2012
Erdős conjectures ("on very flimsy probabilistic grounds") that for every e > 0, a(n) < (log n)^(2+e) infinitely often, while a(n) < (log n)^(2-e) only finitely often. - Charles R Greathouse IV, Mar 11 2015

References

  • S. V. Kotov, The greatest prime factor of a polynomial (in Russian), Mat. Zametki 13 (1973), pp. 515-522.
  • K. Mahler, Über den größten Primteiler spezieller Polynome zweiten Grades, Archiv for mathematik og naturvidenskab 41:6 (1934), pp. 3-26.
  • Georg Pólya, Zur arithmetischen Untersuchung der Polynome, Math. Zeitschrift 1 (1918), pp. 143-148.

Crossrefs

With A037464, the bisections of A076605.
Essentially the same as A069902.
Positions of primes <= p: A085152 (p=5), A085153 (p=7), A252494 (p=11), A252493 (p=13), A252492 (p=17).
Last position of each prime: A002072.

Programs

  • Mathematica
    Table[ Last[ Table[ # [[1]]] & /@ FactorInteger[n^2 - 1]], {n, 3, 160, 2}]
    Table[FactorInteger[n(n+1)][[-1,1]],{n,80}] (* Harvey P. Dale, Sep 28 2021 *)
  • PARI
    gpf(n)=my(f=factor(n)[,1]); f[#f]
    a(n)=if(n<3, n+1, max(gpf(n),gpf(n+1))) \\ Charles R Greathouse IV, Sep 14 2015

Formula

a(n) = Max (A006530(2n), A006530(2n+2)).
Pasten proves that a(n) >> (log log n)^2/(log log log n), see Corollary 1.5. - Charles R Greathouse IV, Oct 14 2024

Extensions

Extended by Robert G. Wilson v, Dec 02 2002

A252493 Numbers n such that n(n+1) is 13-smooth. (Related to the abc conjecture.)

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 20, 21, 24, 25, 26, 27, 32, 35, 39, 44, 48, 49, 54, 55, 63, 64, 65, 77, 80, 90, 98, 99, 104, 120, 125, 143, 168, 175, 195, 224, 242, 324, 350, 351, 363, 384, 440, 539, 624, 675, 728, 1000, 1715, 2079, 2400, 3024, 4095, 4224, 4374, 6655, 9800, 10647, 123200
Offset: 1

Views

Author

M. F. Hasler, Jan 16 2015

Keywords

Comments

Equivalently: Numbers n such that all prime factors of n and n+1 are <= 13, i.e., both are in A080197.
This sequence is complete by a theorem of Stormer, cf. A002071.
This is the 6th row of the table A138180. It has 68=A002071(6)=A145604(1)+...+ A145604(6) terms and ends with A002072(6)=123200. It is the union of all terms in rows 1 through 6 of the table A145605.
Contains A085152, A085153, A252494 as subsequences.

Crossrefs

Programs

  • Maple
    N:= 130000: # to get all entries <= N
    f:= proc(n)
    uses padic;
    evalb(2^ordp(n,2)*3^ordp(n,3)*5^ordp(n,5)*7^ordp(n,7)*11^ordp(n,11)*13^ordp(n,13) = n)
    end proc:
    L:= map(f, [$1..N+1]):
    select(t -> L[t] and L[t+1], [$1..N]); # Robert Israel, Jan 16 2015
  • Mathematica
    Select[Range[123456], FactorInteger[ # (# + 1)][[ -1,1]] <= 13 &]
  • PARI
    for(n=1,123456, vecmax(factor(n++,13)[,1])<17 && vecmax(factor(n--+(n<2),13))<17 && print1(n",")) \\ Skips the next n if n+1 is not 13-smooth: Twice as fast as the naïve version. Instead of vecmax(.)<17 one could use is_A080197().

A252494 Numbers n such that all prime factors of n and n+1 are <= 11. (Related to the abc conjecture.)

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 14, 15, 20, 21, 24, 27, 32, 35, 44, 48, 49, 54, 55, 63, 80, 98, 99, 120, 125, 175, 224, 242, 384, 440, 539, 2400, 3024, 4374, 9800
Offset: 1

Views

Author

M. F. Hasler, Jan 16 2015

Keywords

Comments

This sequence is complete by a theorem of Stormer, cf. A002071.
This is the 5th row of the table A138180. It has 40=A002071(5)=A145604(1)+...+ A145604(5) terms and ends with A002072(5)=9800. It is the union of all terms in rows 1 through 5 of the table A145605.
This is a subsequence of A252493, and contains A085152 and A085153 as subsequences.

Crossrefs

Programs

  • Mathematica
    Select[Range[10000], FactorInteger[ # (# + 1)][[ -1,1]] <= 11 &]
  • PARI
    for(n=1,9e6,vecmax(factor(n++)[,1])<12 && vecmax(factor(n--+(n<2))[,1])<12 && print1(n",")) \\ Skips 2 if n+1 is not 11-smooth: Twice as fast as the naive version.

A275156 The 108 numbers n such that n(n+1) is 17-smooth.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 20, 21, 24, 25, 26, 27, 32, 33, 34, 35, 39, 44, 48, 49, 50, 51, 54, 55, 63, 64, 65, 77, 80, 84, 90, 98, 99, 104, 119, 120, 125, 135, 143, 153, 168, 169, 175, 195, 220, 224, 242, 255, 272, 288, 324, 350, 351, 363, 374, 384, 440, 441, 539, 560, 594, 624, 675, 714, 728, 832, 935, 1000, 1088, 1155, 1224, 1274, 1700, 1715, 2057, 2079, 2400, 2430, 2499, 2600, 3024, 4095, 4224, 4374, 4913, 5831, 6655, 9800, 10647, 12375, 14399, 28560, 31212, 37179, 123200, 194480, 336140
Offset: 1

Views

Author

Jean-François Alcover, Nov 13 2016

Keywords

Comments

This is the 7th row of the table A138180.

References

Crossrefs

Programs

  • Mathematica
    pMax = 17; smoothMax = 10^12; smoothNumbers[p_?PrimeQ, max_] := Module[{a, aa, k, pp, iter}, k = PrimePi[p]; aa = Array[a, k]; pp = Prime[Range[k]]; iter = Table[{a[j], 0, PowerExpand@Log[pp[[j]], max/Times @@ (Take[pp, j - 1]^Take[aa, j - 1])]}, {j, 1, k}]; Table[Times @@ (pp^aa), Sequence @@ iter // Evaluate] // Flatten // Sort]; Select[(Sqrt[1 + 4*smoothNumbers[pMax, smoothMax]] - 1)/2, IntegerQ]
  • PARI
    is(n)=my(t=510510); n*=n+1; while((t=gcd(n,t))>1, n/=t); n==1 \\ Charles R Greathouse IV, Nov 13 2016

A275164 The 167 numbers n such that n(n+1) is 19-smooth.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 24, 25, 26, 27, 32, 33, 34, 35, 38, 39, 44, 48, 49, 50, 51, 54, 55, 56, 63, 64, 65, 75, 76, 77, 80, 84, 90, 95, 98, 99, 104, 119, 120, 125, 132, 135, 143, 152, 153, 168, 169, 170, 175, 189, 195, 208, 209, 220, 224, 242, 255, 272, 285, 288, 323, 324, 342, 350, 351, 360, 363, 374, 384, 399, 440
Offset: 1

Views

Author

Jean-François Alcover, Nov 14 2016

Keywords

Comments

See A002071.
The full list of 167 terms is given in the b-file (this is the 8th row of the table A138180).

Crossrefs

Programs

  • Mathematica
    pMax = 19; smoothMax = 10^15; smoothNumbers[p_?PrimeQ, max_] := Module[{a, aa, k, pp, iter}, k = PrimePi[p]; aa = Array[a, k]; pp = Prime[Range[k]]; iter = Table[{a[j], 0, PowerExpand@Log[pp[[j]], max/Times @@ (Take[pp, j - 1]^Take[aa, j - 1])]}, {j, 1, k}]; Table[Times @@ (pp^aa), Sequence @@ iter // Evaluate] // Flatten // Sort]; Select[(Sqrt[1 + 4*smoothNumbers[pMax, smoothMax]] - 1)/2, IntegerQ]
Showing 1-7 of 7 results.