cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A225395 Replace each prime number with its rank in the recursive prime factorization of n.

This page as a plain text file.
%I A225395 #27 Jan 01 2025 22:11:40
%S A225395 1,1,2,1,3,2,4,1,2,3,5,2,6,4,6,1,7,2,8,3,8,5,9,2,3,6,4,4,10,6,11,1,10,
%T A225395 7,12,2,12,8,12,3,13,8,14,5,6,9,15,2,4,3,14,6,16,4,15,4,16,10,17,6,18,
%U A225395 11,8,1,18,10,19,7,18,12,20,2,21,12,6,8,20,12,22,3,2,13,23,8,21,14,20,5,24,6,24,9,22,15,24,2,25,4,10,3,26,14,27,6,24,16,28
%N A225395 Replace each prime number with its rank in the recursive prime factorization of n.
%C A225395 a(A000040(n)) = n, hence all natural numbers appear in this sequence.
%C A225395 a(2n) = n.
%C A225395 It appears that a(35) = 12 is the only instance where a composite index yields a larger value than any smaller index. Checked to 10^7. - _Charles R Greathouse IV_, Jul 30 2016
%H A225395 Paul Tek, <a href="/A225395/b225395.txt">Table of n, a(n) for n = 1..10000</a>
%H A225395 Paul Tek, <a href="/A225395/a225395.txt">Perl program for this sequence</a>
%F A225395 Multiplicative, with a(prime(i)^j) = i^a(j).
%F A225395 a(n) = prod(A049084(A027748(k))^a(A124010(k)): k=1..A001221(n)). - _Reinhard Zumkeller_, May 10 2013
%e A225395 The number 9967 is the 1228th prime number.
%e A225395 Hence a(9967) = 1228.
%e A225395 The recursive prime factorization of 31250 is 2*5^(2*3).
%e A225395 The numbers 2, 3 and 5 are respectively the 1st, 2nd and 3rd prime numbers.
%e A225395 Hence a(31250) = a(2*5^(2*3)) = 1*3^(1*2) = 9.
%t A225395 a[1] = 1; a[p_?PrimeQ] := a[p] = PrimePi[p]; a[n_] := a[n] = Times @@ (PrimePi[#[[1]]]^a[#[[2]]]& /@ FactorInteger[n]); Table[a[n], {n, 1, 100}] (* _Jean-François Alcover_, Jun 07 2013 *)
%o A225395 (Perl) # See Tek link.
%o A225395 (Haskell)
%o A225395 a225395 n = product $ zipWith (^)
%o A225395     (map a049084 $ a027748_row n) (map a225395 $ a124010_row n)
%o A225395 -- _Reinhard Zumkeller_, May 10 2013
%o A225395 (PARI) a(n)=if(n<3, return(1)); my(f=factor(n)); prod(i=1,#f~, primepi(f[i,1])^a(f[i,2])) \\ _Charles R Greathouse IV_, Jul 30 2016
%Y A225395 Cf. A000040, A000026, A156061.
%K A225395 nonn,mult,nice
%O A225395 1,3
%A A225395 _Paul Tek_, May 06 2013