(PARI)
almost_prime_count(N,k)=如果(k==1,返回(primepi(N)));(f(m,p,k,j=0)=my(c=0,s=sqrtnint(N\m,k));如果(k==2,对于素数(q=p,s,c+=primepi(N\(m*q))-j);j+=1),对于素数(q=p,s,c+=f(m*q,q,k-1,j);j+=1));c) ;f(1,2,k);
a(n)=如果(n==0,1,almost_prime_count(8^n,n))\\丹尼尔·苏图2023年7月10日