OFFSET
1,1
COMMENTS
No terms are divisible by 2 or 3; no terms are in A000961. - Robert Israel, Oct 10 2023
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
with(numtheory):
q:= n-> (f-> ormap(d-> irem(d-1, f)>0, divisors(n)))(min(factorset(n))-1):
select(q, [$2..600])[]; # Alois P. Heinz, Feb 12 2021
MATHEMATICA
Select[Range[2, 600], Function[{d, k}, AnyTrue[d, Mod[#, k] != 0 &]] @@ {Divisors[#] - 1, FactorInteger[#][[1, 1]] - 1} &] (* Michael De Vlieger, Feb 12 2021 *)
PROG
(MATLAB)
n=300; % gives all terms of the sequence not exceeding n
A=[];
for i=2:n
lpf=2;
while mod(i, lpf)~=0
lpf=lpf+1;
end
for d=1:i
if mod(i, d)==0 && mod(d-1, lpf-1)~=0
A=[A i];
break
end
end
end
CROSSREFS
KEYWORD
nonn
AUTHOR
Maxim Karimov, Jan 02 2021
STATUS
approved