%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