%I #2 Sep 24 2013 00:41:43
%S 3,5,743,59,312408461,13669,10920223,85691034670497287,
%T 795398356293665458781,14301636907,4280339305565030602944575375011
%N Primes in A136133 (integer log of harmonic numbers sequence).
%t sopfr = Function[x, Plus @@ Map[Times @@ # &, FactorInteger[x]]]; Select[sopfr /@ HarmonicNumber[Range[100]], PrimeQ]
%Y Cf. A136133, A136121, A136122, A001008, A002805.
%K nonn
%O 0,1
%A _Carlos Alves_, Dec 16 2007