A052551
Expansion of 1/((1 - x)*(1 - 2*x^2)).
Original entry on oeis.org
1, 1, 3, 3, 7, 7, 15, 15, 31, 31, 63, 63, 127, 127, 255, 255, 511, 511, 1023, 1023, 2047, 2047, 4095, 4095, 8191, 8191, 16383, 16383, 32767, 32767, 65535, 65535, 131071, 131071, 262143, 262143, 524287, 524287, 1048575, 1048575, 2097151, 2097151
Offset: 0
encyclopedia(AT)pommard.inria.fr, Jan 25 2000
- S. Wolfram, A New Kind of Science, Wolfram Media, 2002; p. 170.
- Vincenzo Librandi, Table of n, a(n) for n = 0..2000
- INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 488
- Donatella Merlini and Massimo Nocentini, Algebraic Generating Functions for Languages Avoiding Riordan Patterns, Journal of Integer Sequences, Vol. 21 (2018), Article 18.1.3.
- N. J. A. Sloane, On the Number of ON Cells in Cellular Automata, arXiv:1503.01168 [math.CO], 2015.
- Eric Weisstein's World of Mathematics, Elementary Cellular Automaton
- S. Wolfram, A New Kind of Science
- Wolfram Research, Wolfram Atlas of Simple Programs
- Index entries for sequences related to cellular automata
- Index to 2D 5-Neighbor Cellular Automata
- Index to Elementary Cellular Automata
- Index entries for linear recurrences with constant coefficients, signature (1,2,-2).
Column 2 (offset by two) of
A304972.
-
Flat(List([1..21],n->[2^n-1,2^n-1])); # Muniru A Asiru, Oct 16 2018
-
[2^Floor(n/2)-1: n in [2..50]]; // Vincenzo Librandi, Aug 16 2011
-
spec := [S,{S=Prod(Sequence(Prod(Z,Union(Z,Z))),Sequence(Z))},unlabeled]: seq(combstruct[count](spec,size=n), n=0..20);
-
Table[StirlingS2[Floor[n/2] + 2, 2], {n, 0, 50}] (* Robert A. Russell, Dec 20 2017 *)
Drop[LinearRecurrence[{1, 2, -2}, {0, 1, 1}, 50], 1] (* Robert A. Russell, Oct 14 2018 *)
CoefficientList[Series[1/((1-x)*(1-2*x^2)), {x, 0, 50}], x] (* Stefano Spezia, Oct 16 2018 *)
2^(1+Floor[(Range[0,50])/2])-1 (* Federico Provvedi, Nov 22 2018 *)
((-1)^#(Sqrt[2]-1)+Sqrt[2]+1)2^((#-1)/2)-1&@Range[0, 50] (* Federico Provvedi, Nov 23 2018 *)
-
x='x+O('x^50); Vec(1/((1-x)*(1-2*x^2))) \\ Altug Alkan, Mar 19 2018
-
[2^(floor(n/2)) -1 for n in (2..50)] # G. C. Greubel, Mar 04 2019
A152175
Triangle read by rows: T(n,k) is the number of k-block partitions of an n-set up to rotations.
Original entry on oeis.org
1, 1, 1, 1, 1, 1, 1, 3, 2, 1, 1, 3, 5, 2, 1, 1, 7, 18, 13, 3, 1, 1, 9, 43, 50, 20, 3, 1, 1, 19, 126, 221, 136, 36, 4, 1, 1, 29, 339, 866, 773, 296, 52, 4, 1, 1, 55, 946, 3437, 4281, 2303, 596, 78, 5, 1, 1, 93, 2591, 13250, 22430, 16317, 5817, 1080, 105, 5, 1, 1, 179, 7254, 51075, 115100, 110462, 52376, 13299, 1873, 147, 6, 1
Offset: 1
Triangle begins with T(1,1):
1;
1, 1;
1, 1, 1;
1, 3, 2, 1;
1, 3, 5, 2, 1;
1, 7, 18, 13, 3, 1;
1, 9, 43, 50, 20, 3, 1;
1, 19, 126, 221, 136, 36, 4, 1;
1, 29, 339, 866, 773, 296, 52, 4, 1;
1, 55, 946, 3437, 4281, 2303, 596, 78, 5, 1;
1, 93, 2591, 13250, 22430, 16317, 5817, 1080, 105 , 5, 1;
1, 179, 7254, 51075, 115100, 110462, 52376, 13299, 1873, 147, 6, 1;
1, 315, 20125, 194810, 577577, 717024, 439648, 146124, 27654, 3025, 187, 6, 1;
...
For T(4,2)=3, the set partitions are AAAB, AABB, and ABAB.
For T(4,3)=2, the set partitions are AABC and ABAC.
- M. R. Nester (1999). Mathematical investigations of some plant interaction designs. PhD Thesis. University of Queensland, Brisbane, Australia. [See A056391 for pdf file of Chap. 2]
A(1,n,k) in formula is the Stirling subset number
A008277.
-
(* Using recursion formula from Gilbert and Riordan:*)
Adn[d_, n_] := Adn[d, n] = Which[0==n, 1, 1==n, DivisorSum[d, x^# &],
1==d, Sum[StirlingS2[n, k] x^k, {k, 0, n}],
True, Expand[Adn[d, 1] Adn[d, n-1] + D[Adn[d, n - 1], x] x]];
Table[CoefficientList[DivisorSum[n, EulerPhi[#] Adn[#, n/#] &]/(x n), x],
{n, 1, 10}] // Flatten (* Robert A. Russell, Feb 23 2018 *)
Adnk[d_,n_,k_] := Adnk[d,n,k] = If[n>0 && k>0, Adnk[d,n-1,k]k + DivisorSum[d,Adnk[d,n-1,k-#] &], Boole[n==0 && k==0]]
Table[DivisorSum[n,EulerPhi[#]Adnk[#,n/#,k]&]/n,{n,1,12},{k,1,n}] // Flatten (* Robert A. Russell, Oct 16 2018 *)
-
\\ see A056391 for Polya enumeration functions
T(n,k) = NonequivalentStructsExactly(CyclicPerms(n), k); \\ Andrew Howroyd, Oct 14 2017
-
R(n) = {Mat(Col([Vecrev(p/y, n) | p<-Vec(intformal(sum(m=1, n, eulerphi(m) * subst(serlaplace(-1 + exp(sumdiv(m, d, y^d*(exp(d*x + O(x*x^(n\m)))-1)/d))), x, x^m))/x))]))}
{ my(A=R(12)); for(n=1, #A, print(A[n, 1..n])) } \\ Andrew Howroyd, Sep 20 2019
A059053
Number of chiral pairs of necklaces with n beads and two colors (color complements being equivalent); i.e., turning the necklace over neither leaves it unchanged nor simply swaps the colors.
Original entry on oeis.org
0, 0, 0, 0, 0, 0, 0, 1, 2, 7, 12, 31, 58, 126, 234, 484, 906, 1800, 3402, 6643, 12624, 24458, 46686, 90157, 172810, 333498, 641340, 1238671, 2388852, 4620006, 8932032, 17302033, 33522698, 65042526, 126258960, 245361172, 477091232
Offset: 0
For a(7) = 1, the chiral pair is AAABABB-AAABBAB.
For a(8) = 2, the chiral pairs are AAAABABB-AAAABBAB and AAABAABB-AAABBAAB.
-
Prepend[Table[DivisorSum[n, EulerPhi[#] StirlingS2[n/# + If[Divisible[#,2],1,0], 2] &] / (2n) - StirlingS2[1+Floor[n/2],2] / 2, {n, 1, 40}],0] (* Robert A. Russell, Oct 02 2018 *)
-
a(n) = {if(n<1, 0, (sumdiv(n, k, eulerphi(2*k) * 2^(n/k)) / (2*n) - 2^(n\2))/2)}; \\ Andrew Howroyd, Nov 03 2019
A294791
Triangle read by rows, 1 <= k <= n: T(n,k) = non-isomorphic colorings of a toroidal n X k grid using exactly two colors under translational symmetry and swappable colors.
Original entry on oeis.org
0, 1, 4, 1, 7, 31, 3, 23, 179, 2107, 3, 55, 1095, 26271, 671103, 7, 189, 7327, 350063, 17896831, 954459519, 9, 595, 49939, 4794087, 490853415, 52357746895, 5744387279871, 19, 2101, 349715, 67115111, 13743921631, 2932032057731, 643371380132743, 144115188277194943, 29, 7315, 2485591, 954444607, 390937468407, 166799988703927, 73201365371896619
Offset: 1
For the 2 X 2 grid and two colors we find T(2,2) = 4:
+---+ +---+ +---+ +---+
|X| | |X| | |X|X| |X| |
+-+-+ +-+-+ +-+-+ +-+-+
| | | | |X| | | | |X| |
+-+-+ +-+-+ +-+-+ +-+-+
- F. Harary and E. Palmer, Graphical Enumeration, Academic Press, 1973.
A056303
Number of primitive (period n) n-bead necklace structures using exactly two different colored beads.
Original entry on oeis.org
0, 1, 1, 2, 3, 5, 9, 16, 28, 51, 93, 170, 315, 585, 1091, 2048, 3855, 7280, 13797, 26214, 49929, 95325, 182361, 349520, 671088, 1290555, 2485504, 4793490, 9256395, 17895679, 34636833, 67108864, 130150493, 252645135, 490853403, 954437120, 1857283155
Offset: 1
- M. R. Nester (1999). Mathematical investigations of some plant interaction designs. PhD Thesis. University of Queensland, Brisbane, Australia. [See A056391 for pdf file of Chap. 2]
-
vector(100, n, sumdiv(n, d, (d%2)*(moebius(d)*2^(n/d)))/(2*n)-!(n-1)) \\ Altug Alkan, Oct 26 2015
-
from sympy import divisors, mobius
def a000048(n): return 1 if n<1 else sum([mobius(d)*2**(n/d) for d in divisors(n) if d%2 == 1])/(2*n)
def a(n): return a000048(n) - 0**(n - 1) # Indranil Ghosh, Apr 28 2017
A056357
Number of bracelet structures using exactly two different colored beads.
Original entry on oeis.org
0, 1, 1, 3, 3, 7, 8, 17, 22, 43, 62, 121, 189, 361, 611, 1161, 2055, 3913, 7154, 13647, 25481, 48733, 92204, 176905, 337593, 649531, 1246862, 2405235, 4636389, 8964799, 17334800, 33588233, 65108061, 126390031, 245492243, 477353375, 928772649, 1808676325
Offset: 1
- M. R. Nester (1999). Mathematical investigations of some plant interaction designs. PhD Thesis. University of Queensland, Brisbane, Australia. [See A056391 for pdf file of Chap. 2]
-
with(numtheory);
rho:=n->(3+(-1)^n)/2;
f:=n->2^((n+rho(n))/2-2) + (1/(4*n))*(add(phi(d)*rho(d)*2^(n/d), d in divisors(n))) - 1;
# N. J. A. Sloane, Jul 13 2011
-
a(n) = {if(n<1, 0, 2^(n\2-1) - 1 + sumdiv(n, k, eulerphi(2*k) * 2^(n/k)) / (4*n))}; \\ Andrew Howroyd, Oct 24 2019
A187767
Number of bicolored cyclic patterns n X n.
Original entry on oeis.org
0, 2, 3, 10, 15, 35, 63, 138, 255, 527, 1023, 2083, 4095, 8255, 16383, 32906, 65535, 131327, 262143, 524815, 1048575, 2098175, 4194303, 8390691, 16777215, 33558527, 67108863, 134225983, 268435455, 536887295, 1073741823, 2147516554, 4294967295, 8590000127, 17179869183
Offset: 1
a(4)=10 is represented below. See Links for more examples.
. 1000 0100 0010 0001 0101 1010 1001 0110 1100 0011
. 0100 0001 0100 0001 0101 0101 1100 1100 0011 0011
. 0010 0100 1000 0001 0101 1010 0110 1001 1100 0011
. 0001 0001 0001 0001 0101 0101 0011 0011 0011 0011
The number of patterns made of vertical stripes only is
A056295(n).
-
cyPatt[n_]:=Block[{b,c},c[v_,q_:1]:=Table[RotateLeft[v,i q],{i,n}]; b=Union[(First@Union[c@#,c[1-#]])& /@ IntegerDigits[Range[2^n/2-1], 2,n]]; Union@Flatten[Table[c[e,j],{j,n},{e,b}],1]];
(*count*) a[n_] := Length@cyPatt@n; Print["Seq = ",a/@Range[12]];
(*show*) showP[p_] := GraphicsGrid@Partition[ArrayPlot/@p,8,8,1,Null];
showP[cyPatt[6]]
-
b(n)=sumdiv(n,d,(d%2)*(moebius(d)*2^(n/d)))/(2*n);
a(n)=sumdiv(n,d,d*b(d)) - 1; \\ Andrew Howroyd, Jun 02 2017
A327734
Number of n-bead bracelet structures using exactly two different colored beads that are not self-equivalent under either a nonzero rotation or reversal (turning over bracelet).
Original entry on oeis.org
0, 0, 0, 0, 0, 0, 1, 2, 7, 12, 31, 57, 126, 232, 484, 900, 1800, 3388, 6643, 12594, 24457, 46624, 90157, 172680, 333498, 641088, 1238664, 2388357, 4620006, 8931052
Offset: 1
Case n=7: There are 7 necklace structures using exactly two colors. See example in A056295. Of these, 5 are achiral and the other two a chiral pair: AAABABB/AAABBAB. a(7) = 1, since this sequence considers chiral pairs as the same bracelet structure.
Showing 1-8 of 8 results.
Comments