A147643 Numbers A associated with the records of the merit function of the ABC conjecture admitting only C which are powers of 23.
7, 17, 162, 2
Offset: 1
Extensions
Edited by M. F. Hasler, Jan 16 2015
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.
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
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().
a = {{1, 1}}; aa = {1}; bb = {}; rr = {}; Do[logmax = 0; k = 2^x; w = Floor[(k - 1)/2]; Do[m = FactorInteger[n (k - n)]; rad = 1; Do[rad = rad m[[s]][[1]], {s, 1, Length[m]}]; log = Log[k]/Log[rad]; If[log > logmax, bmin = k - n; amax = n; logmax = log; r = rad], {n, 1, w, 2}]; Print[{x, amax}]; AppendTo[aa, amax]; AppendTo[bb, bmin]; AppendTo[rr, r]; AppendTo[a, {x, logmax}], {x, 2, 15}]; aa (* Artur Jasinski with assistance of M. F. Hasler *)
A143700(n) = {my(b=1, m=2^n-b); forstep(a=3, 2^(n-1), 2, A007947(a)*A007947(2^n-a)A007947((2^n-a)*b=a)); b; } \\ M. F. Hasler, Nov 13 2008
logmax = 0; aa = {}; bb = {}; cc = {}; dd = {}; ee = {}; ff = {}; gg \ = {}; Do[min = 10^100; max = 0; ile = 0; Do[If[GCD[m, n, n - m] == 1, ile = ile + 1; s = m n (n - m); k = FactorInteger[s]; g = 1; Do[g = g k[[p]][[1]], {p, 1, Length[k]}]; If[g > max, max = g; mmax = m]; If[g < min, min = g; mmin = m]], {m, 1, n - 1}]; AppendTo[aa, min]; AppendTo[bb, max]; AppendTo[cc, mmax]; AppendTo[dd, mmin]; AppendTo[gg, ile]; If[(Log[n]/Log[min]) > logmax, logmax = (Log[n]/Log[min]); AppendTo[ee, {N[logmax], n, mmin, min, mmax, max}]; Print[{N[logmax], n, mmin, min, mmax, max}]; AppendTo[ff, n]], {n, 2, 129}]; bb (* Artur Jasinski *) Table[Max[Times @@ FactorInteger[#][[All, 1]] & /@ ((m = Range[1, n - 1])*(n - m)*n)], {n, 2, 46}] (* Ivan Neretin, May 21 2015 *)
logmax = 0; aa = {}; bb = {}; cc = {}; dd = {}; ee = {}; ff = {}; gg \ = {}; Do[min = 10^100; max = 0; ile = 0; Do[If[GCD[m, n, n - m] == 1, ile = ile + 1; s = m n (n - m); k = FactorInteger[s]; g = 1; Do[g = g k[[p]][[1]], {p, 1, Length[k]}]; If[g > max, max = g; mmax = m]; If[g < min, min = g; mmin = m]], {m, 1, n - 1}]; AppendTo[aa, min]; AppendTo[bb, max]; AppendTo[cc, mmax]; AppendTo[dd, mmin]; AppendTo[gg, ile]; If[(Log[n]/Log[min]) > logmax, logmax = (Log[n]/Log[min]); AppendTo[ee, {N[logmax], n, mmin, min, mmax, max}]; Print[{N[logmax], n, mmin, min, mmax, max}]; AppendTo[ff, n]], {n, 2, 129}]; cc (* Artur Jasinski *)
(A,B,C) = (1,5,6) defines the first record, L=0.5268.. (A,B,C)=(1,11,12) reaches L=0.5931.. (A,B,C) = (1,17,18) reaches L=0.6249.. The first C-number selected from A033845 that does not generate a new record is 72.
Digits := 120 : A007947 := proc(n) local f,p; f := ifactors(n)[2] ; mul( op(1,p),p=f) ; end: L := proc(A,B,C) local rad; rad := A007947(A*B*C) ; evalf(log(C)/log(rad)) ; end: isA033845 := proc(n) local f,p; f := ifactors(n)[2] ; for p in f do if not op(1,p) in {2,3} then RETURN(false) ; fi; od: RETURN( (n mod 2 = 0 ) and (n mod 3 = 0 ) ) ; end: crek := -1 : for C from 3 do if isA033845(C) then for A from 1 to C/2 do B := C-A ; if gcd(A,B) = 1 then l := L(A,B,C) ; if l > crek then print(C) ; crek := l ; fi; fi; od: fi; od: # R. J. Mathar, Aug 24 2009
The case C=2 does not create a valid (A,B,C) triple, so C=4 is the first case, which sets a first record L=0.7737 with (A,B,C)=(1,3,4). The next admitted case, C=8, sets a new record L=0.7879 with (A,B,C)=(1,7,8), and so do (A,B,C)=(1,15,16) with L=0.8151. For C=32, we consider the largest L possible for A<B<C, which is (A,B,C)=(5,27,32) with L=1.0189. The value L=0.839 from (A,B,C)=(1,31,32) at the same C is smaller and discarded.
Digits := 120 : A007947 := proc(n) local f, p; f := ifactors(n)[2] ; mul( op(1, p), p=f) ; end: L := proc(A, B, C) local rad; rad := A007947(A*B*C) ; evalf(log(C)/log(rad)) ; end: crek := -1 : for x from 2 do C := 2^x ; for A from 1 to C/2 do B := C-A ; if gcd(A, B) = 1 then l := L(A, B, C) ; if l > crek then print(C) ; crek := l ; fi; fi; od: od: # R. J. Mathar, Aug 28 2009
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.
Comments