A259936 -id:A259936 - cp's OEIS frontend

A316441 a(n) = Sum (-1)^k where the sum is over all factorizations of n into factors > 1 and k is the number of factors.

1, -1, -1, 0, -1, 0, -1, -1, 0, 0, -1, 0, -1, 0, 0, 1, -1, 0, -1, 0, 0, 0, -1, 1, 0, 0, -1, 0, -1, 1, -1, -1, 0, 0, 0, 1, -1, 0, 0, 1, -1, 1, -1, 0, 0, 0, -1, 0, 0, 0, 0, 0, -1, 1, 0, 1, 0, 0, -1, 1, -1, 0, 0, 1, 0, 1, -1, 0, 0, 1, -1, 0, -1, 0, 0, 0, 0, 1, -1

Offset: 1

Gus Wiseman, Jul 03 2018

sign

First term greater than 1 in absolute value is a(256) = 2.

			The factorizations of 24 are (2*2*2*3), (2*2*6), (2*3*4), (2*12), (3*8), (4*6), (24); so a(24) = 1 - 2 + 3 - 1 = 1.

Antti Karttunen, Table of n, a(n) for n = 1..65537
Index entries for sequences computed from exponents in factorization of n

Cf. A001222, A001055, A045778, A114592, A162247, A190938, A259936, A281116, A303386.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Sum[(-1)^Length[f],{f,facs[n]}],{n,200}]

A316441(n, m=n, k=0) = if(1==n, (-1)^k, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A316441(n/d, d, k+1))); (s)); \\ Antti Karttunen, Sep 08 2018, after Michael B. Porter's code for A001055

Dirichlet g.f.: Product_{n > 1} 1/(1 + 1/n^s).

Secondary offset added by Antti Karttunen, Sep 08 2018

A294617 Number of ways to choose a set partition of a strict integer partition of n.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 10, 12, 17, 24, 44, 51, 76, 98, 138, 217, 272, 366, 493, 654, 848, 1284, 1560, 2115, 2718, 3610, 4550, 6024, 8230, 10296, 13354, 17144, 21926, 27903, 35556, 44644, 59959, 73456, 94109, 117735, 150078, 185800, 235719, 290818, 365334, 467923

Offset: 0

Gus Wiseman, Nov 05 2017

nonn

From Gus Wiseman, Sep 17 2024: (Start)
Also the number of strict integer compositions of n whose leaders, obtained by splitting into maximal increasing subsequences and taking the first term of each, are decreasing. For example, the strict composition (3,6,2,1,4) has maximal increasing subsequences ((3,6),(2),(1,4)), with leaders (3,2,1), so is counted under a(16). The a(0) = 1 through a(7) = 12 compositions are:
  ()  (1)  (2)  (3)    (4)    (5)    (6)      (7)
                (1,2)  (1,3)  (1,4)  (1,5)    (1,6)
                (2,1)  (3,1)  (2,3)  (2,4)    (2,5)
                              (3,2)  (4,2)    (3,4)
                              (4,1)  (5,1)    (4,3)
                                     (1,2,3)  (5,2)
                                     (2,1,3)  (6,1)
                                     (2,3,1)  (1,2,4)
                                     (3,1,2)  (2,1,4)
                                     (3,2,1)  (2,4,1)
                                              (4,1,2)
                                              (4,2,1)
(End)

			The a(6) = 10 set partitions are: {{6}}, {{1},{5}}, {{5,1}}, {{2},{4}}, {{4,2}}, {{1},{2},{3}}, {{1},{3,2}}, {{2,1},{3}}, {{3,1},{2}}, {{3,2,1}}.

Alois P. Heinz, Table of n, a(n) for n = 0..5000
Gus Wiseman, Sequences counting and ranking compositions by their leaders (for six types of runs).

Cf. A000009, A000110, A063834, A258466, A259936, A279375, A279785, A279790.
Row sums of A330460 and of A330759.
This is a strict case of A374689, weak version A189076.
A011782 counts compositions, strict A032020.
A238130, A238279, A333755 count compositions by number of runs.
Cf. A188920, A238343, A274174, A373949, A374518, A376263.
Cf. A374631, A374632, A374634, A374635, A374686, A374688.

b:= proc(n, i, t) option remember; `if`(n>i*(i+1)/2, 0,
      `if`(n=0, combinat[bell](t), b(n, i-1, t)+
      `if`(i>n, 0, b(n-i, min(n-i, i-1), t+1))))
    end:
a:= n-> b(n$2, 0):
seq(a(n), n=0..50);  # Alois P. Heinz, Nov 07 2017

Table[Total[BellB[Length[#]]&/@Select[IntegerPartitions[n],UnsameQ@@#&]],{n,25}]
(* Second program: *)
b[n_, i_, t_] := b[n, i, t] = If[n > i (i + 1)/2, 0, If[n == 0, BellB[t], b[n, i - 1, t] + If[i > n, 0, b[n - i, Min[n - i, i - 1], t + 1]]]];
a[n_] := b[n, n, 0];
a /@ Range[0, 50] (* Jean-François Alcover, May 10 2021, after Alois P. Heinz *)

A279375(n) <= a(n) <= A279790(n).

G.f.: Sum_{k>=0} Bell(k) * x^(k*(k + 1)/2) / Product_{j=1..k} (1 - x^j). - Ilya Gutkovskiy, Jan 28 2020

A305149 Number of factorizations of n whose distinct factors are pairwise indivisible and greater than 1.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 2, 2, 2, 1, 3, 1, 2, 2, 3, 1, 3, 1, 3, 2, 2, 1, 4, 2, 2, 2, 3, 1, 5, 1, 2, 2, 2, 2, 6, 1, 2, 2, 4, 1, 5, 1, 3, 3, 2, 1, 5, 2, 3, 2, 3, 1, 4, 2, 4, 2, 2, 1, 8, 1, 2, 3, 4, 2, 5, 1, 3, 2, 5, 1, 6, 1, 2, 3, 3, 2, 5, 1, 5, 3, 2, 1, 8, 2, 2, 2, 4, 1, 8, 2, 3, 2, 2, 2, 6, 1, 3, 3, 6, 1, 5, 1, 4, 5

Offset: 1

Gus Wiseman, May 26 2018

nonn

			The a(60) = 8 factorizations are (2*2*3*5), (2*2*15), (3*4*5), (3*20), (4*15), (5*12), (6*10), (60). Missing from this list are (2*3*10), (2*5*6), (2*30).

Antti Karttunen, Table of n, a(n) for n = 1..10000
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..100000

Cf. A001055, A001970, A007716, A034444, A045778, A259936, A281116, A285572, A302242, A303386, A303431, A305001, A305148, A305150.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[Select[facs[n],Select[Tuples[Union[#],2],UnsameQ@@#&&Divisible@@#&]=={}&]],{n,100}]

pairwise_indivisible(v) = { for(i=1,#v,for(j=i+1,#v,if(!(v[j]%v[i]),return(0)))); (1); };
A305149(n, m=n, facs=List([])) = if(1==n, pairwise_indivisible(Set(facs)), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A305149(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Oct 08 2018

More terms from Antti Karttunen, Oct 08 2018

A340608 The number of prime factors of n (A001222) is relatively prime to the maximum prime index of n (A061395).

Original entry on oeis.org

2, 3, 4, 5, 7, 8, 10, 11, 12, 13, 15, 16, 17, 18, 19, 22, 23, 25, 27, 28, 29, 31, 32, 33, 34, 37, 40, 41, 42, 43, 44, 46, 47, 48, 51, 53, 55, 59, 60, 61, 62, 63, 64, 66, 67, 68, 69, 70, 71, 72, 73, 76, 77, 79, 80, 82, 83, 85, 88, 89, 90, 93, 94, 97, 98, 99

Offset: 1

Gus Wiseman, Jan 27 2021

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 sequence of terms together with their prime indices begins:
     2: {1}          22: {1,5}          44: {1,1,5}
     3: {2}          23: {9}            46: {1,9}
     4: {1,1}        25: {3,3}          47: {15}
     5: {3}          27: {2,2,2}        48: {1,1,1,1,2}
     7: {4}          28: {1,1,4}        51: {2,7}
     8: {1,1,1}      29: {10}           53: {16}
    10: {1,3}        31: {11}           55: {3,5}
    11: {5}          32: {1,1,1,1,1}    59: {17}
    12: {1,1,2}      33: {2,5}          60: {1,1,2,3}
    13: {6}          34: {1,7}          61: {18}
    15: {2,3}        37: {12}           62: {1,11}
    16: {1,1,1,1}    40: {1,1,1,3}      63: {2,2,4}
    17: {7}          41: {13}           64: {1,1,1,1,1,1}
    18: {1,2,2}      42: {1,2,4}        66: {1,2,5}
    19: {8}          43: {14}           67: {19}

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

Note: Heinz numbers are given in parentheses below.
These are the Heinz numbers of the partitions counted by A200750.
The case of equality is A047993 (A106529).
The divisible instead of coprime version is A168659 (A340609).
The dividing instead of coprime version is A168659 (A340610), with strict case A340828 (A340856).
A001222 counts prime factors.
A006141 counts partitions whose length equals their minimum (A324522).
A051424 counts singleton or pairwise coprime partitions (A302569).
A056239 adds up prime indices.
A061395 selects the maximum prime index.
A067538 counts partitions whose length divides their sum (A316413).
A067538 counts partitions whose maximum divides their sum (A326836).
A112798 lists the prime indices of each positive integer.
A259936 counts singleton or pairwise coprime factorizations.
A325134 = A001222 + A061395.
A326845 = A056239 * A061395.
A326849 counts partitions whose sum divides length times maximum (A326848).
A327516 counts pairwise coprime partitions (A302696).
Cf. A003114, A039900, A096401, A244990/A244991, A340606, A340653, A340691, A340787/A340788.

Select[Range[100],GCD[PrimeOmega[#],PrimePi[FactorInteger[#][[-1,1]]]]==1&]

A304717 Number of connected strict integer partitions of n with pairwise indivisible parts.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 2, 2, 1, 3, 2, 4, 3, 5, 2, 5, 4, 6, 3, 7, 6, 9, 5, 9, 8, 13, 10, 15, 9, 15, 13, 18, 14, 22, 21, 26, 19, 29, 24, 36, 31, 40, 35, 45, 38, 54, 55, 59, 55, 70, 69, 84, 74, 89, 86, 107, 103, 119, 115, 143, 143, 159

Offset: 1

Gus Wiseman, May 17 2018

nonn

Given a finite set S of positive integers greater than one, let G(S) be the simple labeled graph with vertex set S and edges between any two vertices with a common divisor. For example, G({6,14,15,35}) is a 4-cycle. A multiset S is said to be connected if G(S) is a connected graph.

			The a(34) = 13 connected strict integer partitions with pairwise indivisible parts are (34), (18,16), (20,14), (22,12), (24,10), (26,8), (28,6), (30,4), (14,12,8), (15,10,9), (20,8,6), (14,10,6,4), (15,9,6,4). Their corresponding multiset multisystems (see A112798, A302242) are the following.
         (34): {{1,7}}
       (30 4): {{1,2,3},{1,1}}
       (28 6): {{1,1,4},{1,2}}
       (26 8): {{1,6},{1,1,1}}
      (24 10): {{1,1,1,2},{1,3}}
      (22 12): {{1,5},{1,1,2}}
      (20 14): {{1,1,3},{1,4}}
     (20 8 6): {{1,1,3},{1,1,1},{1,2}}
      (18 16): {{1,2,2},{1,1,1,1}}
    (15 10 9): {{2,3},{1,3},{2,2}}
   (15 9 6 4): {{2,3},{2,2},{1,2},{1,1}}
    (14 12 8): {{1,4},{1,1,2},{1,1,1}}
  (14 10 6 4): {{1,4},{1,3},{1,2},{1,1}}

Cf. A000009, A003963, A006126, A048143, A051424, A054921, A076078, A259936, A281116, A285572, A285573, A286518, A286520, A293993, A302242, A303362, A304714, A304716.

zsm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[Less@@#,GCD@@s[[#]]]>1&]},If[c==={},s,zsm[Union[Append[Delete[s,List/@c[[1]]],LCM@@s[[c[[1]]]]]]]]];
Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&Length[zsm[#]]===1&&Select[Tuples[#,2],UnsameQ@@#&&Divisible@@#&]==={}&]],{n,30}]

A305150 Number of factorizations of n into distinct, pairwise indivisible factors greater than 1.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, May 26 2018

nonn

			The a(60) = 6 factorizations are (3 * 4 * 5), (3 * 20), (4 * 15), (5 * 12), (6 * 10), (60). Missing from this list are (2 * 3 * 10), (2 * 5 * 6), (2 * 30).

Antti Karttunen, Table of n, a(n) for n = 1..16384
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..100000

Cf. A001055, A001970, A007716, A045778, A162247, A259936, A275024, A285572, A281113, A281116, A303386, A303431, A305001, A305148, A305149, A305253.

facs[n_] := If[n <= 1, {{}}, Join@@Table[Map[Prepend[#, d] &, Select[facs[n/d], Min@@ # >= d &]], {d, Rest[Divisors[n]]}]]; Table[Length[Select[facs[n], UnsameQ@@ # && Select[Tuples[Union[#], 2], UnsameQ@@ # && Divisible@@ # &] == {} &]], {n, 100}]

A305150(n, m=n, facs=List([])) = if(1==n, 1, my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m)&&factorback(apply(x -> (x%d),Vec(facs))), newfacs = List(facs); listput(newfacs,d); s += A305150(n/d, d-1, newfacs))); (s)); \\ Antti Karttunen, Dec 06 2018

a(n) <= A045778(n) <= A001055(n). - Antti Karttunen, Dec 06 2018

More terms from Antti Karttunen, Dec 06 2018

A317624 Number of integer partitions of n where all parts are > 1 and whose LCM is n.

Original entry on oeis.org

0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 3, 1, 5, 1, 1, 1, 17, 1, 1, 1, 7, 1, 60, 1, 1, 1, 1, 1, 76, 1, 1, 1, 55, 1, 105, 1, 11, 10, 1, 1, 187, 1, 6, 1, 13, 1, 30, 1, 111, 1, 1, 1, 5043, 1, 1, 15, 1, 1, 230, 1, 17, 1, 242, 1, 4173, 1, 1, 12, 19, 1

Offset: 0

Gus Wiseman, Aug 01 2018

nonn

			The a(20) = 5 partitions are (20), (10,4,4,2), (10,4,2,2,2), (5,5,4,4,2), (5,5,4,2,2,2).
The a(45) = 10 partitions:
  (45),
  (15,15,9,3,3), (15,9,9,9,3),
  (15,9,9,3,3,3,3), (15,9,5,5,5,3,3), (9,9,9,5,5,5,3),
  (15,9,3,3,3,3,3,3,3), (9,9,5,5,5,3,3,3,3), (9,5,5,5,5,5,5,3,3),
  (9,5,5,5,3,3,3,3,3,3,3).
From _David A. Corneth_, Sep 08 2018: (Start)
Let sum(t) denote the sum of elements of a tuple t. The tuples t with distinct divisors of 45 that have lcm(t) = 45 and sum(t) <= 45 are {(45) and (3, 9, 15), (3, 5, 9, 15), (3, 5, 9), (5, 9), (9, 15), (5, 9, 15)}. For each such tuple t, find the number of partitions of 45 - s(t) into distinct parts of t.
For the tuple (45), there is 1 partition of 45 - 45 = 0 into parts with 45. That is: {()}.
For the tuple (3, 9, 15), there are 4 partitions of 45 - (3 + 9 + 15) = 18 into parts with 3, 9 and 15. They are {(3, 15), (9, 9), (3, 3, 3, 9), (3, 3, 3, 3, 3, 3)}.
For the tuple (3, 5, 9), there are 4 partitions of 45 - (3 + 5 + 9) = 28 into parts with 3, 5 and 9; they are {(5, 5, 9, 9), (3, 3, 3, 5, 5, 9), (3, 5, 5, 5, 5, 5), (3, 3, 3, 3, 3, 3, 5, 5)}.
For the tuple (3, 5, 9, 15), there is 1 partition of 45 - (3 + 5 + 9 + 15) = 13 into parts with 3, 5, 9 and 15. That is (3, 5, 5).
The other tuples, (5, 9), (9, 15), and (5, 9, 15); they give no extra tuples. That's because there is no solution to the Diophantine equation for 5x + 9y = 45 - (5 + 9), corresponding to the tuple (5, 9) with nonnegative x, y.
That also excludes (9, 15); if there is a solution for that, there would also be a solution for (5, 9). This could whittle down the number of seeds even further. Similarly, (5, 9, 15) gives no solution.
Therefore a(45) = 1 + 4 + 4 + 1 = 10.
(End)
In general, there are A318670(n) (<= A069626(n)) such seed sets of divisors where to start extending the partition from. (See the second PARI program which uses subroutine toplevel_starting_sets.) - _Antti Karttunen_, Sep 08 2018

Antti Karttunen, Table of n, a(n) for n = 0..359
Index entries for sequences related to lcm's

Cf. A000837, A018818, A066874, A067538, A069626, A074761, A074971, A143773, A259936, A281116, A285572, A290103, A305566, A316429, A316431, A316432, A316433, A318670.

Table[Length[Select[IntegerPartitions[n],And[Min@@#>=2,LCM@@#==n]&]],{n,30}]

strong_divisors_reversed(n) = vecsort(select(x -> (x>1), divisors(n)), , 4);
partitions_into_lcm(orgn,n,parts,from=1,m=1) = if(!n,(m==orgn),my(k = #parts, s=0); for(i=from,k,if(parts[i]<=n, s += partitions_into_lcm(orgn,n-parts[i],parts,i,lcm(m,parts[i])))); (s));
A317624(n) = if(n<=1,0,partitions_into_lcm(n,n,strong_divisors_reversed(n))); \\ Antti Karttunen, Sep 07 2018

strong_divisors_reversed(n) = vecsort(select(x -> (x>1), divisors(n)), , 4);
partitions_into(n,parts,from=1) = if(!n,1, if(#parts==from, (0==(n%parts[from])), my(s=0); for(i=from,#parts,if(parts[i]<=n, s += partitions_into(n-parts[i],parts,i))); (s)));
toplevel_starting_sets(orgn,n,parts,from=1,ss=List([])) = { my(k = #parts, s=0, newss); if(lcm(Vec(ss))==orgn,s += partitions_into(n,ss)); for(i=from,k,if(parts[i]<=n, newss = List(ss); listput(newss,parts[i]); s += toplevel_starting_sets(orgn,n-parts[i],parts,i+1,newss))); (s) };
A317624(n) = if(n<=1,0,toplevel_starting_sets(n,n,strong_divisors_reversed(n))); \\ Antti Karttunen, Sep 08-10 2018

A327517 Number of factorizations of n that are empty or have at least two factors, all of which are > 1 and pairwise coprime.

Original entry on oeis.org

1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 4, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 0, 4, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 4, 0, 1, 1, 0, 1, 4, 0, 1, 1, 4, 0, 1, 0, 1, 1, 1, 1, 4, 0, 1, 0, 1, 0, 4, 1, 1, 1

Offset: 1

Gus Wiseman, Sep 19 2019

nonn

Gus Wiseman, Sequences counting and encoding certain classes of multisets

See link for additional cross-references.

Cf. A302569, A302696, A304711, A305079, A327076, A327392.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[Select[facs[n],#=={}||CoprimeQ@@#&]],{n,100}]

a(n > 1) = A259936(n) - 1 = A000110(A001221(n)) - 1.

A284118 Sum of nonprime squarefree divisors of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 7, 1, 1, 1, 11, 1, 7, 1, 15, 16, 1, 1, 7, 1, 11, 22, 23, 1, 7, 1, 27, 1, 15, 1, 62, 1, 1, 34, 35, 36, 7, 1, 39, 40, 11, 1, 84, 1, 23, 16, 47, 1, 7, 1, 11, 52, 27, 1, 7, 56, 15, 58, 59, 1, 62, 1, 63, 22, 1, 66, 128, 1, 35, 70, 130, 1, 7, 1, 75, 16, 39, 78, 150, 1, 11, 1, 83, 1, 84, 86, 87, 88, 23, 1, 62

Offset: 1

Ilya Gutkovskiy, Mar 20 2017

nonn

			a(30) = 62 because 30 has 8 divisors {1, 2, 3, 5, 6, 10, 15, 30} among which 5 are nonprime squarefree {1, 6, 10, 15, 30} therefore 1 + 6 + 10 + 15 + 30 = 62.

G. C. Greubel, Table of n, a(n) for n = 1..5000
Index entries for sequences related to sums of divisors.

Cf. A000469, A001221, A008472, A008683, A048250, A259936.

nmax = 90; Rest[CoefficientList[Series[x/(1 - x) + Sum[Sign[PrimeNu[k] - 1] MoebiusMu[k]^2 k x^k/(1 - x^k), {k, 2, nmax}], {x, 0, nmax}], x]]
Table[Total[Select[Divisors[n], #1 == 1 || (SquareFreeQ[#1] && PrimeNu[#1] > 1) &]], {n, 90}]

Vec((x/(1 - x)) + sum(k=2, 90, sign(omega(k) - 1) * moebius(k)^2 * k * x^k/(1 - x^k)) + O(x^91)) \\ Indranil Ghosh, Mar 21 2017

from sympy import divisors
from sympy.ntheory.factor_ import core
def a(n): return sum([i for i in divisors(n) if core(i)==i and isprime(i)==0]) # Indranil Ghosh, Mar 21 2017

G.f.: x/(1 - x) + Sum_{k>=2} sgn(omega(k)-1)*mu(k)^2*k*x^k/(1 - x^k), where omega(k) is the number of distinct primes dividing k (A001221) and mu(k) is the Moebius function (A008683).
a(n) = Sum_{d|n, d nonprime squarefree} d.
a(n) = 1 if n is a prime power.
a(n) = A048250(n) - A008472(n). - Amiram Eldar, Jun 05 2025

A354870 Number of nonprime squarefree divisors of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 5, 1, 1, 2, 2, 2, 2, 1, 2, 2, 2, 1, 5, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 5, 1, 2, 2, 1, 2, 5, 1, 2, 2, 5, 1, 2, 1, 2, 2, 2, 2, 5, 1, 2, 1, 2, 1, 5, 2, 2, 2, 2, 1, 5, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 1, 5, 1, 2, 5

Offset: 1

Antti Karttunen, Jun 11 2022

nonn

Number of terms of A000469 that divide n.

Antti Karttunen, Table of n, a(n) for n = 1..65537
Index entries for sequences computed from exponents in factorization of n

Inverse Möbius transform of A354819.
Cf. A000469, A046523.
Differs from A259936 for the first time at n=210, where a(210) = 12, while A259936(210) = 15.
Cf. A034444, A001221.

with(numtheory): seq(2^nops(factorset(n)) - nops(factorset(n)), n=1..80); # Ridouane Oudra, Dec 07 2023

a[n_] := DivisorSum[n, 1 &, ! PrimeQ[#] && SquareFreeQ[#] &]; Array[a, 100] (* Amiram Eldar, Jun 11 2022 *)

A354819(n) = ((1!=bigomega(n))&&issquarefree(n));
A354870(n) = sumdiv(n,d,A354819(d));

a(n) = Sum_{d|n} A354819(d).
For all n >= 1, a(n) = a(A046523(n)).
a(n) = A034444(n) - A001221(n). - Ridouane Oudra, Dec 07 2023

cp's OEIS Frontend

A316441 a(n) = Sum (-1)^k where the sum is over all factorizations of n into factors > 1 and k is the number of factors.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A294617 Number of ways to choose a set partition of a strict integer partition of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A305149 Number of factorizations of n whose distinct factors are pairwise indivisible and greater than 1.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A340608 The number of prime factors of n (A001222) is relatively prime to the maximum prime index of n (A061395).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A304717 Number of connected strict integer partitions of n with pairwise indivisible parts.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A305150 Number of factorizations of n into distinct, pairwise indivisible factors greater than 1.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A317624 Number of integer partitions of n where all parts are > 1 and whose LCM is n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

A327517 Number of factorizations of n that are empty or have at least two factors, all of which are > 1 and pairwise coprime.