login
If n is the k-th prime power > 1, a(n) = k, otherwise a(n) = 0.
3

%I #14 Jun 17 2021 09:13:53

%S 0,1,2,3,4,0,5,6,7,0,8,0,9,0,0,10,11,0,12,0,0,0,13,0,14,0,15,0,16,0,

%T 17,18,0,0,0,0,19,0,0,0,20,0,21,0,0,0,22,0,23,0,0,0,24,0,0,0,0,0,25,0,

%U 26,0,0,27,0,0,28,0,0,0,29,0,30,0,0,0,0,0,31,0,32,0,33,0,0,0,0,0,34,0,0,0,0,0,0,0,35

%N If n is the k-th prime power > 1, a(n) = k, otherwise a(n) = 0.

%C This is a ("corrected") variant of A095874, which uses the list of "powers of primes" A000961 instead of prime powers A246655. - _M. F. Hasler_, Jun 16 2021

%H Antti Karttunen, <a href="/A322981/b322981.txt">Table of n, a(n) for n = 1..16384</a>

%F a(n) = A010055(n) * A025528(n) = A069513(n) * A025528(n).

%F a(n) = A025528(A069513(n)*n), when assuming that A025528(0) = 0.

%F a(A000961(1+n)) = n for all n >= 1.

%o (PARI)

%o up_to = 16384;

%o A322981list(up_to) = { my(v=vector(up_to), k=0); for(n=1,up_to,if(isprimepower(n),k++; v[n] = k, v[n] = 0)); (v); };

%o v322981 = A322981list(up_to);

%o A322981(n) = v322981[n];

%o (PARI) A322981(n)=if(isprimepower(n),sum(i=1,exponent(n),primepi(sqrtnint(n,i)))) \\ _M. F. Hasler_, Jun 16 2021

%Y Cf. A000961 (powers of primes, including 1), A010055, A025528, A069513, A095874 (analog based on A000961), A246655 (prime powers > 1).

%Y Cf. A049084.

%K nonn

%O 1,3

%A _Antti Karttunen_, Jan 01 2019