login
Primes from twin primes by taking the factorial of each digit and adding them up.
0

%I #13 Jul 07 2023 09:54:39

%S 11,80651,67,251,368881,80651,23,1451,251,23,10093,61,362927,363269,

%T 10093,372973,263,10111,59,362941,107,363773,10139,107,364513,10139,

%U 725827,263,365173,364513,12967,737287,10093,10111,374041,776401,162727,725827,725827,368647,10093,263,363269,257,363157,80657,443543,73,444971,727211,735853,735853,19,29,109,363271,20173,17,257,10099

%N Primes from twin primes by taking the factorial of each digit and adding them up.

%H Dario Alpern, <a href="https://www.alpertron.com.ar/ECM.HTM">Factorization using the Elliptic Curve Method</a>

%H Chris Caldwell, <a href="https://t5k.org/lists/small/100ktwins.txt">The First 100,000 Twin Primes </a>

%e 11 from 101 and 103: 1! + 0! + 1! + 1! + 0! + 3! = 11; Similarly 80651 from 281 and 283; 368881 from 569 and 571.

%t Select[Total[Flatten[IntegerDigits/@#]!]&/@Select[Partition[ Prime[ Range[2500]],2,1],Last[#]-First[#]==2&],PrimeQ] (* _Harvey P. Dale_, Jun 21 2011 *)

%Y Cf. A001359

%K base,nonn

%O 1,1

%A _Parthasarathy Nambi_, Oct 01 2009

%E More terms from _Harvey P. Dale_, Jun 21 2011