A124010 -id:A124010 - cp's OEIS frontend

A046897 Sum of divisors of n that are not divisible by 4.

1, 3, 4, 3, 6, 12, 8, 3, 13, 18, 12, 12, 14, 24, 24, 3, 18, 39, 20, 18, 32, 36, 24, 12, 31, 42, 40, 24, 30, 72, 32, 3, 48, 54, 48, 39, 38, 60, 56, 18, 42, 96, 44, 36, 78, 72, 48, 12, 57, 93, 72, 42, 54, 120, 72, 24, 80, 90, 60, 72, 62, 96, 104, 3, 84, 144, 68, 54, 96, 144, 72

Offset: 1

N. J. A. Sloane

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

The o.g.f. is (theta_3(0,x)^4 - 1)/8, see the Hardy reference, eqs. 9.2.1, 9.2.3 and 9.2.4 on p. 133 for Sum' m*u_m. Also Hardy-Wright, p. 314. See also the Somos, Jan 25 2008 formula below. - Wolfdieter Lang, Dec 11 2016

			G.f. = q + 3*q^2 + 4*q^3 + 3*q^4 + 6*q^5 + 12*q^6 + 8*q^7 + 3*q^8 + 13*q^9 + ...

J. M. Borwein, D. H. Bailey and R. Girgensohn, Experimentation in Mathematics, A K Peters, Ltd., Natick, MA, 2004. x+357 pp. See p. 194.
G. H. Hardy, Ramanujan: twelve lectures on subjects suggested by his life and work, AMS Chelsea Publishing, Providence, Rhode Island 2002, p. 133.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Clarendon Press, Oxford, Fifth edition, 1979, p. 314.
P. A. MacMahon, Combinatory Analysis, Cambridge Univ. Press, London and New York, Vol. 1, 1915 and Vol. 2, 1916; see vol. 2, p 31, Article 273.
C. J. Moreno and S. S. Wagstaff, Jr., Sums of Squares of Integers, Chapman & Hall, 2006.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Michael Somos, Introduction to Ramanujan theta functions.
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions.

Cf. A000203, A000118, A051731, A069733, A027748, A124010, A190621, A000593 (not divis. by 2), A046913 (not divis. by 3), A116073 (not divis. by 5).

a046897 1 = 1
a046897 n = product $ zipWith
            (\p e -> if p == 2 then 3 else div (p ^ (e + 1) - 1) (p - 1))
            (a027748_row n) (a124010_row n)
-- Reinhard Zumkeller, Aug 12 2015

A := Basis( ModularForms( Gamma0(4), 2), 72); B := (A[1] - 1)/8 + A[2]; B; /* Michael Somos, Dec 30 2014 */

A046897 := proc(n) if n mod 4 = 0 then numtheory[sigma](n)-4*numtheory[sigma](n/4) ; else numtheory[sigma](n) ; end if; end proc: # R. J. Mathar, Mar 23 2011

a[n_] := Sum[ Boole[ !Divisible[d, 4]]*d, {d, Divisors[n]}]; Table[ a[n], {n, 1, 71}] (* Jean-François Alcover, Dec 12 2011 *)
DivisorSum[#1, # &, Mod[#, 4] != 0 &] & /@ Range[71] (* Jayanta Basu, Jun 30 2013 *)
a[ n_] := SeriesCoefficient[ (EllipticTheta[ 3, 0, q]^4 - 1) / 8, {q, 0, n}]; (* Michael Somos, Dec 30 2014 *)
f[2, e_] := 3; f[p_, e_] := (p^(e+1)-1)/(p-1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 15 2020 *)

{a(n) = if( n<1, 0, sumdiv(n, d, if(d%4, d)))};

a(n) = (-1)^(n+1)*Sum_{d divides n} (-1)^(n/d+d)*d. Multiplicative with a(2^e) = 3, a(p^e) = (p^(e+1)-1)/(p-1) for an odd prime p. - Vladeta Jovovic, Sep 10 2002 [For a proof of the multiplicative property, see for example Moreno and Wagstaff, p. 33. - N. J. A. Sloane, Nov 09 2016]
G.f.: Sum_{k>0} x^k/(1+(-x)^k)^2, or Sum_{k>0} k*x^k/(1+(-x)^k). - Vladeta Jovovic, Dec 16 2002
Expansion of (1 - phi(q)^4) / 8 in powers of q where phi() is a Ramanujan theta function. - Michael Somos, Jan 25 2008
Equals inverse Mobius transform of A190621. - Gary W. Adamson, Jul 03 2008
A000118(n) = 8*a(n) for all n>0.
Dirichlet g.f.: (1 - 4^(1-s)) * zeta(s) * zeta(s-1). - Michael Somos, Oct 21 2015
L.g.f.: log(Product_{k>=1} (1 - x^(4*k))/(1 - x^k)) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, Mar 14 2018
From Peter Bala, Dec 19 2021: (Start)
Logarithmic g.f.: Sum_{n >= 1} a(n)*x^n/n = Sum_{n >= 1} x^n*(1 + x^n + x^(2*n))/( n*(1 - x^(4*n)) )
G.f.: Sum_{n >= 1} x^n*(x^(6*n) + 2*x^(5*n) + 3*x^(4*n) + 3*x^(2*n) + 2*x^n + 1)/(1 - x^(4*n))^2. (End)
Sum_{k=1..n} a(k) ~ (Pi^2/16) * n^2. - Amiram Eldar, Oct 04 2022

A209061 Exponentially squarefree numbers.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69

Offset: 1

Reinhard Zumkeller, Mar 13 2012

nonn

Numbers having only squarefree exponents in their canonical prime factorization.
According to the formula of Theorem 3 [Toth], the density of the exponentially squarefree numbers is 0.9559230158619... (A262276). - Peter J. C. Moses and Vladimir Shevelev, Sep 10 2015
From Vladimir Shevelev, Sep 24 2015: (Start)
A generalization. Let S be a finite or infinite increasing integer sequence s=s(n), s(0)=0.
Let us call a positive number N an exponentially S-number, if all exponents in its prime power factorization are in the sequence S.
Let {u(n)} be the characteristic function of S. Then, for the density h=h(S) of the exponentially S-numbers, we have the representations
h(S) = Product_{prime p} Sum_{j in S} (p-1)/p^(j+1) = Product_{p} (1 + Sum_{j>=1} (u(j) - u(j-1))/p^j). In particular, if S = {0,1}, then the exponentially S-numbers are squarefree numbers; if S consists of 0 and {2^k}A138302%20(see%20%5BShevelev%5D,%202007);%20if%20S%20consists%20of%200%20and%20squarefree%20numbers,%20then%20u(n)=%7Cmu(n)%7C,%20where%20mu(n)%20is%20the%20M%C3%B6bius%20function%20(A008683),%20we%20obtain%20the%20density%20h%20of%20the%20exponentially%20squarefree%20numbers%20(cf.%20Toth's%20link,%20Theorem%203);%20the%20calculation%20of%20h%20with%20a%20very%20high%20degree%20of%20accuracy%20belongs%20to%20_Juan%20Arias-de-Reyna">{k>=0}, then the exponentially S-numbers form A138302 (see [Shevelev], 2007); if S consists of 0 and squarefree numbers, then u(n)=|mu(n)|, where mu(n) is the Möbius function (A008683), we obtain the density h of the exponentially squarefree numbers (cf. Toth's link, Theorem 3); the calculation of h with a very high degree of accuracy belongs to _Juan Arias-de-Reyna (A262276). Note that if S contains 1, then h(S) >= 1/zeta(2) = 6/Pi^2; otherwise h(S) = 0. Indeed, in the latter case, the density of the sequence of exponentially S-numbers does not exceed the density of A001694, which equals 0. (End)
The term "exponentially squarefree number" was apparently coined by Subbarao (1972). - Amiram Eldar, May 28 2025

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Xiaodong Cao and Wenguang Zhai, Some arithmetic functions involving exponential divisors, JIS 13 (2010), Article 10.3.7.
Y.-F. S. Petermann, Arithmetical functions involving exponential divisors: note on two papers by L. Toth, Ann. Univ. Sci. Budapest, Sect. Comp. 32 (2010) 143-149.
Vladimir Shevelev, Compact integers and factorials, Acta Arithmetica 126:3 (2007), pp. 195-236.
Vladimir Shevelev, Exponentially S-numbers, arXiv:1510.05914 [math.NT], 2015.
Vladimir Shevelev, Set of all densities of exponentially S-numbers, arXiv:1511.03860 [math.NT], 2015.
Vladimir Shevelev, A fast computation of density of exponentially S-numbers, arXiv:1602.04244 [math.NT], 2016.
Vladimir Shevelev, S-exponential numbers, Acta Arithmetica, Vol. 175(2016), 385-395.
H. M. Stark, On the asymptotic density of the k-free integers, Proc. Amer. Soc. 17 (1966), 1211-1214.
M. V. Subbarao, On some arithmetic convolutions, in: A. A. Gioia and D. L. Goldsmith (eds.), The Theory of Arithmetic Functions: Proceedings of the Conference at Western Michigan University, April 29 - May 1, 1971, Berlin, Heidelberg: Springer Berlin Heidelberg, 2006, pp. 247-271; alternative link.
László Tóth, On certain arithmetic functions involving exponential divisors, II., Annales Univ. Sci. Budapest., Sect. Comp., 27 (2007), 155-166 and arXiv:0708.3557 [math.NT], 2007-2009.

Complement of A130897.
A005117, A004709, and A046100 are subsequences.
Cf. A001221, A001694, A008683, A008966, A036537, A115063, A138302, A197680, A262276, A262675, A268335, A270428.

a209061 n = a209061_list !! (n-1)
a209061_list = filter
   (all (== 1) . map (a008966 . fromIntegral) . a124010_row) [1..]

Select[Range@ 69, Times @@ Boole@ Map[SquareFreeQ, Last /@ FactorInteger@ #] > 0 &] (* Michael De Vlieger, Sep 07 2015 *)

is(n)=my(f=factor(n)[,2]); for(i=1,#f,if(!issquarefree(f[i]), return(0))); 1 \\ Charles R Greathouse IV, Sep 02 2015

A166234(a(n)) <> 0.
Product_{k=1..A001221(n)} A008966(A124010(n,k)) = 1.
One can prove that the principal term of Toth's asymptotics for the density of this sequence (cf. Toth's link, Theorem 3) equals also Product_{prime p}(Sum_{j in S}(p-1)/p^{j+1})*x, where S is the set of 0 and squarefree numbers. The remainder term O(x^(0.2+t)), where t>0 is arbitrarily small, was obtained by L. Toth while assuming the Riemann Hypothesis. - Vladimir Shevelev, Sep 12 2015

A355536 Irregular triangle read by rows where row n lists the differences between adjacent prime indices of n; if n is prime, row n is empty.

Original entry on oeis.org

0, 1, 0, 0, 0, 2, 0, 1, 3, 1, 0, 0, 0, 1, 0, 0, 2, 2, 4, 0, 0, 1, 0, 5, 0, 0, 0, 3, 1, 1, 0, 0, 0, 0, 3, 6, 1, 0, 1, 0, 7, 4, 0, 0, 2, 1, 2, 0, 4, 0, 1, 8, 0, 0, 0, 1, 0, 2, 0, 5, 0, 5, 1, 0, 0, 2, 0, 0, 3, 6, 9, 0, 1, 1, 10, 0, 2, 0, 0, 0, 0, 0, 3, 1, 3, 0, 6

Offset: 2

Gus Wiseman, Jul 12 2022

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
The version where zero is prepended to the prime indices is A287352.
One could argue that row n = 1 is empty, but adding it changes only the offset, not the data.

			Triangle begins (showing n, prime indices, differences*):
   2:    (1)       .
   3:    (2)       .
   4:   (1,1)      0
   5:    (3)       .
   6:   (1,2)      1
   7:    (4)       .
   8:  (1,1,1)    0 0
   9:   (2,2)      0
  10:   (1,3)      2
  11:    (5)       .
  12:  (1,1,2)    0 1
  13:    (6)       .
  14:   (1,4)      3
  15:   (2,3)      1
  16: (1,1,1,1)  0 0 0

Gus Wiseman, Sequences counting and ranking integer partitions by the differences of their successive parts.

Row-lengths are A001222 minus one.
The prime indices are A112798, sum A056239.
Row-sums are A243055.
Constant rows have indices A325328.
The Heinz numbers of the rows plus one are A325352.
Strict rows have indices A325368.
Row minima are A355524.
Row maxima are A286470, also A355526.
An adjusted version is A358169, reverse A355534.
Cf. A066312, A124010, A129654, A243056, A287352, A325394, A355523, A355528, A355531.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Differences[primeMS[n]],{n,2,100}]

A332282 Numbers whose unsorted prime signature is not unimodal.

Original entry on oeis.org

300, 588, 600, 980, 1176, 1200, 1452, 1500, 1960, 2028, 2100, 2205, 2352, 2400, 2420, 2904, 2940, 3000, 3300, 3380, 3388, 3468, 3900, 3920, 4056, 4116, 4200, 4332, 4410, 4704, 4732, 4800, 4840, 5100, 5445, 5700, 5780, 5808, 5880, 6000, 6348, 6468, 6600, 6615

Offset: 1

Gus Wiseman, Feb 19 2020

nonn

The unsorted prime signature of a positive integer (row n of A124010) is the sequence of exponents it is prime factorization.
A sequence of positive integers is unimodal if it is the concatenation of a weakly increasing and a weakly decreasing sequence.
Also Heinz numbers of integer partitions with non-unimodal run-lengths. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

			The sequence of terms together with their prime indices begins:
   300: {1,1,2,3,3}
   588: {1,1,2,4,4}
   600: {1,1,1,2,3,3}
   980: {1,1,3,4,4}
  1176: {1,1,1,2,4,4}
  1200: {1,1,1,1,2,3,3}
  1452: {1,1,2,5,5}
  1500: {1,1,2,3,3,3}
  1960: {1,1,1,3,4,4}
  2028: {1,1,2,6,6}
  2100: {1,1,2,3,3,4}
  2205: {2,2,3,4,4}
  2352: {1,1,1,1,2,4,4}
  2400: {1,1,1,1,1,2,3,3}
  2420: {1,1,3,5,5}
  2904: {1,1,1,2,5,5}
  2940: {1,1,2,3,4,4}
  3000: {1,1,1,2,3,3,3}
  3300: {1,1,2,3,3,5}
  3380: {1,1,3,6,6}

MathWorld, Unimodal Sequence

The opposite version is A332642.
These are the Heinz numbers of the partitions counted by A332281.
Non-unimodal permutations are A059204.
Non-unimodal compositions are A115981.
Non-unimodal normal sequences are A328509.
Cf. A001523, A007052, A056239, A112798, A124010, A227038, A332280, A332284, A332286, A332639, A332643, A332671.

unimodQ[q_]:=Or[Length[q]==1,If[q[[1]]<=q[[2]],unimodQ[Rest[q]],OrderedQ[Reverse[q]]]];
Select[Range[1000],!unimodQ[Last/@FactorInteger[#]]&]

A359178 Numbers with a unique smallest prime exponent.

Original entry on oeis.org

2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 16, 17, 18, 19, 20, 23, 24, 25, 27, 28, 29, 31, 32, 37, 40, 41, 43, 44, 45, 47, 48, 49, 50, 52, 53, 54, 56, 59, 61, 63, 64, 67, 68, 71, 72, 73, 75, 76, 79, 80, 81, 83, 88, 89, 92, 96, 97, 98, 99, 101, 103, 104, 107, 108, 109, 112, 113, 116, 117

Offset: 1

Jens Ahlström, Jan 08 2023

nonn

180 is the smallest number with a unique smallest prime exponent that is not a member of A130091.

			2 = 2^1 is a term since it has 1 as a unique smallest exponent.
6 = 2^1 * 3^1 is not a term since it has two primes with the same smallest exponent.
180 = 2^2 * 3^2 * 5^1 is a term since it has 1 as a unique smallest exponent.

Jens Ahlström, Table of n, a(n) for n = 1..9999
Wikipedia, Mode (statistics).

For parts instead of multiplicities we have A247180, counted by A002865.
For greatest instead of smallest we have A356862, counted by A362608.
The complement is A362606, counted by A362609.
Partitions of this type are counted by A362610.
These are the positions of 1's in A362613, for modes A362611.
A001221 counts prime exponents and A001222 adds them up.
A027746 lists prime factors, A112798 indices, A124010 exponents.
Cf. A102750, A130091, A327473, A327476, A359908, A362605, A362615.

q[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, Count[e, Min[e]] == 1]; Select[Range[2, 200], q] (* Amiram Eldar, Jan 08 2023 *)

isok(n) = if (n>1, my(f=factor(n), e = vecmin(f[,2])); #select(x->(x==e), f[,2], 1) == 1); \\ Michel Marcus, Jan 27 2023

from sympy import factorint
def ok(k):
  c = sorted(factorint(k).values())
  return len(c) == 1 or c[0] != c[1]
print([k for k in range(2, 118) if ok(k)])

from itertools import count, islice
from sympy import factorint
def A359178_gen(startvalue=2): # generator of terms >= startvalue
    return filter(lambda n:(f:=list(factorint(n).values())).count(min(f))==1,count(max(startvalue,2)))
A359178_list = list(islice(A359178_gen(),20)) # Chai Wah Wu, Feb 08 2023

A055229 Greatest common divisor of largest square dividing n and squarefree part of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 3, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1

Offset: 1

Labos Elemer, Jun 21 2000

Record values occur at cubes of squarefree numbers: a(A062838(n)) = A005117(n) and a(m) < A005117(n) for m < A062838(n). - Reinhard Zumkeller, Apr 09 2010

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A000188, A007913, A007947, A008833.

Cf. A027748, A005117, A062838, A060476, A124010, A220218.

a055229 n = product $ zipWith (^) ps (map (flip mod 2) es) where
   (ps, es) = unzip $
              filter ((> 1) . snd) $ zip (a027748_row n) (a124010_row n)
-- Reinhard Zumkeller, Oct 27 2015

a[n_] := With[{sf = Times @@ Power @@@ ({#[[1]], Mod[#[[2]], 2]}& /@ FactorInteger[n])}, GCD[sf, n/sf]]; Table[a[n], {n, 1, 105}] (* Jean-François Alcover, Feb 05 2014 *)

a(n)=my(c=core(n));gcd(c,n/c) \\ Charles R Greathouse IV, Nov 20 2012

a(n) = gcd[A008833(n), A007913(n)].
Multiplicative with a(p^e)=1 for even e, a(p)=1, a(p^e)=p for odd e>1. - Vladeta Jovovic, Apr 30 2002
A220218(a(n)) = 1; A060476(a(n)) > 1 for n > 1. - Reinhard Zumkeller, Nov 30 2015
a(n) = core(n)*rad(n/core(n))/rad(n), where core = A007913 and rad = A007947. - Conjecture by Velin Yanev, proof by David J. Seal, Sep 19 2017
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} ((p^3 + p^2 + p - 1)/(p^2 * (p + 1))) = 1.2249749939341923764... . - Amiram Eldar, Oct 08 2022

A076467 Perfect powers m^k where m is a positive integer and k >= 3.

Original entry on oeis.org

1, 8, 16, 27, 32, 64, 81, 125, 128, 216, 243, 256, 343, 512, 625, 729, 1000, 1024, 1296, 1331, 1728, 2048, 2187, 2197, 2401, 2744, 3125, 3375, 4096, 4913, 5832, 6561, 6859, 7776, 8000, 8192, 9261, 10000, 10648, 12167, 13824, 14641, 15625, 16384, 16807

Offset: 1

Robert G. Wilson v, Oct 14 2002

nonn

If p|n with p prime then p^3|n.

Alois P. Heinz, Table of n, a(n) for n = 1..10000 (terms n = 1..250 from Reinhard Zumkeller)

Subsequence of A036966.

Cf. A001597, A076468, A076469, A076470.

a076467 n = a076467_list !! (n-1)
a076467_list = 1 : filter ((> 2) . foldl1 gcd . a124010_row) [2..]
-- Reinhard Zumkeller, Apr 13 2012

import qualified Data.Set as Set (null)
import Data.Set (empty, insert, deleteFindMin)
a076467 n = a076467_list !! (n-1)
a076467_list = 1 : f [2..] empty where
   f xs'@(x:xs) s | Set.null s || m > x ^ 3 = f xs $ insert (x ^ 3, x) s
                  | m == x ^ 3  = f xs s
                  | otherwise = m : f xs' (insert (m * b, b) s')
                  where ((m, b), s') = deleteFindMin s
-- Reinhard Zumkeller, Jun 18 2013

N:= 10^5: # to get all terms <= N
S:= {1, seq(seq(m^k, m = 2 .. floor(N^(1/k))),k=3..ilog2(N))}:
sort(convert(S,list)); # Robert Israel, Sep 30 2015

a = {1}; Do[ If[ Apply[ GCD, Last[ Transpose[ FactorInteger[n]]]] > 2, a = Append[a, n]; Print[n]], {n, 2, 17575}]; a
(* Second program: *)
n = 10^5; Join[{1}, Table[m^k, {k, 3, Floor[Log[2, n]]}, {m, 2, Floor[n^(1/k)]}] // Flatten // Union] (* Jean-François Alcover, Feb 13 2018, after Robert Israel *)

is(n)=ispower(n)>2||n==1 \\ Charles R Greathouse IV, Sep 03 2015, edited for n=1 by M. F. Hasler, May 26 2018

A076467(lim)={my(L=List(1),lim2=logint(lim,2),m,k);for(k=3,lim2, for(m=2,sqrtnint(lim,k),listput(L, m^k);));listsort(L,1);L}
b076467(lim)={my(L=A076467(lim)); for(i=1,#L,print(i ," ",L[i]));} \\ Anatoly E. Voevudko, Sep 29 2015, edited by M. F. Hasler, May 25 2018

A076467_vec(LIM,S=List(1))={for(x=2,sqrtnint(LIM,3),for(k=3, logint(LIM, x), listput(S, x^k))); Set(S)} \\ M. F. Hasler, May 25 2018

from sympy import mobius, integer_nthroot
def A076467(n):
    def f(x): return int(n-1+x-integer_nthroot(x,4)[0]+sum(mobius(k)*(integer_nthroot(x,k)[0]+integer_nthroot(x,k<<1)[0]-2) for k in range(3,x.bit_length())))
    kmin, kmax = 1,2
    while f(kmax) >= kmax:
        kmax <<= 1
    while True:
        kmid = kmax+kmin>>1
        if f(kmid) < kmid:
            kmax = kmid
        else:
            kmin = kmid
        if kmax-kmin <= 1:
            break
    return kmax # Chai Wah Wu, Aug 14 2024

For n > 1: GCD(exponents in prime factorization of a(n)) > 2, cf. A124010. - Reinhard Zumkeller, Apr 13 2012

Sum_{n>=1} 1/a(n) = 2 - zeta(2) + Sum_{k>=2} mu(k)*(2 - zeta(k) - zeta(2*k)) = 1.3300056287... - Amiram Eldar, Jul 02 2022

Edited by Robert Israel, Sep 30 2015

A000379 Numbers where total number of 1-bits in the exponents of their prime factorization is even; a 2-way classification of integers: complement of A000028.

Original entry on oeis.org

1, 6, 8, 10, 12, 14, 15, 18, 20, 21, 22, 26, 27, 28, 32, 33, 34, 35, 36, 38, 39, 44, 45, 46, 48, 50, 51, 52, 55, 57, 58, 62, 63, 64, 65, 68, 69, 74, 75, 76, 77, 80, 82, 85, 86, 87, 91, 92, 93, 94, 95, 98, 99, 100, 106, 111, 112, 115, 116, 117, 118, 119, 120, 122, 123, 124, 125, 129

Offset: 1

N. J. A. Sloane

This sequence and A000028 (its complement) give the unique solution to the problem of splitting the positive integers into two classes in such a way that products of pairs of distinct elements from either class occur with the same multiplicities [Lambek and Moser]. Cf. A000069, A001969.
See A000028 for precise definition, Maple program, etc.
The sequence contains products of even number of distinct terms of A050376. - Vladimir Shevelev, May 04 2010
From Vladimir Shevelev, Oct 28 2013: (Start)
Numbers m such that the infinitary Möbius function (A064179) of m equals 1. (This follows from the definition of A064179.)
A number m is in the sequence iff the number k = k(m) of terms of A050376 that divide m with odd maximal exponent is even (see example).
(End)
Numbers k for which A064547(k) [or equally, A268386(k)] is even. Numbers k for which A010060(A268387(k)) = 0. - Antti Karttunen, Feb 09 2016
The sequence is closed under the commutative binary operation A059897(.,.). As integers are self-inverse under A059897(.,.), it therefore forms a subgroup of the positive integers considered as a group under A059897(.,.). Specifically (expanding on the comment above dated May 04 2010) it is the subgroup of even length words in A050376, which is the group's lexicographically earliest ordered minimal set of generators. A000028, the set of odd length words in A050376, is its complementary coset. - Peter Munn, Nov 01 2019
From Amiram Eldar, Oct 02 2024: (Start)
Numbers whose number of infinitary divisors (A037445) is a square.
Numbers whose exponentially odious part (A367514) has an even number of distinct prime factors, i.e., numbers k such that A092248(A367514(k)) = 0. (End)

			If m = 120, then the maximal exponent of 2 that divides 120 is 3, for 3 it is 1, for 4 it is 1, for 5 it is 1. Thus k(120) = 4 and 120 is a term. - _Vladimir Shevelev_, Oct 28 2013

Joe Roberts, Lure of the Integers, Math. Assoc. America, 1992, p. 22.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

N. J. A. Sloane, Table of n, a(n) for n = 1..10000
J. Lambek and L. Moser, On some two way classifications of integers, Canad. Math. Bull. 2 (1959), 85-89.
Eric Weisstein's World of Mathematics, Group.
Wikipedia, Generating set of a group.

Subsequences: A030229, A238748, A262675, A268390.
Subsequence of A268388 (apart from the initial 1).
Complement: A000028.
Sequences used in definitions of this sequence: A133008, A050376, A059897, A064179, A064547, A124010 (prime exponents), A268386, A268387, A010060.
Other 2-way classifications: A000069/A001969 (to which A000120 and A010060 are relevant), A000201/A001950.
This is different from A123240 (e.g., does not contain 180). The first difference occurs already at n=31, where A123240(31) = 60, a value which does not occur here, as a(31+1) = 62. The same is true with respect to A131181, as A131181(31) = 60.
Cf. A030231, A037445, A092248, A367514.

a000379 n = a000379_list !! (n-1)
a000379_list = filter (even . sum . map a000120 . a124010_row) [1..]
-- Reinhard Zumkeller, Oct 05 2011

Select[ Range[130], EvenQ[ Count[ Flatten[ IntegerDigits[#, 2]& /@ Transpose[ FactorInteger[#]][[2]]], 1]]&] // Prepend[#, 1]& (* Jean-François Alcover, Apr 11 2013, after Harvey P. Dale *)

is(n)=my(f=factor(n)[,2]); sum(i=1,#f,hammingweight(f[i]))%2==0 \\ Charles R Greathouse IV, Aug 31 2013
(Scheme, two variants)
(define A000379 (MATCHING-POS 1 1 (COMPOSE even? A064547)))
(define A000379 (MATCHING-POS 1 1 (lambda (n) (even? (A000120 (A268387 n))))))
;; Both require also my IntSeq-library. - Antti Karttunen, Feb 09 2016

Edited by N. J. A. Sloane, Dec 20 2007, to restore the original definition.

A100716 Numbers k such that p^p divides k for some prime p.

Original entry on oeis.org

4, 8, 12, 16, 20, 24, 27, 28, 32, 36, 40, 44, 48, 52, 54, 56, 60, 64, 68, 72, 76, 80, 81, 84, 88, 92, 96, 100, 104, 108, 112, 116, 120, 124, 128, 132, 135, 136, 140, 144, 148, 152, 156, 160, 162, 164, 168, 172, 176, 180, 184, 188, 189, 192, 196, 200, 204, 208, 212

Offset: 1

Leroy Quet, Dec 10 2004

nonn

Complement of A048103; A129251(a(n)) > 0; A051674 is a subsequence; A129254 = (terms a(k) such that a(k+1)=a(k)+1). - Reinhard Zumkeller, Apr 07 2007

A027748(a(n),k) <= A124010(a(n),k) for some k<=A001221(a(n)). - Reinhard Zumkeller, Apr 28 2012

			54 is included because 3^3 divides 54.

Reinhard Zumkeller, Table of n, a(n) for n = 1..1000

Complement: A048103.
Positions of nonzeros in A129251.
Cf. A100717, A129150, A129152.
Cf. A054744.
Cf. A051674 (a subsequence).
Subsequence of A276079 from which it differs for the first time at n=175, where a(175) = 628, while A276079(175) = 625, a value missing from here.

a100716 n = a100716_list !! (n-1)
a100716_list = filter (\x -> or $
   zipWith (<=) (a027748_row x) (map toInteger $ a124010_row x)) [1..]
-- Reinhard Zumkeller, Apr 28 2012
(Scheme, with Antti Karttunen's IntSeq-library)
(define A100716 (NONZERO-POS 1 1 A129251))
;; Antti Karttunen, Aug 18 2016

fQ[n_] := Union[ Table[ #[[1]] <= #[[2]]] & /@ FactorInteger[n]][[ -1]] == True; Select[ Range[2, 215], fQ[ # ] &] (* Robert G. Wilson v, Dec 14 2004 *)
f[n_] := Module[{aux=FactorInteger[n]}, Last@Union@Table[aux[[i,1]] <=  aux[[i,2]], {i,Length[aux]}] == True]; Select[Range[2,215], f] (* José María Grau Ribas, Jan 25 2012 *)
Rest@ Select[Range@ 216, Times @@ Boole@ Map[First@ # > Last@ # &, FactorInteger@ #] == 0 &] (* Michael De Vlieger, Aug 19 2016 *)

is(n)=forprime(p=2,default(primelimit),if(n%p^p==0,return(1));if(p^p>n,return(0))) \\ Charles R Greathouse IV, Jan 24 2012

from itertools import count, islice
from sympy import factorint
def A100716_gen(startvalue=1): # generator of terms >= startvalue
    return filter(lambda n:any(map(lambda d:d[1]>=d[0],factorint(n).items())),count(max(startvalue,1)))
A100716_list = list(islice(A100716_gen(),30)) # Chai Wah Wu, Jan 05 2023

a(n) ~ k*n with k = 1/(1 - Product(1 - p^-p)) = 3.5969959469... where the product is over all primes p. - Charles R Greathouse IV, Jan 24 2012

More terms from Robert G. Wilson v, Dec 14 2004

A112526 Characteristic function for powerful numbers.

Original entry on oeis.org

1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0

Offset: 1

Franklin T. Adams-Watters, Sep 09 2005

A signed multiplicative variant is defined by b(n) = a(n)*mu(n) with mu = A008683, such that b(p^e)=0 if e=1 and b(p^e)= -1 if e>1. This has Dirichlet series Sum_{n>=1} b(n)/n = A005596 and Sum_{n>=1} b(n)/n^2 = A065471. - R. J. Mathar, Apr 04 2011

			a(72) = 1 because 72 = 2^3*3^2 has all exponents > 1.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Powerful Number.
Index entries for characteristic functions.

Differs from characteristic function of perfect powers A075802 at Achilles numbers A052486.
Cf. A001694 (powerful numbers), A124010, A001221, A027746.
Cf. A008683, A008966, A005596, A065471, A082695, A063524.
Cf. A002117, A078434, A005361.

a112526 1 = 1
a112526 n = fromEnum $ (> 1) $ minimum $ a124010_row n
-- Reinhard Zumkeller, Jun 03 2015, Sep 16 2011

cfpn[n_]:=If[n==1||Min[Transpose[FactorInteger[n]][[2]]]>1,1,0]; Array[ cfpn,120] (* Harvey P. Dale, Jul 17 2012 *)

for(n=1, 100, print1(direuler(p=2, n, (1+X^3)/(1-X^2))[n], ", ")) \\ Vaclav Kotesovec, Jul 15 2022

a(n) = ispowerful(n); \\ Amiram Eldar, Jul 02 2025

from sympy import factorint
def A112526(n): return int(all(e>1 for e in factorint(n).values())) # Chai Wah Wu, Sep 15 2024

Multiplicative with a(p^e) = 1 - 0^(e-1), e > 0 and p prime.
Dirichlet g.f.: zeta(2*s)*zeta(3*s)/zeta(6*s), e.g., A082695 at s=1.
a(n) * A008966(n) = A063524(n). - Reinhard Zumkeller, Sep 16 2011
a(n) = {m: Min{A124010(m,k): k=1..A001221(m)} > 1}. - Reinhard Zumkeller, Jun 03 2015
Sum_{k=1..n} a(k) ~ zeta(3/2)*sqrt(n)/zeta(3) + 6*zeta(2/3)*n^(1/3)/Pi^2. - Vaclav Kotesovec, Feb 08 2019
a(n) = Sum_{d|n} A005361(d)*A008683(n/d). - Ridouane Oudra, Jul 03 2025

cp's OEIS Frontend

A046897 Sum of divisors of n that are not divisible by 4.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Formula

A209061 Exponentially squarefree numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Formula

A355536 Irregular triangle read by rows where row n lists the differences between adjacent prime indices of n; if n is prime, row n is empty.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A332282 Numbers whose unsorted prime signature is not unimodal.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A359178 Numbers with a unique smallest prime exponent.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Python

A055229 Greatest common divisor of largest square dividing n and squarefree part of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Formula

A076467 Perfect powers m^k where m is a positive integer and k >= 3.

Views

Author

Keywords

Comments

Links