A086571 a(1)=1, for n>1: a(n) = smallest divisor of n occurring earlier at most once.

1, 1, 3, 2, 5, 2, 7, 4, 3, 5, 11, 4, 13, 7, 15, 8, 17, 6, 19, 10, 21, 11, 23, 6, 25, 13, 9, 14, 29, 10, 31, 8, 33, 17, 35, 9, 37, 19, 39, 20, 41, 14, 43, 22, 15, 23, 47, 12, 49, 25, 51, 26, 53, 18, 55, 28, 57, 29, 59, 12, 61, 31, 21, 16, 65, 22, 67, 34, 69, 35, 71, 18, 73

Offset: 1

Author

Reinhard Zumkeller, Jul 22 2003

Keywords

• nonn

Comments

a(A086939(n)) = a(A086938(n)) = n. - Reinhard Zumkeller, Jul 24 2003

Examples

			a(p) = p for primes p>2;
divisor set of n=12: {1,2,3,4,6,12}, divisors occurring <12: 1=a(1)=a(2), 2=a(4)=a(6), 3=a(3)=a(9) and 4=a(8); as 4 occurs only once a(12)=4.
		

Links

• Ivan Neretin, Table of n, a(n) for n = 1..10000

Crossrefs

Cf. A086572, A086938, A086939.

Programs

Mathematica

a = {1}; twice = {}; Do[nxt = Complement[Divisors[n], twice][[1]]; If[MemberQ[a, nxt], AppendTo[twice, nxt]]; AppendTo[a, nxt], {n, 2, 10}]; a (* Ivan Neretin, May 28 2015 *)