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-10 of 22 results. Next

A300858 a(n) = A243823(n) - A243822(n).

Original entry on oeis.org

0, 0, 0, 0, 0, -1, 0, 1, 1, -1, 0, -1, 0, 1, 2, 4, 0, -1, 0, 3, 4, 3, 0, 3, 3, 5, 6, 7, 0, -5, 0, 11, 6, 7, 6, 6, 0, 9, 8, 11, 0, 1, 0, 13, 12, 13, 0, 13, 5, 13, 12, 17, 0, 13, 10, 19, 14, 19, 0, 5, 0, 21, 18, 26, 12, 11, 0, 23, 18, 15, 0, 25, 0, 25, 24, 27
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 is the difference between the latter and the former species of nondivisors in the cototient of n.
Since A045763(n) = A243822(n) + A243823(n), this sequence examines the balance of the two components among nondivisors in the cototient of n.
For positive n < 6 and for p prime, a(n) = a(p) = 0, thus a(A046022(n)) = 0.
For prime powers p^e, a(p^e) = A243823(p^e), since A243822(p^e) = 0, thus a(n) = A243823(n) for n in A000961.
Value of a(n) is generally nonnegative; a(n) is negative for n = {6, 10, 12, 18, 30}; a(30) = -5, but a(n) = -1 for the rest of the aforementioned numbers. These five numbers are a subset of A295523.

Examples

			a(6) = -1 since the only nondivisor in the cototient of 6 is 4, which divides 6^e with e > 1; therefore 0 - 1 = -1.
a(8) = 1 since the only nondivisor in the cototient of 8 is 6, and 6 does not divide 8^e with e > 1, therefore 1 - 0 = 1.
Some values of a(n) and related sequences:
   n  a(n) A243823(n) A243822(n)    A272619(n)       A272618(n)
  -------------------------------------------------------------
   1   0          0          0      -                -
   2   0          0          0      -                -
   3   0          0          0      -                -
   4   0          0          0      -                -
   5   0          0          0      -                -
   6  -1          0          1      -                {4}
   7   0          0          0      -                -
   8   1          1          0      {6}              -
   9   1          1          0      {6}              -
  10  -1          1          2      {6}              {4,8}
  11   0          0          0      -                -
  12  -1          1          2      {10}             {8,9}
  13   0          0          0      -                -
  14   1          3          2      {6,10,12}        {4,8}
  15   2          3          1      {6,10,12}        {9}
  16   4          4          0      {6,10,12,14}     -
  17   0          0          0      -                -
  18  -1          3          4      {10,14,15}       {4,8,12,16}
  19   0          0          0      -                -
  20   3          5          2      {6,12,14,15,18}  {8,16}
  ...
		

Crossrefs

Programs

  • Mathematica
    f[n_] := Count[Range@ n, _?(PowerMod[n, Floor@ Log2@ n, #] == 0 &)]; Array[#1 - #3 + 1 - 2 #2 + #4 & @@ {#, f@ #, EulerPhi@ #, DivisorSigma[0, #]} &, 76]
  • PARI
    a(n) = 1 + n + numdiv(n) - eulerphi(n) - 2*sum(k=1, n, if(gcd(n,k)-1, 0, moebius(k)*(n\k))); \\ Michel Marcus, Mar 17 2018

Formula

a(n) = 1 + n - A000010(n) - 2*A010846(n) + A000005(n).

A292867 Indices of records in A243823.

Original entry on oeis.org

1, 8, 14, 16, 20, 22, 26, 28, 32, 38, 40, 44, 46, 48, 50, 52, 54, 56, 58, 62, 64, 68, 72, 78, 80, 86, 88, 92, 94, 96, 100, 108, 114, 122, 124, 126, 130, 132, 138, 144, 156, 160, 162, 174, 186, 192, 204, 216, 222, 228, 234, 240, 246, 252, 258, 264, 270, 276, 282
Offset: 1

Views

Author

Michael De Vlieger, Oct 02 2017

Keywords

Comments

Except for A292867(1) = 1, all terms are even.
Some conjectures:
1. The only prime powers p^e in this sequence are {8, 16, 32, 64}.
2. Squarefree terms m appear throughout. (There are 261 squarefree values among the first 1261 terms.)
3. Terms that set records for omega(m) are 1, followed by 2^e, with 3 <= e <= 6, then 2^e * 3 with 6 <= e <= 8, then 2^7 * A002110(k) with k >= 1.
4. Primorials A002110(n) for n >= 6 appear in this sequence. The first primorials in m are terms 6 through 8 of A002110 (i.e., 30030, 510510, 9699690) at n = 419, 774, 1258, respectively.
5. Outside of a(n) with 2 <= n <= 21 and n = {29, 30}, all terms of A244052 are also in this sequence. This observation applies to the smallest 104 terms of A244052.
6. For very large n, all terms are also in A244052. For small n, few terms of A244052 appear and are separated by many other numbers. Since numbers m in A244052 are products of k primes, many of which are the smallest primes, phi is minimized and A010846(m) becomes infinitesimal in comparison to m. Therefore A243823(m) is tantamount to the cototient of m. The size of n required to observe this agreement between this sequence and A244052 is unknown.

Examples

			8 is in the sequence since it is the first number n such that A243823(n) > 0. 14 appears immediately after 8 since A243823(14) = 3, and 3 is greater than the values that precede it.
Table of indices a(n) of records b(n) in A243823 = n - phi(n) - A010846(n) + 1:
.
   n  a(n) b(n) phi(a(n)) A010846(a(n))
  -------------------------------------
   1    1    0         1         1
   2    8    1         4         4
   3   14    3         6         6
   4   16    4         8         5
   5   20    5         8         8
   6   22    6        10         7
   7   26    8        12         7
   8   28    9        12         8
   9   32   11        16         6
  10   38   13        18         8
  11   40   14        16        11
  12   44   16        20         9
  13   46   17        22         8
  14   48   18        16        15
  15   50   19        20        12
  16   52   20        24         9
  17   54   21        18        16
  18   56   22        24        11
  19   58   23        28         8
  20   62   25        30         8
  ...
		

Crossrefs

Programs

A292868 Records in A243823.

Original entry on oeis.org

0, 1, 3, 4, 5, 6, 8, 9, 11, 13, 14, 16, 17, 18, 19, 20, 21, 22, 23, 25, 26, 27, 31, 32, 35, 36, 38, 39, 40, 45, 46, 52, 53, 54, 55, 58, 60, 62, 68, 74, 78, 79, 85, 90, 98, 104, 108, 119, 120, 123, 126, 128, 135, 136, 142, 144, 146, 153, 158, 164, 166, 171, 176
Offset: 1

Views

Author

Michael De Vlieger, Oct 02 2017

Keywords

Comments

See A292867 for comments, linked tables, and conjectures. - Michael De Vlieger, Nov 17 2017

Examples

			First 20 terms and their positions b(n) in A243823:
   n  a(n) b(n)
  -------------
   1    0    1
   2    1    8
   3    3   14
   4    4   16
   5    5   20
   6    6   22
   7    8   26
   8    9   28
   9   11   32
  10   13   38
  11   14   40
  12   16   44
  13   17   46
  14   18   48
  15   19   50
  16   20   52
  17   21   54
  18   22   56
  19   23   58
  20   25   62
		

Crossrefs

Programs

A294575 Numbers n such that 2 * A243823(n) > n.

Original entry on oeis.org

144, 162, 174, 186, 192, 198, 200, 204, 216, 220, 222, 228, 230, 234, 238, 240, 246, 250, 252, 258, 260, 264, 266, 270, 276, 280, 282, 288, 290, 294, 300, 306, 308, 310, 312, 318, 320, 322, 324, 330, 336, 340, 342, 348, 350, 354, 360, 364, 366, 370, 372, 374
Offset: 1

Views

Author

Michael De Vlieger, Nov 17 2017

Keywords

Comments

Consider A243823(n), which is the number of m < n that are products of at least one prime p | n and at least one prime q that does not divide n. These numbers m in the cototient of n do not divide a power of n. This sequence lists numbers n where such numbers m are predominant.

Examples

			A272619(10) = (6), A243823(10) = 1, so 10 is not in the sequence. A243823(144) = 74, which is greater than half of all numbers less than or equal to 144, thus 144 is the first term.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[10^3], Function[n, 2 (n - (Count[Range@ n, _?(PowerMod[n, Floor@ Log2@ n, #] == 0 &)] + EulerPhi[n] - 1)) > n]]

A294576 Odd numbers n such that 2 * A243823(n) > n.

Original entry on oeis.org

945, 1155, 1365, 1575, 1785, 1995, 2145, 2205, 2415, 2625, 2805, 2835, 3003, 3045, 3135, 3255, 3315, 3465, 3675, 3705, 3795, 3885, 4095, 4125, 4305, 4455, 4485, 4515, 4725, 4785, 4845, 4935, 5115, 5145, 5355, 5445, 5565, 5655, 5775, 5865, 5985
Offset: 1

Views

Author

Michael De Vlieger, Nov 17 2017

Keywords

Comments

Odd terms in A294575.
2^k - 1 is in the sequence for k = 12, 20, 24, 28, 30, 36, 40, 48, ...

Crossrefs

Programs

  • Mathematica
    Select[Range[1, 6000, 2], Function[n, 2 (n - (Count[Range@ n, _?(PowerMod[n, Floor@ Log2@ n, #] == 0 &)] + EulerPhi[n] - 1)) > n]]

A295221 Numbers k such that 2*A243823(k) = k.

Original entry on oeis.org

156, 190, 224, 286, 352, 416, 544, 578, 608, 736, 928, 992, 1184, 1312, 1376, 1504, 1696, 1888, 1952, 33555776, 33557824, 33558208, 33558464, 33558592, 33559616, 33560768, 33560896, 33562304, 33562432, 33563456, 33564992, 33567808, 33568448, 33568576, 33569216
Offset: 1

Views

Author

Michael De Vlieger, Nov 17 2017

Keywords

Comments

Observations:
1. There is a large gap between a(19) and a(20).
2. Products 2^5 * prime(i), with 3 <= i <= 17, are in the sequence.
3. Products 2^6 * prime(j), with 43391 <= j <= 82025, are in the sequence.
4. a(1) = 2^2 * 3 * 13, and terms 190, 286, and 578 are even, but do not follow the pattern of 2^h*p prime.

Examples

			a(1) = 156 since 2 * (A010846(156) + A000010(156) - 1) = 2 * (31 + 48 - 1) = 2 * 78 = 156.
		

Crossrefs

Programs

A295523 Nonprime numbers n such that A243822(n) >= A243823(n).

Original entry on oeis.org

1, 4, 6, 10, 12, 18, 30
Offset: 1

Views

Author

Michael De Vlieger, Nov 23 2017

Keywords

Comments

Consider numbers m that are nondivisors in the cototient of n, listed in row n of A133995 and counted by A045763(n). This sequence lists numbers n for which there are more m such that m | n^e with e >= 0 than there are m that are products of at least one prime divisor p of n and one nondivisor prime q. The former species of m are "semidivisors" listed in row n of A272618 and counted by A243822(n), while the latter are "semitotatives" listed in row n of A272619 and counted by A243823(n). These two species constitute the only species of nondivisors in the cototient of n.
Primes p have no nondivisors in the cototient, i.e., A045763(p) = 0, therefore A243822(p) and A243823(p) also are 0. The equality of these latter two sequences is trivial in the case of primes.
Prime powers p^e except for p^e = 4 have A243823(p^e) > A243822(p^e), since A243822(p^e) = 0. All powers p^k with 0 <= k <= e divide p^e.
The sequence is finite because there exist a lot more nondivisor primes q than p | n as n increases. Therefore there are more numbers m in row n of A272619 than there are in row n of A272618, since the former are products p*q and the latter are products only of p.

Examples

			1 is in the sequence because it is not prime and there are no nondivisors in the cototient, therefore A243822(1) = A243823(1) = 0.
4 is in the sequence because it is the very smallest composite; nondivisors in the cototient of n are composite and since 4 | 4, both A243822(4) and A243823(4) = 0.
6 is in the sequence because it is the only number for which A243822(6) = 1 but A243823(6) = 0. A272618(6) = 4; 4 | 6^2.
10 is in the sequence because it has 2 semidivisors 4 | 10^2 and 8 | 10^3, while only 1 semitotative 6 = 2 * 3.
14 is not in the sequence since it has 2 semidivisors (4 and 8) but 3 semitotatives (6, 10, and 12).
List of terms n followed by row n of A272618 and A272619:
   1, {}, {}
   4, {}, {}
   6, {4}, {}
  10, {4,8}, {6}
  12, {8,9}, {10}
  18, {4,8,12,16}, {10,14,15}
  30, {4,8,9,12,16,18,20,24,25,27}, {14,21,22,26,28}
		

Crossrefs

Programs

  • Mathematica
    Select[Range@ 30, Function[n, And[! PrimeQ@ n, #2 - #1 >= n - (#2 + #3 - 1)] & @@ {DivisorSigma[0, n], Count[Range@ n, _?(PowerMod[n, #, #] == 0 &)], EulerPhi@ n}]]

A243822 Number of k < n such that rad(k) | n but k does not divide n, where rad = A007947.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 2, 0, 2, 1, 0, 0, 4, 0, 2, 1, 3, 0, 3, 0, 3, 0, 2, 0, 10, 0, 0, 2, 4, 1, 5, 0, 4, 2, 3, 0, 11, 0, 3, 2, 4, 0, 5, 0, 6, 2, 3, 0, 8, 1, 3, 2, 4, 0, 14, 0, 4, 2, 0, 1, 14, 0, 4, 2, 12, 0, 6, 0, 5, 3, 4, 1, 15, 0, 4, 0, 5, 0, 16, 1, 5, 3, 3, 0, 20, 1, 4, 3, 5, 1, 8, 0, 7, 2, 6
Offset: 1

Views

Author

Michael De Vlieger, Jun 11 2014

Keywords

Comments

Former name: number of "semidivisors" of n, numbers m < n that do not divide n but divide n^e for some integer e > 1. See ACM Inroads paper.

Examples

			From _Michael De Vlieger_, Aug 11 2024: (Start)
Let S(n) = row n of A162306 and let D(n) = row n of A027750.a(2) = 0 since S(2) \ D(2) = {1, 2} \ {1, 2} is null.
a(10) = 2 since S(10) \ D(10) = {1, 2, 4, 5, 8, 10} \ {1, 2, 5, 10} = {4, 8}.a(16) = 0 since S(16) \ D(16) = {1, 2, 4, 8, 16} \ {1, 2, 4, 8, 16} is null, etc.Table of a(n) and S(n) \ D(n):
   n  a(n)  row n of A272618.
  ---------------------------
   6    1   {4}
  10    2   {4, 8}
  12    2   {8, 9}
  14    2   {4, 8}
  15    1   {9}
  18    4   {4, 8, 12*, 16}
  20    2   {8, 16}
  21    1   {9}
  22    3   {4, 8, 16}
  24    3   {9, 16, 18*}
  26    3   {4, 8, 16}
  28    2   {8, 16}
  30   10   {4, 8, 9, 12, 16, 18, 20, 24, 25, 27}
Terms in A272618 marked with an asterisk are counted by A355432. All other terms are counted by A361235. (End)
		

Crossrefs

Programs

Formula

a(n) = A010846(n) - A000005(n) = card({row n of A162306} \ {row n of A027750}).
a(n) = A045763(n) - A243823(n).
a(n) = (Sum_{1<=k<=n, gcd(n,k)=1} mu(k)*floor(n/k)) - tau(n). - Michael De Vlieger, May 10 2016, after Benoit Cloitre at A010846.
From Michael De Vlieger, Aug 11 2024" (Start)
a(n) = 0 for n in A000961, a(n) > 0 for n in A024619.
a(n) = A051953(n) - A000005(n) + 1 = n - A000010(n) - A000005(n) - A243823(n) + 1.
a(n) = A355432(n) + A361235(n).
a(n) = A355432(n) for n in A360768.
a(n) = A361235(n) for n not in A360768.
a(n) = number of terms in row n of A272618.
a(n) = sum of row n of A304570. (End)

Extensions

New name from David James Sycamore, Aug 11 2024

A272618 Irregular array read by rows: n-th row contains (in ascending order) the nondivisors 1 <= k < n such that all the prime divisors p of k also divide n.

Original entry on oeis.org

0, 0, 0, 0, 0, 4, 0, 0, 0, 4, 8, 0, 8, 9, 0, 4, 8, 9, 0, 0, 4, 8, 12, 16, 0, 8, 16, 9, 4, 8, 16, 0, 9, 16, 18, 0, 4, 8, 16, 0, 8, 16, 0, 4, 8, 9, 12, 16, 18, 20, 24, 25, 27, 0, 0, 9, 27, 4, 8, 16, 32, 25, 8, 16, 24, 27, 32, 0, 4, 8, 16, 32, 9, 27, 16, 25, 32, 0, 4, 8, 9, 12, 16, 18, 24, 27, 28, 32
Offset: 1

Views

Author

Michael De Vlieger, May 03 2016

Keywords

Comments

The k are the "semidivisors" or nondivisor regular numbers of n as counted by A243822(n).
All nonzero terms k are composite and pertain to composite rows n. This is because prime k must either divide or be coprime to n, and k = 1 is both a divisor of and coprime to n.
Row n for prime p contains zero, since numbers 1 <= k < p must either divide or be coprime to prime p.
Row n for prime powers p^e contains zero, since there is only one prime divisor p of p^e and every power 1 <= m <= e of p divides p^e.
Row n = 4 is a special case of composite n that contains zero. This is because 4 is the smallest composite number; there are no composites k < n.
Thus rows n for composite n > 4 contain at least 1 nonzero value.
In base n, 1/a(n) has a terminating expansion with at least 2 places.

Examples

			For n = 12, the numbers 1 <= k < n such that the prime divisors p of k also divide n are {2, 3, 4, 6, 8, 9}; {2, 3, 4, 6} divide n = 12, thus row n = 12 is {8, 9}.
n: k
1: 0
2: 0
3: 0
4: 0
5: 0
6: 4
7: 0
8: 0
9: 0
10: 4 8
11: 0
12: 8 9
13: 0
14: 4 8
15: 9
16: 0
17: 0
18: 4 8 12 16
19: 0
20: 8 16
		

References

  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, pp. 144-145, Theorem 136.

Crossrefs

Union of A027750 and nonzero terms of a(n) = A162306, thus A000005(n) + A243822(n) = A010846(n).
The union of nonzero terms of a(n) and A272619 = A133995, thus A243822(n) + A243823(n) = A045763(n).

Programs

  • Mathematica
    Table[With[{r = First /@ FactorInteger@ n}, Select[Range@ n,
    And[SubsetQ[r, Map[First, FactorInteger@ #]], ! Divisible[n, #]] &]], {n, 30}] /. {} -> 0 // Flatten (* Michael De Vlieger, May 03 2016 *)

A272619 Irregular array read by rows: n-th row contains (in ascending order) the numbers 1 <= k < n such that at least one prime divisor p of k also divides n and at least one prime divisor q of k is coprime to n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 6, 6, 6, 0, 10, 0, 6, 10, 12, 6, 10, 12, 6, 10, 12, 14, 0, 10, 14, 15, 0, 6, 12, 14, 15, 18, 6, 12, 14, 15, 18, 6, 10, 12, 14, 18, 20, 0, 10, 14, 15, 20, 21, 22, 10, 15, 20, 6, 10, 12, 14, 18, 20, 22, 24, 6, 12, 15, 18, 21, 24, 6, 10, 12, 18, 20, 21, 22, 24, 26, 0, 14, 21, 22
Offset: 1

Views

Author

Michael De Vlieger, May 03 2016

Keywords

Comments

The k are the "semitotatives" of n as counted by A243823(n).
All nonzero terms k are composite and pertain to composite rows n. This is because prime k must either divide or be coprime to n, and k = 1 is both a divisor of and coprime to n. Further, the terms k must have at least two distinct prime divisors p and q.
Row n for prime p contains zero, since numbers 1 <= k < p must either divide or be coprime to prime p.
Row n for prime powers p^e contains all the numbers k in the corresponding row of A133995. There is only one prime divisor p of p^e and every power 1 <= m <= e of p divides p^e, thus none of the terms of the corresponding row of A133995 are in A272618(n).
Rows n = 4 and 6 are special cases of composite n that contains zero. 4 is the smallest composite number; there are no composites k < n. 6 has the prime divisors 2 and 3, thus 5 is the smallest prime coprime to 6; the product of the minimum prime divisor and minimum prime coprime to 6 is 10, which exceeds 6 and falls outside the considered range. The situation is not so for composite n > 6. Thus rows n for composite n > 6 contain at least 1 nonzero value.
The smallest k of row n = A096014(n) < n, i.e., those values of A096014(n) pertaining to composite n > 6, a product of the smallest prime divisor p of n and the smallest prime q coprime to n. The smallest k of n are even squarefree semiprimes since 2 either divides n or is coprime to n and k is by definition a number with at least two distinct primes. The smallest k = 2p for p^2 sets record values for A096014(n) when we ignore values pertaining to prime n, n = 4, and n = 6.
In base n, 1/a(n) has a mixed recurrent expansion.

Examples

			For n = 12, the numbers 1 <= k < n such that the prime divisors p of k also divide n are {2, 3, 4, 6, 8, 9}; {2, 3, 4, 6} divide n = 12, thus row n = 12 is {8, 9}.
n:   k
1:   0
2:   0
3:   0
4:   0
5:   0
6:   0
7:   0
8:   6
9:   6
10:  6
11:  0
12: 10
13:  0
14:  6 10 12
15:  6 10 12
16:  6 10 12 14
17:  0
18: 10 14 15
19:  0
20:  6 12 14 15 18
		

References

  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, pp. 144-5, Theorem 136.

Crossrefs

The union of nonzero terms of a(n) and A272618 = A133995, thus A243822(n) + A243823(n) = A045763(n).

Programs

  • Mathematica
    Table[With[{r = First /@ FactorInteger@ n}, Select[Range@ n, Function[m, And[! SubsetQ[r, First /@ FactorInteger@ m], 1 < GCD[m, n] < n]]]], {n, 30}] /. {} -> {0} // Flatten (* Michael De Vlieger, May 03 2016 *)
Showing 1-10 of 22 results. Next