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.

Previous Showing 21-30 of 72 results. Next

A225244 Number of partitions of n into squarefree divisors of n.

Original entry on oeis.org

1, 1, 2, 2, 3, 2, 8, 2, 5, 4, 11, 2, 27, 2, 14, 14, 9, 2, 64, 2, 40, 18, 20, 2, 125, 6, 23, 10, 53, 2, 742, 2, 17, 26, 29, 26, 343, 2, 32, 30, 195, 2, 1654, 2, 79, 136, 38, 2, 729, 8, 341, 38, 92, 2, 1000, 38, 265, 42, 47, 2, 14188, 2, 50, 184, 33, 44, 5257, 2
Offset: 0

Views

Author

Reinhard Zumkeller, May 05 2013

Keywords

Comments

a(n) <= A018818(n);
a(n) = A018818(n) iff n is squarefree: a(A005117(n)) = A018818(A005117(n));
a(A000040(n)) = 2.

Examples

			a(8) = #{2+2+2+2, 2+2+2+1+1, 2+2+1+1+1+1, 2+6x1, 8x1} = 5;
a(9) = #{3+3+3, 3+3+1+1+1, 3+1+1+1+1+1+1, 9x1} = 4;
a(10) = #{10, 5+5, 5+2+2+1, 5+2+1+1+1, 5+5x1, 2+2+2+2+2, 2+2+2+2+1+1, 2+2+2+1+1+1+1, 2+2+6x1, 2+8x1, 10x1} = 11;
a(11) = #{11, 1+1+1+1+1+1+1+1+1+1+1} = 2;
a(12) = #{6+6, 6+3+3, 6+3+2+1, 6+3+1+1+1, 6+2+2+2, 6+2+2+1+1, 6+2+1+1+1+1, 6+6x1, 3+3+3+3, 3+3+3+2+1, 3+3+3+1+1+1, 3+3+2+2+2, 3+3+2+2+1+1, 3+3+2+4x1, 3+3+6x1, 3+2+2+2+2+1, 3+2+2+2+1+1+1, 3+2+2+5x1, 3+2+7x1, 3+8x1, 2+2+2+2+2+2, 2+2+2+2+2+1+1, 2+2+2+2+1+1+1+1, 2+2+2+6x1, 2+2+8x1, 2+10x1, 12x1} = 27;
a(13) = #{11, 1+1+1+1+1+1+1+1+1+1+1+1+1} = 2;
a(14) = #{14, 7+7, 7+2+2+2+1, 7+2+2+1+1+1, 7+2+5x1, 7+7x1, 7x2, 6x2+1+1, 5x2+1+1+1+1, 4x2+6x1, 2+2+2+8x1, 2+2+10x1, 2+12x1, 14x1} = 14;
a(15) = #{15, 5+5+5, 5+5+3+1+1, 5+5+5x1, 5+3+3+3+1, 5+3+3+1+1+1+1, 5+3+7x1, 5+10x1, 3+3+3+3+3, 3+3+3+3+1+1+1, 3+3+3+6x1, 3+3+9x1, 3+12x1, 15x1} = 14.
		

Crossrefs

Programs

  • Haskell
    a225244 n = p (a206778_row n) n where
       p _          0 = 1
       p []         _ = 0
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
  • Maple
    with(numtheory):
    a:= proc(n) local b, l; l:= sort([select(issqrfree, divisors(n))[]]):
          b:= proc(m, i) option remember; `if`(m=0 or i=1, 1,
                `if`(i<1, 0, b(m, i-1)+`if`(l[i]>m, 0, b(m-l[i], i))))
              end; forget(b):
          b(n, nops(l))
        end:
    seq(a(n), n=0..100); # Alois P. Heinz, Feb 05 2014
  • Mathematica
    a[0] = 1; a[n_] := Module[{b, l}, l = Select[Divisors[n], SquareFreeQ]; b[m_, i_] := b[m, i] = If[m == 0 || i == 1, 1, If[i < 1, 0, b[m, i - 1] + If[l[[i]] > m, 0, b[m - l[[i]], i]]]]; b[n, Length[l]]]; Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Oct 27 2015, after Alois P. Heinz *)

Formula

a(n) = [x^n] Product_{d|n, mu(d) != 0} 1/(1 - x^d), where mu() is the Moebius function (A008683). - Ilya Gutkovskiy, Jul 26 2017

A320778 Inverse Euler transform of the Euler totient function phi = A000010.

Original entry on oeis.org

1, 1, 0, 1, 0, 2, -3, 4, -4, 4, -9, 14, -19, 30, -42, 50, -76, 128, -194, 286, -412, 598, -909, 1386, -2100, 3178, -4763, 7122, -10758, 16414, -25061, 38056, -57643, 87568, -133436, 203618, -311128, 475536, -726355, 1109718, -1697766, 2601166, -3987903, 6114666
Offset: 0

Views

Author

Gus Wiseman, Oct 22 2018

Keywords

Comments

The Euler transform of a sequence q is the sequence of coefficients of x^n, n > 0, in the expansion of Product_{n > 0} 1/(1 - x^n)^q(n). The constant term 1 is sometimes taken to be the zeroth part of the Euler transform.

Crossrefs

Number theoretical functions: A000005, A000010, A000203, A001055, A001221, A001222, A008683, A010054.
Inverse Euler transforms: A059966, A320767, A320776, A320777, A320779, A320780, A320781, A320782.

Programs

  • Maple
    # The function EulerInvTransform is defined in A358451.
    a := EulerInvTransform(n -> ifelse(n=0, 1, NumberTheory:-Totient(n))):
    seq(a(n), n = 0..43); # Peter Luschny, Nov 21 2022
  • Mathematica
    EulerInvTransform[{}]={};EulerInvTransform[seq_]:=Module[{final={}},For[i=1,i<=Length[seq],i++,AppendTo[final,i*seq[[i]]-Sum[final[[d]]*seq[[i-d]],{d,i-1}]]];
    Table[Sum[MoebiusMu[i/d]*final[[d]],{d,Divisors[i]}]/i,{i,Length[seq]}]];
    EulerInvTransform[Array[EulerPhi,30]]

A368422 Number of non-isomorphic set multipartitions of weight n satisfying a strict version of the axiom of choice.

Original entry on oeis.org

1, 1, 2, 4, 9, 18, 43, 95, 233, 569
Offset: 0

Views

Author

Gus Wiseman, Dec 26 2023

Keywords

Comments

A set multipartition is a finite multiset of finite nonempty sets. The weight of a set multipartition is the sum of cardinalities of its elements. Weight is generally not the same as number of vertices.
The axiom of choice says that, given any sequence of nonempty sets, it is possible to choose a sequence containing an element from each. In the strict version, the elements of this sequence must be distinct, meaning none is chosen more than once.

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(5) = 18 set multipartitions:
  {{1}}  {{1,2}}    {{1,2,3}}      {{1,2,3,4}}        {{1,2,3,4,5}}
         {{1},{2}}  {{1},{2,3}}    {{1,2},{1,2}}      {{1},{2,3,4,5}}
                    {{2},{1,2}}    {{1},{2,3,4}}      {{1,2},{3,4,5}}
                    {{1},{2},{3}}  {{1,2},{3,4}}      {{1,4},{2,3,4}}
                                   {{1,3},{2,3}}      {{2,3},{1,2,3}}
                                   {{3},{1,2,3}}      {{4},{1,2,3,4}}
                                   {{1},{2},{3,4}}    {{1},{2,3},{2,3}}
                                   {{1},{3},{2,3}}    {{1},{2},{3,4,5}}
                                   {{1},{2},{3},{4}}  {{1},{2,3},{4,5}}
                                                      {{1},{2,4},{3,4}}
                                                      {{1},{4},{2,3,4}}
                                                      {{2},{1,3},{2,3}}
                                                      {{2},{3},{1,2,3}}
                                                      {{3},{1,3},{2,3}}
                                                      {{4},{1,2},{3,4}}
                                                      {{1},{2},{3},{4,5}}
                                                      {{1},{2},{4},{3,4}}
                                                      {{1},{2},{3},{4},{5}}
		

Crossrefs

The case of unlabeled graphs is A134964, complement A140637.
Set multipartitions have ranks A302478, cf. A073576.
The case of labeled graphs is A133686, complement A367867.
The complement without repeats is A368094 connected A368409.
Without repeats we have A368095, connected A368410.
The complement allowing repeats is A368097, ranks A355529.
Allowing repeated elements gives A368098, ranks A368100.
Factorizations of this type are counted by A368414, complement A368413.
The complement is counted by A368421.
A000110 counts set partitions, non-isomorphic A000041.
A003465 counts covering set-systems, unlabeled A055621.
A007716 counts non-isomorphic multiset partitions, connected A007718.
A058891 counts set-systems, unlabeled A000612, connected A323818.
A283877 counts non-isomorphic set-systems, connected A300913.

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mpm[n_]:=Join@@Table[Union[Sort[Sort /@ (#/.x_Integer:>s[[x]])]&/@sps[Range[n]]],{s,Flatten[MapIndexed[Table[#2,{#1}]&,#]]& /@ IntegerPartitions[n]}];
    brute[m_]:=First[Sort[Table[Sort[Sort /@ (m/.Rule@@@Table[{i,p[[i]]},{i,Length[p]}])], {p,Permutations[Union@@m]}]]];
    Table[Length[Union[brute /@ Select[mpm[n],And@@UnsameQ@@@#&&Select[Tuples[#], UnsameQ@@#&]!={}&]]],{n,0,6}]

A379310 Number of nonsquarefree prime indices of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 2, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Dec 27 2024

Keywords

Comments

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.

Examples

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

Crossrefs

Positions of first appearances are A000420.
Positions of zero are A302478, counted by A073576 (strict A087188).
No squarefree parts: A379307, counted by A114374 (strict A256012).
One squarefree part: 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:
- A330944 nonprime, see A000586, A000607, A076610, A330945.
- A379311 old prime, see A204389, A320629, A379312-A379315.

Programs

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

Formula

Totally additive with a(prime(k)) = A107078(k) = 1 - A008966(k).

A320767 Inverse Euler transform applied once to {1,-1,0,0,0,...}, twice to {1,0,0,0,0,...}, or three times to {1,1,1,1,1,...}.

Original entry on oeis.org

1, 1, -2, 1, -1, 2, -3, 4, -5, 8, -13, 18, -25, 40, -62, 90, -135, 210, -324, 492, -750, 1164, -1809, 2786, -4305, 6710, -10460, 16264, -25350, 39650, -62057, 97108, -152145, 238818, -375165, 589520, -927200, 1459960, -2300346, 3626200, -5720274, 9030450
Offset: 0

Views

Author

Gus Wiseman, Oct 20 2018

Keywords

Comments

The Euler transform of a sequence q is the sequence of coefficients of x^n, n > 0, in the expansion of Product_{n > 0} 1/(1 - x^n)^q(n). The constant term 1 is sometimes taken to be the zeroth part of the Euler transform.

Crossrefs

Programs

  • Mathematica
    EulerInvTransform[{}]={};EulerInvTransform[seq_]:=Module[{final={}},For[i=1,i<=Length[seq],i++,AppendTo[final,i*seq[[i]]-Sum[final[[d]]*seq[[i-d]],{d,i-1}]]];
    Table[Sum[MoebiusMu[i/d]*final[[d]],{d,Divisors[i]}]/i,{i,Length[seq]}]];
    Nest[EulerInvTransform,Array[DiscreteDelta,50,0],2]

A368421 Number of non-isomorphic set multipartitions of weight n contradicting a strict version of the axiom of choice.

Original entry on oeis.org

0, 0, 1, 2, 7, 16, 47, 116, 325, 861
Offset: 0

Views

Author

Gus Wiseman, Dec 26 2023

Keywords

Comments

A set multipartition is a finite multiset of finite nonempty sets. The weight of a set multipartition is the sum of cardinalities of its elements. Weight is generally not the same as number of vertices.
The axiom of choice says that, given any sequence of nonempty sets Y, it is possible to choose a sequence containing an element from each. In the strict version, the elements of this sequence must be distinct, meaning none is chosen more than once.

Examples

			Non-isomorphic representatives of the a(2) = 1 through a(5) = 16 set multipartitions:
  {{1},{1}}  {{1},{1},{1}}  {{1},{1},{2,3}}    {{1},{1},{2,3,4}}
             {{1},{2},{2}}  {{1},{2},{1,2}}    {{2},{1,2},{1,2}}
                            {{2},{2},{1,2}}    {{3},{3},{1,2,3}}
                            {{1},{1},{1},{1}}  {{1},{1},{1},{2,3}}
                            {{1},{1},{2},{2}}  {{1},{1},{3},{2,3}}
                            {{1},{2},{2},{2}}  {{1},{2},{2},{1,2}}
                            {{1},{2},{3},{3}}  {{1},{2},{2},{3,4}}
                                               {{1},{2},{3},{2,3}}
                                               {{1},{3},{3},{2,3}}
                                               {{2},{2},{2},{1,2}}
                                               {{1},{1},{1},{1},{1}}
                                               {{1},{1},{2},{2},{2}}
                                               {{1},{2},{2},{2},{2}}
                                               {{1},{2},{2},{3},{3}}
                                               {{1},{2},{3},{3},{3}}
                                               {{1},{2},{3},{4},{4}}
		

Crossrefs

The case of unlabeled graphs is A140637, complement A134964.
Set multipartitions have ranks A302478, cf. A073576.
The case of labeled graphs is A367867, complement A133686.
With distinct edges we have A368094 connected A368409.
The complement with distinct edges is A368095, connected A368410.
Allowing repeated elements gives A368097, ranks A355529.
The complement allowing repeats is A368098, ranks A368100.
Factorizations of this type are counted by A368413, complement A368414.
The complement is counted by A368422.
A000110 counts set partitions, non-isomorphic A000041.
A003465 counts covering set-systems, unlabeled A055621.
A007716 counts non-isomorphic multiset partitions, connected A007718.
A058891 counts set-systems, unlabeled A000612, connected A323818.
A283877 counts non-isomorphic set-systems, connected A300913.

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mpm[n_]:=Join@@Table[Union[Sort[Sort /@ (#/.x_Integer:>s[[x]])]&/@sps[Range[n]]],{s,Flatten[MapIndexed[Table[#2,{#1}]&,#]]& /@ IntegerPartitions[n]}];
    brute[m_]:=First[Sort[Table[Sort[Sort /@ (m/.Rule@@@Table[{i,p[[i]]},{i,Length[p]}])], {p,Permutations[Union@@m]}]]];
    Table[Length[Union[brute /@ Select[mpm[n],And@@UnsameQ@@@#&&Select[Tuples[#], UnsameQ@@#&]=={}&]]],{n,0,6}]

A377041 First term of the n-th differences of the squarefree numbers. Inverse zero-based binomial transform of A005117.

Original entry on oeis.org

1, 1, 0, 1, -3, 6, -8, 3, 22, -92, 252, -578, 1189, -2255, 3991, -6617, 10245, -14626, 18666, -19635, 12104, 13090, -69122, 171478, -332718, 552138, -798629, 982514, -901485, 116219, 2351842, -8715135, 23856206, -57926011, 130281064, -273804584, 535390333
Offset: 0

Views

Author

Gus Wiseman, Oct 18 2024

Keywords

Crossrefs

The version for primes is A007442, noncomposites A030016, composites A377036.
This is the first column of A377038.
For nonsquarefree numbers we have A377049.
For prime-powers we have A377054.
A000040 lists the primes, differences A001223, seconds A036263.
A005117 lists the squarefree numbers, complement A013929 (differences A078147).
A073576 counts integer partitions into squarefree numbers, factorizations A050320.
A377042 gives first position of 0 in each row of A377038.

Programs

  • Mathematica
    q=Select[Range[100],SquareFreeQ];
    t=Table[Sum[(-1)^(j-k)*Binomial[j,k]*q[[1+k]],{k,0,j}],{j,0,Length[q]/2}]

Formula

The inverse zero-based binomial transform of a sequence (q(0), q(1), q(2), ...) is the sequence p given by:
p(j) = sum_{k=0..j} (-1)^(j-k) binomial(j,k) q(k)

A322526 Number of integer partitions of n whose product of parts is a squarefree number.

Original entry on oeis.org

1, 1, 2, 3, 3, 5, 6, 8, 9, 10, 13, 15, 17, 21, 24, 27, 30, 36, 41, 46, 51, 57, 65, 73, 82, 90, 101, 109, 121, 134, 150, 164, 177, 193, 214, 232, 253, 278, 300, 324, 351, 386, 419, 452, 484, 521, 563, 610, 658, 706, 758, 809, 868, 938, 1006, 1071, 1140, 1220, 1307
Offset: 0

Views

Author

Gus Wiseman, Dec 14 2018

Keywords

Comments

The parts of such a partition must also be squarefree and distinct except for any number of 1's.

Examples

			The a(8) = 9 partitions are (53), (71), (521), (611), (5111), (32111), (311111), (2111111), (11111111). Missing from this list are (8), (62), (44), (431), (422), (4211), (41111), (332), (3311), (3221), (2222), (22211), (221111).
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],SquareFreeQ[Times@@#]&]],{n,30}]

A376655 Sorted positions of first appearances in the second differences of consecutive squarefree numbers (A005117).

Original entry on oeis.org

1, 2, 3, 5, 6, 30, 61, 150, 514, 1025, 5153, 13390, 13391, 131964, 502651, 664312, 4387185, 5392318, 20613826
Offset: 1

Views

Author

Gus Wiseman, Oct 07 2024

Keywords

Comments

Warning: Do not confuse with A246655 (prime-powers exclusive).

Examples

			The squarefree numbers (A005117) are:
  1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 26, 29, 30, 31, 33, ...
with first differences (A076259):
  1, 1, 2, 1, 1, 3, 1, 2, 1, 1, 2, 2, 2, 1, 1, 3, 3, 1, 1, 2, 1, 1, 2, 1, 1, 2, ...
with first differences (A376590):
  0, 1, -1, 0, 2, -2, 1, -1, 0, 1, 0, 0, -1, 0, 2, 0, -2, 0, 1, -1, 0, 1, -1, 0, ...
with sorted first appearances at (A376655):
  1, 2, 3, 5, 6, 30, 61, 150, 514, 1025, 5153, 13390, 13391, ...
		

Crossrefs

For first differences we had A376311 (first appearances in A076259).
These are the sorted positions of first appearances in A376590.
For prime-powers instead of squarefree numbers we have A376653/A376654.
For primes instead of squarefree numbers we have A376656.
A000040 lists the prime numbers, differences A001223.
A005117 lists squarefree numbers, complement A013929 (differences A078147).
A073576 counts integer partitions into squarefree numbers, factorizations A050320.
For second differences: A036263 (prime), A073445 (composite), A376559 (perfect-power), A376562 (non-perfect-power), A376593 (nonsquarefree), A376596 (prime-power inclusive), A376599 (non-prime-power inclusive).
For squarefree: A376591 (inflections and undulations), A376592 (nonzero curvature).

Programs

  • Mathematica
    q=Differences[Select[Range[1000],SquareFreeQ],2];
    Select[Range[Length[q]],!MemberQ[Take[q,#-1],q[[#]]]&]

Extensions

a(14)-a(19) from Chai Wah Wu, Oct 07 2024

A320776 Inverse Euler transform of the number of prime factors (with multiplicity) function A001222.

Original entry on oeis.org

1, 0, 1, 1, 1, 0, -1, -1, 0, 1, 0, -1, -1, -1, 1, 3, 3, -2, -5, -4, 0, 7, 7, 0, -9, -10, 2, 15, 15, -3, -27, -30, 3, 46, 51, 1, -71, -91, -7, 117, 157, 23, -194, -265, -57, 318, 465, 111, -536, -821, -230, 893, 1456, 505, -1485, -2559, -1036, 2433, 4483, 2022
Offset: 0

Views

Author

Gus Wiseman, Oct 22 2018

Keywords

Comments

The Euler transform of a sequence q is the sequence of coefficients of x^n, n > 0, in the expansion of Product_{n > 0} 1/(1 - x^n)^q(n). The constant term 1 is sometimes taken to be the zeroth part of the Euler transform.

Crossrefs

Number theoretical functions: A000005, A000010, A000203, A001055, A001221, A001222, A008683, A010054.
Inverse Euler transforms: A059966, A320767, A320777, A320778, A320779, A320780, A320781, A320782.

Programs

  • Maple
    # The function EulerInvTransform is defined in A358451.
    a := EulerInvTransform(n -> ifelse(n=0, 1, NumberTheory:-NumberOfPrimeFactors(n))):
    seq(a(n), n = 0..59); # Peter Luschny, Nov 21 2022
  • Mathematica
    EulerInvTransform[{}]={};EulerInvTransform[seq_]:=Module[{final={}},For[i=1,i<=Length[seq],i++,AppendTo[final,i*seq[[i]]-Sum[final[[d]]*seq[[i-d]],{d,i-1}]]];
    Table[Sum[MoebiusMu[i/d]*final[[d]],{d,Divisors[i]}]/i,{i,Length[seq]}]];
    EulerInvTransform[Array[PrimeOmega,100]]
Previous Showing 21-30 of 72 results. Next