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 23 results. Next

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 *)

A133995 Irregular array read by rows: n-th row contains (in numerical order) the positive integers <= n which are neither divisors of n nor are coprime to n. A 0 is put into row n if there are no such integers.

Original entry on oeis.org

0, 0, 0, 0, 0, 4, 0, 6, 6, 4, 6, 8, 0, 8, 9, 10, 0, 4, 6, 8, 10, 12, 6, 9, 10, 12, 6, 10, 12, 14, 0, 4, 8, 10, 12, 14, 15, 16, 0, 6, 8, 12, 14, 15, 16, 18, 6, 9, 12, 14, 15, 18, 4, 6, 8, 10, 12, 14, 16, 18, 20, 0, 9, 10, 14, 15, 16, 18, 20, 21, 22, 10, 15, 20, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 6, 12, 15, 18, 21, 24
Offset: 1

Views

Author

Leroy Quet, Oct 01 2007

Keywords

Comments

Row n has length A264441(n).
The number of nonzero entries in row n is A045763(n).
Row n has a 0 if every positive integer <= n is coprime to n or divides n.
From Michael De Vlieger, Aug 19 2017: (Start)
When row n is not empty (and here represented by 0), the terms of row n are composite, since primes p < n must either divide or be coprime to n and the empty product 1 both divides and is coprime to all numbers. For the following, let prime p divide n and prime q be coprime to n.
Row n is empty for n < 8 except n = 6.
There are two distinct species of term m of row n. The first are nondivisor regular numbers g in A272618(n) that divide some integer power e > 1 of n. In other words, these numbers are products of primes p that also divide n and no primes q that are coprime to n, yet g itself does not divide n. Prime powers n = p^k cannot have numbers g in A272618(n) since they have only one distinct prime divisor p; all regular numbers g = p^e with 0 <= e <= k divide p^k. The smallest n = 6 for which there is a number in A272618. The number 4 is the smallest composite and is equal to n = 4 thus must divide it; 4 is coprime to 5. The number 4 is neither coprime to nor a divisor of 6.
The second are numbers h in A272619(n) that are products of at least one prime p that divides n and one prime q that is coprime to n.
The smallest n = 8 for which there is a number in A272619 is 8; the number 6 is the product of the smallest two distinct primes. 6 divides 6 and is coprime to 7. The number 6 is neither coprime to nor a divisor of the prime power 8; 4 divides 8 and does not appear in a(8).
There can be no other species since primes p <= n divide n and q < n are coprime to n, and products of primes q exclusive of any p are coprime to n.
As a consequence of these two species, rows 1 <= n <= 5 and n = 7 are empty and thus have 0 in row n.
(End)

Examples

			The divisors of 12 are: 1,2,3,4,6,12. The positive integers which are <= 12 and are coprime to 12 are: 1,5,7,11. So row 12 contains the positive integers <= 12 which are in neither of these two lists: 8,9,10.
The irregular triangle T(n, k) begins:
n\k 1  2  3  4  5  6  7 ...
1:  0
2:  0
3:  0
4:  0
5:  0
6:  4
7:  0
8:  6
9:  6
10: 4  6  8
11: 0
12: 8  9 10
13: 0
14: 4  6  8 10 12
15: 6  9 10 12
16: 6 10 12 14
17: 0
18: 4  8 10 12 14 15 16
19: 0
20: 6  8 12 14 15 16 18
... formatted by _Wolfdieter Lang_, Jan 16 2016
		

Crossrefs

Programs

  • Maple
    row:= proc(n) local r;
       r:= remove(t -> member(igcd(t, n), [1, t]), [$1..n]):
       if r = [] then 0 else op(r) fi
    end proc:
    A:= [seq](row(n), n=1..30); # Robert Israel, Jan 19 2016
  • Mathematica
    Table[Select[Range@ n, Nor[Divisible[n, #], CoprimeQ[n, #]] &] /. {} -> {0}, {n, 27}] // Flatten (* Michael De Vlieger, Aug 19 2017 *)

Formula

a(n) = union(A272618(n), A272619(n)). - Michael De Vlieger, Aug 19 2017

Extensions

More terms from Alvin Hoover Belt, Jan 21 2008
Edited by Wolfdieter Lang, Jan 16 2016

A299990 a(n) = A243822(n) - A000005(n).

Original entry on oeis.org

-1, -2, -2, -3, -2, -3, -2, -4, -3, -2, -2, -4, -2, -2, -3, -5, -2, -2, -2, -4, -3, -1, -2, -5, -3, -1, -4, -4, -2, 2, -2, -6, -2, 0, -3, -4, -2, 0, -2, -5, -2, 3, -2, -3, -4, 0, -2, -5, -3, 0, -2, -3, -2, 0, -3, -5, -2, 0, -2, 2, -2, 0, -4, -7, -3, 6, -2, -2
Offset: 1

Views

Author

Michael De Vlieger, Feb 25 2018

Keywords

Comments

Since A010846(n) = A000005(n) + A243822(n), this sequence examines the balance of the two components among "regular" numbers.
Value of a(n) is generally less frequently negative as n increases.
a(1) = -1.
For primes p, a(p) = -2 since 1 | p and the cototient is restricted to the divisor p.
For perfect prime powers p^e, a(p^e) = -(e + 1), since all m < p^e in the cototient of p^e that do not have a prime factor q coprime to p^e are powers p^k with 1 < p^k <= p^e; all such p^k divide p^e.
Generally for n with A001221(n) = 1, a(n) = -1 * A000005(n), since the cototient is restricted to divisors, and in the case of p^e > 4, divisors and numbers in A272619(p^e) not counted by A010846(p^e).
For m >= 3, a(A002110(m)) is positive.
For m >= 9, a(A244052(m)) is positive.

Examples

			a(6) = -3 since 6 has 4 divisors, and 4 | 6^2; A243822(6) = 1 and A000005(6) = 4; 1 - 4 = -3. Alternatively, A010846(6) = 5; 5 - 2*4 = -3.
a(30) = 2 since 30 has 8 divisors and the numbers {4, 8, 9, 12, 16, 18, 20, 24, 25, 27} divide 30^e with e > 1; A243822(30) = 10 and A000005(30) = 8; 10 - 8 = 2. Alternatively, A010846(30) = 18; 18 - 2*8 = 2.
Some values of a(n) and related sequences:
   n  a(n) A010846(n) A243822(n) A000005(n) A272618(n)
  ----------------------------------------------------
   1   -1          1          0          1  0
   2   -2          2          0          2  0
   3   -2          2          0          2  0
   4   -3          3          0          3  0
   5   -2          2          0          2  0
   6   -3          5          1          4  {4}
   7   -2          2          0          2  0
   8   -4          4          0          4  0
   9   -3          3          0          3  0
  10   -2          6          2          4  {4,8}
  11   -2          2          0          2  0
  12   -4          8          2          6  {8,9}
  ...
  30    2         18         10          8  {4,8,9,12,16,18,20,24,25,27}
  ...
  34    0          8          4          4  {4,8,16,32}
  ...
		

Crossrefs

Programs

  • Mathematica
    Table[Count[Range[n], _?(PowerMod[n, Floor@ Log2@ n, #] == 0 &)] - 2 DivisorSigma[0, n], {n, 68}]

Formula

a(n) = A010846(n) - 2*A000005(n).

A360543 a(n) = number of numbers k < n, gcd(k, n) > 1, such that omega(k) > omega(n) and rad(n) | rad(k), where omega(n) = A001221(n) and rad(n) = A007947(n).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 6, 0, 0, 0, 0, 11, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 2, 5, 1, 0, 0, 0, 2, 0, 1, 0, 0, 0, 0, 0, 0, 1, 26, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 2, 0, 0, 0, 0, 3, 23, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 7, 0, 3, 1, 4
Offset: 1

Views

Author

Michael De Vlieger, Mar 06 2023

Keywords

Examples

			a(4) = 0 since k = 1..3 are prime powers.
a(8) = 1 since only k = 6 is such that p = 3, q = 5, but gcd(6, 10) = 2.
a(9) = 1 since the following satisfies definition: {6},
a(16) = 4, i.e., {6, 10, 12, 14},
a(25) = 3, i.e., {10, 15, 20},
a(27) = 6, i.e., {6, 12, 15, 18, 21, 24},
a(32) = 11, i.e., {6, 10, 12, 14, 18, 20, 22, 24, 26, 28, 30},
a(36) = 1, i.e., {30},
a(40) = 1, i.e., {30},
a(45) = 1, i.e., {30}, etc.
		

Crossrefs

Programs

  • Mathematica
    nn = 120; rad[n_] := rad[n] = Times @@ FactorInteger[n][[All, 1]]; c = Select[Range[4, nn], CompositeQ]; Table[Function[{q, r}, Count[TakeWhile[c, # <= n &], _?(And[PrimeNu[#] > q, Divisible[rad[#], r]] &)]] @@ {PrimeNu[n], rad[n]}, {n, nn}]

Formula

a(n) = A243823(n) - A360480(n).
a(n) = A045763(n) - A243822(n) - A360480(n).
a(n) = A051953(n) - A000005(n) - A243822(n) - A360480(n).
a(n) = A051953(n) - A010846(n) - A360480(n).
a(n) = A243823(n) = A045763(n) for n in A246547.
For prime power n = p^e, n > 1, a(n) = p^(e-1) - e.
For n in A360765, a(n) > 0.

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

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).

A360224 Number of k < n such that gcd(k, n) > 1, gcd(n^2-1, k) = 1, and rad(k) does not divide n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 3, 0, 0, 0, 2, 0, 3, 0, 2, 1, 4, 0, 5, 0, 5, 2, 1, 0, 4, 0, 6, 2, 6, 0, 12, 0, 5, 3, 7, 0, 14, 0, 5, 2, 10, 0, 11, 0, 4, 5, 13, 0, 19, 0, 12, 7, 7, 1, 13, 0, 14, 3, 11, 0, 31, 0, 13, 9, 8, 2, 21, 0, 19, 7, 21, 0, 18, 2, 13, 9, 22, 0, 21, 1, 16, 10, 16
Offset: 1

Views

Author

Michael De Vlieger, May 19 2023

Keywords

Comments

Number of terms in row n of A272619 that are coprime to (n-1)*(n+1).

Examples

			Let S(n) = row n of A272619.
a(p) = 0 since S(p) is empty.
a(4) = 0 since S(4) is empty.
a(6) = 0 since S(6) is empty.
a(8) = 0 since S(8) = {6}, but gcd(6,(8+1)) = 3.
a(10) = 0 since S(10) = {6}, but gcd(6,(10-1)) = 3.
a(12) = 1 since S(12) = {10}, and gcd(10,143) = 1.
a(16) = 1 since S(16) = {6, 10, 12, 14}, but only 14 is such that gcd(14, 255) = 1.
a(18) = 3 since S(18) = {10, 14, 15}, and none of these share a prime factor with 323.
a(20) = 0 since S(20) = {6, 12, 14, 15, 18}, and all of these share a factor with 21.
		

Crossrefs

Programs

  • Mathematica
    Table[Count[Range[k], _?(Nor[CoprimeQ[#, k], GCD[k^2 - 1, #] > 1, Divisible[k, Times @@ FactorInteger[#][[All, 1]]]] &)], {k, 120}]
  • PARI
    rad(n) = factorback(factorint(n)[, 1]); \\ A007947
    a(n) = sum(k=1, n-1, (gcd(k,n)>1) && (gcd(n^2-1, k) == 1) && (n % rad(k))); \\ Michel Marcus, May 20 2023

Formula

a(n) <= A243822(n).

A360480 a(n) = number of numbers k < n, with gcd(k, n) > 1, such that there is at least one prime divisor p | k that does not divide n, and at least one prime divisor q | n that does not divide k.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 3, 3, 0, 0, 3, 0, 5, 5, 6, 0, 6, 0, 8, 0, 9, 0, 5, 0, 0, 8, 11, 7, 10, 0, 13, 10, 13, 0, 12, 0, 16, 13, 17, 0, 16, 0, 18, 14, 20, 0, 19, 11, 21, 16, 23, 0, 19, 0, 25, 19, 0, 13, 25, 0, 27, 20, 27, 0, 27, 0, 30, 25, 31, 13, 32, 0, 32, 0, 34, 0, 33, 17, 36, 25, 37
Offset: 1

Views

Author

Michael De Vlieger, Feb 28 2023

Keywords

Comments

a(n) = 0 for prime powers, since the definition implies omega(n) >= 2.
a(n) is the cardinality of k in the cototient of n such that rad(k) mod rad(n) != 0 and rad(n) mod rad(k) != 0.

Examples

			a(6) = 0 since k = 1..5 are prime powers.
a(10) = 1 since only k = 6 is such that p = 3, q = 5, but gcd(6, 10) = 2.
a(14) = 3 since for k = 6, we have p = 3, q = 7, for k = 10, we have p = 5, q = 7, and for k = 12, we have the same situation as we have with 6, while 6, 10, 12, and 14 are even.
Table listing k <= n counted by row n = 10..22 of this sequence:
a(10) = 1:  6  .  .  .  .
a(11) = 0:  .  .  .  .  .  .
a(12) = 1:  .  .  .  .  10 .  .
a(13) = 0:  .  .  .  .  .  .  .  .
a(14) = 3:  6  .  .  .  10 .  12 .  .
a(15) = 3:  6  .  .  .  10 .  12 .  .  .
a(16) = 0:  .  .  .  .  .  .  .  .  .  .  .
a(17) = 0:  .  .  .  .  .  .  .  .  .  .  .  .
a(18) = 3:  .  .  .  .  10 .  .  .  14 15 .  .  .
a(19) = 0:  .  .  .  .  .  .  .  .  .  .  .  .  .  .
a(20) = 5:  6  .  .  .  .  .  12 .  14 15 .  .  18 .  .
a(21) = 5:  6  .  .  .  .  .  12 .  14 15 .  .  18 .  .  .
a(22) = 6:  6  .  .  .  10 .  12 .  14 .  .  .  18 .  20 .  .
		

Crossrefs

Programs

  • Mathematica
    rad[x_] := rad[x] = Times @@ FactorInteger[x][[All, 1]]; Table[k = rad[n]; Count[Range[n], _?(Nor[CoprimeQ[#1, n], Divisible[#2, k], Divisible[k, #2]] & @@ {#, rad[#]} &)], {n, 88}]

Formula

a(n) <= A243823(n).
a(n) = | { k < n : gcd(k,n) > 1 AND (rad(k)|rad(n) NOR rad(n)|rad(k)) } |.

A381497 a(n) = sum of numbers k < n such that 1 < gcd(k,n) and rad(k) != rad(n), where rad = A007947.

Original entry on oeis.org

0, 0, 0, 0, 0, 9, 0, 6, 6, 25, 0, 36, 0, 49, 45, 42, 0, 81, 0, 100, 84, 121, 0, 144, 45, 169, 96, 196, 0, 315, 0, 210, 198, 289, 175, 354, 0, 361, 273, 430, 0, 609, 0, 484, 435, 529, 0, 648, 140, 655, 459, 676, 0, 801, 385, 826, 570, 841, 0, 1260, 0, 961, 798
Offset: 1

Views

Author

Michael De Vlieger, Mar 02 2025

Keywords

Comments

Analogous to A066760(n), the sum of row n of A133995, and A381499(n), sum of row n of A272619.

Examples

			Table of n and a(n) for select n, showing prime power decomposition of both and row n of A381094:
   n   Factor(n) a(n)  Factor(a(n))  Row n of A381094
  -------------------------------------------------------------------
   6   2 * 3       9   3^2           {2,3,4}
   8   2^3         6   2 * 3         {6}
   9   3^2         6   2 * 3         {6}
  10   2 * 5      25   5^2           {2,4,5,6,8}
  12   2^2 * 3    36   2^2 * 3^2     {2,3,4,8,9,10}
  14   2 * 7      49   7^2           {2,4,6,7,8,10,12}
  15   3 * 5      45   3^2 * 5       {3,5,6,9,10,12}
  16   2^4        42   2 * 3 * 7     {6,10,12,14}
  18   2 * 3^2    81   3^4           {2,3,4,8,9,10,14,15,16}
  20   2^2 * 5   100   2^2 * 5^2     {2,4,5,6,8,12,14,15,16,18}
  21   3 * 7      84   2^2 * 3 * 7   {3,6,7,9,12,14,15,18}
  22   2 * 11    121   11^2          {2,4,6,8,10,11,12,14,16,18,20}
  24   2^3 * 3   144   2^4 * 3^2     {2,3,4,8,9,10,14,15,16,20,21,22}
a(6) = (2+4) + (3) = 9,
a(n) = 6 for n in {8, 9} since 6 is the only number less than n that shares a factor with n but does not have the same squarefree kernel as n.
a(10) = (2+4+6+8) + (5) = 25.
a(12) = (2+4+8+10) + (3+9) = 36.
a(14) = (2+4+6+8+10+12) + (7) = 49.
a(15) = (3+6+9+12) + (5+10) = 45.
a(16) = (6+10+12+14) = 42, etc.
		

Crossrefs

Programs

  • Mathematica
    rad[x_] := rad[x] = Times @@ FactorInteger[x][[All, 1]]; Table[r = rad[n]; Total@ Select[Range[n], Nor[CoprimeQ[#, n], rad[#] == r] &], {n, 120}]

Formula

a(n) is the sum of row n of A381094.
a(n) = 0 for prime n and n = 4.
a(n) = A067392(n) - A381498(n).
Showing 1-10 of 23 results. Next