A070321 -id:A070321 - cp's OEIS frontend

A365683 The largest exponentially squarefree divisor of n.

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 8, 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, 24, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68

Offset: 1

Amiram Eldar, Sep 15 2023

First differs from A058035 at n = 32.

The number of these divisors is A365680(n) and their sum is A365682(n).

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A070321, A209061, A365680, A365682.

f[p_, e_] := Module[{k = e}, While[! SquareFreeQ[k], k--]; p^k]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]

s(e) = {my(k = e); while(!issquarefree(k), k--); k;};
a(n) = {my(f = factor(n)); prod(i = 1, #f~, f[i,1]^s(f[i,2]));}

Multiplicative with a(p^e) = p^A070321(e).
a(n) <= n, with equality if and only if n is exponentially squarefree number (A209061).
Sum_{k=1..n} a(k) ~ c*n^2, where c = 0.487850776747... = (1/2) * Product_{p prime} (1 + Sum_{k>=1} (p^f(k) - p^(f(k-1)+1))/p^(2*k)), f(k) = A070321(k) and f(0) = 0.

A378087 First-differences of A067535 (least positive integer >= n that is squarefree).

Original entry on oeis.org

1, 1, 2, 0, 1, 1, 3, 0, 0, 1, 2, 0, 1, 1, 2, 0, 2, 0, 2, 0, 1, 1, 3, 0, 0, 3, 0, 0, 1, 1, 2, 0, 1, 1, 2, 0, 1, 1, 2, 0, 1, 1, 3, 0, 0, 1, 4, 0, 0, 0, 2, 0, 2, 0, 2, 0, 1, 1, 2, 0, 1, 3, 0, 0, 1, 1, 2, 0, 1, 1, 2, 0, 1, 3, 0, 0, 1, 1, 3, 0, 0, 1, 2, 0, 1, 1, 2

Offset: 1

Gus Wiseman, Nov 20 2024

nonn

Does this contain all nonnegative integers? The positions of first appearances begin: 4, 1, 3, 7, 47, 241, 843, 22019, 217069, ...

Ones are A007674.
Zeros are A013929, complement A005117.
Positions of first appearances are A020754 (except first term) = A045882 - 1.
First-differences of A067535.
Twos are A280892.
For prime-powers we have A377780, differences of A000015.
The nonsquarefree opposite is A378036, differences of A378033.
The restriction to primes + 1 is A378037 (opposite A378038), differences of A112926.
For nonsquarefree numbers we have A378039, see A377783, A377784, A378040.
The opposite is A378085, differences of A070321.
A000040 lists the primes, differences A001223, seconds A036263.
A005117 lists the squarefree numbers.
A013929 lists the nonsquarefree numbers, differences A078147, seconds A376593.
A061398 counts squarefree numbers between primes, zeros A068360.
A061399 counts nonsquarefree numbers between primes, zeros A068361.
Cf. A005117, A007675, A068781, A073247, A073248, A378084.
Cf. A013928, A053797, A053806, A072284, A120327, A224363.

Differences[Table[NestWhile[#+1&,n,#>1&&!SquareFreeQ[#]&],{n,100}]]

A379306 Number of squarefree prime indices of n.

Original entry on oeis.org

0, 1, 1, 2, 1, 2, 0, 3, 2, 2, 1, 3, 1, 1, 2, 4, 1, 3, 0, 3, 1, 2, 0, 4, 2, 2, 3, 2, 1, 3, 1, 5, 2, 2, 1, 4, 0, 1, 2, 4, 1, 2, 1, 3, 3, 1, 1, 5, 0, 3, 2, 3, 0, 4, 2, 3, 1, 2, 1, 4, 0, 2, 2, 6, 2, 3, 1, 3, 1, 2, 0, 5, 1, 1, 3, 2, 1, 3, 1, 5, 4, 2, 1, 3, 2, 2, 2

Offset: 1

Gus Wiseman, Dec 25 2024

nonn

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 prime indices of 39 are {2,6}, so a(39) = 2.
The prime indices of 70 are {1,3,4}, so a(70) = 2.
The prime indices of 98 are {1,4,4}, so a(98) = 1.
The prime indices of 294 are {1,2,4,4}, a(294) = 2.
The prime indices of 1911 are {2,4,4,6}, so a(1911) = 2.
The prime indices of 2548 are {1,1,4,4,6}, so a(2548) = 3.

Positions of first appearances are A000079.
Positions of zero are A379307, counted by A114374 (strict A256012).
Positions of one are A379316, counted by A379308 (strict A379309).
A000040 lists the primes, differences A001223.
A005117 lists the squarefree numbers, differences A076259.
A008966 is the characteristic function for the squarefree numbers.
A013929 lists the nonsquarefree numbers, differences A078147.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A061398 counts squarefree numbers between primes, zeros A068360.
A377038 gives k-th differences of squarefree numbers.
Other counts of prime indices:
- A087436 postpositive, see A038550.
- A257991 odd, see A000041, A000070, A066207, A349158.
- A257992 even, see A000009, A038348, A066208, A379317.
- A257994 prime, see A002095, A096258, A320628, A331386, A331915, A379304, A379305.
- A330944 nonprime, see A000586, A000607, A076610, A330945.
- A379300 composite, see A023895, A034891, A036497, A302540, A379301.
- A379310 nonsquarefree, see A302478.
- A379311 old prime, see A204389, A320629, A379312-A379315.
Cf. A000720, A013928, A057627, A068361, A070321, A071403, A072284, A112925, A112929, A120327, A377430, A378086.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Select[prix[n],SquareFreeQ]],{n,100}]

Totally additive with a(prime(k)) = A008966(k).

A081217 Greatest squarefree number not exceeding n-th prime power.

Original entry on oeis.org

1, 2, 3, 3, 5, 7, 7, 7, 11, 13, 15, 17, 19, 23, 23, 26, 29, 31, 31, 37, 41, 43, 47, 47, 53, 59, 61, 62, 67, 71, 73, 79, 79, 83, 89, 97, 101, 103, 107, 109, 113, 119, 123, 127, 127, 131, 137, 139, 149, 151, 157, 163, 167, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223

Offset: 1

Reinhard Zumkeller, Mar 10 2003

nonn

a(n) = A070321(A000961(n)).

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A081218.

Cf. A156898, A156899, A156900, A005117, A000961

n:= 0:
for m from 1 while n <= 100 do
  if numtheory:-issqrfree(m) then s:= m fi;
  if nops(numtheory:-factorset(m)) <= 1 then n:= n+1; A[n]:= s fi
od:
seq(A[i],i=1..100); # Robert Israel, Jan 08 2018

nmax = 100; coeff = 5;
A000961 = Join[{1}, Select[Range[2, coeff*nmax],PrimePowerQ]];
a070321[n_] := For[k = n, True, k--, If[SquareFreeQ[k], Return[k]]];
a[n_] := If[n > Length[A000961], Print["error: increase insufficient coeff ", coeff], a070321[A000961[[n]] ] ];
Array[a, nmax] (* Jean-François Alcover, Mar 08 2019 *)

A373415 Maximum of the n-th maximal run of squarefree numbers.

Original entry on oeis.org

3, 7, 11, 15, 17, 19, 23, 26, 31, 35, 39, 43, 47, 51, 53, 55, 59, 62, 67, 71, 74, 79, 83, 87, 89, 91, 95, 97, 103, 107, 111, 115, 119, 123, 127, 131, 134, 139, 143, 146, 149, 151, 155, 159, 161, 163, 167, 170, 174, 179, 183, 187, 191, 195, 197, 199, 203, 206

Offset: 1

Gus Wiseman, Jun 05 2024

nonn

The minimum is given by A072284.
A run of a sequence (in this case A005117) is an interval of positions at which consecutive terms differ by one.
Consists of all squarefree numbers k such that k + 1 is not squarefree.

			Row-maxima of:
   1   2   3
   5   6   7
  10  11
  13  14  15
  17
  19
  21  22  23
  26
  29  30  31
  33  34  35
  37  38  39
  41  42  43
  46  47
  51
  53
  55
  57  58  59

Gus Wiseman, Four statistics for runs and antiruns of prime, nonprime, squarefree, and nonsquarefree numbers

Functional neighbors: A006093, A007674, A067774, A072284, A120992, A373413.
A005117 lists the squarefree numbers, first differences A076259.
A013929 lists the nonsquarefree numbers, first differences A078147.
Cf. A049093, A049094, A061398, A069010, A077641, A077643, A112925, A112926, A143658, A372683, A373125, A373126.

Last/@Split[Select[Range[100],SquareFreeQ],#1+1==#2&]//Most

a(n) = A070321(A072284(n+1) - 1).

A378358 Least non-perfect-power >= n.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Nov 24 2024

nonn

Perfect-powers (A001597) are 1 and numbers with a proper integer root, complement A007916.

The version for prime-powers is A000015, for non-prime-powers A378372.
The union is A007916, complement A001597.
The version for nonsquarefree numbers is A067535, negative A120327 (subtract A378369).
The version for composite numbers is A113646.
The version for prime numbers is A159477.
The run-lengths are A375706.
Terms appearing only once are A375738, multiple times A375703.
The version for perfect-powers is A377468.
Subtracting from n gives A378357.
The opposite version is A378363, for perfect-powers A081676.
A000040 lists the primes, differences A001223.
A000961 lists the powers of primes, differences A057820.
A001597 lists the perfect-powers, differences A053289.
A007916 lists the non-perfect-powers, differences A375706.
A069623 counts perfect-powers <= n.
A076411 counts perfect-powers < n.
A131605 lists perfect-powers that are not prime-powers.
A377432 counts perfect-powers between primes, zeros A377436.
Cf. A007918, A014210, A031218, A045542, A052410, A065514, A070321, A076412, A188951, A216765, A345531, A378033.

perpowQ[n_]:=n==1||GCD@@FactorInteger[n][[All,2]]>1;
Table[NestWhile[#+1&,n,perpowQ[#]&],{n,100}]

from sympy import mobius, integer_nthroot
def A378358(n):
    def f(x): return int(1-sum(mobius(k)*(integer_nthroot(x,k)[0]-1) for k in range(2,x.bit_length())))
    a = max(1,n-f(n-1))
    m, k = a, f(a)+a
    while m != k: m, k = k, f(k)+a
    return m # Chai Wah Wu, Nov 26 2024

from sympy import perfect_power
def A378358(n): return n if n>1 and perfect_power(n)==False else n+1 if perfect_power(n+1)==False else n+2 # Chai Wah Wu, Nov 27 2024

a(n) = n - A378357(n).

A377782 First-differences of A031218(n) = greatest number <= n that is 1 or a prime-power.

Original entry on oeis.org

1, 1, 1, 1, 0, 2, 1, 1, 0, 2, 0, 2, 0, 0, 3, 1, 0, 2, 0, 0, 0, 4, 0, 2, 0, 2, 0, 2, 0, 2, 1, 0, 0, 0, 0, 5, 0, 0, 0, 4, 0, 2, 0, 0, 0, 4, 0, 2, 0, 0, 0, 4, 0, 0, 0, 0, 0, 6, 0, 2, 0, 0, 3, 0, 0, 3, 0, 0, 0, 4, 0, 2, 0, 0, 0, 0, 0, 6, 0, 2, 0, 2, 0, 0, 0, 0, 0

Offset: 1

Gus Wiseman, Nov 16 2024

nonn

Note 1 is a power of a prime (A000961) but not a prime-power (A246655).

Positions of 1 are A006549.
Positions of 0 are A080765 = A024619 - 1, complement A181062 = A000961 - 1.
Positions of 2 are A120432 (except initial terms).
Sorted positions of first appearances appear to include A167236 - 1.
Positions of terms > 1 are A373677.
The restriction to primes minus 1 is A377289.
Below, A (B) indicates that A is the first-differences of B:
- This sequence is A377782 (A031218), which has restriction to primes A065514 (A377781).
- The opposite is A377780 (A000015), restriction A377703 (A345531).
- For nonsquarefree we have A378036 (A378033), opposite A378039 (A120327).
- For squarefree we have A378085 (A112925), restriction A378038 (A070321).
A000040 lists the primes, differences A001223.
A000961 and A246655 list prime-powers, differences A057820.
A024619 lists the non-prime-powers, differences A375735, seconds A376599.
A361102 lists the non-powers of primes, differences A375708.
A378034 gives differences of A378032 (restriction of A378033).
Prime-powers between primes: A053607, A080101, A366833, A377057, A377286, A377287.
Cf. A053289, A343249, A375706, A377281, A377282.

Differences[Table[NestWhile[#-1&,n,#>1&&!PrimePowerQ[#]&],{n,100}]]

A378363 Greatest number <= n that is 1 or not a perfect-power.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Nov 24 2024

Perfect-powers (A001597) are 1 and numbers with a proper integer root, complement A007916.

			In the non-perfect-powers ... 5, 6, 7, 10, 11 ... the greatest term <= 8 is 7, so a(8) = 7.

The union is A007916, complement A001597.
The version for prime numbers is A007917 or A151799, opposite A159477.
The version for prime-powers is A031218, opposite A000015.
The version for squarefree numbers is A067535, opposite A070321.
The version for perfect-powers is A081676, opposite A377468.
The version for composite numbers is A179278, opposite A113646.
Terms appearing multiple times are A375704, opposite A375703.
The run-lengths are A375706.
Terms appearing only once are A375739, opposite A375738.
The version for nonsquarefree numbers is A378033, opposite A120327.
The opposite version is A378358.
Subtracting n gives A378364, opposite A378357.
The version for non-prime-powers is A378367 (subtracted A378371), opposite A378372.
A000040 lists the primes, differences A001223.
A000961 lists the powers of primes, differences A057820.
A001597 lists the perfect-powers, differences A053289.
A007916 lists the non-perfect-powers, differences A375706.
A069623 counts perfect-powers <= n.
A076411 counts perfect-powers < n.
A131605 lists perfect-powers that are not prime-powers.
A377432 counts perfect-powers between primes, zeros A377436.
Cf. A007918, A014210, A014234, A045542, A052410, A065514, A076412, A162966, A188951, A345531.

perpowQ[n_]:=n==1||GCD@@FactorInteger[n][[All,2]]>1;
Table[NestWhile[#-1&,n,#>1&&perpowQ[#]&],{n,100}]

from sympy import mobius, integer_nthroot
def A378363(n):
    def f(x): return int(1-sum(mobius(k)*(integer_nthroot(x,k)[0]-1) for k in range(2,x.bit_length())))
    a = n-f(n)
    m, k = a, f(a)+a
    while m != k: m, k = k, f(k)+a
    return m # Chai Wah Wu, Nov 26 2024

A378367 Greatest non prime power <= n, allowing 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 6, 6, 6, 6, 10, 10, 12, 12, 14, 15, 15, 15, 18, 18, 20, 21, 22, 22, 24, 24, 26, 26, 28, 28, 30, 30, 30, 33, 34, 35, 36, 36, 38, 39, 40, 40, 42, 42, 44, 45, 46, 46, 48, 48, 50, 51, 52, 52, 54, 55, 56, 57, 58, 58, 60, 60, 62, 63, 63, 65, 66, 66

Offset: 1

Gus Wiseman, Nov 29 2024

Non prime powers allowing 1 (A361102) are numbers that are not a prime power (A246655), namely 1, 6, 10, 12, 14, 15, 18, 20, 21, 22, 24, ...

			The greatest non prime power <= 7 is 6, so a(7) = 6.

Sequences obtained by subtracting each term from n are placed in parentheses below.
For prime we have A007917 (A064722).
For nonprime we have A179278 (A010051 almost).
For perfect power we have A081676 (A069584).
For squarefree we have A070321.
For nonsquarefree we have A378033.
For non perfect power we have A378363.
The opposite is A378372, subtracting n A378371.
For prime power we have A031218 (A276781 - 1).
Subtracting from n gives (A378366).
A000015 gives the least prime power >= n (A378370).
A000040 lists the primes, differences A001223.
A000961 and A246655 list the prime powers, differences A057820.
A024619 and A361102 list the non prime powers, differences A375708 and A375735.
A151800 gives the least prime > n (A013632), weak version A007918 (A007920).
Prime powers between primes: A053607, A080101, A304521, A366833, A377057.
Cf. A007916, A065514, A113646, A345531 (A377281), A377051, A377054, A377282, A377468 (A074984), A378358, A378457.
Cf. A356068.

Table[NestWhile[#-1&,n,PrimePowerQ[#]&],{n,100}]

a(n) = n - A378366(n).

a(n) = A361102(A356068(n)). - Ridouane Oudra, Aug 22 2025

A081210 In prime factorization of n replace each prime power p^e with the greatest squarefree number <= p^e.

Original entry on oeis.org

1, 2, 3, 3, 5, 6, 7, 7, 7, 10, 11, 9, 13, 14, 15, 15, 17, 14, 19, 15, 21, 22, 23, 21, 23, 26, 26, 21, 29, 30, 31, 31, 33, 34, 35, 21, 37, 38, 39, 35, 41, 42, 43, 33, 35, 46, 47, 45, 47, 46, 51, 39, 53, 52, 55, 49, 57, 58, 59, 45, 61, 62, 49, 62, 65, 66, 67, 51, 69, 70, 71, 49, 73

Offset: 1

Reinhard Zumkeller, Mar 10 2003

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

Cf. A005117, A070321, A081212, A081213, A081214.

A081210 := proc(n)
    local a,pe;
    a :=1 ;
    for pe in ifactors(n)[2] do
        a := a*A070321(op(1,pe)^op(2,pe)) ;
    end do:
    a ;
end proc:
seq(A081210(n),n=1..100) ; # R. J. Mathar, May 25 2023

gsf[n_] := For[k = n, True, k--, If[ SquareFreeQ[k], Return[k]]]; a[n_] := Times @@ gsf /@ Power @@@ FactorInteger[n]; Table[a[n], {n, 1, 80}] (* Jean-François Alcover, Mar 27 2013 *)

a(n) = {my(f = factor(n)); prod(i = 1, #f~, if(f[i,2] == 1, f[i,1], my(k = f[i,1]^f[i,2]); while(!issquarefree(k), k--); k));} \\ Amiram Eldar, Jun 09 2025

Multiplicative with a(p^e) = A070321(p^e), p prime.
a(n) <= n and a(n) = n iff n is squarefree (A005117).
A081211(n) = a(a(n)), see A081212, A081213 and A081214 for iterations until a fixed point is reached.

cp's OEIS Frontend

A365683 The largest exponentially squarefree divisor of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A378087 First-differences of A067535 (least positive integer >= n that is squarefree).

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

A379306 Number of squarefree prime indices of n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A081217 Greatest squarefree number not exceeding n-th prime power.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

A373415 Maximum of the n-th maximal run of squarefree numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A378358 Least non-perfect-power >= n.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

Python

Python

Formula

A377782 First-differences of A031218(n) = greatest number <= n that is 1 or a prime-power.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

A378363 Greatest number <= n that is 1 or not a perfect-power.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Python