|
%I #33 May 15 2023 08:43:43
%S 3,4,5,6,7,4,3,5,11,6,13,7,3,4,17,6,19,4,3,11,23,6,5,13,3,4,29,10,31,
%T 4,3,17,5,6,37,19,3,4,41,6,43,4,3,23,47,6,7,5,3,4,53,6,5,4,3,29,59,10,
%U 61,31,3,4,5,6,67,4,3,5,71,6,73,37,3,4,7,6,79
%N Smallest divisor of n that is greater than the smallest prime not dividing n (A053669(n)).
%C Sequence is similar to A126800 but differs for the first time at n = 30 and thereafter at n = 30k, where k = 3, 5, 7, 9, 11, 13, 14, 15 ... The generating function for k is not known.
%C The sequence starts at n = 3 as there are no qualifying divisors for n = 1 or n = 2.
%H K. Spage, <a href="/A245689/b245689.txt">Table of n, a(n) for n = 3..1000</a>
%F a(n) = n if n is an odd prime.
%e For n = 30 the smallest prime not dividing n is 7 and the smallest divisor of 30 that is greater than 7 is 10, so a(30) = 10.
%p a:= proc(n)
%p uses numtheory;
%p local F,p;
%p if n::odd then p:= 2
%p else
%p F:= map(pi,factorset(n));
%p p:= ithprime(min(map(`+`,F,1) minus F));
%p fi;
%p min(select(`>`,divisors(n),p));
%p end proc:
%p seq(a(n),n=3..100); # _Robert Israel_, Jul 31 2014
%t A053669[n_] := Module[{p}, For[p = 2, True, p = NextPrime[p], If[CoprimeQ[n, p], Return[p]]]];
%t A245689[n_] := SelectFirst[Divisors[n], # > A053669[n]&];
%t Table[A245689[n], {n, 3, 100}] (* _Jean-François Alcover_, May 15 2023 *)
%o (PARI) A053669(n)={forprime(p=2, ,if(n%p, return(p)))}
%o A245689(n) ={my(c=A053669(n)+1);while(n%c,c++);c}
%Y Cf. A053669, A245690.
%K nonn
%O 3,1
%A _K. Spage_, Jul 29 2014
|
