A030632 -id:A030632 - cp's OEIS frontend

A065985 Numbers k such that d(k) / 2 is prime, where d(k) = number of divisors of k.

6, 8, 10, 12, 14, 15, 18, 20, 21, 22, 26, 27, 28, 32, 33, 34, 35, 38, 39, 44, 45, 46, 48, 50, 51, 52, 55, 57, 58, 62, 63, 65, 68, 69, 74, 75, 76, 77, 80, 82, 85, 86, 87, 91, 92, 93, 94, 95, 98, 99, 106, 111, 112, 115, 116, 117, 118, 119, 122, 123, 124, 125, 129, 133, 134

Offset: 1

Joseph L. Pe, Dec 10 2001

nonn

Numbers whose sorted prime signature (A118914) is either of the form {2*p-1} or {1, p-1}, where p is a prime. Equivalently, disjoint union of numbers of the form q^(2*p-1) where p and q are primes, and numbers of the form r * q^(p-1), where p, q and r are primes and r != q. - Amiram Eldar, Sep 09 2024

Harry J. Smith, Table of n, a(n) for n = 1..1000

Cf. A000005, A118914.
Numbers with exactly 2*p divisors: A030513 (p=2), A030515 (p=3), A030628 \ {1} (p=5), A030632 (p=7), A137485 (p=11), A137489 (p=13), A175744 (p=17), A175747 (p=19).
Subsequences: A006881, A030078, A050997, A054753, A095990, A096156, A138031, A178739, A179665, A189987.

Select[Range[1, 1000], PrimeQ[DivisorSigma[0, # ] / 2] == True &]

n=0; for (m=1, 10^9, f=numdiv(m)/2; if (frac(f)==0 && isprime(f), write("b065985.txt", n++, " ", m); if (n==1000, return))) \\ Harry J. Smith, Nov 05 2009

is(n)=n=numdiv(n)/2; denominator(n)==1 && isprime(n) \\ Charles R Greathouse IV, Oct 15 2015

A350767 a(1)=1. Thereafter, a(n+1) is the least unused number k such that either d(j(n)) properly divides d(k) or d(k) properly divides d(j(n)), where j(n) = a(n)+1 and d is the divisor counting function A000005.

Original entry on oeis.org

1, 6, 8, 12, 10, 14, 2, 15, 48, 18, 20, 3, 28, 21, 5, 7, 11, 4, 22, 24, 32, 13, 17, 9, 19, 23, 26, 29, 27, 25, 30, 33, 31, 37, 40, 34, 41, 35, 49, 43, 47, 16, 38, 42, 39, 46, 44, 53, 51, 59, 45, 54, 56, 60, 50, 61, 66, 52, 55, 57, 67, 71, 58, 62, 72, 63, 192, 65

Offset: 1

David James Sycamore, Jan 14 2022

nonn

If d(j(n)) is prime p then d(a(n+1)) must be properly divisible by p. In practice the proper divisor for computation of a(n+1) toggles between d(j(n)) and d(k).
Conjecture: This is a permutation of the positive integers. Numbers with the same number (tau) of divisors appear in their natural orders (e.g., primes, semiprimes, squares).
The plot, after the first few terms, resolves itself into points tightly packed on and around a straight line of slope 1, with exceptional points appearing as significant upward or downward "spikes".
When d(j(n)) is prime p appearing for the first time in the sequence J = {d(j(a(n)), n>=1}, then a(n+1) is the smallest number with 2p divisors, which produces a significantly large upward spike above the straight line (6, 12, 48, 192, 3072, 12288, ...).
When d(j(a(n)) is 2p, seen for the first time in J, then a(n+1) is the smallest number with p divisors, which produces a large downward spike, below the straight line (2, 4, 16, 64, 1024, 4096, ...).
The sequence of fixed points starts: 1, 46, 69, 74, 110, 140, 142, 152, 154, 178, ... apparently becoming denser as n increases.

			a(1)=1, so j(1)=2, d(j(1))=2, a prime, so we need the smallest unused k such that d(k) is properly divisible by 2, hence a(2)=6.
a(2)=6, j(2)=4, d(j(2))=3, a prime so we need the smallest unused k such that d(k) is properly divisible by 3, hence a(3)=8.

Michael De Vlieger, Table of n, a(n) for n = 1..10001
Michael De Vlieger, Scatterplot of a(n), n = 1..256, color coded to show primes in red, evens in blue, 2 in magenta, odd composites in black. A gold highlight indicates terms such that (a(n)+1) | a(n+1). Outlying points are annotated.
Michael De Vlieger, Log-log scatterplot of a(n), n = 1..2^16. Terms in A3680(p) such that p is prime appear annotated in red, and those in A061286 appear in blue.

Cf. A000040, A001248, A030514, A054753, A030516, A030626, A085986, A178739, A030629, A030630, A030631, A030632, A030633, A030634, A005179, A061286, A003680.

Block[{c, d, j, k, u, nn}, nn = 120; j = u = 2; c[] = 0; c[1] = 1; Do[d[i] = DivisorSigma[0, i], {i, 2^(Ceiling@ Log2[nn] + 3)}]; {1}~Join~Reap[Do[k = u; While[Nand[Or[Divisible[d[j], d[k]], Divisible[d[k], d[j]]], d[j] != d[k], c[k] == 0], k++]; Sow[k]; c[k] = i; j = k + 1; If[k == u, While[c[u] != 0, u++]], {i, nn}] ][[-1, -1]] ] (* _Michael De Vlieger, Jan 14 2022 *)

isok(k, last, set) = if (!setsearch(set, k), my(ndk=numdiv(k), ndl=numdiv(last+1)); (ndl != ndk) && ((!(ndk % ndl)) || (!(ndl % ndk))));
lista(nn) = {my(last = 1, list = List(last), set = Set(list)); for (n=2, nn, my(k=1); while (!isok(k, last, set), k++); listput(list, k); set = Set(list); last = k;); Vec(list);} \\ Michel Marcus, Jan 15 2022

More terms from Michael De Vlieger, Jan 14 2022

cp's OEIS Frontend

A065985 Numbers k such that d(k) / 2 is prime, where d(k) = number of divisors of k.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

PARI

A350767 a(1)=1. Thereafter, a(n+1) is the least unused number k such that either d(j(n)) properly divides d(k) or d(k) properly divides d(j(n)), where j(n) = a(n)+1 and d is the divisor counting function A000005.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions