A028893 Duplicate of A029561.
1595, 6785, 53867, 67727, 102377, 296003, 740027, 961877, 998867, 1048817
Offset: 1
This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
a(1) = 35 because this is the first squarefree composite number n such that at least one integer b except 0 exists such that for every prime factor p of n applies that p+b divides n+b (-3): 35 = 5*7 and 2, 4 both divide 32.
fQ[n_] := Block[{c = -1, fi = FactorInteger@ n, k, lmt, p}, If[Times @@ (Last@# & /@ fi) == 1 < Plus @@ (Last@# & /@ fi), p = First@# & /@ fi; k = -fi[[1, 1]] + 1; lmt = Abs[(n - fi[[-1, 1]]^2)/fi[[-1, 1]]]; While[k < lmt, If[ Union[ IntegerQ@# & /@ ((n + k)/(p + k))] == {True}, c++; If[c > 0, Goto [fini]]]; k++]]; Label[fini]; c > 0]; Select[ Range@ 2000, fQ] (* Robert G. Wilson v, Dec 05 2015 *)
for(n=2,1000000, if(!isprime(n), if(issquarefree(n), f=factor(n); k=0; for(b=-(f[1, 1]-1),n, c=0; for(i=1, #f[, 1], if((n+b)%(f[i, 1]+b)>0, c++)); if(c==0, if(!b==0, k++))); if(k>0, print1(n,", ")))))
Comments