login
a(1) = 1, a(2) = 1, a(n) = n * LeastPrimeFactor(n-1) / LeastPrimeFactor(n)
5

%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