cp's OEIS Frontend

A225704 Composite squarefree numbers n such that p(i)-4 divides n+4, where p(i) are the prime factors of n.

%I A225704 #14 May 17 2013 14:03:21
%S A225704 6,10,14,15,30,35,66,266,455,806,4154,4686,6665,10370,16646,22781,
%T A225704 31146,36305,72086,205871,246506,473711,570011,653666,733586,900581,
%U A225704 904046,1422410,1941971,1969565,2023010,2807255,2821269,3009821,3043274,3355271,3880301
%N A225704 Composite squarefree numbers n such that p(i)-4 divides n+4, where p(i) are the prime factors of n.
%e A225704 Prime factors of 205871 are 29, 31 and 229. We have that (205871+4)/(29-4) = 8235, (205871+4)/(31-4) = 7625 and (205871+4)/(229-4) = 915.
%p A225704 with(numtheory); A225704:=proc(i,j) local c, d, n, ok, p, t;
%p A225704 for n from 1 to i do if not isprime(n) then p:=ifactors(n)[2]; ok:=1;
%p A225704 for d from 1 to nops(p) do if p[d][2]>1 or p[d][1]=j then ok:=0; break; fi;
%p A225704 if  not type((n+j)/(p[d][1]-j),integer) then ok:=0; break; fi; od;
%p A225704 if ok=1 then print(n); fi; fi; od; end: A225704(10^9,4);
%t A225704 t = {}; n = 0; len = -2; While[len <= 262, n++; {p, e} = Transpose[FactorInteger[n]]; If[Length[p] > 1 && Union[e] == {1} && Union[Mod[n + 4, p - 4]] == {0}, AppendTo[t, n]; len = len + Length[IntegerDigits[n]] + 2]]; t (* _T. D. Noe_, May 17 2013 *)
%Y A225704 Cf. A208728, A225702, A225703, A225705-A225720.
%K A225704 nonn
%O A225704 1,1
%A A225704 _Paolo P. Lava_, May 13 2013
%E A225704 Extended by _T. D. Noe_, May 17 2013