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
A002071
Number of pairs of consecutive integers x, x+1 such that all prime factors of both x and x+1 are at most the n-th prime.
Original entry on oeis.org
1, 4, 10, 23, 40, 68, 108, 167, 241, 345, 482, 653, 869, 1153, 1502, 1930, 2454, 3106, 3896
Offset: 1
- N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
- N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
- Lucas A. Brown, stormer.py.
- E. F. Ecklund and R. B. Eggleton, Prime factors of consecutive integers, Amer. Math. Monthly, 79 (1972), 1082-1089.
- D. Eppstein, Smooth pairs.
- D. Eppstein, Python program
- D. H. Lehmer, On a problem of Størmer, Ill. J. Math., 8 (1964), 57-69.
- C. Stormer, Quelques théorèmes sur l'équation de Pell x^2 - Dy^2 = +-1 et leurs applications, Skrifter Videnskabs-selskabet (Christiania), Mat.-Naturv (1897). Kl. I (2).
- Wikipedia, Størmer's theorem
- OEIS Index entries for sequences related to the abc conjecture
Cf.
A138180 (triangle of x values for each n).
-
(* 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]; a[n_] := Module[{sn, cnt}, sn = smoothNumbers[Prime[n], xMaxima[[n]]+1]; cnt = 0; Do[If[sn[[i]]+1 == sn[[i+1]], cnt++], {i, 1, Length[sn]-1}]; cnt]; Table[an = a[n]; Print["a(", n, ") = ", an]; an, {n, 1, 16}] (* Jean-François Alcover, Nov 10 2016 *)
A002072 = {1, 8, 80, 4374, 9800, 123200, 336140, 11859210, 11859210};
Table[Length[Select[Table[Max[FactorInteger[x], FactorInteger[x + 1]], {x, A002072[[n]]}], # <= Prime[n] &]], {n, 7}] (* Robert Price, Oct 29 2018 *)
-
A002071(n)=[1,4,10,23,40,68,108,167,241,345,482,653,869,1153,1502][n] \\ "practical" solution. - M. F. Hasler, Jan 16 2015
-
A002071(n,b=A002072,c=1,p=prime(n))={for(k=2,b(n),vecmax(factor(k++,p)[,1])<=p && vecmax(factor(k--+(k<2),p)[,1])<=p && c++); c} \\ b can be any upper bound for A002072, e.g., n->10^n should work, too. - 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
A145604
Number of pairs of consecutive integers x, x+1 such that both are prime(n)-smooth but both are not prime(n-1)-smooth.
Original entry on oeis.org
1, 3, 6, 13, 17, 28, 40, 59, 74, 104, 137, 171, 216, 284, 349, 428, 524, 652, 790
Offset: 1
A145606
Largest number x such that x and x+1 are prime(n)-smooth but not prime(n-1)-smooth.
Original entry on oeis.org
1, 8, 80, 4374, 9800, 123200, 336140, 11859210, 5142500, 177182720, 1611308699, 3463199999, 63927525375, 421138799639, 1109496723125, 1453579866024, 20628591204480, 31887350832896, 12820120234375, 119089041053696, 2286831727304144, 9591468737351909375, 17451620110781856, 166055401586083680, 49956990469100000, 4108258965739505499, 19316158377073923834000, 386539843111191224
Offset: 1
Terms a(16) onward by
Andrey V. Kulsha, Aug 10 2011, according to Jim White's computations
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.
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
-
Select[Range[10000], FactorInteger[ # (# + 1)][[ -1,1]] <= 11 &]
-
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.
A193945
Integers n such that for all i > n the largest prime factor of i(i+1)(i+2)(i+3) exceeds the largest prime factor of n(n+1)(n+2)(n+3).
Original entry on oeis.org
1, 3, 7, 9, 63, 168, 322, 1518, 1680, 10878, 17575, 70224, 97524, 224846, 612360, 15473807, 61011223, 108384553, 137308589, 999208250, 1497820158, 22377473780, 238092971076, 316419292600, 23931257472314
Offset: 1
Showing 1-10 of 13 results.
Comments