OFFSET
1,1
LINKS
Pierre CAMI, Table of n, a(n) for n = 1..4000
EXAMPLE
Sequence starts with 2,2 as (2*2^1-1)*2*2^1-1=11 is prime and (1*2^1-1)*1*2^1-1 is unity.
(1*2^2-1)*1*2^2-1=11 prime so the second pair is 1,2.
PROG
PFGW & SCRIPTIFY
SCRIPT
DIM n, 0
DIM j
DIM k
DIM ss
DIMS tt
OPENFILEOUT myf, a(n).txt
LABEL a
SET n, n+1
IF n>2000 THEN END
SET ss, 2
LABEL b
SET ss, ss+1
SET j, 1
LABEL c
SET j, j+1
SET k, ss-j
IF k<1 THEN GOTO b
SETS tt, %d, %d, %d\,; n; k; j
PRP (k*j^n-1)*k*j^n-1, tt
IF ISPRP THEN GOTO d
GOTO c
LABEL d
WRITE myf, tt
GOTO a
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Oct 14 2013
STATUS
approved