cp's OEIS Frontend

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.

A259562 Numbers n such that the sum of the distinct prime factors of prime(n)-1 and prime(n+1)-1 are the same.

This page as a plain text file.
%I A259562 #20 Jul 01 2015 19:24:02
%S A259562 2,414,556,3962,4972,6151,6521,8440,8665,13769,13909,15576,16696,
%T A259562 17176,19926,20630,21541,27090,30822,62118,65349,74014,94203,98600,
%U A259562 101231,103058,108333,112332,136036,142714,145588,147150,160730,162366,169137,194681,200837
%N A259562 Numbers n such that the sum of the distinct prime factors of prime(n)-1 and prime(n+1)-1 are the same.
%C A259562 Although there are more terms than A105403 so far, these numbers are still fairly uncommon.
%C A259562 Is this sequence infinite?
%C A259562 It would follow from the generalized Bunyakovsky conjecture that, e.g., there are infinitely many primes p such that p+2, p+12, p+14, 6*p^2+84*p+1 and 6*p^2+84*p+145 are all prime, and there are no primes between 6*p^2+84*p+1 and 6*p^2+84*p+145. If so, then the sequence is infinite, because it contains n where prime(n) = 6*p^2+84*p+1, with prime(n)-1 having distinct prime factors 2,3,p,p+14 and prime(n+1) having distinct prime factors 2,3,p+2,p+12. - _Robert Israel_, Jun 30 2015
%e A259562 The prime factors of prime(414)-1 are 2,3,5,5,19 and the prime factors of prime(415)-1 are 2,2,2,3,7,17. The sum of the distinct entries in each of these lists is 29.
%p A259562 Primes:= select(isprime,[2,seq(2*i+1,i=1..10^6)]):
%p A259562 spf:= map(p -> convert(numtheory:-factorset(p-1),`+`), Primes):
%p A259562 select(t -> spf[t+1]=spf[t], [$1..nops(Primes)-1]);
%t A259562 Select[Range@ 250000, Total[First /@ FactorInteger[Prime@ # - 1]] == Total[First /@ FactorInteger[Prime[# + 1] - 1]] &] (* _Michael De Vlieger_, Jul 01 2015 *)
%o A259562 (PARI) spf(n) = {my(f=factor(n)); sum(k=1, #f~, f[k,1]);}
%o A259562 lista(nn) = {forprime(p=2, nn, if (spf(p-1)==spf(nextprime(p+1)-1), print1(primepi(p), ", ")););} \\ _Michel Marcus_, Jun 30 2015
%Y A259562 Cf. A008472, A105403, A259564.
%K A259562 nonn
%O A259562 1,1
%A A259562 _Otis Tweneboah_, _Pratik Koirala_, _Eugene Fiorini_, _Nathan Fox_, Jun 30 2015
%E A259562 More terms from _Alois P. Heinz_, Jun 30 2015