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 A280387 #28 Feb 19 2019 11:10:51 %S A280387 4,8,9,16,21,25,27,32,36,45,49,64,81,87,91,99,121,125,128,144,169,196, %T A280387 217,243,256,289,325,343,361,400,417,481,512,529,559,625,685,697,703, %U A280387 721,729,745,749,775,801,841,925,931,961,1024,1156,1157,1261,1331 %N A280387 Composite numbers n such that sum of proper divisors of n divides sum of proper divisors of n^n. %C A280387 Terms are 2^2, 2^3, 3^2, 2^4, 3*7, 5^2, 3^3, 2^5, 2^2*3^2, 3^2*5, 7^2, 2^6, 3^4, 3*29, 7*13, 3^2*11, 11^2, 5^3, ... %C A280387 Terms that are not Duffinian numbers are 45, 87, 91, 99, 196, 703, 745, 775, 801, 931, ... %e A280387 Composite number 21 is a term because (sigma(21) - 21) = 11 divides (sigma(21^21) - 21^21) = 4381940263463668467705506011 %t A280387 Select[Range[10^3], And[CompositeQ@ #, Divisible @@ Map[DivisorSigma[1, #] - # &, {#^#, #}]] &] (* _Michael De Vlieger_, Jan 02 2017 *) %o A280387 (PARI) is(n) = !isprime(n) && (sigma(n^n)-n^n)%(sigma(n)-n)==0; %Y A280387 Cf. A000203, A001065, A003624, A275123. %K A280387 nonn,less %O A280387 1,1 %A A280387 _Altug Alkan_, Jan 01 2017 %E A280387 More terms from _Amiram Eldar_, Feb 19 2019