（PARI）ProdNzD（x）={p=1；while（x>9，d=x-10*（x\10）；if（d，p*=d）；x\=10）；return（p*x）}iscube（x）={f=因子（x）~；for（i=1，length（f），if（t=f[2，i]%3，return 2.txt“，n++，”“，m）；如果（n==1000，中断））}\\哈里·史密斯2009年8月6日
（PARI）ok（k）={my（d=数字（k））；ispower（prod（i=1，#d，if（d[i]，d[i]，1）），3）}\\安德鲁·霍罗伊德，2018年6月29日
