(PARI)f(n)=总和(n,d,如果(gcd(d,n/d)==1,d))-n;
isok6(n)=ifrer(f(f(f(f(n))))==n,E,0);
isok3(n)=ifrer(f(f(n)))==n,E,0);
isok2(n)=ifrer(f(f(n))==n,E,0);
isok1(n)=ifrer(f(n)==n,E,0);
isok(n)=isok6(n)&&!isok1(n)&&!isok2(n)&&!isok3(n);
(PARI)
A063919号(n) =我的(f=系数(n));prod(i=1,#f~,f[i,1]^f[i、2]+1)-n