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 466 results. Next

A367902 Number of sets of nonempty subsets of {1..n} satisfying a strict version of the axiom of choice.

Original entry on oeis.org

1, 2, 7, 61, 1771, 187223, 70038280, 90111497503, 397783376192189
Offset: 0

Views

Author

Gus Wiseman, Dec 05 2023

Keywords

Comments

The axiom of choice says that, given any set of nonempty sets Y, it is possible to choose a set containing an element from each. The strict version requires this set to have the same cardinality as Y, meaning no element is chosen more than once.

Examples

			The a(2) = 7 set-systems:
  {}
  {{1}}
  {{2}}
  {{1,2}}
  {{1},{2}}
  {{1},{1,2}}
  {{2},{1,2}}
		

Crossrefs

The version for simple graphs is A133686, covering A367869.
The version without singletons is A367770.
The complement allowing empty edges is A367901.
The complement is A367903, without singletons A367769, ranks A367907.
For a unique choice we have A367904, ranks A367908.
These set-systems have ranks A367906.
A000372 counts antichains, covering A006126, nonempty A014466.
A003465 counts covering set-systems, unlabeled A055621.
A058891 counts set-systems, unlabeled A000612.
A059201 counts covering T_0 set-systems.
A323818 counts covering connected set-systems.
A326031 gives weight of the set-system with BII-number n.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Subsets[Range[n]]], Select[Tuples[#],UnsameQ@@#&]!={}&]],{n,0,3}]

Formula

a(n) = A370636(2^n-1). - Alois P. Heinz, Mar 09 2024

Extensions

a(6)-a(8) from Christian Sievers, Jul 25 2024

A320655 Number of factorizations of n into semiprimes. Number of multiset partitions of the multiset of prime factors of n, into pairs.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 18 2018

Keywords

Comments

The characteristic function of nonzero terms is A065043. - R. J. Mathar, Jan 18 2021

Examples

			The a(900) = 5 factorizations into semiprimes:
  900 = (4*9*25)
  900 = (4*15*15)
  900 = (6*6*25)
  900 = (6*10*15)
  900 = (9*10*10)
The a(900) = 5 multiset partitions into pairs:
  {{1,1},{2,2},{3,3}}
  {{1,1},{2,3},{2,3}}
  {{1,2},{1,2},{3,3}}
  {{1,2},{1,3},{2,3}}
  {{2,2},{1,3},{1,3}}
		

Crossrefs

Programs

  • Mathematica
    semfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[semfacs[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],PrimeOmega[#]==2&]}]];
    Table[Length[semfacs[n]],{n,100}]
  • PARI
    A320655(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((2==bigomega(d)&&(d<=m)), s += A320655(n/d, d))); (s)); \\ Antti Karttunen, Dec 06 2020

Extensions

Data section extended up to 105 terms by Antti Karttunen, Dec 06 2020

A120733 Number of matrices with nonnegative integer entries and without zero rows or columns such that sum of all entries is equal to n.

Original entry on oeis.org

1, 1, 5, 33, 281, 2961, 37277, 546193, 9132865, 171634161, 3581539973, 82171451025, 2055919433081, 55710251353953, 1625385528173693, 50800411296363617, 1693351638586070209, 59966271207156833313, 2248276994650395873861, 88969158875611127548481
Offset: 0

Views

Author

Vladeta Jovovic, Aug 18 2006, Aug 21 2006

Keywords

Comments

The number of such matrices up to rows/columns permutations are given in A007716.
Dimensions of the graded components of the Hopf algebra MQSym (Matrix quasi-symmetric functions). - Jean-Yves Thibon (jyt(AT)univ-mlv.fr), Oct 23 2006
From Kyle Petersen, Aug 10 2016: (Start)
Number of cells in the two-sided Coxeter complex of the symmetric group. Inclusion of faces corresponds to refinement of matrices, see Section 6 of Petersen paper. The number of cells in the type B analog is given by A275787.
Also known as "two-way contingency tables" in the Diaconis-Gangolli reference. (End)

Examples

			a(2) = 5:
[1 0]   [0 1]   [1]   [1 1]   [2]
[0 1]   [1 0]   [1]
From _Gus Wiseman_, Nov 14 2018: (Start)
The a(3) = 33 matrices:
  [3][21][12][111]
.
  [2][20][11][11][110][101][1][10][10][100][02][011][01][01][010][001]
  [1][01][10][01][001][010][2][11][02][011][10][100][20][11][101][110]
.
  [1][10][10][10][100][100][01][01][010][01][010][001][001]
  [1][10][01][01][010][001][10][10][100][01][001][100][010]
  [1][01][10][01][001][010][10][01][001][10][100][010][100]
(End)
		

Crossrefs

Row sums of A261781.

Programs

  • Maple
    t1 := M -> add( add( add( (-1)^(n-j)*binomial(n, j)*((1-x)^(-j)-1)^m, j=0..n), n=0..M), m=0..M); s := series(t1(20),x,20); gfun[seriestolist](%); # N. J. A. Sloane, Jan 14 2009
  • Mathematica
    a[n_] := Sum[2^(-2-r-s)*Binomial[n+r*s-1, n], {r, 0, Infinity}, {s, 0, Infinity}]; Table[Print[an = a[n]]; an, {n, 0, 19}] (* Jean-François Alcover, May 15 2012, after Vladeta Jovovic *)
    Flatten[{1,Table[1/n!*Sum[(-1)^(n-k)*StirlingS1[n,k]*Sum[m!*StirlingS2[k, m],{m,k}]^2,{k,n}],{n,20}]}] (* Vaclav Kotesovec, May 07 2014 *)
    multsubs[set_,k_]:=If[k==0,{{}},Join@@Table[Prepend[#,set[[i]]]&/@multsubs[Drop[set,i-1],k-1],{i,Length[set]}]]; Table[Length[Select[multsubs[Tuples[Range[n],2],n],And[Union[First/@#]==Range[Max@@First/@#],Union[Last/@#]==Range[Max@@Last/@#]]&]],{n,5}] (* Gus Wiseman, Nov 14 2018 *)

Formula

a(n) = (1/n!)*Sum_{k=0..n} (-1)^(n-k)*Stirling1(n,k)*A000670(k)^2.
G.f.: Sum_{m>=0,n>=0} Sum_{j=0..n} (-1)^(n-j)*C(n,j)*((1-x)^(-j)-1)^m.
a(n) = Sum_{r>=0,s>=0} binomial(r*s+n-1,n)/2^(r+s+2).
G.f.: Sum_{n>=0} 1/(2-(1-x)^(-n))/2^(n+1). - Vladeta Jovovic, Oct 30 2006
a(n) ~ 2^(log(2)/2-2) * n! / (log(2))^(2*n+2). - Vaclav Kotesovec, May 07 2014

Extensions

More terms from N. J. A. Sloane, Jan 14 2009

A320656 Number of factorizations of n into squarefree semiprimes. Number of multiset partitions of the multiset of prime factors of n, into strict pairs.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 18 2018

Keywords

Examples

			The a(4620) = 6 factorizations into squarefree semiprimes:
  4620 = (6*10*77)
  4620 = (6*14*55)
  4620 = (6*22*35)
  4620 = (10*14*33)
  4620 = (10*21*22)
  4620 = (14*15*22)
The a(4620) = 6 multiset partitions into strict pairs:
  {{1,2},{1,3},{4,5}}
  {{1,2},{1,4},{3,5}}
  {{1,2},{1,5},{3,4}}
  {{1,3},{1,4},{2,5}}
  {{1,3},{2,4},{1,5}}
  {{1,4},{2,3},{1,5}}
The a(69300) = 10 factorizations into squarefree semiprimes:
  69300 = (6*6*35*55)
  69300 = (6*10*15*77)
  69300 = (6*10*21*55)
  69300 = (6*10*33*35)
  69300 = (6*14*15*55)
  69300 = (6*15*22*35)
  69300 = (10*10*21*33)
  69300 = (10*14*15*33)
  69300 = (10*15*21*22)
  69300 = (14*15*15*22)
The a(69300) = 10 multiset partitions into strict pairs:
  {{1,2},{1,2},{3,4},{3,5}}
  {{1,2},{1,3},{2,3},{4,5}}
  {{1,2},{1,3},{2,4},{3,5}}
  {{1,2},{1,3},{2,5},{3,4}}
  {{1,2},{1,4},{2,3},{3,5}}
  {{1,2},{2,3},{1,5},{3,4}}
  {{1,3},{1,3},{2,4},{2,5}}
  {{1,3},{1,4},{2,3},{2,5}}
  {{1,3},{2,3},{2,4},{1,5}}
  {{1,4},{2,3},{2,3},{1,5}}.
The a(210) = 3 factorizations into squarefree semiprimes: 210 = (6*35) = (10*21) = (14*15). - _Antti Karttunen_, Nov 02 2022
		

Crossrefs

Programs

  • Mathematica
    bepfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[bepfacs[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],SquareFreeQ[#]&&PrimeOmega[#]==2&]}]];
    Table[Length[bepfacs[n]],{n,100}]
  • PARI
    A320656(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m)&&issquarefree(d)&&2==bigomega(d), s += A320656(n/d, d))); (s)); \\ Antti Karttunen, Nov 02 2022

Formula

a(A002110(n)) = A123023(n). - Antti Karttunen, Nov 02 2022

Extensions

Data section extended up to a(120) and the secondary offset added by Antti Karttunen, Nov 02 2022

A306017 Number of non-isomorphic multiset partitions of weight n in which all parts have the same size.

Original entry on oeis.org

1, 1, 4, 6, 17, 14, 66, 30, 189, 222, 550, 112, 4696, 202, 5612, 30914, 63219, 594, 453125, 980, 3602695, 5914580, 1169348, 2510, 299083307, 232988061, 23248212, 2669116433, 14829762423, 9130, 170677509317, 13684, 1724710753084, 2199418340875, 14184712185, 38316098104262
Offset: 0

Views

Author

Gus Wiseman, Jun 17 2018

Keywords

Comments

A multiset partition of weight n is a finite multiset of finite nonempty multisets whose sizes sum to n.
Number of distinct nonnegative integer matrices with all row sums equal and total sum n up to row and column permutations. - Andrew Howroyd, Sep 05 2018
From Gus Wiseman, Oct 11 2018: (Start)
Also the number of non-isomorphic multiset partitions of weight n in which each vertex appears the same number of times. For n = 4, non-isomorphic representatives of these 17 multiset partitions are:
{{1,1,1,1}}
{{1,1,2,2}}
{{1,2,3,4}}
{{1},{1,1,1}}
{{1},{1,2,2}}
{{1},{2,3,4}}
{{1,1},{1,1}}
{{1,1},{2,2}}
{{1,2},{1,2}}
{{1,2},{3,4}}
{{1},{1},{1,1}}
{{1},{1},{2,2}}
{{1},{2},{1,2}}
{{1},{2},{3,4}}
{{1},{1},{1},{1}}
{{1},{1},{2},{2}}
{{1},{2},{3},{4}}
(End)

Examples

			Non-isomorphic representatives of the a(4) = 17 multiset partitions:
  {{1,1,1,1}}
  {{1,1,2,2}}
  {{1,2,2,2}}
  {{1,2,3,3}}
  {{1,2,3,4}}
  {{1,1},{1,1}}
  {{1,1},{2,2}}
  {{1,2},{1,2}}
  {{1,2},{2,2}}
  {{1,2},{3,3}}
  {{1,2},{3,4}}
  {{1,3},{2,3}}
  {{1},{1},{1},{1}}
  {{1},{1},{2},{2}}
  {{1},{2},{2},{2}}
  {{1},{2},{3},{3}}
  {{1},{2},{3},{4}}
		

Crossrefs

Programs

  • Mathematica
    permcount[v_List] := Module[{m = 1, s = 0, k = 0, t}, For[i = 1, i <= Length[v], i++, t = v[[i]]; k = If[i > 1 && t == v[[i - 1]], k + 1, 1]; m *= t*k; s += t]; s!/m];
    K[q_List, t_, k_] := SeriesCoefficient[1/Product[g = GCD[t, q[[j]]]; (1 - x^(q[[j]]/g))^g, {j, 1, Length[q]}], {x, 0, k}];
    RowSumMats[n_, m_, k_] := Module[{s = 0}, Do[s += permcount[q]* SeriesCoefficient[Exp[Sum[K[q, t, k]/t*x^t, {t, 1, n}]], {x, 0, n}], {q, IntegerPartitions[m]}]; s/m!];
    a[n_] := a[n] = If[n==0, 1, If[PrimeQ[n], 2 PartitionsP[n], Sum[ RowSumMats[ n/d, n, d], {d, Divisors[n]}]]];
    Table[Print[n, " ", a[n]]; a[n], {n, 0, 35}] (* Jean-François Alcover, Nov 07 2019, after Andrew Howroyd *)
  • PARI
    \\ See A318951 for RowSumMats.
    a(n)={sumdiv(n,d,RowSumMats(n/d,n,d))} \\ Andrew Howroyd, Sep 05 2018

Formula

For p prime, a(p) = 2*A000041(p).
a(n) = Sum_{d|n} A331485(n/d, d). - Andrew Howroyd, Feb 09 2020

Extensions

Terms a(11) and beyond from Andrew Howroyd, Sep 05 2018

A318284 Number of multiset partitions of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 2, 2, 3, 4, 5, 5, 9, 7, 7, 11, 11, 12, 16, 15, 15, 26, 22, 21, 29, 19, 30, 36, 31, 30, 66, 38, 42, 52, 56, 52, 47, 45, 57, 92, 77, 67, 77, 74, 101, 98, 135, 64, 137, 97, 176, 135, 109, 109, 118, 105, 231, 249, 97, 141, 181, 139, 297, 198, 385, 195, 269
Offset: 1

Views

Author

Gus Wiseman, Aug 23 2018

Keywords

Examples

			The a(12) = 11 multiset partitions of {1,1,2,3}:
  {{1,1,2,3}}
  {{1},{1,2,3}}
  {{2},{1,1,3}}
  {{3},{1,1,2}}
  {{1,1},{2,3}}
  {{1,2},{1,3}}
  {{1},{1},{2,3}}
  {{1},{2},{1,3}}
  {{1},{3},{1,2}}
  {{2},{3},{1,1}}
  {{1},{1},{2},{3}}
		

Crossrefs

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[facs[Times@@Prime/@nrmptn[n]]],{n,60}]
  • PARI
    permcount(v) = {my(m=1,s=0,k=0,t); for(i=1,#v,t=v[i]; k=if(i>1&&t==v[i-1],k+1,1); m*=t*k;s+=t); s!/m}
    sig(n)={my(f=factor(n)); concat(vector(#f~, i, vector(f[i,2], j, primepi(f[i,1]))))}
    count(sig)={my(n=vecsum(sig), A=O(x*x^vecmax(sig)), s=0); forpart(p=n, my(q=1/prod(i=1, #p, 1 - x^p[i] + A)); s+=prod(i=1, #sig, polcoef(q,sig[i]))*permcount(p)); s/n!}
    a(n)={if(n==1, 1, my(s=sig(n)); if(#s==1, numbpart(s[1]), count(sig(n))))} \\ Andrew Howroyd, Dec 10 2018

Formula

a(n) = A001055(A181821(n)).
a(prime(n)^k) = A219727(n,k). - Andrew Howroyd, Dec 10 2018

A300913 Number of non-isomorphic connected set-systems of weight n.

Original entry on oeis.org

1, 1, 1, 2, 4, 7, 18, 37, 96, 239, 658, 1810, 5358, 16057, 50373, 161811, 536964, 1826151, 6380481, 22822280, 83587920, 312954111, 1197178941, 4674642341, 18620255306, 75606404857, 312763294254, 1317356836235, 5646694922172, 24618969819915, 109125629486233, 491554330852608
Offset: 0

Views

Author

Gus Wiseman, Jun 19 2018

Keywords

Comments

The weight of a set-system is the sum of cardinalities of the sets. Weight is generally not the same as number of vertices.

Examples

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

Crossrefs

Programs

Formula

Inverse Euler transform of A283877.

Extensions

a(11)-a(31) from Jean-François Alcover, Nov 07 2019

A305551 Number of partitions of partitions of n where all partitions have the same sum.

Original entry on oeis.org

1, 1, 3, 4, 9, 8, 22, 16, 43, 41, 77, 57, 201, 102, 264, 282, 564, 298, 1175, 491, 1878, 1509, 2611, 1256, 7872, 2421, 7602, 8026, 16304, 4566, 38434, 6843, 48308, 41078, 56582, 28228, 221115, 21638, 146331, 208142, 453017, 44584, 844773, 63262, 1034193, 1296708
Offset: 0

Views

Author

Gus Wiseman, Jun 20 2018

Keywords

Examples

			The a(4) = 9 partitions of partitions where all partitions have the same sum:
(4), (31), (22), (211), (1111),
(2)(2), (2)(11), (11)(11),
(1)(1)(1)(1).
		

Crossrefs

Programs

  • Mathematica
    Table[Sum[Binomial[PartitionsP[n/k]+k-1,k],{k,Divisors[n]}],{n,60}]
  • PARI
    a(n)={if(n<1, n==0, sumdiv(n, d, binomial(numbpart(n/d) + d - 1, d)))} \\ Andrew Howroyd, Jun 22 2018

Formula

a(n) = Sum_{d|n} binomial(A000041(n/d) + d - 1, d).

A275024 Total weight of the n-th twice-prime-factored multiset partition.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Nov 12 2016

Keywords

Comments

A multiset partition is a finite multiset of finite nonempty multisets of positive integers. The n-th twice-prime-factored multiset partition is constructed by factoring n into prime numbers and then factoring each prime index plus 1 into prime numbers. This produces a unique multiset of multisets of prime numbers which can then be normalized (see example) to produce each possible multiset partition as n ranges over all positive integers.

Examples

			The sequence of multiset partitions begins:
(), ((1)), ((2)), ((1)(1)), ((11)), ((1)(2)), ((3)),
((1)(1)(1)), ((2)(2)), ((1)(11)), ((12)), ((1)(1)(2)),
((4)), ((1)(3)), ((2)(11)), ((1)(1)(1)(1)), ((111)),
((1)(2)(2)), ((22)), ((1)(1)(11)), ((2)(3)), ((1)(12)),
((13)), ((1)(1)(1)(2)), ((11)(11)), ((1)(4)), ((2)(2)(2)),
((1)(1)(3)), ((5)), ((1)(2)(11)), ((112)), ((1)(1)(1)(1)(1)),
((2)(12)), ((1)(111)), ((3)(11)), ((1)(1)(2)(2)), ((6)), ...
		

Crossrefs

Programs

  • Mathematica
    Table[Total[Cases[FactorInteger[n],{p_,k_}:>k*PrimeOmega[PrimePi[p]+1]]],{n,1,100}]

Formula

If prime(k) has weight equal to the number of prime factors (counting multiplicity) of k+1, then a(n) is the sum of weights over all prime factors (counting multiplicity) of n.

A302478 Products of prime numbers of squarefree index.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 15, 16, 17, 18, 20, 22, 24, 25, 26, 27, 29, 30, 31, 32, 33, 34, 36, 39, 40, 41, 43, 44, 45, 47, 48, 50, 51, 52, 54, 55, 58, 59, 60, 62, 64, 65, 66, 67, 68, 72, 73, 75, 78, 79, 80, 81, 82, 83, 85, 86, 87, 88, 90, 93, 94
Offset: 1

Views

Author

Gus Wiseman, Apr 08 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n.

Examples

			Entry A302242 describes a correspondence between positive integers and multiset multisystems. In this case it gives the following sequence of set multisystems.
01:  {}
02:  {{}}
03:  {{1}}
04:  {{},{}}
05:  {{2}}
06:  {{},{1}}
08:  {{},{},{}}
09:  {{1},{1}}
10:  {{},{2}}
11:  {{3}}
12:  {{},{},{1}}
13:  {{1,2}}
15:  {{1},{2}}
16:  {{},{},{},{}}
17:  {{4}}
18:  {{},{1},{1}}
20:  {{},{},{2}}
22:  {{},{3}}
24:  {{},{},{},{1}}
25:  {{2},{2}}
26:  {{},{1,2}}
27:  {{1},{1},{1}}
29:  {{1,3}}
30:  {{},{1},{2}}
31:  {{5}}
32:  {{},{},{},{},{}}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],Or[#===1,And@@SquareFreeQ/@PrimePi/@FactorInteger[#][[All,1]]]&]
  • PARI
    ok(n)={!#select(p->!issquarefree(primepi(p)), factor(n)[,1])} \\ Andrew Howroyd, Aug 26 2018
Previous Showing 21-30 of 466 results. Next