Anonymous user
Nice primes: Difference between revisions
→{{header|PARI/GP}}
Line 1,160:
=={{header|PARI/GP}}==
<lang PARI/GP>nicePrimes( s, e ) = { local( m );
forprime( p = s, e,
m = p; \\
while( m > 9, \\ m == m mod 9
m = sumdigits( m ) ); \\
if( isprime( m ),
print1( i, " " ) ) );
}</lang>
or
<lang PARI/GP>apply( s -> s, select( s -> isprime( s % 9 ), primes( [500, 1000] )))</lang>
| |||