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.

A066112 Numbers k such that sigma_4(k)/sigma_2(k) is an integer but not a prime.

This page as a plain text file.
%I A066112 #24 Feb 23 2025 02:40:26
%S A066112 1,16,36,48,49,64,81,100,121,144,162,180,196,225,245,256,324,361,400,
%T A066112 432,441,484,500,529,576,605,625,648,676,729,784,841,900,931,980,1024,
%U A066112 1089,1156,1200,1225,1280,1296,1369,1444,1521,1600,1620,1681,1764,1805
%N A066112 Numbers k such that sigma_4(k)/sigma_2(k) is an integer but not a prime.
%H A066112 Amiram Eldar, <a href="/A066112/b066112.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from Harry J. Smith)
%e A066112 The sequence includes squares, twice squares (such as 162 and 648), and other numbers (such as 48 and 180). The sigma_4/sigma_2 quotients usually have more than one distinct prime factor. Exception: sigma_4(48)/sigma_2(48) = 5732210/3410 = 1681 = 41^2.
%t A066112 Do[s=DivisorSigma[4, n]; z=DivisorSigma[2, n]; If[IntegerQ[s/z]&&!PrimeQ[s/z], Print[n]], {n, 1, 10000}]
%o A066112 (PARI) isok(k) = { my(f=sigma(k, 4)/sigma(k, 2)); !frac(f) && !isprime(f) } \\ _Harry J. Smith_, Feb 01 2010
%Y A066112 Cf. A001157, A001159, A046871, A066109, A066110, A066111.
%K A066112 nonn
%O A066112 1,2
%A A066112 _Labos Elemer_, Dec 06 2001
%E A066112 Edited by _Jon E. Schoenfield_, Dec 24 2016