OFFSET
1,1
COMMENTS
Subsequence of A072443.
LINKS
David A. Corneth, Table of n, a(n) for n = 1..10408 (terms <= 2*10^10)
EXAMPLE
101556 = 156*651 = 273*372; 2096640 = 1092*1920 = 1365*1536.
PROG
(PARI) {for(n=100, 6500000, my(k=logint(n, 100), f=divisors(n), v=[]); for(h=1, matsize(f)[2], if(10^k<f[h]&& f[h]<10^(k+1), v=concat(v, f[h]))); my(b=matsize(v)[2]); if(b>1, my(w=[]); for(i=1, b, my(s=[], a=v[i]); while(a>0, my(d=divrem(a, 10)); a=d[1]; s=concat(d[2], s)); w=concat(w, [vecsort(s)])); my(c=0); for(i=1, b-1, for(j=i+1, b, if(c<2&&w[i]==w[j], if(v[i]*v[j]==n, if(c==1, print1(n, ", "); c=2, c=1)))))))}
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Klaus Brockhaus, Nov 14 2002
EXTENSIONS
Offset changed to 1 by and more terms from David A. Corneth, Sep 08 2024
STATUS
approved