%I #10 May 20 2015 10:46:27
%S 1,1,2,6,2,15,2,28,6,15,2,66,2,91,10,24,2,153,2,190,14,33,2,276,10,65,
%T 18,42,2,435,2,496,22,51,14,90,2,703,26,60,2,861,2,946,30,69,2,1128,
%U 14,175,34,78,2,1431,22,140,38,87,2,1770,2,1891,42,96,26,165,2,2278,46,105
%N a(1) = 1, a(2) = 1, a(n) = n * LeastPrimeFactor(n-1) / LeastPrimeFactor(n)
%H Ivan Neretin, <a href="/A126287/b126287.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = A032742(n)*A020639(n-1), for n>2. - _Michel Marcus_, May 20 2015
%e a(6) = 6 * LeastPrimeFactor(5) / LeastPrimeFactor(6) = 6 * 5 / 2 = 15
%p N:= 1000: # to get a(1) to a(n)
%p a[1]:= 1: a[2]:= 1:
%p b:= 2:
%p for n from 3 to N do
%p c:= min(numtheory:-factorset(n));
%p a[n]:= n*b/c;
%p b:= c;
%p od:
%p seq(a[n],n=1..N); # _Robert Israel_, May 20 2015
%t a[1] := 1; a[2] := 1; a[n_] := n*FactorInteger[n - 1][[1, 1]]/FactorInteger[n][[1, 1]]; Table[a[n], {n, 70}] (* _Ivan Neretin_, May 20 2015 *)
%Y Cf. A126286, A126288, A126289, A020639.
%K easy,nonn
%O 1,3
%A _Lior Manor_, Dec 25 2006