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

A305829 Factor n into distinct Fermi-Dirac primes (A050376), normalize by replacing every instance of the k-th Fermi-Dirac prime with k, then multiply everything together.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jun 10 2018

Keywords

Comments

Let f(n) = A050376(n) be the n-th Fermi-Dirac prime. Every positive integer n has a unique factorization of the form n = f(s_1)*...*f(s_k) where the s_i are strictly increasing positive integers. Then a(n) = s_1 * ... * s_k.
Multiplicative because for coprime m and n the Fermi-Dirac factorizations of m and n are disjoint and their union is the Fermi-Dirac factorization of m * n. - Andrew Howroyd, Aug 02 2018

Crossrefs

Programs

  • Mathematica
    nn=100;
    FDfactor[n_]:=If[n===1,{},Sort[Join@@Cases[FactorInteger[n],{p_,k_}:>Power[p,Cases[Position[IntegerDigits[k,2]//Reverse,1],{m_}->2^(m-1)]]]]];
    FDprimeList=Array[FDfactor,nn,1,Union];FDrules=MapIndexed[(#1->#2[[1]])&,FDprimeList];
    Table[Times@@(FDfactor[n]/.FDrules),{n,nn}]
  • PARI
    \\ here isfd is membership test for A050376.
    isfd(n)={my(e=isprimepower(n)); e && e == 1<Andrew Howroyd, Aug 02 2018

A305832 Number of connected components of the n-th FDH set-system.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 2, 2, 2, 1, 1, 2, 1, 3, 1, 3, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 2, 2, 3, 2, 2, 1, 3, 1, 2, 2, 2, 1, 3, 1, 2, 2, 3, 1, 2, 1, 2, 2, 2, 2, 2, 1, 2, 1, 2, 1, 3, 2, 2, 1
Offset: 1

Views

Author

Gus Wiseman, Jun 10 2018

Keywords

Comments

Let f(n) = A050376(n) be the n-th Fermi-Dirac prime. Every positive integer n has a unique factorization of the form n = f(s_1)*...*f(s_k) where the s_i are strictly increasing positive integers. The n-th FDH set-system is obtained by repeating this factorization on each index s_i.

Examples

			Let f = A050376. The FD-factorization of 765 is 5*9*17 or f(4)*f(6)*f(10) = f(4)*f(2*3)*f(2*5) with connected components {{{4}},{{2,3},{2,5}}}, so a(765) = 2.
		

Crossrefs

Programs

  • Mathematica
    FDfactor[n_]:=If[n===1,{},Sort[Join@@Cases[FactorInteger[n],{p_,k_}:>Power[p,Cases[Position[IntegerDigits[k,2]//Reverse,1],{m_}->2^(m-1)]]]]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>1]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    nn=100;FDprimeList=Array[FDfactor,nn,1,Union];FDrules=MapIndexed[(#1->#2[[1]])&,FDprimeList];
    Table[Length[csm[FDfactor[#]/.FDrules&/@(FDfactor[n]/.FDrules)]],{n,nn}]

A316268 FDH numbers of connected strict integer partitions.

Original entry on oeis.org

2, 3, 4, 5, 7, 9, 11, 13, 15, 16, 17, 19, 23, 25, 27, 29, 31, 36, 37, 39, 41, 43, 45, 47, 49, 51, 53, 59, 61, 64, 65, 67, 69, 71, 73, 79, 81, 83, 85, 87, 89, 92, 97, 101, 103, 107, 108, 109, 111, 113, 115, 117, 119, 121, 124, 127, 129, 131, 135, 137, 139, 144
Offset: 1

Views

Author

Gus Wiseman, Jun 28 2018

Keywords

Comments

Let f(n) = A050376(n) be the n-th Fermi-Dirac prime. The FDH number of a strict integer partition (y_1,...,y_k) is f(y_1)*...*f(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 set or strict partition S is said to be connected if G(S) is a connected graph.

Examples

			Sequence of connected strict integer partitions begins (1), (2), (3), (4), (5), (6), (7), (8), (4,2), (9), (10), (11), (12), (13), (6,2).
		

Crossrefs

Programs

  • Mathematica
    nn=100;
    FDfactor[n_]:=If[n===1,{},Sort[Join@@Cases[FactorInteger[n],{p_,k_}:>Power[p,Cases[Position[IntegerDigits[k,2]//Reverse,1],{m_}->2^(m-1)]]]]];
    FDrules=MapIndexed[(#1->#2[[1]])&,Array[FDfactor,nn,1,Union]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>1]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[nn],Length[csm[primeMS/@(FDfactor[#]/.FDrules)]]==1&]
Showing 1-3 of 3 results.