2番目以降の四つ子素数

計算時間を贅沢に

> library(gmp)
> n = 30*0:10000
> n = n[isprime(n+11) > 0 & isprime(n+13) > 0 & isprime(n+17) > 0 & isprime(n+19) > 0]
> t(sapply(n, "+", c(11, 13, 17, 19)))
        [,1]   [,2]   [,3]   [,4]
 [1,]     11     13     17     19
 [2,]    101    103    107    109
 [3,]    191    193    197    199
 [4,]    821    823    827    829
 [5,]   1481   1483   1487   1489
 [6,]   1871   1873   1877   1879
 [7,]   2081   2083   2087   2089
 [8,]   3251   3253   3257   3259
 [9,]   3461   3463   3467   3469
[10,]   5651   5653   5657   5659
：