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-10 of 21 results. Next

A001037 Number of degree-n irreducible polynomials over GF(2); number of n-bead necklaces with beads of 2 colors when turning over is not allowed and with primitive period n; number of binary Lyndon words of length n.

Original entry on oeis.org

1, 2, 1, 2, 3, 6, 9, 18, 30, 56, 99, 186, 335, 630, 1161, 2182, 4080, 7710, 14532, 27594, 52377, 99858, 190557, 364722, 698870, 1342176, 2580795, 4971008, 9586395, 18512790, 35790267, 69273666, 134215680, 260300986, 505286415, 981706806, 1908866960, 3714566310, 7233615333, 14096302710, 27487764474
Offset: 0

Views

Author

Keywords

Comments

Also dimensions of free Lie algebras - see A059966, which is essentially the same sequence.
This sequence also represents the number N of cycles of length L in a digraph under x^2 seen modulo a Mersenne prime M_q=2^q-1. This number does not depend on q and L is any divisor of q-1. See Theorem 5 and Corollary 3 of the Shallit and Vasiga paper: N=sum(eulerphi(d)/order(d,2)) where d is a divisor of 2^(q-1)-1 such that order(d,2)=L. - Tony Reix, Nov 17 2005
Except for a(0) = 1, Bau-Sen Du's [1985/2007] Table 1, p. 6, has this sequence as the 7th (rightmost) column. Other columns of the table include (but are not identified as) A006206-A006208. - Jonathan Vos Post, Jun 18 2007
"Number of binary Lyndon words" means: number of binary strings inequivalent modulo rotation (cyclic permutation) of the digits and not having a period smaller than n. This provides a link to A103314, since these strings correspond to the inequivalent zero-sum subsets of U_m (m-th roots of unity) obtained by taking the union of U_n (n|m) with 0 or more U_d (n | d, d | m) multiplied by some power of exp(i 2Pi/n) to make them mutually disjoint. (But not all zero-sum subsets of U_m are of that form.) - M. F. Hasler, Jan 14 2007
Also the number of dynamical cycles of period n of a threshold Boolean automata network which is a quasi-minimal positive circuit of size a multiple of n and which is updated in parallel. - Mathilde Noual (mathilde.noual(AT)ens-lyon.fr), Feb 25 2009
Also, the number of periodic points with (minimal) period n in the iteration of the tent map f(x):=2min{x,1-x} on the unit interval. - Pietro Majer, Sep 22 2009
Number of distinct cycles of minimal period n in a shift dynamical system associated with a totally disconnected hyperbolic iterated function system (see Barnsley link). - Michel Marcus, Oct 06 2013
From Jean-Christophe Hervé, Oct 26 2014: (Start)
For n > 0, a(n) is also the number of orbits of size n of the transform associated to the Kolakoski sequence A000002 (and this is true for any map with 2^n periodic points of period n). The Kolakoski transform changes a sequence of 1's and 2's by the sequence of the lengths of its runs. The Kolakoski sequence is one of the two fixed points of this transform, the other being the same sequence without the initial term. A025142 and A025143 are the periodic points of the orbit of size 2. A027375(n) = n*a(n) gives the number of periodic points of minimal period n.
For n > 1, this sequence is equal to A059966 and to A060477, and for n = 1, a(1) = A059966(1)+1 = A060477(1)-1; this because the n-th term of all 3 sequences is equal to (1/n)*sum_{d|n} mu(n/d)*(2^d+e), with e = -1/0/1 for resp. A059966/this sequence/A060477, and sum_{d|n} mu(n/d) equals 1 for n = 1 and 0 for all n > 1. (End)
Warning: A000031 and A001037 are easily confused, since they have similar formulas.
From Petros Hadjicostas, Jul 14 2020: (Start)
Following Kam Cheong Au (2020), let d(w,N) be the dimension of the Q-span of weight w and level N of colored multiple zeta values (CMZV). Here Q are the rational numbers.
Deligne's bound says that d(w,N) <= D(w,N), where 1 + Sum_{w >= 1} D(w,N)*t^w = (1 - a*t + b*t^2)^(-1) when N >= 3, where a = phi(N)/2 + omega(N) and b = omega(N) - 1 (with omega(N) = A001221(N) being the number of distinct primes of N).
For N = 3, a = phi(3)/2 + omega(3) = 2/2 + 1 = 2 and b = omega(3) - 1 = 0. It follows that D(w, N=3) = A000079(w) = 2^w.
For some reason, Kam Cheong Au (2020) assumes Deligne's bound is tight, i.e., d(w,N) = D(w,N). He sets Sum_{w >= 1} c(w,N)*t^w = log(1 + Sum_{w >= 1} d(w,N)*t^w) = log(1 + Sum_{w >= 1} D(w,N)*t^w) = -log(1 - a*t + b*t^2) for N >= 3.
For N = 3, we get that c(w, N=3) = A000079(w)/w = 2^w/w.
He defines d*(w,N) = Sum_{k | w} (mu(k)/k)*c(w/k,N) to be the "number of primitive constants of weight w and level N". (Using the terminology of A113788, we may perhaps call d*(w,N) the number of irreducible colored multiple zeta values at weight w and level N.)
Using standard techniques of the theory of g.f.'s, we can prove that Sum_{w >= 1} d*(w,N)*t^w = Sum_{s >= 1} (mu(s)/s) Sum_{k >= 1} c(k,N)*(t^s)^k = -Sum_{s >= 1} (mu(s)/s)*log(1 - a*t^s + b*t^(2*s)).
For N = 3, we saw that a = 2 and b = 0, and hence d*(w, N=3) = a(w) = Sum_{k | w} (mu(k)/k) * 2^(w/k) / (w/k) = (1/w) * Sum_{k | w} mu(k) * 2^(w/k) for w >= 1. See Table 1 on p. 6 in Kam Cheong Au (2020). (End)

Examples

			Binary strings (Lyndon words, cf. A102659):
a(0) = 1 = #{ "" },
a(1) = 2 = #{ "0", "1" },
a(2) = 1 = #{ "01" },
a(3) = 2 = #{ "001", "011" },
a(4) = 3 = #{ "0001", "0011", "0111" },
a(5) = 6 = #{ "00001", "00011", "00101", "00111", "01011", "01111" }.
		

References

  • Michael F. Barnsley, Fractals Everywhere, Academic Press, San Diego, 1988, page 171, Lemma 3.
  • E. R. Berlekamp, Algebraic Coding Theory, McGraw-Hill, NY, 1968, p. 84.
  • E. L. Blanton, Jr., S. P. Hurd and J. S. McCranie. On the digraph defined by squaring mod m, when m has primitive roots. Congr. Numer. 82 (1991), 167-177.
  • P. J. Freyd and A. Scedrov, Categories, Allegories, North-Holland, Amsterdam, 1990. See 1.925.
  • M. Lothaire, Combinatorics on Words, Addison-Wesley, Reading, MA, 1983, pp. 65, 79.
  • Robert M. May, "Simple mathematical models with very complicated dynamics." Nature, Vol. 261, June 10, 1976, pp. 459-467; reprinted in The Theory of Chaotic Attractors, pp. 85-93. Springer, New York, NY, 2004. The sequences listed in Table 2 are A000079, A027375, A000031, A001037, A000048, A051841. - N. J. A. Sloane, Mar 17 2019
  • Guy Melançon, Factorizing infinite words using Maple, MapleTech Journal, vol. 4, no. 1, 1997, pp. 34-42, esp. p. 36.
  • 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]
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence in entries N0046 and N0287).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Column 2 of A074650.
Row sums of A051168, which gives the number of Lyndon words with fixed number of zeros and ones.
Euler transform is A000079.
See A058943 and A102569 for initial terms. See also A058947, A011260, A059966.
Irreducible over GF(2), GF(3), GF(4), GF(5), GF(7): A058943, A058944, A058948, A058945, A058946. Primitive irreducible over GF(2), GF(3), GF(4), GF(5), GF(7): A058947, A058949, A058952, A058950, A058951.
Cf. A000031 (n-bead necklaces but may have period dividing n), A014580, A046211, A046209, A006206-A006208, A038063, A060477, A103314.
See also A102659 for the list of binary Lyndon words themselves.

Programs

  • Haskell
    a001037 0 = 1
    a001037 n = (sum $ map (\d -> (a000079 d) * a008683 (n `div` d)) $
                           a027750_row n) `div` n
    -- Reinhard Zumkeller, Feb 01 2013
    
  • Maple
    with(numtheory): A001037 := proc(n) local a,d; if n = 0 then RETURN(1); else a := 0: for d in divisors(n) do a := a+mobius(n/d)*2^d; od: RETURN(a/n); fi; end;
  • Mathematica
    f[n_] := Block[{d = Divisors@ n}, Plus @@ (MoebiusMu[n/d]*2^d/n)]; Array[f, 32]
  • PARI
    A001037(n)=if(n>1,sumdiv(n,d,moebius(d)*2^(n/d))/n,n+1) \\ Edited by M. F. Hasler, Jan 11 2016
    
  • PARI
    {a(n)=polcoeff(1-sum(k=1,n,moebius(k)/k*log(1-2*x^k+x*O(x^n))),n)} \\ Paul D. Hanna, Oct 13 2010
    
  • PARI
    a(n)=if(n>1,my(s);forstep(i=2^n+1,2^(n+1),2,s+=polisirreducible(Mod(1,2) * Pol(binary(i))));s,n+1) \\ Charles R Greathouse IV, Jan 26 2012
    
  • Python
    from sympy import divisors, mobius
    def a(n): return sum(mobius(d) * 2**(n//d) for d in divisors(n))/n if n>1 else n + 1 # Indranil Ghosh, Apr 26 2017

Formula

For n >= 1:
a(n) = (1/n)*Sum_{d | n} mu(n/d)*2^d.
A000031(n) = Sum_{d | n} a(d).
2^n = Sum_{d | n} d*a(d).
a(n) = A027375(n)/n.
a(n) = A000048(n) + A051841(n).
For n > 1, a(n) = A059966(n) = A060477(n).
G.f.: 1 - Sum_{n >= 1} moebius(n)*log(1 - 2*x^n)/n, where moebius(n) = A008683(n). - Paul D. Hanna, Oct 13 2010
From Richard L. Ollerton, May 10 2021: (Start)
For n >= 1:
a(n) = (1/n)*Sum_{k=1..n} mu(gcd(n,k))*2^(n/gcd(n,k))/phi(n/gcd(n,k)).
a(n) = (1/n)*Sum_{k=1..n} mu(n/gcd(n,k))*2^gcd(n,k)/phi(n/gcd(n,k)). (End)
a(n) ~ 2^n / n. - Vaclav Kotesovec, Aug 11 2021

Extensions

Revised by N. J. A. Sloane, Jun 10 2012

A107754 Number of subsets of the n-th roots of unity that sum to 1.

Original entry on oeis.org

1, 1, 1, 2, 1, 6, 1, 8, 4, 18, 1, 60, 1, 66, 20, 128, 1, 600, 1, 612, 68, 1026, 1, 6000, 16, 4098, 256, 8580, 1, 95226, 1, 32768
Offset: 1

Views

Author

T. D. Noe, May 23 2005

Keywords

Crossrefs

Cf. A103314 (number of subsets of the n-th roots of unity summing to zero) and A108417 (number of subsets of the n-th roots of unity summing to the absolute value of 1).

Programs

  • Mathematica
    << DiscreteMath`Combinatorica`; f[n_] := Plus @@ Table[ Count[ KSubsets[ Range[n], k], q_List /; Chop[ -1 + Plus @@ (E^((2.*Pi*I*q)/n))] === 0], {k, 0, n}]; Table[ f[n], {n, 24}] (* Robert G. Wilson v, Jun 03 2005 *)

Formula

For prime p, a(p^i) = 2^(p^(i-1)-1).

A262181 a(n) = total number of convex equilateral n-gons with corner angles of m*Pi/n (0 < m <= n).

Original entry on oeis.org

1, 2, 1, 11, 1, 42, 64, 202, 1, 1557, 1, 5539, 32298, 30666, 1, 405200, 1, 1035642
Offset: 3

Views

Author

Stuart E Anderson, Sep 14 2015

Keywords

Comments

An n-gon is a polygon with n corners and n sides, each of which is a straight line segment joining two corners. An n-gon or polygon P is said to be a simple polygon (or a Jordan polygon) if the only points of the plane belonging to two polygon edges of P are the polygon vertices of P. Such a polygon has a well-defined interior and exterior. Simple polygons are topologically equivalent to a disk, hence zero angles are not allowed; allowable angles are m*Pi/n (where m and n are integers and 0 < m <= n). An n-gon is convex if it contains all the diagonal segments connecting any pair of its points. A convex polygon is sometimes strictly defined as a polygon with all its interior angles less than Pi. We use the less strict definition where every internal or interior angle is less than or equal to Pi, that is, straight angles are permitted.
Conjecture: There is only one convex equilateral n-gon for prime n.

Examples

			For n = 3 there is one convex n-gon, the equilateral triangle, with m angle factors (3 3 3); so a(3) = 1.
For n = 4 there are two convex n-gons, the square and a rhombus, with respective m angle factors (2 2 2 2) and (1 3 1 3); so a(4) = 2.
For n = 5, there is the regular pentagon, m factors (3 3 3 3 3); so a(5) = 1.
For n = 6 there are 11 convex n-gons; here are the m factors:(1 5 6 1 5 6), (1 6 5 1 6 5), (2 4 6 2 4 6), (2 5 5 2 5 5), (2 6 2 6 2 6), (2 6 4 2 6 4), (3 3 6 3 3 6), (3 4 5 3 4 5), (3 5 3 5 3 5), (3 5 4 3 5 4), (4 4 4 4 4 4); so a(6) = 11.
		

Crossrefs

A262244 for concave polygons with corner angles of m*Pi/n (0 < m < 2n), where m and n are integers.

Formula

a(n) = A292355(n) for n prime or twice prime. - Andrew Howroyd, Sep 14 2017
a(n) = -(1+(-1)^n)/2 + (1/(2*n))*(A321415(n) - binomial(3*n-1, n) + Sum_{d|n} phi(n/d) * binomial(3*d-1, d)). - Andrew Howroyd, Nov 09 2018

Extensions

a(10) corrected and a(12)-a(17) from Andrew Howroyd, Sep 14 2017
a(18)-a(20) from Andrew Howroyd, Nov 09 2018

A103306 Triangle read by rows: T(n,k) = number of k-subsets of the n-th roots of 1 that add to zero (0 <= k <= n).

Original entry on oeis.org

1, 1, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 2, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 3, 2, 3, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 4, 0, 6, 0, 4, 0, 1, 1, 0, 0, 3, 0, 0, 3, 0, 0, 1, 1, 0, 5, 0, 10, 2, 10, 0, 5, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 6, 4, 15, 12, 24, 12, 15, 4, 6, 0, 1, 1, 0, 0, 0, 0
Offset: 0

Views

Author

Wouter Meeussen, Mar 11 2005

Keywords

Comments

Observe that T(n,k) = binomial(n,k) (mod n). Because the sum of the n n-th roots of unity is 0 for n>1, each row is symmetric for n>1. Hence only k=0..floor(n/2) need to be computed. - T. D. Noe, Jan 16 2008

Examples

			Triangle begins:
{1},
{1, 0},
{1, 0, 1},
{1, 0, 0, 1},
{1, 0, 2, 0, 1},
{1, 0, 0, 0, 0, 1},
{1, 0, 3, 2, 3, 0, 1},
{1, 0, 0, 0, 0, 0, 0, 1},
{1, 0, 4, 0, 6, 0, 4, 0, 1},
{1, 0, 0, 3, 0, 0, 3, 0, 0, 1},
T(10,4)=10, counting {1,2,6,7}, {1,3,6,8}, {1,4,6,9}, {1,5,6,10}, {2,3,7,8}, {2,4,7,9}, {2,5,7,10}, {3,4,8,9}, {3,5,8,10}, {4,5,9,10}.
		

Crossrefs

Row sums give A103314.

Programs

  • Mathematica
    < n/2 >= 1}}, Count[Subsets[Range[n], {k}], subset_/;PossibleZeroQ[ExpToTrig[Sum[Exp[2*Pi*I*m/n], {m, subset}]]]]]; Table[T[n, k], {n, 0, 12}, {k, 0, n}] // TableForm (* David M. Zimmerman, Sep 23 2020 *)

A110981 a(n) = the number of aperiodic subsets S of the n-th roots of 1 with zero sum (i.e., there is no r different from 1 such that r*S=S).

Original entry on oeis.org

1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 24, 0, 6, 0, 0, 0, 236, 0, 0, 0, 18, 0, 3768, 0, 0, 0, 0, 0, 20384, 0, 0, 0, 7188, 0, 227784, 0, 186, 480, 0, 0, 1732448, 0, 237600, 0, 630, 0, 16028160, 0, 306684, 0, 0, 0, 341521732, 0, 0, 4896, 0, 0, 1417919208
Offset: 1

Views

Author

Max Alekseyev, Jan 20 2008

Keywords

Comments

We count these subsets only modulo rotations (multiplication by a nontrivial root of unity).
A103314(n) = a(n)*n + 2^n - A001037(n)*n. Note that as soon as a(n)=0, we have simply A103314(n) = 2^n - A001037(n)*n. This makes it especially interesting to study those n for which a(n)=0. Quite surprisingly, it appears that the sequence of such n coincides with A102466.
From Max Alekseyev, Jan 31 2008: (Start)
Every subset of the set U(n) = { 1=r^0, r^1, ..., r^(n-1) } of the n-th power roots of 1 (where r is a fixed primitive root) defines a binary word w of the length n where the j-th bit is 1 iff the root r^j is included in the subset.
If d is the period of w with respect to cyclic rotations (thus d|n) then the periodic part of w uniquely defines some binary Lyndon word of the length d (see A001037). In turn, each binary Lyndon word of the length d, where d
The binary Lyndon words of the length n are different in this respect: only some of them correspond to n distinct zero-sum subsets of U(n) while the others do not correspond to such subsets at all. A110981(n) gives the number of binary Lyndon words of the length n that correspond to zero-sum subsets of U(n). (End)

Crossrefs

Formula

a(n) = A001037(n) - A107847(n) ( = A001037(n) - (2^n - A103314(n))/n ). - M. F. Hasler, Jan 31 2008

Extensions

Additional comments from M. F. Hasler, Jan 31 2008

A164896 Number of subsets (up to cyclic shifts) of the n-th roots of 1 with zero sum.

Original entry on oeis.org

1, 2, 2, 3, 2, 5, 2, 6, 4, 9, 2, 19, 2, 21, 10, 36, 2, 94, 2, 117, 22, 189, 2, 618, 8, 633, 60, 1203, 2, 6069, 2, 4116, 190, 7713, 26, 35324, 2, 27597, 634, 59706, 2, 328835, 2, 190935, 2728, 364725, 2, 2435780, 20, 1579884, 7714, 2582061, 2, 21013770, 194, 9894294, 27598, 18512793, 2, 377367015, 2, 69273669, 104832, 134219796, 638, 1678410951
Offset: 1

Author

Joerg Arndt, Aug 30 2009

Keywords

Comments

Cyclic shifts correspond to multiplication by a root of unity.
a(n)=2 for n prime, corresponding to the empty and the full subset. - Joerg Arndt, Jun 10 2011

Examples

			a(6) = 5 because these subsets add to zero: (left: as bitstring, right: subset)
  ......  (empty sum)
  ..1..1  0 3
  .1.1.1  0 2 4
  .11.11  0 1 3 4
  111111  0 1 2 3 4 5 (all roots of unity)
		

Crossrefs

Cf. A066656, A103314, A110981 (counts subsets with bitstrings being Lyndon words).

Formula

a(n) = A110981(n) + Sum_{d|n,dA001037(d) = A110981(n) + A000031(n) - A001037(n). - Max Alekseyev, Apr 08 2013
a(n) = A110981(n) + A066656(n). - Andrew Howroyd, Mar 22 2023

Extensions

a(32)-a(39) from Joerg Arndt, Jun 10 2011
Terms a(40) onward from Max Alekseyev, Apr 08 2013

A321414 Array read by antidiagonals: T(n,k) is the number of n element multisets of the 2k-th roots of unity with zero sum.

Original entry on oeis.org

0, 0, 1, 0, 2, 0, 0, 3, 0, 1, 0, 4, 2, 3, 0, 0, 5, 0, 6, 0, 1, 0, 6, 0, 10, 6, 4, 0, 0, 7, 4, 15, 0, 12, 0, 1, 0, 8, 0, 21, 2, 20, 12, 5, 0, 0, 9, 0, 28, 24, 35, 0, 21, 0, 1, 0, 10, 6, 36, 0, 64, 10, 35, 22, 6, 0, 0, 11, 0, 45, 0, 84, 84, 70, 0, 33, 0, 1
Offset: 1

Author

Andrew Howroyd, Nov 08 2018

Keywords

Comments

Equivalently, the number of closed convex paths of length n whose steps are the 2k-th roots of unity up to translation. For even n, there will be k paths of zero area consisting of n/2 steps in one direction followed by n/2 steps in the opposite direction.

Examples

			Array begins:
  =========================================================
  n\k| 1  2  3  4   5   6   7    8    9   10   11    12
  ---|-----------------------------------------------------
   1 | 0  0  0  0   0   0   0    0    0    0    0     0 ...
   2 | 1  2  3  4   5   6   7    8    9   10   11    12 ...
   3 | 0  0  2  0   0   4   0    0    6    0    0     8 ...
   4 | 1  3  6 10  15  21  28   36   45   55   66    78 ...
   5 | 0  0  6  0   2  24   0    0   54    4    0    96 ...
   6 | 1  4 12 20  35  64  84  120  183  220  286   396 ...
   7 | 0  0 12  0  10  84   2    0  270   40    0   624 ...
   8 | 1  5 21 35  70 174 210  330  657  715 1001  1749 ...
   9 | 0  0 22  0  30 236  14    0 1028  220    0  3000 ...
  10 | 1  6 33 56 128 420 462  792 2097 2010 3003  6864 ...
  11 | 0  0 36  0  70 576  56    0 3312  880    2 11976 ...
  12 | 1  7 50 84 220 926 924 1716 6039 5085 8008 24216 ...
  ...
T(5, 3) = 6 because there are 6 rotations of the following figure:
       o---o
      /     \
     o---o---o
.
T(6, 3) = 12 because there are 4 basic shapes illustrated below which with rotations and reflections give 3 + 2 + 1 + 6 = 12 convex paths.
                        o        o---o     o---o
                       / \      /     \     \   \
    o===o===o===o     o   o    o       o     o   o
                     /     \    \     /       \   \
                    o---o---o    o---o         o---o
		

Crossrefs

Main diagonal is A321415.
Columns include A053090(n+3), A321416, A321417, A321419.

Programs

  • PARI
    \\ only supports k with at most one odd prime factor.
    T(n, k)={my(r=valuation(k, 2), p); polcoef(if(k>>r == 1, 1/(1-x^2)^k + O(x*x^n), if(isprimepower(k>>r, &p), ((2/(1 - x^p) - 1)/(1 - x^2 + O(x*x^n))^p)^(k/p), error("Cannot handle k=", k) )), n)}

Formula

G.f. of column k = 2^r: 1/(1 - x^2)^k - 1.
G.f. of column k = 2^r*p^e: ((2/(1 - x^p) - 1)/(1 - x^2)^p)^(k/p) - 1 for odd prime p.

A107753 Number of primitive subsets of the n-th roots of unity summing to zero.

Original entry on oeis.org

1, 2, 2, 3, 2, 6, 2, 5, 4, 8, 2, 11, 2, 10, 9, 9, 2, 16, 2, 15, 11, 14, 2, 21, 6, 16, 10, 19, 2, 212, 2, 17, 15, 20, 13, 31, 2, 22, 17, 29, 2
Offset: 1

Author

T. D. Noe, May 23 2005

Keywords

Comments

A primitive subset has no nonempty proper subset whose members sum to zero. Note that a(30) is the first term for which the formulas do not apply. For n=30, there are 1,0,15,10,0,5,30,60,60,30 primitive subsets of size 0,1,2,...,9.

Crossrefs

Cf. A103314 (number of subsets of the n-th roots of unity summing to zero), A107754 (number of subsets of the n-th roots of unity summing to one).

Formula

For primes p and q, if n = p^i, then a(n)=1+n/p; if n=p^i q^j, then a(n)=1+n/p+n/q.

A107847 Related to sums of the n-th roots of unity: sums in a circular wedge (excluding the origin).

Original entry on oeis.org

1, 1, 2, 3, 6, 9, 18, 30, 56, 99, 186, 333, 630, 1161, 2182, 4080, 7710, 14508, 27594, 52371, 99858, 190557, 364722, 698634, 1342176, 2580795, 4971008, 9586377, 18512790, 35786499, 69273666, 134215680, 260300986, 505286415, 981706806
Offset: 1

Author

T. D. Noe, May 25 2005

Keywords

Comments

Consider the 2^n sums formed from all the subsets of the n-th roots of unity. The number A103314(n) tells how many of these sums are zero. The remaining sums fall into n wedges centered at the origin. The number a(n) gives the number of sums that fall into each wedge. Interestingly, a(n) coincides with A059966(n) when n is either p^k or pq for primes p and q.

Crossrefs

Cf. A103314 (number of subsets of the n-th roots of unity summing to zero), A107848 (number of subsets of the n-th roots of unity summing to a real number).
Cf. also A110981.

Formula

a(n) = (2^n - A103314(n))/n.
a(n) = A001037(n) - A110981(n). - Max Alekseyev, Jan 14 2008

A107848 Number of subsets of the n-th roots of unity summing to a real number.

Original entry on oeis.org

2, 4, 4, 8, 8, 24, 16, 48, 40, 144, 64, 336, 128, 864, 416, 1728, 512, 8304, 1024, 10656, 4032, 31104, 4096, 116256, 9248, 186624, 40000, 374976, 32768, 3537024, 65536, 2239488
Offset: 1

Author

T. D. Noe, May 25 2005

Keywords

Crossrefs

Cf. A103314 (number of subsets of the n-th roots of unity summing to zero).

Programs

  • Mathematica
    Needs["DiscreteMath`Combinatorica`"]; Table[Plus@@Table[Count[(KSubsets[Range[n], k]), q_List/;Im[Chop[Plus@@(E^(2.*Pi*I*q/n))]]==0], {k, 0, n}], {n, 20}]

Formula

For prime n, a(n)=2^((n+1)/2).
Showing 1-10 of 21 results. Next