OFFSET
1,1
COMMENTS
PROG
(PARI)
{ N = 10^5 + 2;
default(seriesprecision, N);
V = Vec( eta(q) * eta(q^3) * eta(q^5) * eta(q^15) );
forprime(p=2, N, if( V[p]%p == 0, print1(p, ", ") ) );
} \\ Joerg Arndt, Sep 11 2016
(PARI) \\ Much slower than the above, but maybe useful for isolated values
is(n)=if(!isprime(n), return(0)); my(s, t, y='y); for(x=1, n, s+=#polrootsmod(y^2+x*y+y-x^3-x^2+10*x+10, n); if(s>n, return(0))); s==n \\ Charles R Greathouse IV, Sep 12 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Sep 10 2016
EXTENSIONS
More terms from Joerg Arndt, Sep 11 2016
STATUS
approved