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

A300860 Indices of records in A300858.

Original entry on oeis.org

1, 8, 15, 16, 26, 27, 28, 32, 44, 52, 56, 62, 64, 76, 80, 88, 96, 100, 104, 112, 122, 124, 128, 144, 160, 176, 184, 192, 200, 216, 246, 248, 250, 256, 272, 276, 282, 288, 318, 320, 324, 348, 354, 366, 372, 384, 414, 426, 432, 468, 474, 486, 516, 522, 528, 534
Offset: 1

Views

Author

Michael De Vlieger, Mar 14 2018

Keywords

Comments

Consider numbers in the cototient of n, listed in row n of A121998. For composite n > 4, there are nondivisors m in the cototient, listed in row n of A133995. Of these m, there are two species. The first are m that divide n^e with integer e > 1, while the last do not divide n^e. These are listed in row n of A272618 and A272619, and counted by A243822(n) and A243823(n), respectively. This sequence lists the record setters in the sequence A300858(n), which is a function that represents the difference between the latter and the former species of nondivisors in the cototient of n.
Odd terms m < 36,000,000: {1, 15, 27}.
Smallest term m with A001221(m) = {0, 1, 2, ..., 8} = {1, 8, 15, 246, 2010, 9870, 30030, 510510, 9699690} (the last 3 terms are in A002110).
Smallest term m with A001222(m) = {0, 2, 3, ..., 12} = {1, 15, 8, 16, 32, 64, 128, 256, 768, 1536, 7680, 53760, 3843840} (includes 2^e with 3 <= e <= 8). Note, A300858(p) for p prime = 0.

Examples

			8 is in the sequence because A300858(n) for n < 8 is negative or 0 after A300858(1) = 0. A300858(8) = A243823(8) - A243822(8) = 1 - 0 = 1. Within the cototient of 8 there is one nondivisor (6) and it does not divide 8^e for integer e. (All prime powers m have A243822(m) = 0 and for m > 4, A243823(m) is positive.)
15 is in the sequence because -1 <= A300858(n) <= 1 for n < 15. A300858(15) = 2. Within the cototient of 15 there are 4 nondivisors; of these 3 (i.e., {6, 10, 12}) do not divide 15^e for integer e, but 9 | 15^2. Therefore 3 - 1 = 2 and 2 exceeds all values A300858(n) for n < 15.
		

Crossrefs

Programs

  • Mathematica
    f[n_] := Count[Range@ n, _?(PowerMod[n, Floor@ Log2@ n, #] == 0 &)]; With[{s = Array[#1 - #3 + 1 - 2 #2 + #4 & @@ {#, f@ #, EulerPhi@ #, DivisorSigma[0, #]} &, 550]},  Map[FirstPosition[s, #][[1]] &, Union@ FoldList[Max, s]] ]
  • PARI
    a300858(n) = 1 + n + numdiv(n) - eulerphi(n) - 2*sum(k=1, n, if(gcd(n, k)-1, 0, moebius(k)*(n\k))) \\ after Michel Marcus
    r=-1; for(i=1, oo, if(a300858(i) > r, print1(i, ", "); r=a300858(i))) \\ Felix Fröhlich, Mar 30 2018

A300861 Records in A300858.

Original entry on oeis.org

0, 1, 2, 4, 5, 6, 7, 11, 13, 17, 19, 21, 26, 27, 31, 35, 37, 40, 43, 47, 49, 51, 57, 66, 73, 79, 81, 93, 95, 109, 111, 113, 119, 120, 127, 129, 133, 153, 155, 163, 172, 173, 177, 185, 189, 211, 213, 223, 245, 247, 253, 271, 277, 279, 283, 289, 301, 303, 309, 336
Offset: 1

Views

Author

Michael De Vlieger, Mar 28 2018

Keywords

Comments

Consider numbers in the cototient of n, listed in row n of A121998. For composite n > 4, there are nondivisors m in the cototient, listed in row n of A133995. Of these m, there are two species. The first are m that divide n^e with integer e > 1, while the last do not divide n^e. These are listed in row n of A272618 and A272619, and counted by A243822(n) and A243823(n), respectively. This sequence lists the records in A300858, which is a function that represents the difference between the latter and the former species of nondivisors in the cototient of n.

Examples

			0 is the first term since A300858(1) = 0. A300858 is 0 or negative for n < 8.
A300858(8) = A243823(8) - A243822(8) = 1 - 0 = 1. Within the cototient of 8 there is one nondivisor (6) and it does not divide 8^e for integer e. (All prime powers m have A243822(m) = 0 and for m > 4, A243823(m) is positive.) Therefore 1 is the next term. Between 8 and 15, -1 <= A300858(n) <= 1.
A300858(15) = 2. Within the cototient of 15 there are 4 nondivisors; of these 3 (i.e., {6, 10, 12}) do not divide 15^e for integer e, but 9 | 15^2. Therefore 3 - 1 = 2 and 2 exceeds all values A300858(n) for n < 15, and appears after 1.
		

Crossrefs

Programs

  • Mathematica
    f[n_] := Count[Range@ n, _?(PowerMod[n, Floor@ Log2@ n, #] == 0 &)]; Union@ FoldList[Max, Array[#1 - #3 + 1 - 2 #2 + #4 & @@ {#, f@ #, EulerPhi@ #, DivisorSigma[0, #]} &, 600]]
  • PARI
    a300858(n) = 1 + n + numdiv(n) - eulerphi(n) - 2*sum(k=1, n, if(gcd(n, k)-1, 0, moebius(k)*(n\k))) \\ after Michel Marcus in A300858
    r=-1; for(x=1, oo, if(a300858(x) > r, r=a300858(x); print1(r, ", "))) \\ Felix Fröhlich, Mar 30 2018

A300859 Where records occur in A045763.

Original entry on oeis.org

1, 6, 10, 14, 18, 22, 26, 30, 36, 38, 42, 50, 54, 60, 66, 78, 84, 90, 102, 114, 120, 126, 132, 138, 150, 168, 174, 180, 186, 198, 204, 210, 234, 240, 246, 252, 258, 264, 270, 294, 300, 318, 330, 360, 378, 390, 420, 450, 462, 480, 504, 510, 540, 546, 570, 600
Offset: 1

Views

Author

Michael De Vlieger, Mar 15 2018

Keywords

Comments

The cototient of n consists of numbers 1 < m <= n that are not coprime to n, i.e., gcd(m,n) > 1. These numbers have at least one prime divisor p that also divides n. The cototient of n contains the divisors d of n; the remaining nondivisors in the cototient of n are listed in A133995. The counting function of A133995 is A045763(n). There are two species of numbers in the nondivisor-cototient of n: those in row n of A272618, of which A243822(n) is counting function, and those in row n of A272619, of which A243823(n) is the counting function. The former species divides n^e for integer e > 1, while the latter does not divide any integer power of n.
A045763(p) = 0 for p prime, therefore there are no primes in a(n).
Except for prime terms (i.e., 2), A002110 is a subset as primorials minimize the totient function. The divisor counting function is increasingly vanishingly small compared to the totient function for A002110(i) as i increases, and A002110(i) for 1 < i <= 9 is observed in a(n).
Conjectures based on 1255 terms of a(n) < 36,000,000:
1. There are no prime powers p^e > 1 in a(n), i.e., the intersection of a(n) and A000961 is {1}.
2. A293555 is a subset of A300859. Numbers that have a lot of nondivisors m | n^e with e > 1 (i.e., in row n of A272618 and counted by A243822(n)) tend to reduce the totient and increasingly have fewer divisors than highly composite numbers, widening the nondivisor-cototient.
3. A300156 is a subset of A300859. Numbers that have more nondivisors m | n^e with e > 1 (i.e., in row n of A272618 and counted by A243822(n)) than divisors tend to reduce the totient and have fewer divisors than highly composite numbers (i.e., those n in A002182), widening the nondivisor-cototient.
Increasingly many terms k in A262867 also appear in a(n) as k increases. A292867 lists record-setters in A243823, which is the counting function of one of the two species of nondivisors in the cototient of n.

Examples

			6 is in the sequence because there is 1 nondivisor in the cototient of 6 (i.e., 4), and that total exceeds 0 for all smaller positive numbers.
10 follows 6 because there are 3 nondivisors in the cototient (4, 6, 8), and this exceeds the total of 1 for n = 6, 8, and 9.
		

References

  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, p. 352 (sixth edition), see Theorem 327.

Crossrefs

Programs

  • Mathematica
    With[{s = Array[1 + # - EulerPhi@ # - DivisorSigma[0, #] &, 10^3]}, Map[FirstPosition[s, #][[1]] &, Union@ FoldList[Max, s]]]
    DeleteDuplicates[Table[{n,n+1-DivisorSigma[0,n]-EulerPhi[n]},{n,600}],GreaterEqual[#1 [[2]],#2 [[2]]]&][[;;,1]] (* Harvey P. Dale, Mar 29 2023 *)

Formula

a(n) = 1 + n - A000010(n) - A000005(n).
Showing 1-3 of 3 results.