login
a(n) = smallest prime which is > smallest prime dividing n and is coprime to n.
1

%I #17 Feb 18 2018 15:31:22

%S 2,3,5,3,7,5,11,3,5,3,13,5,17,3,7,3,19,5,23,3,5,3,29,5,7,3,5,3,31,7,

%T 37,3,5,3,11,5,41,3,5,3,43,5,47,3,7,3,53,5,11,3,5,3,59,5,7,3,5,3,61,7,

%U 67,3,5,3,7,5,71,3,5,3,73,5,79,3,7

%N a(n) = smallest prime which is > smallest prime dividing n and is coprime to n.

%H Harvey P. Dale, <a href="/A117369/b117369.txt">Table of n, a(n) for n = 1..1000</a>

%e a(6) = 5 because 5 is the smallest prime which is both greater than the smallest prime dividing 6, which is 2 and is coprime to 6.

%t a[1] := 2; a[n_] := Module[{}, k = PrimePi[FactorInteger[n][[1, 1]]]; k++; While[Not[GCD[Prime[k], n] == 1 ], k++ ]; Prime[k]]; Table[a[i], {i, 1, 80}] (* _Stefan Steinerberger_ and Patrick Hanslmaier, Jun 03 2007 *)

%t spdn[n_]:=Module[{s=FactorInteger[n][[1,1]],p},p=NextPrime[s];While[ !CoprimeQ[ p,n],p=NextPrime[p]];p]; Array[spdn,80] (* _Harvey P. Dale_, Feb 18 2018 *)

%Y Cf. A079068, A117367.

%K nonn

%O 1,1

%A _Leroy Quet_, Mar 10 2006

%E More terms from _Stefan Steinerberger_ and Patrick Hanslmaier, Jun 03 2007