|
%I #17 Sep 14 2015 11:41:48
%S 2,2,2,3,3,2,2,3,2,5,2,3,2,3,5,3,2,2,5,3,2,7,2,5,2,3,7,3,2,5,2,3,5,2,
%T 7,2,3,3,7,2,3,2,11,5,2,5,2,3,7,2,3,2,3,5,11,2,3,2,7,5,2,11,3,2,5,7,2,
%U 3,13,2,5,3,13,3,11,2,7,2,5,3,2,2,7,3,5,2,13,7,2,3,5,3,2,11,3,17,2,3
%N Lesser prime factor of n-th product of two distinct primes.
%C a(n)*A070647(n) = A006881(n); a(n) < A070647(n);
%C a(n) = A020639(A006881(n)).
%H Reinhard Zumkeller, <a href="/A096916/b096916.txt">Table of n, a(n) for n = 1..10000</a>
%t f[n_]:=Last/@FactorInteger[n]=={1,1};f1[n_]:=Min[First/@FactorInteger[n]];f2[n_]:=Max[First/@FactorInteger[n]];lst={};Do[If[f[n],AppendTo[lst,f1[n]]],{n,0,6!}];lst (* _Vladimir Joseph Stephan Orlovsky_, Apr 10 2010 *)
%o (Haskell)
%o a096916 = a020639 . a006881 -- _Reinhard Zumkeller_, Sep 23 2011
%o (PARI) go(x)=my(v=List()); forprime(p=2, sqrtint(x\1), forprime(q=p+1, x\p, listput(v, [p*q,p]))); apply(v->v[2], vecsort(Vec(v),1)) \\ _Charles R Greathouse IV_, Sep 14 2015
%Y Cf. A084126, A195758.
%K nonn,look
%O 1,1
%A _Reinhard Zumkeller_, Jul 15 2004
|
