|
PROG
|
(PARI) a087436(n) = bigomega (n >> valuation (n, 2));
a359640(maxp) = {my(k=2, pp=5); forprime (p=7, maxp, my(mi=oo, ma=0); if (p-pp>2, for (j=pp+1, p-1, my(mo=a087436(j)); if (mo<k, mi=ma=0; break); mi=min(mo, mi); ma=max(mo, ma)); if (mi==k && ma==k, print1(pp, ", "); k++)); pp=p)};
a359640(3*10^8)
|