RosettaCodeData/Task/Digital-root-Multiplicative-digital-root/PARI-GP/digital-root-multiplicative-digital-root.pari
2016-12-05 22:15:40 +01:00

3 lines
220 B
Text

a(n)=my(i);while(n>9,n=factorback(digits(n));i++);[i,n];
apply(a, [123321, 7739, 893, 899998])
v=vector(10,i,[]); forstep(n=0,oo,1, t=a(n)[2]+1; if(#v[t]<5,v[t]=concat(v[t],n); if(vecmin(apply(length,v))>4, return(v))))