defnum_factors(n,p,rpt):results=[]foriinrange(1,int(n**0.5)+1):ifn%i==0:ifn==p:rpt.append(i)rpt.append(int(n/i))results.append(i)results.append(int(n/i))ifn==p:rpt=sorted(set(rpt))rpt=rpt[1:]results=sorted(set(results))results=results[1:]k=set(results).intersection(set(rpt))# print kiflen(k)==0:numbers.append(n)returnnumbers#################################################################p=10000r_p_to=[]# relatively prime tonumbers=[]forjinrange(p,0,-1):l=num_factors(j,p,r_p_to)printsorted(l)print'Number of positive integers relatively prime to %d: %d'%(p,len(l))print'Sum of positive integers relatively prime to %d: %d'%(p,sum(l))

$</code> ... <code>$</code>...<code>."> Easy Math Editor

`*italics*`

or`_italics_`

italics`**bold**`

or`__bold__`

boldNote: you must add a full line of space before and after lists for them to show up correctlyparagraph 1

paragraph 2

`[example link](https://brilliant.org)`

`> This is a quote`

Remember to wrap math in $</span> ... <span>$ or $</span> ... <span>$ to ensure proper formatting.`2 \times 3`

`2^{34}`

`a_{i-1}`

`\frac{2}{3}`

`\sqrt{2}`

`\sum_{i=1}^3`

`\sin \theta`

`\boxed{123}`

## Comments

Sort by:

TopNewestAre you testing features out or something?

Log in to reply