N-smooth numbers: Difference between revisions

m
Line 2,366:
S7(1691)
S7(10^6)
my(p=509,n=30000,c=20); [p, n, c, upto(n+c,primes([0,p]))[n..n+c-1] ]
</syntaxhighlight>
{{out}}
<pre>
<pre>%1 = [1, 2, 3, 4, 5, 6, 8, 9, 10, 12, 15, 16, 18, 20, 24, 25, 27, 30, 32, 36]
%2 = 3317760
%3 = 4157409948433216829957008507500000000
%4 = [509, 30000, 20, [62601, 62602, 62604, 62607, 62608, 62609, 62611, 62618, 62620, 62622, 62624, 62625, 62626, 62628, 62629, 62634, 62640, 62643, 62645, 62646]]
</pre>
 
15

edits