login
Totally multiplicative sequence with a(p) = 31.
1

%I #13 Sep 08 2022 08:45:48

%S 1,31,31,961,31,961,31,29791,961,961,31,29791,31,961,961,923521,31,

%T 29791,31,29791,961,961,31,923521,961,961,29791,29791,31,29791,31,

%U 28629151,961,961,961,923521,31,961,961,923521,31,29791,31,29791,29791

%N Totally multiplicative sequence with a(p) = 31.

%H G. C. Greubel, <a href="/A165852/b165852.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = A009975(A001222(n)) = 31^bigomega(n) = 31^A001222(n).

%t 31^PrimeOmega[Range[100]] (* _G. C. Greubel_, Apr 13 2016 *)

%o (PARI) a(n) = 31^bigomega(n); \\ _Michel Marcus_, Apr 14 2016

%o (Magma) [n eq 1 select 1 else 31^(&+[p[2]: p in Factorization(n)]): n in [1..100]]; // _Vincenzo Librandi_, Apr 14 2016

%Y Cf. A001222, A009975.

%K nonn,mult

%O 1,2

%A _Jaroslav Krizek_, Sep 28 2009