login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

a(n) is the smallest b > 1 such that when c is equal to any of the first n composites the congruence b^(c-1) == 1 (mod c) is satisfied, i.e., smallest b larger than 1 such that any member of the set of smallest n composites is a base-b Fermat pseudoprime.
0

%I #14 Jul 28 2019 10:26:26

%S 5,13,25,73,361,361,2521,2521,5041,5041,5041,5041,55441,55441,277201,

%T 3603601,10810801,10810801,10810801,21621601,21621601,367567201,

%U 367567201,367567201

%N a(n) is the smallest b > 1 such that when c is equal to any of the first n composites the congruence b^(c-1) == 1 (mod c) is satisfied, i.e., smallest b larger than 1 such that any member of the set of smallest n composites is a base-b Fermat pseudoprime.

%e For n = 4: The four smallest composites are 4, 6, 8, 9 and for those four values of c the congruence b^(c-1) == 1 (mod c) is satisfied with b = 73. Since 73 is the smallest such value of b > 1, a(4) = 73.

%o (PARI) composites(n) = my(v=[]); forcomposite(c=1, , v=concat(v, [c]); if(#v >= n, return(v)))

%o a(n) = my(cp=composites(n)); for(b=2, oo, for(k=1, #cp, if(Mod(b, cp[k])^(cp[k]-1)!=1, break, if(k==#cp, return(b)))))

%Y Cf. A242742, A256236.

%K nonn,more

%O 1,1

%A _Felix Fröhlich_, Jul 27 2019