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

A289509 Numbers k such that the gcd of the indices j for which the j-th prime prime(j) divides k is 1.

Original entry on oeis.org

2, 4, 6, 8, 10, 12, 14, 15, 16, 18, 20, 22, 24, 26, 28, 30, 32, 33, 34, 35, 36, 38, 40, 42, 44, 45, 46, 48, 50, 51, 52, 54, 55, 56, 58, 60, 62, 64, 66, 68, 69, 70, 72, 74, 75, 76, 77, 78, 80, 82, 84, 85, 86, 88, 90, 92, 93, 94, 95, 96, 98, 99, 100, 102, 104
Offset: 1

Views

Author

Christopher J. Smyth, Jul 11 2017

Keywords

Comments

Any integer k in the sequence encodes (by 'Heinz encoding' cf. A056239) a multiset of integers whose gcd is 1, namely the multiset containing r_j copies of j if k factors as Product_j prime(j)^{r_j} with gcd_j j = 1.
Clearly the sequence contains all even numbers and no odd primes or odd prime powers. It also clearly contains all numbers that are divisible by consecutive primes.
The sequence is the list of those k such that A289508(k) = 1.
It is also the list of those k such that A289506(k) = A289507(k).
Heinz numbers of integer partitions with relatively prime parts, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - Gus Wiseman, Apr 13 2018

Examples

			6 is a term because 6 = p_1*p_2 and gcd(1,2) = 1.
From _Gus Wiseman_, Apr 13 2018: (Start)
Sequence of integer partitions with relatively prime parts begins:
02 : (1)
04 : (11)
06 : (21)
08 : (111)
10 : (31)
12 : (211)
14 : (41)
15 : (32)
16 : (1111)
18 : (221)
20 : (311)
22 : (51)
24 : (2111)
26 : (61)
28 : (411)
30 : (321)
32 : (11111)
33 : (52)
34 : (71)
35 : (43)
36 : (2211)
38 : (81)
40 : (3111)
(End)
		

Crossrefs

Programs

  • Maple
    p:=1:for ind to 10000 do p:=nextprime(p);primeindex[p]:=ind;od:
    out:=[]:for n from 2 to 100 do m:=[];f:=ifactors(n)[2];g:=0;
    for k to nops(f) do mk:=primeindex[f[k][1]];m:=[op(m),mk];
    g:=gcd(g,mk);od; if g=1 then out:=[op(out),n];fi;od:out;
  • Mathematica
    Select[Range[200],GCD@@PrimePi/@FactorInteger[#][[All,1]]===1&] (* Gus Wiseman, Apr 13 2018 *)
  • PARI
    isok(n) = my(f=factor(n)); gcd(apply(x->primepi(x), f[,1])) == 1; \\ Michel Marcus, Jul 19 2017
    
  • Python
    from sympy import gcd, primepi, primefactors
    def ok(n): return gcd([primepi(p) for p in primefactors(n)]) == 1
    print([n for n in range(1, 151) if ok(n)]) # Indranil Ghosh, Aug 06 2017

A047966 a(n) = Sum_{ d divides n } q(d), where q(d) = A000009 = number of partitions of d into distinct parts.

Original entry on oeis.org

1, 2, 3, 4, 4, 8, 6, 10, 11, 15, 13, 25, 19, 29, 33, 42, 39, 62, 55, 81, 84, 103, 105, 153, 146, 185, 203, 253, 257, 344, 341, 432, 463, 552, 594, 747, 761, 920, 1003, 1200, 1261, 1537, 1611, 1921, 2089, 2410, 2591, 3095, 3270, 3815, 4138, 4769, 5121, 5972, 6394, 7367, 7974, 9066, 9793, 11305, 12077, 13736, 14940
Offset: 1

Views

Author

Keywords

Comments

Number of partitions of n such that every part occurs with the same multiplicity. - Vladeta Jovovic, Oct 22 2004
Christopher and Christober call such partitions uniform. - Gus Wiseman, Apr 16 2018
Equals inverse Mobius transform (A051731) * A000009, where the latter begins (1, 1, 2, 2, 3, 4, 5, 6, 8, ...). - Gary W. Adamson, Jun 08 2009

Examples

			The a(6) = 8 uniform partitions are (6), (51), (42), (33), (321), (222), (2211), (111111). - _Gus Wiseman_, Apr 16 2018
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    b:= proc(n) option remember; `if`(n=0, 1, add(add(
         `if`(d::odd, d, 0), d=divisors(j))*b(n-j), j=1..n)/n)
        end:
    a:= n-> add(b(d), d=divisors(n)):
    seq(a(n), n=1..100);  # Alois P. Heinz, Jul 11 2016
  • Mathematica
    b[n_] := b[n] = If[n==0, 1, Sum[DivisorSum[j, If[OddQ[#], #, 0]&]*b[n-j], {j, 1, n}]/n]; a[n_] := DivisorSum[n, b]; Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Dec 06 2016 after Alois P. Heinz *)
    Table[DivisorSum[n,PartitionsQ],{n,20}] (* Gus Wiseman, Apr 16 2018 *)
  • PARI
    N = 66; q='q+O('q^N);
    D(q)=eta(q^2)/eta(q); \\ A000009
    Vec( sum(e=1,N,D(q^e)-1) ) \\ Joerg Arndt, Mar 27 2014

Formula

G.f.: Sum_{k>0} (-1+Product_{i>0} (1+z^(k*i))). - Vladeta Jovovic, Jun 22 2003
G.f.: Sum_{k>=1} q(k)*x^k/(1 - x^k), where q() = A000009. - Ilya Gutkovskiy, Jun 20 2018
a(n) ~ exp(Pi*sqrt(n/3)) / (4*3^(1/4)*n^(3/4)). - Vaclav Kotesovec, Aug 27 2018

A302696 Numbers whose prime indices (with repetition) are pairwise coprime. Nonprime Heinz numbers of integer partitions with pairwise coprime parts.

Original entry on oeis.org

1, 2, 4, 6, 8, 10, 12, 14, 15, 16, 20, 22, 24, 26, 28, 30, 32, 33, 34, 35, 38, 40, 44, 46, 48, 51, 52, 55, 56, 58, 60, 62, 64, 66, 68, 69, 70, 74, 76, 77, 80, 82, 85, 86, 88, 92, 93, 94, 95, 96, 102, 104, 106, 110, 112, 116, 118, 119, 120, 122, 123, 124, 128, 132
Offset: 1

Views

Author

Gus Wiseman, Apr 11 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. Two or more numbers are coprime if no pair has a common divisor other than 1. A single number is not considered coprime unless it is equal to 1.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
Number 36 = prime(1)*prime(1)*prime(2)*prime(2) is not included in the sequence, because the pair of prime indices {2,2} is not coprime. - Gus Wiseman, Dec 06 2021

Examples

			Sequence of integer partitions with pairwise coprime parts begins: (), (1), (11), (21), (111), (31), (211), (41), (32), (1111), (311), (51), (2111), (61), (411), (321).
Missing from this list are: (2), (3), (4), (22), (5), (6), (7), (221), (8), (42), (9), (33), (222).
		

Crossrefs

Programs

  • Maple
    filter:= proc(n) local F;
       F:= ifactors(n)[2];
       if nops(F)=1 then if F[1][1] = 2 then return true else return false fi fi;
       if ormap(t -> t[2]>1 and t[1] <> 2, F) then return false fi;
       F:= map(t -> numtheory:-pi(t[1]), F);
       ilcm(op(F))=convert(F,`*`)
    end proc:
    select(filter, [$1..200]); # Robert Israel, Sep 10 2020
  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[200],Or[#===1,CoprimeQ@@primeMS[#]]&]
  • PARI
    isA302696(n) = if(isprimepower(n),!(n%2), if(!issquarefree(n>>valuation(n,2)), 0, my(pis=apply(primepi,factor(n)[,1])); (lcm(pis)==factorback(pis)))); \\ Antti Karttunen, Dec 06 2021

Extensions

Clarification (with repetition) added to the definition by Antti Karttunen, Dec 06 2021

A007359 Number of partitions of n into pairwise coprime parts that are >= 2.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 1, 3, 2, 3, 3, 5, 4, 6, 5, 5, 8, 9, 10, 11, 11, 10, 14, 18, 19, 18, 20, 20, 25, 30, 35, 34, 32, 32, 43, 43, 57, 56, 51, 55, 67, 78, 87, 87, 80, 82, 97, 125, 128, 127, 128, 127, 146, 182, 191, 185, 184, 193, 213, 263, 290, 279, 258, 271, 312, 354, 404, 402
Offset: 0

Views

Author

N. J. A. Sloane and Mira Bernstein, following a suggestion from Marc LeBrun, Apr 28 1994

Keywords

Comments

This sequence is of interest for group theory. The partitions counted by a(n) correspond to conjugacy classes of optimal order of the symmetric group of n elements: they have no fixed point, their order is the direct product of their cycle lengths and they are not contained in a subgroup of Sym_p for p < n. A123131 gives the maximum order (LCM) reachable by these partitions.

Examples

			The a(17) = 9 strict partitions into pairwise coprime parts that are greater than 1 are (17), (15,2), (14,3), (13,4), (12,5), (11,6), (10,7), (9,8), (7,5,3,2). - _Gus Wiseman_, Apr 14 2018
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Maple
    with(numtheory):
    b:= proc(n, i, s) option remember; local f;
          if n=0 then 1
        elif i<2 then 0
        else f:= factorset(i);
             b(n, i-1, select(x-> is(x is(x b(n, n, {}):
    seq(a(n), n=0..80);  # Alois P. Heinz, Mar 14 2012
  • Mathematica
    b[n_, i_, s_] := b[n, i, s] = Module[{f}, If[n == 0 || i == 1, 1, If[i<2, 0, f = FactorInteger[i][[All, 1]]; b[n, i-1, Select[s, #Jean-François Alcover, Feb 17 2014, after Alois P. Heinz *)
    Table[Length[Select[IntegerPartitions[n],FreeQ[#,1]&&(Length[#]===1||CoprimeQ@@#)&]],{n,20}] (* Gus Wiseman, Apr 14 2018 *)

Formula

a(n) = A051424(n) - A051424(n-1). - Vladeta Jovovic, Dec 11 2004

Extensions

More precise definition from Vladeta Jovovic, Dec 11 2004
More terms from Pab Ter (pabrlos2(AT)yahoo.com), Nov 13 2005

A302796 Squarefree numbers whose prime indices are relatively prime. Nonprime Heinz numbers of strict integer partitions with relatively prime parts.

Original entry on oeis.org

1, 2, 6, 10, 14, 15, 22, 26, 30, 33, 34, 35, 38, 42, 46, 51, 55, 58, 62, 66, 69, 70, 74, 77, 78, 82, 85, 86, 93, 94, 95, 102, 105, 106, 110, 114, 118, 119, 122, 123, 130, 134, 138, 141, 142, 143, 145, 146, 154, 155, 158, 161, 165, 166, 170, 174, 177, 178, 182
Offset: 1

Views

Author

Gus Wiseman, Apr 13 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. Two or more numbers are relatively prime if they have no common divisor other than 1. A single number is not considered relatively prime unless it is equal to 1.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			Sequence of terms together with their sets of prime indices begins:
01 : {}
02 : {1}
06 : {1,2}
10 : {1,3}
14 : {1,4}
15 : {2,3}
22 : {1,5}
26 : {1,6}
30 : {1,2,3}
33 : {2,5}
34 : {1,7}
35 : {3,4}
38 : {1,8}
42 : {1,2,4}
46 : {1,9}
51 : {2,7}
55 : {3,5}
58 : {1,10}
62 : {1,11}
66 : {1,2,5}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],Or[#===1,SquareFreeQ[#]&&GCD@@PrimePi/@FactorInteger[#][[All,1]]===1]&]
  • PARI
    isok(n) = {if (n == 1, return (1)); if (issquarefree(n), my(f = factor(n)); return (gcd(vector(#f~, k, primepi(f[k,1]))) == 1););} \\ Michel Marcus, Apr 13 2018

A304711 Heinz numbers of integer partitions whose distinct parts are pairwise coprime.

Original entry on oeis.org

2, 4, 6, 8, 10, 12, 14, 15, 16, 18, 20, 22, 24, 26, 28, 30, 32, 33, 34, 35, 36, 38, 40, 44, 45, 46, 48, 50, 51, 52, 54, 55, 56, 58, 60, 62, 64, 66, 68, 69, 70, 72, 74, 75, 76, 77, 80, 82, 85, 86, 88, 90, 92, 93, 94, 95, 96, 98, 99, 100, 102, 104, 106, 108, 110
Offset: 1

Views

Author

Gus Wiseman, May 17 2018

Keywords

Comments

Two parts are coprime if they have no common divisor greater than 1. For partitions of length 1 note that (1) is coprime but (x) is not coprime for x > 1.
First differs from A289509 at a(24) = 44, A289509(24) = 42.

Examples

			Sequence of all partitions whose distinct parts are pairwise coprime begins (1), (11), (21), (111), (31), (211), (41), (32), (1111), (221), (311), (51), (2111), (61), (411), (321), (11111), (52), (71), (43), (2211), (81), (3111).
		

Crossrefs

Programs

  • Mathematica
    Select[Range[200],CoprimeQ@@PrimePi/@FactorInteger[#][[All,1]]&]

A355737 Number of ways to choose a sequence of divisors, one of each prime index of n (with multiplicity), such that the result has no common divisor > 1.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 2, 1, 3, 4, 1, 1, 4, 1, 2, 4, 2, 1, 2, 3, 4, 7, 3, 1, 4, 1, 1, 4, 2, 6, 4, 1, 4, 6, 2, 1, 6, 1, 2, 8, 3, 1, 2, 5, 4, 4, 4, 1, 8, 4, 3, 5, 4, 1, 4, 1, 2, 10, 1, 6, 4, 1, 2, 6, 6, 1, 4, 1, 6, 8, 4, 6, 8, 1, 2, 15, 2, 1, 6, 4, 4
Offset: 1

Views

Author

Gus Wiseman, Jul 17 2022

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 a(2) = 1 through a(18) = 4 choices:
  1  1  11  1  11  1  111  11  11  1  111  1  11  11  1111  1  111
               12          12  13     112     12  13           112
                           21                 14  21           121
                                                  23           122
		

Crossrefs

Dominated by A355731, firsts A355732, primes A355741, prime-powers A355742.
For weakly increasing instead of coprime we have A355735, primes A355745.
Positions of first appearances are A355738.
For strict instead of coprime we have A355739, zeros A355740.
A000005 counts divisors.
A001221 counts distinct prime factors, with sum A001414.
A001222 counts prime factors with multiplicity.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798.
A120383 lists numbers divisible by all of their prime indices.
A289508 gives GCD of prime indices.
A289509 ranks relatively prime partitions, odd A302697, squarefree A302796.
A324850 lists numbers divisible by the product of their prime indices.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Select[Tuples[Divisors/@primeMS[n]],GCD@@#==1&]],{n,100}]

A337485 Number of pairwise coprime integer partitions of n with no 1's, where a singleton is not considered coprime unless it is (1).

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 2, 1, 2, 2, 4, 3, 5, 4, 4, 7, 8, 9, 10, 10, 9, 13, 17, 18, 17, 19, 19, 24, 29, 34, 33, 31, 31, 42, 42, 56, 55, 50, 54, 66, 77, 86, 86, 79, 81, 96, 124, 127, 126, 127, 126, 145, 181, 190, 184, 183, 192, 212, 262, 289, 278, 257, 270, 311
Offset: 0

Views

Author

Gus Wiseman, Sep 21 2020

Keywords

Comments

Such a partition is necessarily strict.
The Heinz numbers of these partitions are the intersection of A005408 (no 1's), A005117 (strict), and A302696 (coprime).

Examples

			The a(n) partitions for n = 5, 7, 12, 13, 16, 17, 18, 19 (A..H = 10..17):
  (3,2)  (4,3)  (7,5)    (7,6)  (9,7)    (9,8)      (B,7)    (A,9)
         (5,2)  (5,4,3)  (8,5)  (B,5)    (A,7)      (D,5)    (B,8)
                (7,3,2)  (9,4)  (D,3)    (B,6)      (7,6,5)  (C,7)
                         (A,3)  (7,5,4)  (C,5)      (8,7,3)  (D,6)
                         (B,2)  (8,5,3)  (D,4)      (9,5,4)  (E,5)
                                (9,5,2)  (E,3)      (9,7,2)  (F,4)
                                (B,3,2)  (F,2)      (B,4,3)  (G,3)
                                         (7,5,3,2)  (B,5,2)  (H,2)
                                                    (D,3,2)  (B,5,3)
                                                             (7,5,4,3)
		

Crossrefs

A005408 intersected with A302696 ranks these partitions.
A007359 considers all singletons to be coprime.
A327516 allows 1's, with non-strict version A305713.
A337452 is the relatively prime instead of pairwise coprime version, with non-strict version A302698.
A337563 is the restriction to partitions of length 3.
A002865 counts partitions with no 1's.
A078374 counts relatively prime strict partitions.
A200976 and A328673 count pairwise non-coprime partitions.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],!MemberQ[#,1]&&CoprimeQ@@#&]],{n,0,30}]

Formula

a(n) = A007359(n) - 1 for n > 1.

A302797 Squarefree numbers whose prime indices are pairwise coprime. Heinz numbers of strict integer partitions with pairwise coprime parts.

Original entry on oeis.org

1, 2, 6, 10, 14, 15, 22, 26, 30, 33, 34, 35, 38, 46, 51, 55, 58, 62, 66, 69, 70, 74, 77, 82, 85, 86, 93, 94, 95, 102, 106, 110, 118, 119, 122, 123, 134, 138, 141, 142, 143, 145, 146, 154, 155, 158, 161, 165, 166, 170, 177, 178, 186, 187, 190, 194, 201, 202
Offset: 1

Views

Author

Gus Wiseman, Apr 13 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. Two or more numbers are coprime if no pair of them has a common divisor other than 1. A single number is not considered coprime unless it is equal to 1.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			Sequence of terms together with their sets of prime indices begins:
01 : {}
02 : {1}
06 : {1,2}
10 : {1,3}
14 : {1,4}
15 : {2,3}
22 : {1,5}
26 : {1,6}
30 : {1,2,3}
33 : {2,5}
34 : {1,7}
35 : {3,4}
38 : {1,8}
46 : {1,9}
51 : {2,7}
55 : {3,5}
58 : {1,10}
62 : {1,11}
66 : {1,2,5}
69 : {2,9}
70 : {1,3,4}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],Or[#===1,SquareFreeQ[#]&&CoprimeQ@@PrimePi/@FactorInteger[#][[All,1]]]&]

A023023 Number of partitions of n into 3 unordered relatively prime parts.

Original entry on oeis.org

1, 1, 2, 2, 4, 4, 6, 6, 10, 8, 14, 12, 16, 16, 24, 18, 30, 24, 32, 30, 44, 32, 50, 42, 54, 48, 70, 48, 80, 64, 80, 72, 96, 72, 114, 90, 112, 96, 140, 96, 154, 120, 144, 132, 184, 128, 196, 150, 192, 168, 234, 162, 240, 192, 240, 210, 290, 192, 310, 240, 288, 256, 336, 240, 374
Offset: 3

Views

Author

Keywords

Examples

			From _Gus Wiseman_, Oct 08 2020: (Start)
The a(3) = 1 through a(13) = 14 triples (A = 10, B = 11):
  111   211   221   321   322   332   432   433   443   543   544
              311   411   331   431   441   532   533   552   553
                          421   521   522   541   542   651   643
                          511   611   531   631   551   732   652
                                      621   721   632   741   661
                                      711   811   641   831   733
                                                  722   921   742
                                                  731   A11   751
                                                  821         832
                                                  911         841
                                                              922
                                                              931
                                                              A21
                                                              B11
(End)
		

Crossrefs

A000741 is the ordered version.
A000837 counts these partitions of any length.
A001399(n-3) does not require relative primality.
A023022 is the 2-part version.
A101271 is the strict case.
A284825 counts the case that is also pairwise non-coprime.
A289509 intersected with A014612 gives the Heinz numbers.
A307719 is the pairwise coprime instead of relatively prime version.
A337599 is the pairwise non-coprime instead of relative prime version.
A008284 counts partitions by sum and length.
A078374 counts relatively prime strict partitions.
A337601 counts 3-part partitions whose distinct parts are pairwise coprime.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n,{3}],GCD@@#==1&]],{n,3,50}] (* Gus Wiseman, Oct 08 2020 *)

Formula

G.f. for the number of partitions of n into m unordered relatively prime parts is Sum(moebius(k)*x^(m*k)/Product(1-x^(i*k), i=1..m), k=1..infinity). - Vladeta Jovovic, Dec 21 2004
a(n) = (n^2/12)*Product_{prime p|n} (1 - 1/p^2) = A007434(n)/12 for n > 3 (proved by Mohamed El Bachraoui). [Jonathan Sondow, May 27 2009]
a(n) = Sum_{k=1..floor(n/3)} Sum_{i=k..floor((n-k)/2)} floor(1/gcd(i,k,n-i-k)). - Wesley Ivan Hurt, Jan 02 2021
Showing 1-10 of 34 results. Next