login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
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. 1
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 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
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.
LINKS
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.
EXAMPLE
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.
MATHEMATICA
nn = 2^10; 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 *)
PROG
(PARI) 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
CROSSREFS
Sequence in context: A100602 A183979 A222175 * A315856 A265360 A182253
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from Michael De Vlieger, Jan 14 2022
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 13:42 EDT 2024. Contains 371254 sequences. (Running on oeis4.)