OFFSET
1,1
PROG
(PARI) for(n=0, 10^6, d=digits(n); p=prod(i=1, #d, d[i]); vp=vecsort(digits(p+n), , 8); vs=vecsort(digits(sumdigits(n)+n), , 8); if(vs==vp&&vp==vecsort(d, , 8)&&vs==vecsort(d, , 8)&&p, print1(n, ", ")))
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Derek Orr, Oct 12 2014
STATUS
approved