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.

Previous Showing 11-20 of 29 results. Next

A138904 Number of rotational symmetries in the binary expansion of a number.

Original entry on oeis.org

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

Views

Author

Max Sills, Apr 03 2008, Apr 04 2008

Keywords

Comments

Mersenne numbers of form (2^n - 1) have n rotational symmetries.
For prime length binary expansions these are the only nontrivial symmetries.
For composite length expansions it seems that when the number of symmetries is nontrivial it is equal to a factor of the length. We're working on an explicit formula.
Discovered in the context of random circulant matrices, examining if there's a correlation between degrees of freedom and number of symmetries in the first row.
When combined with A138954, these two sequences should give a full account of the number of redundant rows in a circulant square matrix with at most two distinct values, where a(n) is the encoding of the first row of the matrix into binary such that value a = 1 and value b = 0.
Discovered on the night of Apr 02, 2008 by Maxwell Sills and Gary Doran.
Conjecture: For binary expansions of length n, there are d(n) distinct values that will show up as symmetries, where d is the divisor function. The symmetry values will be precisely the divisors of n.
Example: for binary expansions of length 12, one sees that d(12) = 6 distinct values show up as symmetries (1, 2, 3, 4, 6, 12).
Conjecture: For numbers whose binary expansion has length n which has proper divisors which are all coprime: There will be only one number of length n with n symmetries. That number is 2^n - 1. For each proper divisor d (excluding 1), you can generate all numbers of length n that have n/d symmetries like so: (2^0 + 2^d + 2^2d ... 2^(n-d)) * a, where 2^(d-1) <= a < (2^d) - 1. The rest of the expansions of length n will have only the trivial symmetry.
Also the number of rotational symmetries of the n-th composition in standard order (graded reverse-lexicographic). This composition (row n of A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of n, prepending 0, taking first differences, and reversing again. - Gus Wiseman, Apr 19 2020
From Gus Wiseman, Apr 19 2020: (Start)
Aperiodic compositions are counted by A000740.
Aperiodic binary words are counted by A027375.
The orderless period of prime indices is A052409.
Numbers whose binary expansion is periodic are A121016.
Periodic compositions are counted by A178472.
Period of binary expansion is A302291.
Compositions by sum and number of distinct rotations are A333941.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Necklaces are A065609.
- Sum is A070939.
- Runs are counted by A124767.
- Strict compositions are A233564.
- Constant compositions are A272919.
- Lyndon compositions are A275692.
- Co-Lyndon compositions are A326774.
- Aperiodic compositions are A328594.
- Reversed co-necklaces are A328595.
- Rotational period is A333632.
- Co-necklaces are A333764.
- Reversed necklaces are A333943.
(End).

Examples

			a(10) = 2 because the binary expansion of 10 is 1010 and it has two rotational symmetries (including identity).
		

Crossrefs

Programs

  • Mathematica
    Table[IntegerLength[n,2]/Length[Union[Array[RotateRight[IntegerDigits[n,2],#]&,IntegerLength[n,2]]]],{n,100}] (* Gus Wiseman, Apr 19 2020 *)

Formula

a(n) = A070939(n)/A302291(n) = A000120(n)/A333632(n). - Gus Wiseman, Apr 19 2020

A329140 Numbers whose prime signature is a periodic word.

Original entry on oeis.org

6, 10, 14, 15, 21, 22, 26, 30, 33, 34, 35, 36, 38, 39, 42, 46, 51, 55, 57, 58, 62, 65, 66, 69, 70, 74, 77, 78, 82, 85, 86, 87, 91, 93, 94, 95, 100, 102, 105, 106, 110, 111, 114, 115, 118, 119, 122, 123, 129, 130, 133, 134, 138, 141, 142, 143, 145, 146, 154
Offset: 1

Views

Author

Gus Wiseman, Nov 09 2019

Keywords

Comments

First differs from A182853 in having 2100 = 2^2 * 3^1 * 5^2 * 7^1.
A number's prime signature (A124010) is the sequence of positive exponents in its prime factorization.
A sequence is aperiodic if its cyclic rotations are all different.

Examples

			The sequence of terms together with their prime signatures begins:
   6: (1,1)
  10: (1,1)
  14: (1,1)
  15: (1,1)
  21: (1,1)
  22: (1,1)
  26: (1,1)
  30: (1,1,1)
  33: (1,1)
  34: (1,1)
  35: (1,1)
  36: (2,2)
  38: (1,1)
  39: (1,1)
  42: (1,1,1)
  46: (1,1)
  51: (1,1)
  55: (1,1)
  57: (1,1)
  58: (1,1)
		

Crossrefs

Complement of A329139.
Periodic compositions are A178472.
Periodic binary words are A152061.
Numbers whose binary expansion is periodic are A121016.
Numbers whose prime signature is a Lyndon word are A329131.
Numbers whose prime signature is a necklace are A329138.

Programs

  • Mathematica
    aperQ[q_]:=Array[RotateRight[q,#1]&,Length[q],1,UnsameQ];
    Select[Range[100],!aperQ[Last/@FactorInteger[#]]&]

A302291 a(n) is the period of the binary expansion of n.

Original entry on oeis.org

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

Views

Author

Rémy Sigrist, Apr 04 2018

Keywords

Comments

Zero is assumed to be represented as 0; otherwise, leading zeros are ignored.
See A302295 for the variant where leading zeros are allowed.

Examples

			The first terms, alongside the binary expansion of n with periodic part in parentheses, are:
  n  a(n)    bin(n)
  -- ----    ------
   0    1    (0)
   1    1    (1)
   2    2    (10)
   3    1    (1)(1)
   4    3    (100)
   5    3    (101)
   6    3    (110)
   7    1    (1)(1)(1)
   8    4    (1000)
   9    4    (1001)
  10    2    (10)(10)
  11    4    (1011)
  12    4    (1100)
  13    4    (1101)
  14    4    (1110)
  15    1    (1)(1)(1)(1)
  16    5    (10000)
  17    5    (10001)
  18    5    (10010)
  19    5    (10011)
  20    5    (10100)
		

Crossrefs

Aperiodic compositions are counted by A000740.
Aperiodic binary words are counted by A027375.
The orderless period of prime indices is A052409.
Numbers whose binary expansion is periodic are A121016.
Periodic compositions are counted by A178472.
Numbers whose prime signature is aperiodic are A329139.
Compositions by number of distinct rotations are A333941.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Necklaces are A065609.
- Sum is A070939.
- Runs are counted by A124767.
- Rotational symmetries are counted by A138904.
- Strict compositions are A233564.
- Constant compositions are A272919.
- Lyndon compositions are A275692.
- Co-Lyndon compositions are A326774.
- Aperiodic compositions are A328594.
- Rotational period is A333632.
- Co-necklaces are A333764.
- Reversed necklaces are A333943.

Programs

  • Mathematica
    Table[If[n==0,1,Length[Union[Array[RotateRight[IntegerDigits[n,2],#]&,IntegerLength[n,2]]]]],{n,0,50}] (* Gus Wiseman, Apr 19 2020 *)
  • PARI
    a(n) = my (l=max(1, #binary(n))); fordiv (l, w, if (#Set(digits(n, 2^w))<=1, return (w)))

Formula

a(n) = A070939(n) / A138904(n).
a(2^n) = n + 1 for any n >= 0.
a(2^n - 1) = 1 for any n >= 0.
a(A020330(n)) = a(n) for any n > 0.

A328671 Numbers whose binary indices are relatively prime and pairwise indivisible.

Original entry on oeis.org

1, 6, 12, 18, 20, 22, 24, 28, 48, 56, 66, 68, 70, 72, 76, 80, 82, 84, 86, 88, 92, 96, 104, 112, 120, 132, 144, 148, 176, 192, 196, 208, 212, 224, 240, 258, 264, 272, 274, 280, 296, 304, 312, 320, 322, 328, 336, 338, 344, 352, 360, 368, 376, 384, 400, 416, 432
Offset: 1

Views

Author

Gus Wiseman, Oct 29 2019

Keywords

Comments

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.

Examples

			The sequence of terms together with their binary expansions and binary indices begins:
    1:         1 ~ {1}
    6:       110 ~ {2,3}
   12:      1100 ~ {3,4}
   18:     10010 ~ {2,5}
   20:     10100 ~ {3,5}
   22:     10110 ~ {2,3,5}
   24:     11000 ~ {4,5}
   28:     11100 ~ {3,4,5}
   48:    110000 ~ {5,6}
   56:    111000 ~ {4,5,6}
   66:   1000010 ~ {2,7}
   68:   1000100 ~ {3,7}
   70:   1000110 ~ {2,3,7}
   72:   1001000 ~ {4,7}
   76:   1001100 ~ {3,4,7}
   80:   1010000 ~ {5,7}
   82:   1010010 ~ {2,5,7}
   84:   1010100 ~ {3,5,7}
   86:   1010110 ~ {2,3,5,7}
   88:   1011000 ~ {4,5,7}
		

Crossrefs

The version for prime indices (instead of binary indices) is A328677.
Numbers whose binary indices are relatively prime are A291166.
Numbers whose distinct prime indices are pairwise indivisible are A316476.
BII-numbers of antichains are A326704.
Relatively prime partitions whose distinct parts are pairwise indivisible are A328676, with strict case A328678.

Programs

  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    stableQ[u_,Q_]:=!Apply[Or,Outer[#1=!=#2&&Q[#1,#2]&,u,u,1],{0,1}];
    Select[Range[100],GCD@@bpe[#]==1&&stableQ[bpe[#],Divisible]&]

Formula

Intersection of A291166 with A326704.

A329134 Numbers whose differences of prime indices are a periodic word.

Original entry on oeis.org

8, 16, 27, 30, 32, 64, 81, 105, 110, 125, 128, 180, 210, 238, 243, 256, 273, 343, 385, 450, 506, 512, 625, 627, 729, 806, 935, 1001, 1024, 1080, 1100, 1131, 1155, 1331, 1394, 1495, 1575, 1729, 1786, 1870, 1887, 2048, 2187, 2197, 2310, 2401, 2431, 2451, 2635
Offset: 1

Views

Author

Gus Wiseman, Nov 09 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.
A sequence is periodic if its cyclic rotations are not all different.

Examples

			The sequence of terms together with their differences of prime indices begins:
     8: (0,0)
    16: (0,0,0)
    27: (0,0)
    30: (1,1)
    32: (0,0,0,0)
    64: (0,0,0,0,0)
    81: (0,0,0)
   105: (1,1)
   110: (2,2)
   125: (0,0)
   128: (0,0,0,0,0,0)
   180: (0,1,0,1)
   210: (1,1,1)
   238: (3,3)
   243: (0,0,0,0)
   256: (0,0,0,0,0,0,0)
   273: (2,2)
   343: (0,0)
   385: (1,1)
   450: (1,0,1,0)
		

Crossrefs

Complement of A329135.
These are the Heinz numbers of the partitions counted by A329144.
Periodic binary words are A152061.
Periodic compositions are A178472.
Numbers whose binary expansion is periodic are A121016.
Numbers whose prime signature is periodic are A329140.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    aperQ[q_]:=Array[RotateRight[q,#1]&,Length[q],1,UnsameQ];
    Select[Range[10000],!aperQ[Differences[primeMS[#]]]&]

A329132 Numbers whose augmented differences of prime indices are a periodic sequence.

Original entry on oeis.org

4, 8, 15, 16, 32, 55, 64, 90, 105, 119, 128, 225, 253, 256, 403, 512, 540, 550, 697, 893, 935, 1024, 1155, 1350, 1357, 1666, 1943, 2048, 2263, 3025, 3071, 3150, 3240, 3375, 3451, 3927, 3977, 4096, 4429, 5123, 5500, 5566, 6731, 7735, 8083, 8100, 8192, 9089
Offset: 1

Views

Author

Gus Wiseman, Nov 06 2019

Keywords

Comments

The augmented differences aug(y) of an integer partition y of length k are given by aug(y)i = y_i - y{i + 1} + 1 if i < k and aug(y)_k = y_k. For example, aug(6,5,5,3,3,3) = (2,1,3,1,1,3).
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.
A sequence is periodic if its cyclic rotations are not all different.

Examples

			The sequence of terms together with their augmented differences of prime indices begins:
     4: (1,1)
     8: (1,1,1)
    15: (2,2)
    16: (1,1,1,1)
    32: (1,1,1,1,1)
    55: (3,3)
    64: (1,1,1,1,1,1)
    90: (2,1,2,1)
   105: (2,2,2)
   119: (4,4)
   128: (1,1,1,1,1,1,1)
   225: (1,2,1,2)
   253: (5,5)
   256: (1,1,1,1,1,1,1,1)
   403: (6,6)
   512: (1,1,1,1,1,1,1,1,1)
   540: (2,1,1,2,1,1)
   550: (3,1,3,1)
   697: (7,7)
   893: (8,8)
		

Crossrefs

Complement of A329133.
These are the Heinz numbers of the partitions counted by A329143.
Periodic binary words are A152061.
Periodic compositions are A178472.
Numbers whose binary expansion is periodic are A121016.
Numbers whose prime signature is periodic are A329140.
Numbers whose differences of prime indices are periodic are A329134.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    aperQ[q_]:=Array[RotateRight[q,#1]&,Length[q],1,UnsameQ];
    aug[y_]:=Table[If[i
    				

A328607 Numbers whose reversed binary expansion, without the most significant digit, is a necklace.

Original entry on oeis.org

0, 1, 2, 3, 4, 6, 7, 8, 12, 14, 15, 16, 24, 26, 28, 30, 31, 32, 48, 52, 56, 58, 60, 62, 63, 64, 96, 100, 104, 106, 108, 112, 116, 118, 120, 122, 124, 126, 127, 128, 192, 200, 208, 212, 216, 220, 224, 228, 232, 234, 236, 240, 244, 246, 248, 250, 252, 254, 255
Offset: 0

Views

Author

Gus Wiseman, Oct 30 2019

Keywords

Comments

Offset is 0 to be consistent with A257250.
A necklace is a finite sequence that is lexicographically minimal among all of its cyclic rotations.

Examples

			The sequence of terms together with their binary expansions and binary indices begins:
    0:        0 ~ {}
    1:        1 ~ {1}
    2:       10 ~ {2}
    3:       11 ~ {1,2}
    4:      100 ~ {3}
    6:      110 ~ {2,3}
    7:      111 ~ {1,2,3}
    8:     1000 ~ {4}
   12:     1100 ~ {3,4}
   14:     1110 ~ {2,3,4}
   15:     1111 ~ {1,2,3,4}
   16:    10000 ~ {5}
   24:    11000 ~ {4,5}
   26:    11010 ~ {2,4,5}
   28:    11100 ~ {3,4,5}
   30:    11110 ~ {2,3,4,5}
   31:    11111 ~ {1,2,3,4,5}
   32:   100000 ~ {6}
   48:   110000 ~ {5,6}
   52:   110100 ~ {3,5,6}
		

Crossrefs

The dual non-reversed version is A257250.
The dual non-reversed version involving all digits is A065609.
The version involving all digits is A328595.
The non-reversed version is A328668.
Binary necklaces are A000031.

Programs

  • Mathematica
    neckQ[q_]:=Array[OrderedQ[{q,RotateRight[q,#]}]&,Length[q]-1,1,And];
    Select[Range[0,100],#<=1||neckQ[Reverse[Rest[IntegerDigits[#,2]]]]&]

A328668 Numbers whose binary expansion, without the most significant digit, is a necklace.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 7, 8, 9, 11, 15, 16, 17, 19, 21, 23, 31, 32, 33, 35, 37, 39, 43, 47, 63, 64, 65, 67, 69, 71, 73, 75, 77, 79, 85, 87, 91, 95, 127, 128, 129, 131, 133, 135, 137, 139, 141, 143, 147, 149, 151, 155, 157, 159, 171, 175, 183, 191, 255, 256, 257
Offset: 0

Views

Author

Gus Wiseman, Oct 26 2019

Keywords

Comments

Offset is 0 to be consistent with A257250.
A necklace is a finite sequence that is lexicographically minimal among all of its cyclic rotations.

Examples

			The sequence of terms together with their binary expansions and binary indices begins:
   0:       0 ~ {}
   1:       1 ~ {1}
   2:      10 ~ {2}
   3:      11 ~ {1,2}
   4:     100 ~ {3}
   5:     101 ~ {1,3}
   7:     111 ~ {1,2,3}
   8:    1000 ~ {4}
   9:    1001 ~ {1,4}
  11:    1011 ~ {1,2,4}
  15:    1111 ~ {1,2,3,4}
  16:   10000 ~ {5}
  17:   10001 ~ {1,5}
  19:   10011 ~ {1,2,5}
  21:   10101 ~ {1,3,5}
  23:   10111 ~ {1,2,3,5}
  31:   11111 ~ {1,2,3,4,5}
  32:  100000 ~ {6}
  33:  100001 ~ {1,6}
  35:  100011 ~ {1,2,6}
		

Crossrefs

The dual version is A257250.
The version involving all digits, taken in reverse, is A328595.
The reversed version is A328607.
Binary necklaces are A000031.
Necklace compositions are A008965.

Programs

  • Mathematica
    neckQ[q_]:=Array[OrderedQ[{q,RotateRight[q,#]}]&,Length[q]-1,1,And];
    Select[Range[0,100],#<=1||neckQ[Rest[IntegerDigits[#,2]]]&]

A329133 Numbers whose augmented differences of prime indices are an aperiodic sequence.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 9, 10, 11, 12, 13, 14, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74
Offset: 1

Views

Author

Gus Wiseman, Nov 09 2019

Keywords

Comments

The augmented differences aug(y) of an integer partition y of length k are given by aug(y)i = y_i - y{i + 1} + 1 if i < k and aug(y)_k = y_k. For example, aug(6,5,5,3,3,3) = (2,1,3,1,1,3).
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.
A finite sequence is aperiodic if its cyclic rotations are all different.

Examples

			The sequence of terms together with their augmented differences of prime indices begins:
    1: ()
    2: (1)
    3: (2)
    5: (3)
    6: (2,1)
    7: (4)
    9: (1,2)
   10: (3,1)
   11: (5)
   12: (2,1,1)
   13: (6)
   14: (4,1)
   17: (7)
   18: (1,2,1)
   19: (8)
   20: (3,1,1)
   21: (3,2)
   22: (5,1)
   23: (9)
   24: (2,1,1,1)
		

Crossrefs

Complement of A329132.
These are the Heinz numbers of the partitions counted by A329136.
Aperiodic binary words are A027375.
Aperiodic compositions are A000740.
Numbers whose binary expansion is aperiodic are A328594.
Numbers whose prime signature is aperiodic are A329139.
Numbers whose differences of prime indices are aperiodic are A329135.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    aperQ[q_]:=Array[RotateRight[q,#1]&,Length[q],1,UnsameQ];
    aug[y_]:=Table[If[i
    				

A333941 Triangle read by rows where T(n,k) is the number of compositions of n with rotational period k.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 0, 2, 2, 0, 0, 3, 2, 3, 0, 0, 2, 4, 6, 4, 0, 0, 4, 6, 9, 8, 5, 0, 0, 2, 6, 15, 20, 15, 6, 0, 0, 4, 8, 24, 32, 35, 18, 7, 0, 0, 3, 10, 27, 56, 70, 54, 28, 8, 0, 0, 4, 12, 42, 84, 125, 120, 84, 32, 9, 0, 0, 2, 10, 45, 120, 210, 252, 210, 120, 45, 10, 0
Offset: 0

Views

Author

Gus Wiseman, Apr 16 2020

Keywords

Comments

A composition of n is a finite sequence of positive integers summing to n.

Examples

			Triangle begins:
   1
   0   1
   0   2   0
   0   2   2   0
   0   3   2   3   0
   0   2   4   6   4   0
   0   4   6   9   8   5   0
   0   2   6  15  20  15   6   0
   0   4   8  24  32  35  18   7   0
   0   3  10  27  56  70  54  28   8   0
   0   4  12  42  84 125 120  84  32   9   0
   0   2  10  45 120 210 252 210 120  45  10   0
   0   6  18  66 168 335 450 462 320 162  50  11   0
Row n = 6 counts the following compositions (empty columns indicated by dots):
  .  (6)       (15)    (114)  (1113)  (11112)  .
     (33)      (24)    (123)  (1122)  (11121)
     (222)     (42)    (132)  (1131)  (11211)
     (111111)  (51)    (141)  (1221)  (12111)
               (1212)  (213)  (1311)  (21111)
               (2121)  (231)  (2112)
                       (312)  (2211)
                       (321)  (3111)
                       (411)
		

Crossrefs

Column k = 1 is A000005.
Row sums are A011782.
Diagonal T(2n,n) is A045630(n).
The strict version is A072574.
A version counting runs is A238279.
Column k = n - 1 is A254667.
Aperiodic compositions are counted by A000740.
Aperiodic binary words are counted by A027375.
The orderless period of prime indices is A052409.
Numbers whose binary expansion is periodic are A121016.
Periodic compositions are counted by A178472.
Period of binary expansion is A302291.
Numbers whose prime signature is aperiodic are A329139.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Necklaces are A065609.
- Sum is A070939.
- Rotational symmetries are counted by A138904.
- Constant compositions are A272919.
- Lyndon compositions are A275692.
- Co-Lyndon compositions are A326774.
- Aperiodic compositions are A328594.
- Rotational period is A333632.
- Co-necklaces are A333764.
- Reversed necklaces are A333943.

Programs

  • Mathematica
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],Function[c,Length[Union[Array[RotateRight[c,#]&,Length[c]]]]==k]]],{n,0,10},{k,0,n}]
  • PARI
    T(n,k)=if(n==0, k==0, sumdiv(n, m, sumdiv(gcd(k,m), d, moebius(d)*binomial(m/d-1, k/d-1)))) \\ Andrew Howroyd, Jan 19 2023

Formula

T(n,k) = Sum_{m|n} Sum_{d|gcd(k,m)} mu(d)*binomial(m/d-1, k/d-1) for n > 0. - Andrew Howroyd, Jan 19 2023
Previous Showing 11-20 of 29 results. Next