You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Con el lenguaje PARI se buscan estos números de forma similar a la<br />
del anterior caso, sustituyendo m=core(n) por m=n/core(n), Los<br />
tenemos publicados en http://oeis.org/A230356<br />
N PC(N) PL(N)<br />
10 1 10<br />
18 9 2<br />
27 9 3<br />
40 4 10<br />
45 9 5<br />
54 9 6<br />
63 9 7<br />
72 36 2<br />
90 9 10<br />
108 36 3<br />
117 9 13<br />
126 9 14<br />
135 9 15<br />
153 9 17<br />
160 16 10<br />
162 81 2<br />
171 9 19<br />
180 36 5<br />
207 9 23<br />
216 36 6<br />
220 4 55<br />
234 9 26<br />
243 81 3<br />
250 25 10<br />
Último ejemplo<br />
Números compuestos que coinciden en su suma de cifras con su<br />
función SOPF (suma de factores primos tomados sin repetición)<br />
22, 94, 105, 114, 136, 140, 160, 166, 202, 222, 234, 250, 265, 274, 346, 355, 361, 382,<br />
424, 438, 445, 454, 516, 517, 526, 532, 562, 634, 702, 706, 712, 732, 812, 913, 915,<br />
922, 1036, 1071, 1086, 1111…<br />
Lo habrás adivinado ya. La hemos publicado en<br />
http://oeis.org/A230357<br />
Código PARI<br />
sopf(n)= { local(f, s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]); return(s) }<br />
digsum(n)={local (d, p); d=0; p=n; while(p, d+=p%10; p=floor(p/10)); return(d)}<br />
{for (n=4,<br />
2*10^3,m=sopf(n);if(digsum(n)==digsum(m)&&mn,write1("final.txt",n,", ")));}<br />
27
Change language