OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..125
EXAMPLE
a(3) = 460317 because the prime divisors of 460317 are 3, 11, 13, 29, 37 =>
(3 + 3) | (460317 + 3) = 460320 = 6*76720;
(11 + 3) | 460320 = 14*32880;
(13 + 3) | 460320 = 16*28770;
(29+3) | 460320 = 32*14385;
(37+3) | 460320 = 40*11508.
MAPLE
with(numtheory):for n from 1 to 23 do:i:=0:for k from 1 to 10^8 while(i=0) do:x:=factorset(k):n1:=nops(x):y:=k+n: j:=0:for m from 1 to n1 do:if n1>=2 and irem(y, x[m]+n)=0 then j:=j+1:else fi:od:if j>4 then i:=1: printf ( "%d %d \n", n, k):else fi:od:od:
MATHEMATICA
numd[n_, k_] := Module[{p=FactorInteger[k][[;; , 1]], c=0}, Do[If[Divisible[n+k, n+p[[i]]], c++], {i, 1, Length[p]}]; c]; a[n_]:=Module[{k=1}, While[numd[n, k] <= 4, k++]; k]; Array[a, 30] (* Amiram Eldar, Sep 09 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Dec 13 2011
STATUS
approved