He de que confesar que el descubrimiento que de que entre los números primos gemelos hay cada cierto número de ellos que, para una misma terminación el resto de los dígitos que forman el número primo, son múltiplo de 3, por ejemplo, para el 59/61, tenemos la secuencia,
6 59 6 61
6 6 59 6 6 61
78 6 6 59 78 6 6 61
3 78 6 6 59 3 78 6 6 61
24 3 78 6 6 59 24 3 78 6 6 61
57 24 3 78 6 6 59 57 24 3 78 6 6 61
318 57 24 3 78 6 6 59 318 57 24 3 78 6 6 61
168 318 57 24 3 78 6 6 59 168 318 57 24 3 78 6 6 61
18 168 318 57 24 3 78 6 6 59 18 168 318 57 24 3 78 6 6 61
348 18 168 318 57 24 3 78 6 6 59 348 18 168 318 57 24 3 78 6 6 61
735 348 18 168 318 57 24 3 78 6 6 59 735 348 18 168 318 57 24 3 78 6 6 61
A partir de este hecho hemos construido un algoritmo que teniendo en cuenta esto nos diera números primos de la misma terminación.
Un algoritmo podría ser (en lenguaje SageMath)
var ('b,a,p,w')
a=1;b=a;
p=59;#puede ser cualquier primo gemelo
while a<99999999:
w = int(str(a) + str(p))
a=a+1;
if is_prime (w) and is_prime (w + 2) :
print(a-1);
print( w,' - ', len(str(w)));
p=w;a=b;
end