（PARI）piptopipm1（n）={本地（x，y）；对于（x=1，n，y=pip（x）^pip（x）-1；打印1（bdiv（y）“，”）；）}
pip（n）={return（质数（质数））}
bdiv（n）={本地（x）；x=ifactor（n）；返回（x[长度（x）]）}
ifactor（n，m=0）={局部（f，j，k，flist）；flist=[]；f=Vec（因子（n，m））；对于（j=1，长度（f[1]），对于（k=1，f[2][j]，flist=concat（flist，f[1][j]））；）；return（flist）}
