OFFSET
1,1
PROG
(PARI) \prime indices such that gd of prime(x)+ k and prime(x+m) + k are equal divpm1(n, m, k) = { local(x, l1, l2, v1, v2); for(x=2, n, v1 = ifactor(prime(x)+ k); v2 = ifactor(prime(x+m)+k); l1 = length(v1); l2 = length(v2); if(v1[l1] == v2[l2], print1(x", ") ) ) } ifactor(n) = \Vector of the prime factors of n { local(f, j, k, flist); flist=[]; f=Vec(factor(n)); for(j=1, length(f[1]), for(k = 1, f[2][j], flist = concat(flist, f[1][j]) ); ); return(flist) }
CROSSREFS
KEYWORD
nonn
AUTHOR
Cino Hilliard, May 01 2005
EXTENSIONS
8 more terms from R. J. Mathar, Feb 14 2008
a(32)-a(45) from Donovan Johnson, Apr 03 2008
STATUS
approved