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.

A328220 Number of strict integer partitions of n with no pair of consecutive parts relatively prime.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 1, 2, 2, 3, 1, 5, 1, 5, 4, 6, 3, 10, 3, 11, 7, 12, 3, 19, 5, 18, 12, 23, 9, 36, 11, 33, 21, 40, 20, 58, 19, 58, 35, 70, 31, 98, 36, 101, 65, 112, 56, 155, 64, 164, 97, 188, 88, 250, 112, 256, 157, 293, 145, 392, 163, 399, 241, 461, 242
Offset: 0

Views

Author

Gus Wiseman, Oct 14 2019

Keywords

Examples

			The a(2) = 1 through a(20) = 11 partitions (A..K = 10..20):
  2  3  4  5  6   7  8   9   A   B  C    D  E    F   G    H    I    J    K
              42     62  63  64     84      86   96  A6   863  A8   964  C8
                             82     93      A4   A5  C4   962  C6   A63  E6
                                    A2      C2   C3  E2        E4        F5
                                    642     842      862       F3        G4
                                                     A42       G2        I2
                                                               864       A64
                                                               963       A82
                                                               A62       C62
                                                               C42       E42
                                                                         8642

Links

Crossrefs

The non-strict case is A328187.
Partitions with all consecutive parts relatively prime are A328172, with strict case A328188.
Strict partitions with relatively prime parts are A078374.
Partitions with no consecutive divisibilities are A328171.
Cf. A000837, A018783, A070211, A289509, A318980, A325545, A328170, A328221, A328335, A328336.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&!MatchQ[#,{_,x_,y_,_}/;GCD[x,y]==1]&]],{n,0,30}]

A328170 Number of integer partitions of n whose parts minus 1 are relatively prime.

Original entry on oeis.org

0, 0, 1, 1, 2, 3, 5, 8, 12, 18, 27, 38, 53, 74, 102, 137, 184, 241, 317, 413, 536, 687, 880, 1112, 1405, 1765, 2215, 2755, 3424, 4229, 5216, 6402, 7847, 9572, 11662, 14148, 17139, 20688, 24940, 29971, 35969, 43044, 51438, 61311, 72985, 86678, 102807, 121675
Offset: 0

Views

Author

Gus Wiseman, Oct 09 2019

Keywords

Comments

A partition is relatively prime if the GCD of its parts is 1. Zeros are ignored when computing GCD, and the empty set has GCD 0.

Examples

			The a(2) = 1 through a(9) = 18 partitions:
  (2)  (21)  (22)   (32)    (42)     (43)      (62)       (54)
             (211)  (221)   (222)    (52)      (332)      (63)
                    (2111)  (321)    (322)     (422)      (72)
                            (2211)   (421)     (431)      (432)
                            (21111)  (2221)    (521)      (522)
                                     (3211)    (2222)     (621)
                                     (22111)   (3221)     (3222)
                                     (211111)  (4211)     (3321)
                                               (22211)    (4221)
                                               (32111)    (4311)
                                               (221111)   (5211)
                                               (2111111)  (22221)
                                                          (32211)
                                                          (42111)
                                                          (222111)
                                                          (321111)
                                                          (2211111)
                                                          (21111111)

Links

Crossrefs

The Heinz numbers of these partitions are given by A328168.
Partitions whose parts are relatively prime are A000837.
Partitions whose parts plus 1 are relatively prime are A318980.
The GCD of the prime indices of n, all minus 1, is A328167(n).
Cf. A007359, A018783, A258409, A289509, A328163, A328164.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],GCD@@(#-1)==1&]],{n,0,30}]
  • PARI
    seq(n)=Vec(sum(d=1, n-1, moebius(d)*(-1/(1-x) + 1/prod(k=0, n\d, 1 - x*x^(k*d) + O(x*x^n)))), -(n+1)) \\ Andrew Howroyd, Oct 17 2019

Formula

G.f.: Sum_{d>=1} mu(d)*(-1/(1-x) + 1/(Prod_{k>=0} 1 - x^(k*d + 1))). - Andrew Howroyd, Oct 17 2019

A328169 GCD of the prime indices of n, all plus 1.

Original entry on oeis.org

0, 2, 3, 2, 4, 1, 5, 2, 3, 2, 6, 1, 7, 1, 1, 2, 8, 1, 9, 2, 1, 2, 10, 1, 4, 1, 3, 1, 11, 1, 12, 2, 3, 2, 1, 1, 13, 1, 1, 2, 14, 1, 15, 2, 1, 2, 16, 1, 5, 2, 1, 1, 17, 1, 2, 1, 3, 1, 18, 1, 19, 2, 1, 2, 1, 1, 20, 2, 1, 1, 21, 1, 22, 1, 1, 1, 1, 1, 23, 2, 3, 2
Offset: 1

Views

Author

Gus Wiseman, Oct 09 2019

Keywords

Comments

Zeros are ignored when computing GCD, and the empty set has GCD 0.
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

			85 has prime indices {3,7}, so a(85) = GCD(4,8) = 4.

Crossrefs

Positions of 0's and 1's are A318981.
Positions of records (first appearances) appear to be A116974.
The GCD of the prime indices of n, all minus 1, is A328167(n).
The LCM of the prime indices of n, all plus 1, is A328219(n).
Partitions whose parts plus 1 are relatively prime are A318980.
Cf. A000837, A056239, A112798, A258409, A289508, A289509, A290103, A328168.

Programs

  • Mathematica
    Table[GCD@@(PrimePi/@First/@If[n==1,{},FactorInteger[n]]+1),{n,100}]

Formula

a(n) = A289508(A003961(n)).

A318981 Numbers whose prime indices plus 1 are relatively prime.

Original entry on oeis.org

1, 6, 12, 14, 15, 18, 21, 24, 26, 28, 30, 35, 36, 38, 39, 42, 45, 48, 51, 52, 54, 56, 58, 60, 63, 65, 66, 69, 70, 72, 74, 75, 76, 77, 78, 84, 86, 87, 90, 91, 95, 96, 98, 102, 104, 105, 106, 108, 111, 112, 114, 116, 117, 119, 120, 122, 123, 126, 130, 132, 133
Offset: 1

Views

Author

Gus Wiseman, Sep 06 2018

Keywords

Comments

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

Examples

			The sequence of integer partitions whose Heinz numbers are in the sequence begins: (), (21), (211), (41), (32), (221), (42), (2111), (61), (411), (321), (43), (2211), (81), (62), (421), (322), (21111).

Crossrefs

Cf. A000837, A007359, A018783, A281116, A289508, A289509, A318978, A318980.

Programs

  • Mathematica
    Select[Range[100],GCD@@(PrimePi/@FactorInteger[#][[All,1]]+1)==1&]

A328219 LCM of the prime indices of n, all plus 1.

Original entry on oeis.org

1, 2, 3, 2, 4, 6, 5, 2, 3, 4, 6, 6, 7, 10, 12, 2, 8, 6, 9, 4, 15, 6, 10, 6, 4, 14, 3, 10, 11, 12, 12, 2, 6, 8, 20, 6, 13, 18, 21, 4, 14, 30, 15, 6, 12, 10, 16, 6, 5, 4, 24, 14, 17, 6, 12, 10, 9, 22, 18, 12, 19, 12, 15, 2, 28, 6, 20, 8, 30, 20, 21, 6, 22, 26
Offset: 1

Views

Author

Gus Wiseman, Oct 16 2019

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.

Crossrefs

Sorted positions of first appearances are A328451.
LCM of prime indices is A290103.
LCM of prime indices minus 1 is A328456.
GCD of prime indices plus 1 is A328169.
Partitions whose parts plus 1 are relatively prime are A318980.
Numbers whose prime indices plus 1 are relatively prime are A318981,
Cf. A003961, A056239, A112798, A258409, A289508, A289509, A328167, A328168.

Programs

  • Mathematica
    Table[If[n==1,1,LCM@@(PrimePi/@First/@FactorInteger[n]+1)],{n,100}]
  • PARI
    A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961
A290103(n) = lcm(apply(p->primepi(p),factor(n)[,1]));
A328219(n) = A290103(A003961(n)); \\ Antti Karttunen, Oct 18 2019

Formula

a(n) = A290103(A003961(n)).
If n = A000040(i_1) * ... * A000040(i_k), then a(n) = lcm(1+i_1,...,1+i_k).

A328451 Sorted positions of first appearances in A328219, where if n = A000040(i_1) * ... * A000040(i_k), then A328219(n) = LCM(1+i_1,...,1+i_k).

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 13, 14, 15, 17, 19, 21, 26, 29, 35, 37, 38, 39, 42, 47, 51, 53, 58, 61, 65, 74, 78, 79, 87, 89, 91, 95, 101, 105, 106, 107, 111, 113, 119, 122, 127, 133, 141, 145, 151, 158, 159, 173, 174, 178, 181, 182, 183, 185, 195, 199, 202, 203, 214, 221
Offset: 1

Views

Author

Gus Wiseman, Oct 17 2019

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.
Indices of 1's in the ordinal transform of A328219. - Antti Karttunen, Oct 18 2019

Examples

			The sequence of terms together with their prime indices begins:
   1: {}
   2: {1}
   3: {2}
   5: {3}
   6: {1,2}
   7: {4}
  13: {6}
  14: {1,4}
  15: {2,3}
  17: {7}
  19: {8}
  21: {2,4}
  26: {1,6}
  29: {10}
  35: {3,4}
  37: {12}
  38: {1,8}
  39: {2,6}
  42: {1,2,4}
  47: {15}

Crossrefs

A subsequence of A005117.
Sorted positions of first appearances in A328219.
The GCD of the prime indices of n, all plus 1, is A328169(n).
The LCM of the prime indices of n, all minus 1, is A328456(n).
Partitions whose parts plus 1 are relatively prime are A318980.
Numbers whose prime indices plus 1 are relatively prime are A318981.
Cf. A001222, A056239, A112798, A289508, A289509, A290103, A328163, A328167, A328168, A328170.

Programs

  • Mathematica
    dav=Table[If[n==1,1,LCM@@(PrimePi/@First/@FactorInteger[n]+1)],{n,100}];
Table[Position[dav,i][[1,1]],{i,dav//.{A___,x_,B___,x_,C___}:>{A,x,B,C}}]
  • PARI
    up_to = 1024;
ordinal_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), pt); for(i=1, length(invec), if(mapisdefined(om,invec[i]), pt = mapget(om, invec[i]), pt = 0); outvec[i] = (1+pt); mapput(om,invec[i],(1+pt))); outvec; };
A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961
A290103(n) = lcm(apply(p->primepi(p),factor(n)[,1]));
A328219(n) = A290103(A003961(n));
vord_trans = ordinal_transform(vector(up_to,n,A328219(n)));
for(n=1,up_to,if(1==vord_trans[n], print1(n,", "))); \\ Antti Karttunen, Oct 18 2019
Showing 1-6 of 6 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.