1,2

a(n) = A000041(n) iff n=1 or n is prime. [corrected by Robert Israel, Sep 04 2018]

a(2*n) = n+1.

Robert Israel, Table of n, a(n) for n = 1..10000

f:= n -> combinat:-numbpart(n, min(numtheory:-factorset(n))):

f(1):= 1:

map(f, [$1..100]); # Robert Israel, Sep 04 2018

Cf. A000041, A020639, A097360.

Sequence in context: A323952 A168140 A088100 * A208525 A209574 A079387

Adjacent sequences: A097356 A097357 A097358 * A097360 A097361 A097362

nonn

Reinhard Zumkeller, Aug 09 2004

approved