%I A098651
%S A098651 1,2,6,15,30,39,78,105,910,7735,18655,98605
%N A098651 Where A098650 reaches a record.
%t A098651 (* first do *) Needs["NumberTheory`NumberTheoryFunctions`"] (* then *)
PrimeFactors[n_] := Flatten[ Table[ # [[1]], {1} ] & /@ FactorInteger[
n ]]; f[n_] := Block[{k = n + 1}, If[EvenQ[k], k++ ]; While[ PrimeQ[k]
|| Union[ PowerMod[ PrimeFactors[n], k - 1, k]] != {1}, k += 2];
k]; a = {1}; b = {9}; Do[ If[ SquareFreeQ[n], c = f[n]; If[c > b[[
-1]], AppendTo[a, n]; AppendTo[b, c]; Print[n]]], {n, 2, 145000}];
a
%Y A098651 Cf. A098650, A098652.
%Y A098651 Sequence in context: A138621 A163061 A033286 this_sequence A087427 A141126
A056520
%Y A098651 Adjacent sequences: A098648 A098649 A098650 this_sequence A098652 A098653
A098654
%K A098651 nonn
%O A098651 1,2
%A A098651 Robert G. Wilson v (rgwv(AT)rgwv.com), Sep 18 2004 17094
|