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
-
Select[Range[10000], FactorInteger[ # (# + 1)][[ -1,1]] <= 5 &] (* T. D. Noe, Mar 03 2008 *)
-
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
A085153
All prime factors of n and n+1 are <= 7. (Related to the abc conjecture.)
Original entry on oeis.org
1, 2, 3, 4, 5, 6, 7, 8, 9, 14, 15, 20, 24, 27, 35, 48, 49, 63, 80, 125, 224, 2400, 4374
Offset: 1
-
Select[Range[10000], FactorInteger[ # (# + 1)][[ -1,1]] <= 7 &] (* T. D. Noe, Mar 03 2008 *)
-
for(n=1,9e6,vecmax(factor(n++)[,1])<8 && vecmax(factor(n--+(n<2))[,1])<8 && print1(n",")) \\ M. F. Hasler, Jan 16 2015
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
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),
...
-
(* 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 *)
-
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
- 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.
- Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
- P. Erdős, Problems and results on number theoretic properties of consecutive integers and related questions, Proceedings of the Fifth Manitoba Conference on Numerical Mathematics (Univ. Manitoba, Winnipeg, Man., 1975), Congress. Numer. XVI , pp. 25-44, Utilitas Math., Winnipeg, Man., 1976.
- M. Keates, On the greatest prime factor of a polynomial (1968), pp. 301-303.
- Hector Pasten, The largest prime factor of n^2+1 and improvements on subexponential ABC, arXiv:2312.03566 [math.NT] (2024)
- A. Schinzel, On two theorems of Gelfond and some of their applications, Acta Arithmetica 13:2 (1967-1968), pp. 177-236.
- Carl Størmer, Quelques théorèmes sur l'équation de Pell x^2 - Dy^2 = +-1 et leurs applications (in French), Skrifter udgivne af Videnskabsselskabet i Christiania: Mathematisk-naturvidenskabelig Klasse (1897).
Last position of each prime:
A002072.
-
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 *)
-
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
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
-
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
-
Select[Range[123456], FactorInteger[ # (# + 1)][[ -1,1]] <= 13 &]
-
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().
A252492
The largest prime factor of n*(n+1) equals 17. (Related to the abc conjecture.)
Original entry on oeis.org
16, 17, 33, 34, 50, 51, 84, 119, 135, 153, 169, 220, 255, 272, 288, 374, 441, 560, 594, 714, 832, 935, 1088, 1155, 1224, 1274, 1700, 2057, 2430, 2499, 2600, 4913, 5831, 12375, 14399, 28560, 31212, 37179, 194480, 336140
Offset: 1
-
Select[Range[345678], FactorInteger[ # (# + 1)][[ -1,1]] == 17 &]
-
for(n=1,9e6,vecmax(factor(n++)[,1])<18 && vecmax(factor(n*n--)[,1])==17 && print1(n",")) \\ Skips 2 if n+1 is not 17-smooth: Twice as fast as the naïve version.
Showing 1-6 of 6 results.
Comments