A371730 - OEIS

%I #10 Apr 10 2024 09:42:07

%S 2,5,9,13,17,25,37,41,49,73,97,109,161,181,289,307,361,433,529,701,

%T 841,881,961,1189,1369,1681,1849,2209,2393,2449,3277,3457,3781,4096,

%U 4159,4841,4861,5741,5851,6049,7921,10201,11449,11881,15049,15401,19601,27026,32768

%N Indices of records in A371729: numbers k with A371729(k) > A371729(m) for all m < k.

%C The corresponding record values are 0, 1, 2, 3, 4, 5, 7, 8, 9, 11, 12, ... (see the link for more values).	

%H Amiram Eldar, <a href="/A371730/b371730.txt">Table of n, a(n) for n = 1..69</a> (terms below 10^6)

%H Amiram Eldar, <a href="/A371730/a371730.txt">Table of n, a(n), A371729(a(n)) for n = 1..69</a>.

%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Pseudoprime">Pseudoprime</a>.

%H <a href="/index/Ps#pseudoprimes">Index entries for sequences related to pseudoprimes</a>.

%t f[k_] := Count[Range[4, k - 1], _?(CompositeQ[#] && PowerMod[k, # - 1, #] == 1 &)]; seq[kmax_] := Module[{s = {}, fm = -1, f1}, Do[f1 = f[k]; If[f1 > fm, fm = f1; AppendTo[s, k]], {k, 2, kmax}]; s]; seq[2500]

%o (PARI) f(n) = {my(c=0); forcomposite(k = 4, n-1, if(Mod(n, k)^(k-1) == 1, c++)); c;}

%o lista(kmax) = {my(fm = -1, f1); for(k = 2, kmax, f1 = f(k); if(f1 > fm, fm = f1; print1(k, ", ")));}

%Y Cf. A371729.

%K nonn

%O 1,1

%A _Amiram Eldar_, Apr 05 2024