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.

Showing 1-10 of 92 results. Next

A302242 Total weight of the n-th multiset multisystem. Totally additive with a(prime(n)) = Omega(n).

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Apr 03 2018

Keywords

Comments

A multiset multisystem is a finite multiset of finite multisets of positive integers. The n-th multiset multisystem is constructed by factoring n into prime numbers and then factoring each prime index into prime numbers and taking their prime indices. This produces a unique multiset multisystem for each n, and every possible multiset multisystem is so constructed as n ranges over all positive integers.

Examples

			Sequence of finite multisets of finite multisets of positive integers begins: (), (()), ((1)), (()()), ((2)), (()(1)), ((11)), (()()()), ((1)(1)), (()(2)), ((3)), (()()(1)), ((12)), (()(11)), ((1)(2)), (()()()()), ((4)), (()(1)(1)), ((111)), (()()(2)).
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    a:= n-> add(add(j[2], j=ifactors(pi(i[1]))[2])*i[2], i=ifactors(n)[2]):
    seq(a(n), n=1..100);  # Alois P. Heinz, Sep 07 2018
  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Total[PrimeOmega/@primeMS[n]],{n,100}]
  • PARI
    a(n,f=factor(n))=sum(i=1,#f~, bigomega(primepi(f[i,1]))*f[i,2]) \\ Charles R Greathouse IV, Nov 10 2021

A281116 Number of factorizations of n>=2 into factors greater than 1 with no common divisor other than 1 (a(1)=0 by convention).

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 2, 0, 1, 1, 0, 0, 2, 0, 2, 1, 1, 0, 3, 0, 1, 0, 2, 0, 4, 0, 0, 1, 1, 1, 5, 0, 1, 1, 3, 0, 4, 0, 2, 2, 1, 0, 5, 0, 2, 1, 2, 0, 3, 1, 3, 1, 1, 0, 8, 0, 1, 2, 0, 1, 4, 0, 2, 1, 4, 0, 9, 0, 1, 2, 2, 1, 4, 0, 5, 0, 1, 0, 8, 1, 1, 1, 3, 0, 8, 1, 2, 1, 1, 1, 7, 0, 2, 2, 5
Offset: 1

Views

Author

Gus Wiseman, Jan 15 2017

Keywords

Comments

Let (e1, e2, ..., ek) be a prime-signature of n (that is, n = p^e1 * q^e2 * ... * r^ek for some primes, p, q, ..., r). Then a(n) is the number of ways of partitioning multiset {e1 x 1, e2 x 2, ..., ek x k} into multisets such that none of the numbers 1 .. k is present in all member multisets of that set partition. - Antti Karttunen, Sep 08 2018

Examples

			a(6)=1:  (2*3)
a(12)=2; (2*2*3)       (3*4)
a(24)=3: (2*2*2*3)     (2*3*4)     (3*8)
a(30)=4: (2*3*5)       (2*15)      (3*10)    (5*6)
a(36)=5: (2*2*3*3)     (2*2*9)     (2*3*6)   (3*3*4)   (4*9)
a(96)=7: (2*2*2*2*2*3) (2*2*2*3*4) (2*2*3*8) (2*3*4*4) (2*3*16) (3*4*8) (3*32).
		

Crossrefs

Programs

  • Mathematica
    postfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[postfacs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[postfacs[n],GCD@@#===1&]],{n,2,100}]
  • PARI
    A281116(n, m=n, facs=List([])) = if(1==n, (1==gcd(Vec(facs))), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A281116(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Sep 08 2018

Extensions

Term a(1) = 0 prepended by Antti Karttunen, Sep 08 2018

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

A303386 Number of aperiodic factorizations of n > 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 2, 3, 1, 4, 1, 4, 2, 2, 1, 7, 1, 2, 2, 4, 1, 5, 1, 6, 2, 2, 2, 7, 1, 2, 2, 7, 1, 5, 1, 4, 4, 2, 1, 12, 1, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 11, 1, 2, 4, 7, 2, 5, 1, 4, 2, 5, 1, 16, 1, 2, 4, 4, 2, 5, 1, 12, 3, 2, 1, 11, 2, 2, 2, 7, 1, 11, 2, 4, 2, 2, 2, 19, 1, 4, 4, 7, 1, 5, 1, 7, 5
Offset: 2

Views

Author

Gus Wiseman, Apr 23 2018

Keywords

Comments

An aperiodic factorization of n is a finite multiset of positive integers greater than 1 whose product is n and whose multiplicities are relatively prime.

Examples

			The a(36) = 7 aperiodic factorizations are (2*2*9), (2*3*6), (2*18), (3*3*4), (3*12), (4*9), and (36). Missing from this list are (2*2*3*3) and (6*6).
		

Crossrefs

Programs

  • Mathematica
    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],GCD@@Length/@Split[#]===1&]],{n,2,100}]
  • PARI
    A001055(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A001055(n/d, d))); (s));
    A052409(n) = { my(k=ispower(n)); if(k, k, n>1); }; \\ From A052409
    A303386(n) = if(1==n,n,my(r); sumdiv(A052409(n),d, ispower(n,d,&r); moebius(d)*A001055(r))); \\ Antti Karttunen, Sep 25 2018

Formula

a(n) = Sum_{d|A052409(n)} mu(d) * A001055(n^(1/d)), where mu = A008683.

Extensions

More terms from Antti Karttunen, Sep 25 2018

A066739 Number of representations of n as a sum of products of positive integers. 1 is not allowed as a factor, unless it is the only factor. Representations which differ only in the order of terms or factors are considered equivalent.

Original entry on oeis.org

1, 1, 2, 3, 6, 8, 14, 19, 32, 44, 67, 91, 139, 186, 269, 362, 518, 687, 960, 1267, 1747, 2294, 3106, 4052, 5449, 7063, 9365, 12092, 15914, 20422, 26639, 34029, 44091, 56076, 72110, 91306, 116808, 147272, 187224, 235201, 297594, 372390, 468844, 584644, 732942
Offset: 0

Views

Author

Naohiro Nomoto, Jan 16 2002

Keywords

Examples

			For n=5, 5 = 4+1 = 2*2+1 = 3+2 = 3+1+1 = 2+2+1 = 2+1+1+1 = 1+1+1+1+1, so a(5) = 8.
For n=8, 8 = 4*2 = 2*2*2 = ... = 4+4 = 2*2+4 = 2*2+2*2 = ...; note that there are 3 ways to factor the terms of 4+4. In general, if a partition contains a number k exactly r times, then the number of ways to factor the k's is the binomial coefficient C(A001055(k)+r-1,r).
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    b:= proc(n, k) option remember;
          `if`(n>k, 0, 1) +`if`(isprime(n), 0,
          add(`if`(d>k, 0, b(n/d, d)), d=divisors(n) minus {1, n}))
        end:
    a:= proc(n) option remember;
          `if`(n=0, 1, add(add(d*b(d, d), d=divisors(j)) *a(n-j), j=1..n)/n)
        end:
    seq(a(n), n=0..60); # Alois P. Heinz, Apr 22 2012
  • Mathematica
    p[ n_, 1 ] := If[ n==1, 1, 0 ]; p[ 1, k_ ] := 1; p[ n_, k_ ] := p[ n, k ]=p[ n, k-1 ]+If[ Mod[ n, k ]==0, p[ n/k, k ], 0 ]; A001055[ n_ ] := p[ n, n ]; a[ n_, 1 ] := 1; a[ 0, k_ ] := 1; a[ n_, k_ ] := If[ k>n, a[ n, n ], a[ n, k ]=a[ n, k-1 ]+Sum[ Binomial[ A001055[ k ]+r-1, r ]a[ n-k*r, k-1 ], {r, 1, Floor[ n/k ]} ] ]; a[ n_ ] := a[ n, n ]; (* p[ n, k ]=number of factorizations of n with factors <= k. a[ n, k ]=number of representations of n as a sum of products of positive integers, with summands <= k *)
    b[n_, k_] := b[n, k] = If[n>k, 0, 1] + If[PrimeQ[n], 0, Sum[If[d>k, 0, b[n/d, d]], {d, Divisors[n] ~Complement~ {1, n}}]]; a[0] = 1; a[n_] := a[n] = If[n == 0, 1, Sum[DivisorSum[j, #*b[#, #]&]*a[n-j], {j, 1, n}]/n]; Table[a[n], {n, 0, 60}] (* Jean-François Alcover, Nov 10 2015, after Alois P. Heinz *)
    facs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#1,d]&)/@Select[facs[n/d],Min@@#1>=d&],{d,Rest[Divisors[n]]}]];
    Table[Length[Union[Sort/@Join@@Table[Tuples[facs/@ptn],{ptn,IntegerPartitions[n]}]]],{n,50}] (* Gus Wiseman, Sep 05 2018 *)
  • Python
    from sympy.core.cache import cacheit
    from sympy import divisors, isprime
    @cacheit
    def b(n, k): return (0 if n>k else 1) + (0 if isprime(n) else sum([0 if d>k else b(n//d, d) for d in divisors(n)[1:-1]]))
    @cacheit
    def a(n): return 1 if n==0 else sum(sum(d*b(d, d) for d in divisors(j))*a(n - j)  for j in range(1, n + 1))//n
    print([a(n) for n in range(61)]) # Indranil Ghosh, Aug 19 2017, after Maple code

Formula

a(n) = Sum_{pi} Product_{m=1..n} binomial(k(m)+A001055(m)-1, k(m)), where pi runs through all partitions k(1) + 2 * k( 2) + ... + n * k(n) = n. a(n)=1/n*Sum_{m=1..n} a(n-m)*b(m), n > 0, a(0)=1, b(m)=Sum_{d|m} d*A001055(d). Euler transform of A001055(n): Product_{m=1..infinity} (1-x^m)^(-A001055(m)). - Vladeta Jovovic, Jan 21 2002

Extensions

Edited by Dean Hickerson, Jan 19 2002

A296119 Number of ways to choose a strict factorization of each factor in a factorization of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 4, 2, 3, 1, 7, 1, 3, 3, 7, 1, 7, 1, 7, 3, 3, 1, 16, 2, 3, 4, 7, 1, 12, 1, 12, 3, 3, 3, 21, 1, 3, 3, 16, 1, 12, 1, 7, 7, 3, 1, 33, 2, 7, 3, 7, 1, 16, 3, 16, 3, 3, 1, 34, 1, 3, 7, 23, 3, 12, 1, 7, 3, 12, 1, 50, 1, 3, 7, 7, 3, 12, 1, 33, 7, 3
Offset: 1

Views

Author

Gus Wiseman, Dec 05 2017

Keywords

Examples

			The a(24) = 16 twice-factorizations:
(2)*(2)*(2)*(3),
(2)*(2)*(2*3), (2)*(2)*(6), (2)*(3)*(4),
(2)*(2*6), (2)*(3*4), (2)*(12), (3)*(2*4), (3)*(8), (4)*(2*3), (4)*(6),
(2*3*4), (2*12), (3*8), (4*6), (24).
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Sum[Times@@(Length[Select[facs[#],UnsameQ@@#&]]&/@fac),{fac,facs[n]}],{n,100}]
  • PARI
    A045778(n, m=n) = ((n<=m) + sumdiv(n, d, if((d>1)&&(d<=m)&&(dA045778(n/d, d-1))));
    A296119(n, m=n) = if(1==n, 1, sumdiv(n, d, if((d>1)&&(d<=m), A045778(d)*A296119(n/d, d)))); \\ Antti Karttunen, Oct 08 2018

Formula

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

A281118 a(1)=1, a(n>1) = number of tree-factorizations of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 4, 2, 2, 1, 6, 1, 2, 2, 12, 1, 6, 1, 6, 2, 2, 1, 20, 2, 2, 4, 6, 1, 8, 1, 32, 2, 2, 2, 28, 1, 2, 2, 20, 1, 8, 1, 6, 6, 2, 1, 76, 2, 6, 2, 6, 1, 20, 2, 20, 2, 2, 1, 38, 1, 2, 6, 112, 2, 8, 1, 6, 2, 8, 1, 116, 1, 2, 6, 6, 2, 8, 1, 76, 12, 2, 1
Offset: 1

Views

Author

Gus Wiseman, Jan 15 2017

Keywords

Comments

A tree-factorization of n>=2 is either (case 1) the number n or (case 2) a sequence of two or more tree-factorizations, one of each part of a weakly increasing factorization of n. These are rooted plane trees and the ordering of branches is important. For example, {{2,2},9}, {2,{2,9}}, {{2,2},{3,3}}, {6,{2,3}}, and {{2,3},6} are distinct tree-factorizations of 36, but {9,{2,2}}, {{2,9},2}, and {{3,3},{2,2}} are not.
a(n) depends only on the prime signature of n. - Andrew Howroyd, Nov 18 2018

Examples

			The a(30)=8 tree-factorizations are 30, 2*15, 2*(3*5), 3*10, 3*(2*5), 5*6, 5*(2*3), 2*3*5.
		

Crossrefs

Programs

  • Mathematica
    postfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[postfacs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    treefacs[n_]:=If[n<=1,{{}},Prepend[Join@@Function[q,Tuples[treefacs/@q]]/@DeleteCases[postfacs[n],{n}],n]];
    Table[Length[treefacs[n]],{n,2,83}]
  • PARI
    seq(n)={my(v=vector(n), w=vector(n)); w[1]=v[1]=1; for(k=2, n, w[k]=v[k]+1; forstep(j=n\k*k, k, -k, my(i=j, e=0); while(i%k==0, i/=k; e++; v[j] += w[k]^e*v[i]))); w} \\ Andrew Howroyd, Nov 18 2018

Formula

a(p^n) = A289501(n) for prime p. - Andrew Howroyd, Nov 18 2018

A302494 Products of distinct primes of squarefree index.

Original entry on oeis.org

1, 2, 3, 5, 6, 10, 11, 13, 15, 17, 22, 26, 29, 30, 31, 33, 34, 39, 41, 43, 47, 51, 55, 58, 59, 62, 65, 66, 67, 73, 78, 79, 82, 83, 85, 86, 87, 93, 94, 101, 102, 109, 110, 113, 118, 123, 127, 129, 130, 134, 137, 139, 141, 143, 145, 146, 149, 155, 157, 158, 163
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 systems.
01: {}
02: {{}}
03: {{1}}
05: {{2}}
06: {{},{1}}
10: {{},{2}}
11: {{3}}
13: {{1,2}}
15: {{1},{2}}
17: {{4}}
22: {{},{3}}
26: {{},{1,2}}
29: {{1,3}}
30: {{},{1},{2}}
31: {{5}}
33: {{1},{3}}
34: {{},{4}}
39: {{1},{1,2}}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],Or[#===1,SquareFreeQ[#]&&And@@SquareFreeQ/@PrimePi/@FactorInteger[#][[All,1]]]&]
  • PARI
    is(n) = if(bigomega(n)!=omega(n), return(0), my(f=factor(n)[, 1]~); for(k=1, #f, if(!issquarefree(primepi(f[k])) && primepi(f[k])!=1, return(0)))); 1 \\ Felix Fröhlich, Apr 10 2018

A302243 Total weight of the n-th twice-odd-factored multiset partition.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Apr 03 2018

Keywords

Comments

A multiset partition is a finite multiset of finite nonempty multisets of positive integers. The n-th twice-odd-factored multiset partition is constructed by factoring 2n + 1 into prime numbers and then factoring each prime index into prime numbers and taking their prime indices.

Examples

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

Crossrefs

Programs

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

Formula

a(n) = A302242(2n + 1).

A295935 Number of twice-factorizations of n where the latter factorizations are constant, i.e., type (P,P,R).

Original entry on oeis.org

1, 1, 1, 3, 1, 2, 1, 5, 3, 2, 1, 5, 1, 2, 2, 12, 1, 5, 1, 5, 2, 2, 1, 10, 3, 2, 5, 5, 1, 5, 1, 18, 2, 2, 2, 15, 1, 2, 2, 10, 1, 5, 1, 5, 5, 2, 1, 22, 3, 5, 2, 5, 1, 10, 2, 10, 2, 2, 1, 13, 1, 2, 5, 40, 2, 5, 1, 5, 2, 5, 1, 28, 1, 2, 5, 5, 2, 5, 1, 22, 12, 2, 1
Offset: 1

Views

Author

Gus Wiseman, Nov 29 2017

Keywords

Comments

a(n) is also the number of ways to choose a perfect divisor of each factor in a factorization of n.

Examples

			The a(24) = 10 twice-factorizations are:
(2)*(2)*(2)*(3), (2)*(3)*(2*2), (3)*(2*2*2)
(2)*(2)*(6), (2*2)*(6),
(2)*(3)*(4),
(2)*(12),
(3)*(8),
(4)*(6),
(24).
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Sum[Product[Length[Divisors[GCD@@FactorInteger[d][[All,2]]]],{d,f}],{f,facs[n]}],{n,100}]

Formula

Dirichlet g.f.: 1/Product_{n > 1}(1 - A089723(n)/n^s).
Showing 1-10 of 92 results. Next