login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

Number of iterations of "tamed variant of arithmetic derivative", A327965 needed to reach 0 from n, or -1 if zero is never reached.
10

%I #10 Oct 01 2019 19:51:54

%S 0,1,2,2,2,2,3,2,3,4,3,2,2,2,5,3,3,2,5,2,4,4,3,2,3,4,4,2,3,2,3,2,3,6,

%T 3,3,4,2,5,2,3,2,3,2,3,3,5,2,3,6,4,3,6,2,3,2,3,4,3,2,3,2,7,4,3,6,3,2,

%U 6,5,3,2,3,2,3,3,3,6,3,2,3,2,3,2,3,4,4,3,4,2,4,3,4,4,7,4,3,2,7,4,4,2,4,2,3,3,3,2,3,2,4,4,4,2,3,3,4,4,3,4,3

%N Number of iterations of "tamed variant of arithmetic derivative", A327965 needed to reach 0 from n, or -1 if zero is never reached.

%C Conjecture: from all n, zero is eventually reached.

%H Antti Karttunen, <a href="/A327966/b327966.txt">Table of n, a(n) for n = 0..90609</a>

%F a(0) = 0; for n > 0, a(n) = 1 + a(A327965(n)).

%F a(p) = 2 for all primes p.

%o (PARI)

%o A003415(n) = {my(fac); if(n<1, 0, fac=factor(n); sum(i=1, matsize(fac)[1], n*fac[i, 2]/fac[i, 1]))}; \\ From A003415

%o A327938(n) = { my(f = factor(n)); for(k=1, #f~, f[k,2] = (f[k,2]%f[k,1])); factorback(f); };

%o A327965(n) = if(n<=1,0,A327938(A003415(n)));

%o A327966(n) = { my(k=0); while(n>0, k++; n = A327965(n)); (k); };

%o \\ Or alternatively, as a recurrence:

%o A327966(n) = if(!n,0,1+A327966(A327965(n)));

%Y Cf. A003415, A256750, A327938, A327965, A327967 (indices of the records).

%K nonn

%O 0,3

%A _Antti Karttunen_, Oct 01 2019