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-4 of 4 results.

A050361 Number of factorizations into distinct prime powers greater than 1.

Original entry on oeis.org

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

Views

Author

Christian G. Bower, Oct 15 1999

Keywords

Comments

a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24 = 2^3*3 and 375 = 3*5^3 both have prime signature (3,1).
The number of unordered factorizations of n into 1 and exponentially odd prime powers, i.e., p^e where p is a prime and e is odd (A246551). - Amiram Eldar, Jun 12 2025

Examples

			From _Gus Wiseman_, Jul 30 2022: (Start)
The A000688(216) = 9 factorizations of 216 into prime powers are:
  (2*2*2*3*3*3)
  (2*2*2*3*9)
  (2*2*2*27)
  (2*3*3*3*4)
  (2*3*4*9)
  (2*4*27)
  (3*3*3*8)
  (3*8*9)
  (8*27)
Of these, the a(216) = 4 strict cases are:
  (2*3*4*9)
  (2*4*27)
  (3*8*9)
  (8*27)
(End)
		

Crossrefs

Cf. A124010.
This is the strict case of A000688.
Positions of 1's are A004709, complement A046099.
The case of primes (instead of prime-powers) is A008966, non-strict A000012.
The non-strict additive version allowing 1's A023893, ranked by A302492.
The non-strict additive version is A023894, ranked by A355743.
The additive version (partitions) is A054685, ranked by A356065.
The additive version allowing 1's is A106244, ranked by A302496.
A001222 counts prime-power divisors.
A005117 lists all squarefree numbers.
A034699 gives maximal prime-power divisor.
A246655 lists all prime-powers (A000961 includes 1), towers A164336.
A296131 counts twice-factorizations of type PQR, non-strict A295935.

Programs

  • Haskell
    a050361 = product . map a000009 . a124010_row
    -- Reinhard Zumkeller, Aug 28 2014
    
  • Maple
    A050361 := proc(n)
        local a,f;
        if n = 1 then
            1;
        else
            a := 1 ;
            for f in ifactors(n)[2] do
                a := a*A000009(op(2,f)) ;
            end do:
        end if;
    end proc: # R. J. Mathar, May 25 2017
  • Mathematica
    Table[Times @@ PartitionsQ[Last /@ FactorInteger[n]], {n, 99}] (* Arkadiusz Wesolowski, Feb 27 2017 *)
  • PARI
    A000009(n,k=(n-!(n%2))) = if(!n,1,my(s=0); while(k >= 1, if(k<=n, s += A000009(n-k,k)); k -= 2); (s));
    A050361(n) = factorback(apply(A000009,factor(n)[,2])); \\ Antti Karttunen, Nov 17 2019

Formula

Dirichlet g.f.: Product_{n is a prime power >1}(1 + 1/n^s).
Multiplicative with a(p^e) = A000009(e).
a(A002110(k))=1.
a(n) = A050362(A101296(n)). - R. J. Mathar, May 26 2017
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} f(1/p) = 1.26020571070524171076..., where f(x) = (1-x) * Product_{k>=1} (1 + x^k). - Amiram Eldar, Oct 03 2023

A356065 Squarefree numbers whose prime indices are all prime-powers.

Original entry on oeis.org

1, 3, 5, 7, 11, 15, 17, 19, 21, 23, 31, 33, 35, 41, 51, 53, 55, 57, 59, 67, 69, 77, 83, 85, 93, 95, 97, 103, 105, 109, 115, 119, 123, 127, 131, 133, 155, 157, 159, 161, 165, 177, 179, 187, 191, 201, 205, 209, 211, 217, 227, 231, 241, 249, 253, 255, 265, 277
Offset: 1

Views

Author

Gus Wiseman, Jul 25 2022

Keywords

Examples

			105 has prime indices {2,3,4}, all three of which are prime-powers, so 105 is in the sequence.
		

Crossrefs

The multiplicative version (factorizations) is A050361, non-strict A000688.
Heinz numbers of the partitions counted by A054685, with 1's A106244, non-strict A023894, non-strict with 1's A023893.
Counting twice-partitions of this type gives A279786, non-strict A279784.
Counting twice-factorizations gives A295935, non-strict A296131.
These are the odd products of distinct elements of A302493.
Allowing prime index 1 gives A302496, non-strict A302492.
The case of primes (instead of prime-powers) is A302590, non-strict A076610.
These are the squarefree positions of 1's in A355741.
This is the squarefree case of A355743, complement A356066.
A001222 counts prime-power divisors.
A005117 lists the squarefree numbers.
A034699 gives maximal prime-power divisor.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
A355742 chooses a prime-power divisor of each prime index.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],SquareFreeQ[#]&&And@@PrimePowerQ/@primeMS[#]&]

Formula

Intersection of A005117 and A355743.

A382215 MM-numbers of multiset partitions into constant blocks with a common sum.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 16, 17, 19, 23, 25, 27, 31, 32, 35, 41, 49, 53, 59, 64, 67, 81, 83, 97, 103, 109, 121, 125, 127, 128, 131, 157, 175, 179, 191, 209, 211, 227, 241, 243, 245, 256, 277, 283, 289, 311, 331, 343, 353, 361, 367, 391, 401, 419, 431, 461
Offset: 1

Views

Author

Gus Wiseman, Mar 21 2025

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, sum A056239. The multiset of multisets with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset of multisets with MM-number 78 is {{},{1},{1,2}}.

Examples

			The terms together with their prime indices of prime indices begin:
   1: {}
   2: {{}}
   3: {{1}}
   4: {{},{}}
   5: {{2}}
   7: {{1,1}}
   8: {{},{},{}}
   9: {{1},{1}}
  11: {{3}}
  16: {{},{},{},{}}
  17: {{4}}
  19: {{1,1,1}}
  23: {{2,2}}
  25: {{2},{2}}
  27: {{1},{1},{1}}
  31: {{5}}
  32: {{},{},{},{},{}}
  35: {{2},{1,1}}
  41: {{6}}
  49: {{1,1},{1,1}}
  53: {{1,1,1,1}}
  59: {{7}}
		

Crossrefs

Twice-partitions of this type are counted by A279789.
For just constant blocks we have A302492, counted by A000688.
For sets of constant multisets we have A302496, counted by A050361.
For just common sums we have A326534, counted by A321455.
Factorizations of this type are counted by A381995.
For strict blocks and distinct sums we have A382201, counted by A381633.
Normal multiset partitions of this type are counted by A382204.
For strict instead of constant blocks we have A382304, counted by A382080.
For sets of constant multisets with distinct sums A382426, counted by A381635.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n], {p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],SameQ@@Total/@prix/@prix[#] && And@@SameQ@@@prix/@prix[#]&]
  • PARI
    is(k) = my(f=factor(k)[, 1]~, k, p, v=vector(#f, i, primepi(f[i]))); for(i=1, #v, k=isprimepower(v[i], &p); if(k||v[i]==1, v[i]=k*primepi(p), return(0))); #Set(v)<2; \\ Jinyuan Wang, Apr 02 2025

Formula

Equals A326534 /\ A302492.

A302498 Numbers that are a power of a prime number whose prime index is itself a power of a prime number.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 16, 17, 19, 23, 25, 27, 31, 32, 41, 49, 53, 59, 64, 67, 81, 83, 97, 103, 109, 121, 125, 127, 128, 131, 157, 179, 191, 211, 227, 241, 243, 256, 277, 283, 289, 311, 331, 343, 353, 361, 367, 401, 419, 431, 461, 509, 512, 529, 547, 563
Offset: 1

Views

Author

Gus Wiseman, Apr 09 2018

Keywords

Comments

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

Examples

			49 is in the sequence because 49 = prime(prime(1)^2)^2.
Entry A302242 describes a correspondence between positive integers and multiset multisystems. In this case it gives the following sequence of constant constant-multiset multisystems.
01: {}
02: {{}}
03: {{1}}
04: {{},{}}
05: {{2}}
07: {{1,1}}
08: {{},{},{}}
09: {{1},{1}}
11: {{3}}
16: {{},{},{},{}}
17: {{4}}
19: {{1,1,1}}
23: {{2,2}}
25: {{2},{2}}
27: {{1},{1},{1}}
31: {{5}}
32: {{},{},{},{},{}}
41: {{6}}
49: {{1,1},{1,1}}
53: {{1,1,1,1}}
59: {{7}}
64: {{},{},{},{},{},{}}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],Or[#===1,PrimePowerQ[#]&&And@@(Or[#===1,PrimePowerQ[#]]&/@PrimePi/@FactorInteger[#][[All,1]])]&]
  • PARI
    ok(n)={my(p); n == 1 || (isprimepower(n, &p) && (p == 2 || isprimepower(primepi(p))))} \\ Andrew Howroyd, Aug 26 2018
Showing 1-4 of 4 results.