%I #47 Jul 12 2018 08:57:35
%S 3,4,3,4,12,6,3,4,3,12,20,24687390,3,72,62,4,20,1102903830,12,
%T 58051620,3,1793172,468,1035844571580,62,882,398,75274140,6,
%U 81206805256038,14,1288005000,78428,93888,37664,24380304369772260,432,3300,21962
%N a(n) is the least x such that the cyclotomic polynomial values Phi(d,x) are prime for all d dividing n.
%C An equivalent formulation is that a(n) is smallest number x such that x^n-1 factors only into its algebraic factors.
%C Many more terms are known, in particular terms at prime indices. Massively composite n are the hardest to find - term 256 alone took a month to find. Contact the author for more terms beyond the gaps.
%C 2 never appears in the sequence because Phi(1,2) = 1, which is irreducible but not prime.
%C a(n) is the smallest number x > 2 such that A001222(x^n-1) = A000005(n). - _Thomas Ordowski_, Jan 31 2018
%C All terms are in A008864. If n is even, a(n) is in A014574. - _Robert Israel_, Jan 31 2018
%H Don Reble, <a href="/A072942/b072942.txt">Table of n, a(n) for n = 1..47</a>
%H Don Reble, <a href="/A072942/a072942.txt">Table of n, a(n) for n = 1..99</a> (with question marks at the unknown entries)
%e a(16)=4 because 4^16-1 = 3*5*17*257*65537, which are the 5 algebraic factors.
%p f:= proc(n) local p, C, x, d;
%p C:= [seq(numtheory:-cyclotomic(d,x), d = numtheory:-divisors(n) minus {1})];
%p p:= 1;
%p do
%p p:= nextprime(p);
%p if andmap(isprime, subs(x=p+1, C)) then return p+1 fi
%p od:
%p end proc:
%p map(f, [$1..29]); # _Robert Israel_, Jan 31 2018
%t Table[With[{d = Divisors@ n}, SelectFirst[Range[10^3], AllTrue[Cyclotomic[d, #], PrimeQ] &]], {n, 11}] (* _Michael De Vlieger_, Jan 31 2018 *)
%o (PARI) for(d=1,17,ds=divisors(d); print("Searching for d|"d":"ds); forprime(p=2,499999,okc=1; for(c=2,length(ds), if(!isprime(subst(polcyclo(ds[c]),x,p+1)),okc=0; break)); if(okc, for(c=1,length(ds), print("Phi("ds[c]","p+1")="subst(polcyclo(ds[c]),x,p+1))); break)))
%o (PARI) isok(n, x) = {fordiv(n, d, if (! isprime(polcyclo(d, x)), return(0));); return(1);}
%o a(n) = {my(x=2); while (! isok(n, x), x++); x;} \\ _Michel Marcus_, Jan 31 2018
%Y Cf. A000005, A001222, A008864, A014574, A070737.
%K nonn
%O 1,1
%A _Phil Carmody_, Aug 12 2002
%E Corrected and extended by _Don Reble_, Feb 03 2014
%E Edited by _N. J. A. Sloane_, Mar 01 2014 at the suggestion of _Phil Carmody_ and _Don Reble_