OFFSET
1,1
COMMENTS
Although there are more terms than A105403 so far, these numbers are still fairly uncommon.
Is this sequence infinite?
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
EXAMPLE
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.
MAPLE
Primes:= select(isprime, [2, seq(2*i+1, i=1..10^6)]):
spf:= map(p -> convert(numtheory:-factorset(p-1), `+`), Primes):
select(t -> spf[t+1]=spf[t], [$1..nops(Primes)-1]);
MATHEMATICA
Select[Range@ 250000, Total[First /@ FactorInteger[Prime@ # - 1]] == Total[First /@ FactorInteger[Prime[# + 1] - 1]] &] (* Michael De Vlieger, Jul 01 2015 *)
PROG
(PARI) spf(n) = {my(f=factor(n)); sum(k=1, #f~, f[k, 1]); }
lista(nn) = {forprime(p=2, nn, if (spf(p-1)==spf(nextprime(p+1)-1), print1(primepi(p), ", ")); ); } \\ Michel Marcus, Jun 30 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from Alois P. Heinz, Jun 30 2015
STATUS
approved