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

A305079 Number of connected components of the integer partition with Heinz number n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, May 24 2018

Keywords

Comments

First differs from |A305052(n)| at a(169) = 1, A305052(169) = 0.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
Given a finite multiset S of positive integers greater than one, let G(S) be the simple labeled graph with vertex set S and edges between any two vertices with a common divisor greater than 1. For example, G({6,14,15,35}) is a 4-cycle. If S is the integer partition with Heinz number n, a(n) is the number of connected components of G(S).

Examples

			The a(315) = 2 connected components of {2,2,3,4} are {{3},{2,2,4}}.
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    zsm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[Less@@#,GCD@@s[[#]]]>1&]},If[c=={},s,zsm[Sort[Append[Delete[s,List/@c[[1]]],LCM@@s[[c[[1]]]]]]]]];
    Table[Length[zsm[primeMS[n]]],{n,100}]
  • PARI
    zero_first_elem_and_connected_elems(ys) = { my(cs = List([ys[1]]), i=1); ys[1] = 0; while(i<=#cs, for(j=2,#ys,if(ys[j]&&(1!=gcd(cs[i],ys[j])), listput(cs,ys[j]); ys[j] = 0)); i++); (ys); };
    A007814(n) = valuation(n,2);
    A000265(n) = (n/2^A007814(n));
    A305079(n) = if(!(n%2),A007814(n)+A305079(A000265(n)), my(cs = apply(p -> primepi(p),factor(n)[,1]~), s=0); while(#cs, cs = select(c -> c, zero_first_elem_and_connected_elems(cs)); s++); (s)); \\ Antti Karttunen, Nov 10 2018

Formula

For all n, k > 0, we have a(2^n * k) = n + a(k).
For all x, y > 0, we have a(x * y) <= a(x) + a(y).
For x, y > 0 strongly coprime, we have a(x * y) = a(x) + a(y). Strongly coprime means every prime index of x is coprime to every prime index of y, where a prime index of n is a number m such that prime(m) divides n.
a(n) = A305501(A064989(n)) + A007814(n). - Antti Karttunen, Nov 10 2018

Extensions

Terms and Mathematica program corrected by Gus Wiseman, Nov 10 2018

A305504 Heinz numbers of integer partitions whose distinct parts plus 1 are connected.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 13, 16, 17, 19, 20, 22, 23, 25, 27, 29, 31, 32, 33, 34, 37, 40, 41, 43, 44, 46, 47, 49, 50, 53, 55, 57, 59, 61, 62, 64, 66, 67, 68, 71, 73, 79, 80, 81, 82, 83, 85, 88, 89, 92, 93, 94, 97, 99, 100, 101, 103, 107, 109, 110, 113, 115
Offset: 1

Views

Author

Gus Wiseman, Jun 03 2018

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
Given a finite set S of positive integers greater than one, let G(S) be the simple labeled graph with vertex set S and edges between any two vertices with a common divisor greater than 1. For example, G({6,14,15,35}) is a 4-cycle. A partition y is said to be connected if G(U(y + 1)) is a connected graph, where U(y + 1) is the set of distinct successors of the parts of y.
This is intended to be a cleaner form of A305078, where the treatment of empty multisets is arbitrary.

Examples

			The sequence of entries together with the corresponding twice-prime-factored multiset partitions (see A275024) begins:
   1: {}
   2: {{1}}
   3: {{2}}
   4: {{1},{1}}
   5: {{1,1}}
   7: {{3}}
   8: {{1},{1},{1}}
   9: {{2},{2}}
  10: {{1},{1,1}}
  11: {{1,2}}
  13: {{4}}
  16: {{1},{1},{1},{1}}
  17: {{1,1,1}}
  19: {{2,2}}
  20: {{1},{1},{1,1}}
  22: {{1},{1,2}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    zsm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[Less@@#,GCD@@s[[#]]]>1&]},If[c=={},s,zsm[Union[Append[Delete[s,List/@c[[1]]],LCM@@s[[c[[1]]]]]]]]];
    Select[Range[300],Length[zsm[primeMS[#]+1]]<=1&]
Showing 1-2 of 2 results.