|
（PARI）正方形（n）={local（s=Set（））；对于（j=1，n，s=setunion（s，Set（j^2%n））；）；返回；}
qnonr（n）={local（s=Set（））；sq=squares（n）；对于（j=0，n-1，if（长度（setintersect（Set（j），sq））==0，s=setunion（s，Set（j））；）；返回（s）；}
a（n）={s=Set（）；qnr=qnonr（n）；对于（j=1，#qnr，if（issquarefree（eval（qnr[j]）），s=setunion\\米歇尔·马库斯2013年7月23日
| 