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

A089723 a(1)=1; for n>1, a(n) gives number of ways to write n as n = x^y, 2 <= x, 1 <= y.

Original entry on oeis.org

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

Views

Author

Naohiro Nomoto, Jan 07 2004

Keywords

Comments

This function depends only on the prime signature of n. - Franklin T. Adams-Watters, Mar 10 2006
a(n) is the number of perfect divisors of n. Perfect divisor of n is divisor d such that d^k = n for some k >= 1. a(n) > 1 for perfect powers n = A001597(m) for m > 2. - Jaroslav Krizek, Jan 23 2010
Also the number of uniform perfect integer partitions of n - 1. An integer partition of n is uniform if all parts appear with the same multiplicity, and perfect if every nonnegative integer up to n is the sum of a unique submultiset. The Heinz numbers of these partitions are given by A326037. The a(16) = 3 partitions are: (8,4,2,1), (4,4,4,1,1,1), (1,1,1,1,1,1,1,1,1,1,1,1,1,1,1). - Gus Wiseman, Jun 07 2019
The record values occur at 1 and at 2^A002182(n) for n > 1. - Amiram Eldar, Nov 06 2020

Examples

			144 = 2^4 * 3^2, gcd(4,2) = 2, d(2) = 2, so a(144) = 2. The representations are 144^1 and 12^2.
From _Friedjof Tellkamp_, Jun 14 2025: (Start)
n:          1, 2, 3, 4, 5, 6, 7, 8, 9, ...
----------------------------------------------------
1st powers: 1, 1, 1, 1, 1, 1, 1, 1, 1, ... (A000012)
Squares:    1, 0, 0, 1, 0, 0, 0, 0, 1, ... (A010052)
Cubes:      1, 0, 0, 0, 0, 0, 0, 1, 0, ... (A010057)
Quartics:   1, 0, 0, 0, 0, 0, 0, 0, 0, ... (A374016)
...
Sum:       oo, 1, 1, 2, 1, 1, 1, 2, 2, ...
a(1)=1:     1, 1, 1, 2, 1, 1, 1, 2, 2, ... (= this sequence). (End)
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    A089723 := proc(n) local t1,t2,g,j;
    if n=1 then 1 else
    t1:=ifactors(n)[2]; t2:=nops(t1); g := t1[1][2];
    for j from 2 to t2 do g:=gcd(g,t1[j][2]); od:
    tau(g); fi; end;
    [seq(A089723(n),n=1..100)]; # N. J. A. Sloane, Nov 10 2016
  • Mathematica
    Table[DivisorSigma[0, GCD @@ FactorInteger[n][[All, 2]]], {n, 100}] (* Gus Wiseman, Jun 12 2017 *)
  • PARI
    a(n) = if (n==1, 1, numdiv(gcd(factor(n)[,2]))); \\ Michel Marcus, Jun 13 2017
    
  • Python
    from math import gcd
    from sympy import factorint, divisor_sigma
    def a(n):
        if n == 1: return 1
        e = list(factorint(n).values())
        g = e[0]
        for ei in e[1:]: g = gcd(g, ei)
        return divisor_sigma(g, 0)
    print([a(n) for n in range(1, 105)]) # Michael S. Branicky, Jul 15 2021

Formula

If n = Product p_i^e_i, a(n) = d(gcd()). - Franklin T. Adams-Watters, Mar 10 2006
Sum_{n=1..m} a(n) = A255165(m) + 1. - Richard R. Forberg, Feb 16 2015
Sum_{n>=2} a(n)/n^s = Sum_{n>=2} 1/(n^s-1) = Sum_{k>=1} (zeta(s*k)-1) for all real s with Re(s) > 1 (Golomb, 1973). - Amiram Eldar, Nov 06 2020
For n > 1, a(n) = Sum_{i=1..floor(n/2)} floor(n^(1/i))-floor((n-1)^(1/i)). - Wesley Ivan Hurt, Dec 08 2020
Sum_{n>=1} (a(n)-1)/n = 1 (Mycielski, 1951). - Amiram Eldar, Jul 15 2021
From Friedjof Tellkamp, Jun 14 2025: (Start)
a(n) = 1 + A259362(n) = 1 + A010052(n) + A010057(n) + A374016(n) + (...), for n > 1.
G.f.: x + Sum_{j>=2, k>=1} x^(j^k). (End)

A326020 Number of complete subsets of {1..n}.

Original entry on oeis.org

1, 2, 3, 4, 6, 9, 15, 27, 50, 95, 185, 365, 724, 1441, 2873, 5735, 11458, 22902, 45789, 91561, 183102, 366180, 732331, 1464626, 2929209, 5858367, 11716674, 23433277, 46866473, 93732852, 187465596, 374931067, 749861989, 1499723808, 2999447418
Offset: 0

Views

Author

Gus Wiseman, Jun 04 2019

Keywords

Comments

A set of positive integers summing to n is complete if every nonnegative integer up to n is the sum of some subset.

Examples

			The a(0) = 1 through a(6) = 15 subsets:
  {}  {}   {}     {}       {}         {}           {}
      {1}  {1}    {1}      {1}        {1}          {1}
           {1,2}  {1,2}    {1,2}      {1,2}        {1,2}
                  {1,2,3}  {1,2,3}    {1,2,3}      {1,2,3}
                           {1,2,4}    {1,2,4}      {1,2,4}
                           {1,2,3,4}  {1,2,3,4}    {1,2,3,4}
                                      {1,2,3,5}    {1,2,3,5}
                                      {1,2,4,5}    {1,2,3,6}
                                      {1,2,3,4,5}  {1,2,4,5}
                                                   {1,2,4,6}
                                                   {1,2,3,4,5}
                                                   {1,2,3,4,6}
                                                   {1,2,3,5,6}
                                                   {1,2,4,5,6}
                                                   {1,2,3,4,5,6}
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]],Union[Plus@@@Subsets[#]]==Range[0,Total[#]]&]],{n,0,10}]

Extensions

a(17)-a(34) from Charlie Neder, Jun 05 2019

A326035 Number of uniform knapsack partitions of n.

Original entry on oeis.org

1, 1, 2, 3, 4, 4, 6, 6, 9, 10, 12, 12, 17, 16, 20, 25, 27, 29, 35, 39, 44, 57, 53, 66, 75, 84, 84, 114, 112, 131, 133, 162, 167, 209, 192, 242, 250, 289, 279, 363, 348, 417, 404, 502, 487, 608, 557, 706, 682, 835, 773, 1004, 922, 1149, 1059, 1344, 1257, 1595
Offset: 0

Views

Author

Gus Wiseman, Jun 04 2019

Keywords

Comments

An integer partition is uniform if all parts appear with the same multiplicity, and knapsack if every distinct submultiset has a different sum.

Examples

			The a(1) = 1 through a(8) = 9 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (21)   (22)    (32)     (33)      (43)       (44)
             (111)  (31)    (41)     (42)      (52)       (53)
                    (1111)  (11111)  (51)      (61)       (62)
                                     (222)     (421)      (71)
                                     (111111)  (1111111)  (521)
                                                          (2222)
                                                          (3311)
                                                          (11111111)
		

Crossrefs

Programs

  • Mathematica
    sums[ptn_]:=sums[ptn]=If[Length[ptn]==1,ptn,Union@@(Join[sums[#],sums[#]+Total[ptn]-Total[#]]&/@Union[Table[Delete[ptn,i],{i,Length[ptn]}]])];
    ks[n_]:=Select[IntegerPartitions[n],Length[sums[Sort[#]]]==Times@@(Length/@Split[#]+1)-1&];
    Table[Length[Select[ks[n],SameQ@@Length/@Split[#]&]],{n,30}]

A326021 Number of complete subsets of {1..n} with maximum n.

Original entry on oeis.org

1, 1, 1, 2, 3, 6, 12, 23, 45, 90, 180, 359, 717, 1432, 2862, 5723, 11444, 22887, 45772, 91541, 183078, 366151, 732295, 1464583, 2929158, 5858307, 11716603, 23433196, 46866379, 93732744, 187465471, 374930922, 749861819, 1499723610
Offset: 1

Views

Author

Gus Wiseman, Jun 04 2019

Keywords

Comments

A set of positive integers summing to n is complete if every nonnegative integer up to n is the sum of some subset.

Examples

			The a(1) = 1 through a(7) = 12 subsets:
  {1}  {1,2}  {1,2,3}  {1,2,4}    {1,2,3,5}    {1,2,3,6}      {1,2,3,7}
                       {1,2,3,4}  {1,2,4,5}    {1,2,4,6}      {1,2,4,7}
                                  {1,2,3,4,5}  {1,2,3,4,6}    {1,2,3,4,7}
                                               {1,2,3,5,6}    {1,2,3,5,7}
                                               {1,2,4,5,6}    {1,2,3,6,7}
                                               {1,2,3,4,5,6}  {1,2,4,5,7}
                                                              {1,2,4,6,7}
                                                              {1,2,3,4,5,7}
                                                              {1,2,3,4,6,7}
                                                              {1,2,3,5,6,7}
                                                              {1,2,4,5,6,7}
                                                              {1,2,3,4,5,6,7}
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]],Max@@#==n&&Union[Plus@@@Subsets[#]]==Range[0,Total[#]]&]],{n,10}]

Extensions

a(18)-a(34) from Charlie Neder, Jun 05 2019

A326022 Number of minimal complete subsets of {1..n} with maximum n.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 2, 4, 8, 8, 8, 10, 14, 25, 40, 49, 62
Offset: 1

Views

Author

Gus Wiseman, Jun 04 2019

Keywords

Comments

A set of positive integers summing to m is complete if every nonnegative integer up to m is the sum of some subset. For example, (1,2,3,6,13) is a complete set because we have:
0 = (empty sum)
1 = 1
2 = 2
3 = 3
4 = 1 + 3
5 = 2 + 3
6 = 6
7 = 6 + 1
8 = 6 + 2
9 = 6 + 3
10 = 1 + 3 + 6
11 = 2 + 3 + 6
12 = 1 + 2 + 3 + 6
and the remaining numbers 13-25 are obtained by adding 13 to each of these.

Examples

			The a(3) = 1 through a(9) = 8 subsets:
  {1,2,3}  {1,2,4}  {1,2,3,5}  {1,2,3,6}  {1,2,3,7}  {1,2,4,8}    {1,2,3,4,9}
                    {1,2,4,5}  {1,2,4,6}  {1,2,4,7}  {1,2,3,5,8}  {1,2,3,5,9}
                                                     {1,2,3,6,8}  {1,2,3,6,9}
                                                     {1,2,3,7,8}  {1,2,3,7,9}
                                                                  {1,2,4,5,9}
                                                                  {1,2,4,6,9}
                                                                  {1,2,4,7,9}
                                                                  {1,2,4,8,9}
		

Crossrefs

Programs

  • Mathematica
    fasmin[y_]:=Complement[y,Union@@Table[Union[s,#]&/@Rest[Subsets[Complement[Union@@y,s]]],{s,y}]];
    Table[Length[fasmin[Select[Subsets[Range[n]],Max@@#==n&&Union[Plus@@@Subsets[#]]==Range[0,Total[#]]&]]],{n,10}]

A326037 Heinz numbers of uniform perfect integer partitions.

Original entry on oeis.org

1, 2, 4, 6, 8, 16, 32, 42, 64, 100, 128, 256, 512, 798, 1024, 2048, 2744, 4096, 8192, 16384, 32768, 42294, 52900, 65536
Offset: 1

Views

Author

Gus Wiseman, Jun 04 2019

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
An integer partition of n is uniform if all parts appear with the same multiplicity, and perfect if every nonnegative integer up to n is the sum of a unique submultiset.
The enumeration of these partitions by sum is given by A089723.

Examples

			The sequence of all uniform perfect integer partitions together with their Heinz numbers begins:
      1: ()
      2: (1)
      4: (11)
      6: (21)
      8: (111)
     16: (1111)
     32: (11111)
     42: (421)
     64: (111111)
    100: (3311)
    128: (1111111)
    256: (11111111)
    512: (111111111)
    798: (8421)
   1024: (1111111111)
   2048: (11111111111)
   2744: (444111)
   4096: (111111111111)
   8192: (1111111111111)
  16384: (11111111111111)
		

Crossrefs

Programs

  • Mathematica
    hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
    Select[Range[1000],SameQ@@Last/@FactorInteger[#]&&Sort[hwt/@Divisors[#]]==Range[0,hwt[#]]&]

Formula

Intersection of A072774 (uniform), A299702 (knapsack), and A325781 (complete).
Showing 1-6 of 6 results.