cp's OEIS Frontend

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.

A180934 Numbers m such that m^k has m divisors for some k >= 1.

This page as a plain text file.
%I A180934 #9 Apr 09 2024 03:26:48
%S A180934 1,2,3,5,7,9,11,13,17,19,23,25,28,29,31,37,40,41,43,45,47,49,53,59,61,
%T A180934 67,71,73,79,81,83,89,97,101,103,107,109,113,121,127,131,137,139,149,
%U A180934 151,153,157,163,167,169,173,179,181,191,193,197,199,211,223,225,227
%N A180934 Numbers m such that m^k has m divisors for some k >= 1.
%C A180934 All primes p are in this sequence, since p^(p-1) has p divisors.
%C A180934 For all odd semiprimes s, s^2 is in this sequence, since s^((s-1)/2) has s divisors.
%H A180934 David W. Wilson, <a href="/A180934/b180934.txt">Table of n, a(n) for n=1..10000</a>
%e A180934 11^10 has 11 divisors, so 11 is in the sequence.
%e A180934 225^7 has 225 divisors, so 225 is in the sequence.
%t A180934 q[n_] := Module[{e = FactorInteger[n][[;; , 2]], k = 1}, While[n > Times @@ (k*e + 1), k++]; n == Times @@ (k*e + 1)]; q[1] = True; Select[Range[250], q] (* _Amiram Eldar_, Apr 09 2024 *)
%Y A180934 A000005(m^k) = m for some k >= 1.
%Y A180934 A180935 gives the corresponding k.
%Y A180934 Cf. A036878, A046315.
%K A180934 nonn
%O A180934 1,2
%A A180934 _David W. Wilson_, Sep 26 2010