（PARI）a=[]；对于（S=1199，对于（x=1，S-1，ispseudoprime（p=x^（1+y=S-x）-y^x）&&a=concat（a，p））；vecsort（a）\\以这种方式计算的列表可能不完整，直到最后一个术语。例如，对于较小的S=x+y，在三个较大的项之后，发现x=3，y=97有一个46位数的素数。
（PARI）m=300；a=[]；对于（x=1，m+5，对于（y=1，m+5，p=x^（y+1）-y^x；如果（p<2^m&&ispseudoprime（p），a=concat（a，p））；a=vecsort（a）\\计算低于2^m的所有项。延斯·克鲁斯·安徒生，2014年8月23日
