A051063 -id:A051063 - cp's OEIS frontend

A276936 Numbers m with at least one distinct prime factor prime(k) such that prime(k)^k divides, but prime(k)^(k+1) does not divide m.

2, 6, 9, 10, 14, 18, 22, 26, 30, 34, 36, 38, 42, 45, 46, 50, 54, 58, 62, 63, 66, 70, 72, 74, 78, 82, 86, 90, 94, 98, 99, 102, 106, 110, 114, 117, 118, 122, 125, 126, 130, 134, 138, 142, 144, 146, 150, 153, 154, 158, 162, 166, 170, 171, 174, 178, 180, 182, 186, 190, 194, 198, 202, 206, 207, 210, 214, 218, 222, 225

Offset: 1

Antti Karttunen, Sep 24 2016

nonn

Numbers m with at least one prime factor such that the exponent of its highest power in m is equal to the index of that prime.

The asymptotic density of this sequence is 1 - Product_{k>=1} (1 - 1/prime(k)^k + 1/prime(k)^(k+1)) = 0.31025035294364447031... - Amiram Eldar, Jan 09 2021

			2 is a member as 2 = prime(1) and as 2^1 divides but 2^2 does not divide 2.
3 is NOT a member as 3 = prime(2) but 3^2 does not divide 3.
4 is NOT a member as 2^2 divides 4.
6 is a member as 2 = prime(1) and 2^1 is a divisor of 6, but 2^2 is not.
9 is a member as 3 = prime(2) and 3^2 divides 9.

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..5500 from Antti Karttunen)

Cf. A109298, A276935.
Intersection with A276078 gives A276937.
Cf. A016825, A051063 (subsequences).
Complement of A325130.

q:= n-> ormap(i-> numtheory[pi](i[1])=i[2], ifactors(n)[2]):
select(q, [$1..225])[];  # Alois P. Heinz, Nov 18 2024

Select[Range[225], AnyTrue[FactorInteger[#], PrimePi[First[#1]] == Last[#1] &] &] (* Amiram Eldar, Jan 09 2021 *)

A373318 Numerator of the asymptotic density of numbers that are unitarily divided by n.

Original entry on oeis.org

1, 1, 2, 1, 4, 1, 6, 1, 2, 1, 10, 1, 12, 3, 8, 1, 16, 1, 18, 1, 4, 5, 22, 1, 4, 3, 2, 3, 28, 2, 30, 1, 20, 4, 24, 1, 36, 9, 8, 1, 40, 1, 42, 5, 8, 11, 46, 1, 6, 1, 32, 3, 52, 1, 8, 3, 4, 7, 58, 1, 60, 15, 4, 1, 48, 5, 66, 2, 44, 6, 70, 1, 72, 9, 8, 9, 60, 2, 78

Offset: 1

Amiram Eldar, Jun 01 2024

Numbers that are unitarily divided by n are numbers k such that n is a unitary divisor of k, or equivalently, numbers of the form m*n, with gcd(m, n) = 1.

			Fractions begin with: 1, 1/4, 2/9, 1/8, 4/25, 1/18, 6/49, 1/16, 2/27, 1/25, 10/121, 1/36, ...
For n = 2, the numbers that are unitarily divided by 2 are the numbers of the form 4*k+2 whose asymptotic density is 1/4. Therefore a(2) = numerator(1/4) = 1.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Unitary Divisor.
Wikipedia, Unitary divisor.

Cf. A000010, A003277, A034444, A064549, A076512, A077610, A109395, A173557, A373319(denominators), A373320.
Cf. A001620, A013661, A306016.
Numbers that are unitarily divided by k: A000027 (k=1), A016825 (k=2), A016051 (k=3), A017113 (k=4), A051062 (k=8), A051063 (k=9).

a[n_] := Numerator[EulerPhi[n]/n^2]; Array[a, 100]

PARI
```
a(n) = numerator(eulerphi(n)/n^2);
```

a(n) = 1 if and only if n is in A090778.
a(n) = A000010(n) if and only if n is a cyclic number (A003277).
Let f(n) = a(n)/A373319(n). Then:
f(n) = A000010(n)/n^2 = A076512(n)/(n*A109395(n)).
f(n) = A173557(n)/A064549(n).
f(n) is multiplicative with f(p^e) = (1 - 1/p)/p^e.
Sum_{k=1..n} f(k) = (log(n) + gamma - zeta'(2)/zeta(2)) / zeta(2), where gamma is Euler's constant (A001620).

cp's OEIS Frontend

A276936 Numbers m with at least one distinct prime factor prime(k) such that prime(k)^k divides, but prime(k)^(k+1) does not divide m.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica