This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A067671 #7 Mar 29 2014 13:58:31 %S A067671 4,16,21,27,36,64,256,288,648,729,1024,1444,1458,1764,1936,2304,3125, %T A067671 4096,4361,5184,6272,7688,8277,9408,11664,16384,18432,19683,22472, %U A067671 22987,26244,28125,29403,31199,41472,43264,59577,65536,67712,73008 %N A067671 The prime factors of n are also prime factors of the decimal encoding (A067599) of the prime factorization of n. %e A067671 21 = 3^1 * 7^1 has prime factors 3,7, which are also prime factors of the corresponding decimal encoding 3171 = 3^1 * 7^1 * 151^1. %t A067671 (*f gives the decimal encoding of the prime factorization of n*) f[n_] := FromDigits[Flatten[IntegerDigits[FactorInteger[n]]]]; (*g gives the list of prime factors of n*) g[n_] := Module[{a, l, t}, a = FactorInteger[n]; l = Length[a]; Table[a[[i]][[1]], {i, 1, l}]]; %t A067671 (*main routine*) j[n] := Module[{l1 = g[n], l2 = g[f[n]]}, (Intersection[l1, l2] == l1)]; Select[Range[2, 10^5], j[ # ] &] %Y A067671 Cf. A067599. %K A067671 nonn,base %O A067671 1,1 %A A067671 _Joseph L. Pe_, Feb 04 2002