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.

A147802 Least m coprime to 3 minimizing A007947(m*(3^n-m)).

Original entry on oeis.org

1, 1, 2, 1, 1, 25, 139, 289, 31, 1, 16096, 49, 7424, 588665, 619115, 83521, 8000000, 1515625, 505620842, 17643776, 244140625, 5443635008
Offset: 1

Views

Author

M. F. Hasler, Nov 13 2008

Keywords

Comments

Related to the abc conjecture: Since m is coprime to 3, it is also coprime to 3^n and thus to 3^n-m. Thus A007947(m*(3^n-m)*3^n) = 3*A007947(m(3^n-m)).

Crossrefs

Cf. A007947, A143700 (analog for 2^n), A147300 (general case).

Programs

  • Mathematica
    rad[n_] := Times @@ FactorInteger[n][[All, 1]];
    a[n_] := MinimalBy[Select[Range[3^n - 2], CoprimeQ[#, 3] &], rad[# (3^n - #)] &][[1]];
    Reap[Do[Print[n, " ", an = a[n]]; Sow[an], {n, 1, 16}]][[2, 1]] (* Jean-François Alcover, Mar 27 2020 *)
  • PARI
    A147802(n, p=3) = {local(b, m=n=p^n); for(a=1, (n-1)\2, a%p || next; A007947(n-a)*A007947(a)A007947((n-a)*b=a)); b; }

Extensions

a(17) from Jean-François Alcover, Mar 28 2020
a(18)-a(22) from Giovanni Resta, Mar 29 2020

A143702 a(n) is the minimal values of A007947((2^n)*m*(2^n-m)).

Original entry on oeis.org

2, 6, 14, 30, 30, 42, 30, 78, 182, 1110, 570, 1830, 6666, 2310, 2534, 5538, 9870, 20010, 141270, 14070, 480090, 155490, 334110, 1794858, 2463270, 2132130, 2349390
Offset: 1

Views

Author

Artur Jasinski, Nov 10 2008

Keywords

Comments

The product of distinct prime divisors of (2^n)*m*(2^n-m) is also called the radical of that number: rad((2^n)*m*(2^n-m)).
For numbers m see A143700.

Crossrefs

Programs

  • Mathematica
    aa = {1}; bb = {1}; rr = {2}; 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, 2*r]; AppendTo[a, {x, logmax}], {x, 2, 15}]; rr (* Artur Jasinski with assistance of M. F. Hasler *)

Extensions

Name changed and a(1) added by Jinyuan Wang, Aug 11 2020

A147803 Least m coprime to 5 minimizing A007947(m*(5^n-m)).

Original entry on oeis.org

1, 1, 4, 49, 128, 9, 36864, 19332, 4508, 121, 2
Offset: 1

Views

Author

M. F. Hasler, Nov 13 2008

Keywords

Comments

The minima are given in A147800.
This is related to the abc conjecture: Since m is coprime to 5, it is also coprime to 5^n and thus to 5^n-m. Thus the squarefree kernel A007947(m*(5^n-m)*5^n) = 5*A007947(m*(5^n-m)).

Crossrefs

Cf. A007947, A147298 (general case), A147800 (value of minima), A143700 (analog for 2^n), A147802 (analog for 3^n), A147300 (analog for any number).

Programs

  • PARI
    A147803(n,p=5) = {my(b, m=n=p^n); for(a=1, n\2, a%p || next; A007947(n-a)*A007947(a)A007947((n-a)*b=a)); b; }

A179017 Odd numbers c such that c*(c^2 - 1)/4 is squarefree.

Original entry on oeis.org

3, 5, 11, 13, 21, 29, 43, 59, 61, 67, 69, 77, 83, 85, 93, 115, 123, 131, 133, 139, 141, 155, 157, 165, 173, 187, 203, 205, 211, 213, 219, 221, 227, 229, 237, 259, 267, 277, 283, 285, 291, 309, 317, 331, 347, 355, 357, 365, 371, 373, 381, 389, 403, 411, 419, 421
Offset: 1

Views

Author

Artur Jasinski, Jun 24 2010

Keywords

Comments

Original title was: "Numbers c such that (c^2-1)c is square free and gcd(c-1,c,c+1)=1", but (c^2-1)c is never squarefree for odd c, and gcd(n,n+1) is always = 1. - M. F. Hasler, Nov 03 2013
These numbers c with distribution a+b=c such that a=(c-1)/2 (see A172186) and b=(c+1)/2 (see A179019) have minimal possible values with function L(a,b,c) = log(c)/log(N(a,b,c)) = log(c)/log((c^2-1)c/4).
This function is minimal orbital in hypothesis (a,b,c).
There are no numbers or distributions which have value L less than log(c)/log((c^2-1)*c/4).
Equivalently, odd squarefree numbers c such that (c^2 - 1)/4 is also squarefree. - Jon E. Schoenfield, Feb 13 2023
The asymptotic density of this sequence is Product_{p prime} (1 - 3/p^2) = A206256 = 0.125486980905... (Tsang, 1985). - Amiram Eldar, Feb 26 2024

Crossrefs

Programs

  • Mathematica
    aa = {}; Do[If[(GCD[x, (x - 1)/2] == 1) && (GCD[x, (x + 1)/2] == 1) && (GCD[(x - 1)/2, (x + 1)/2] == 1), If[SquareFreeQ[(x^2 - 1) x/4], AppendTo[aa, x]]], {x, 2, 1000}]; aa
  • PARI
    forstep(n=3,421,2,issquarefree(n*(n^2-1)/4)&&print1(n",")) \\ M. F. Hasler, Nov 03 2013
    
  • PARI
    is(n)=n%2 && issquarefree(n) && issquarefree(n^2\4) \\ Charles R Greathouse IV, Mar 11 2014

Formula

a(n) = 2*A172186(n) + 1. - Bernard Schott, Mar 06 2023

Extensions

Edited by M. F. Hasler, Nov 03 2013

A147804 Least m coprime to 7 minimizing A007947(m*(7^n-m)).

Original entry on oeis.org

1, 1, 100, 1, 423, 28561, 3072, 124609, 119232
Offset: 1

Views

Author

M. F. Hasler, Nov 13 2008

Keywords

Comments

The minima are given in A147799.
This is related to the abc conjecture: Since m is coprime to 7, it is also coprime to 7^n and thus to 7^n-m. Thus the squarefree kernel A007947(m*(7^n-m)*7^n) = 7*A007947(m(7^n-m)).

Crossrefs

Cf. A007947, A147799 (value of minima), A143700, A147802, A147803 (analog for 2^n, 3^n, 5^n), A147300 (analog for any number).

Programs

  • PARI
    A147804(n,p=7)={my(b, m=3*n=p^n, t); for(a=1, n\2, a%p || next; m>2*(t=A007947(a)) || next; m>(t*=A007947(n-a)) || next; m=t; b=a); b; }

Extensions

Typo in title corrected by M. F. Hasler, Nov 17 2008

A179019 a(n) = (A179017(n)+1)/2.

Original entry on oeis.org

2, 3, 6, 7, 11, 15, 22, 30, 31, 34, 35, 39, 42, 43, 47, 58, 62, 66, 67, 70, 71, 78, 79, 83, 87, 94, 102, 103, 106, 107, 110, 111, 114, 115, 119, 130, 134, 139, 142, 143, 146, 155, 159, 166, 174, 178, 179, 183, 186, 187, 191, 195, 202, 206, 210, 211, 214, 215, 218, 219
Offset: 1

Views

Author

Artur Jasinski, Jun 24 2010

Keywords

Comments

For numbers a and c, see A172186 and A179017. Numbers b are this sequence.
These numbers c, with distribution a+b=c such that a=(c-1)/2 and b=(c+1)/2, have minimal possible values with function L(a,b,c) = log(c)/log(N[a,b,c]) = log(c)/log((c^2-1)c/4).
There exist no numbers or distributions for which L < log(c)/log((c^2-1)c/4). - Artur Jasinski

Crossrefs

Programs

  • Mathematica
    aa = {}; Do[If[(GCD[x, (x - 1)/2] == 1) && (GCD[x, (x + 1)/2] == 1) && (GCD[(x - 1)/2, (x + 1)/2] == 1), If[SquareFreeQ[(x^2 - 1) x/4], AppendTo[aa, (x + 1)/2]]], {x, 2, 1000}]; aa

Formula

a(n) = A179017(n) - A172186(n). - Hugo Pfoertner, Mar 22 2020

A143701 a(n) is the least odd number 2^n - m minimizing A007947(m*(2^n - m)).

Original entry on oeis.org

1, 3, 7, 15, 27, 63, 125, 243, 343, 999, 1805, 3721, 8181, 16335, 32761, 65533, 112847, 190269, 519375, 1046875, 1953125, 4192479, 8385125, 16775019, 24398405, 66976875, 134216625
Offset: 1

Views

Author

Artur Jasinski, Nov 10 2008

Keywords

Comments

a(n) is the smallest odd number such that the product of distinct prime divisors of (2^n)*a(n)*(2^n-a(n)) is the smallest for the range a(n) <= 2^x - a(n) < 2^x.
The product of distinct prime divisors of m*(2^n-m) is also called the radical of that number: rad(m*(2^n-m)).

Crossrefs

Programs

  • Mathematica
    aa = {1}; bb = {1}; 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}]; bb (* Artur Jasinski with assistance of M. F. Hasler *)

Formula

a(n) = 2^n - A143700(n).

Extensions

a(1) added by Jinyuan Wang, Aug 11 2020
Showing 1-7 of 7 results.