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

A034017 Numbers that are primitively represented by x^2 + xy + y^2.

Original entry on oeis.org

0, 1, 3, 7, 13, 19, 21, 31, 37, 39, 43, 49, 57, 61, 67, 73, 79, 91, 93, 97, 103, 109, 111, 127, 129, 133, 139, 147, 151, 157, 163, 169, 181, 183, 193, 199, 201, 211, 217, 219, 223, 229, 237, 241, 247, 259, 271, 273, 277, 283, 291, 301, 307, 309, 313, 327, 331
Offset: 1

Views

Author

Keywords

Comments

Gives the location of the nonzero terms of A000086.
Starting at a(3), a(n)^2 is the ordered semiperimeter of primitive integer Soddyian triangles (see A210484). - Frank M Jackson, Feb 04 2013
A000086(a(n)) > 0; a(n) = A004611(k) or a(n) = 3*A004611(k) for n > 3 and an appropriate k. - Reinhard Zumkeller, Jun 23 2013
The number of structure units in an icosahedral virus is 20*a(n), see Stannard link. - Charles R Greathouse IV, Nov 03 2015
From Wolfdieter Lang, Apr 09 2021: (Start)
The positive definite binary quadratic form F = [1, 1, 1], that is x^2 + x*y + y^2, has discriminant Disc = -3, and class number 1 (see Buell, Examples, p. 19, first line: Delta = -3, h = 1). This reduced form is equivalent to the form [1,-1, 1], but to no other reduced one (see Buell, Theorem 2.4, p. 15).
This form F represents a positive integer k (= a(n)) properly if and only if A002061(j+1) = 2*T(j) + 1 = j^2 + j + 1 == 0 (mod k), for j from {0, 1, ..., k-1}. This congruence determines the representative parallel primitive forms (rpapfs) of discriminant Disc = -3 and representation of a positive integer number k, given by [k, 2*j+1, c(j)], and c(j) is determined from Disc =-3 as c(j) = ((2*j+1)^2 + 3)/(4*k) = (j^2 + j + 1)/k. Each rpapf has a first reduced form, the so-called right neighbor form, namely [1, 1, 1] for k = 1 = a(1) (the already reduced parallel form from j = 0), and [1, -1, 1] for k = a(n), with n >= 2.
Only odd numbers k are eligible for representation, because 2*T(j) + 1, with the triangular numbers T = A000217, is odd. The odd k with at least one solution of the congruence are then the members of the present sequence.
The solutions of the reduced forms F = [1, 1, 1] and F' = [1, -1, 1] representing k are related by type I equivalence because of the first two entries ([a, a, c] == [a, -a, c]), and also by type II equivalence because [a, b, a] == [a, -b, a], for positive b. These transformation matrices are R_I = Matrix([1, -1],[0, 1]) and R_{II} = Matrix([0, -1], [1, 0]), respectively, to obtain the forms with negative second entry from the ones with positive second entry. The corresponding solutions (x, y)^t (t for transposed) are related by the inverse of these matrices.
The table with the A341422(n) solutions j of the congruence given above are given in A343232. (End)
Apparently, also the integers k that can be expressed as a quotient of two terms from A002061. - Martin Becker, Aug 14 2022
For some x, y let a(n) = r, x*(x+y) = s, y*(x+y) = t, x*y = u then (r,s,t,u) is a Pythagorean quadruple such that r^2 = s^2 + t^2 + u^2. - Frank M Jackson, Feb 26 2024

References

  • B. C. Berndt and R. A. Rankin, Ramanujan: Letters and Commentary, see p. 184, AMS, Providence, RI, 1995.
  • D. A. Buell, Binary Quadratic Forms, Springer, 1989, pp. 15, 19.

Crossrefs

Cf. A000217, A002061, A002476, A003136, A007645 (primes), A045611, A045897, A226946 (complement), A045897 (subsequence), A341422, A343232.

Programs

  • Haskell
    a034017 n = a034017_list !! (n-1)
    a034017_list = 0 : filter ((> 0) . a000086) [1..]
    -- Reinhard Zumkeller, Jun 23 2013
  • Maple
    N:= 1000: # to get all terms <= N
    P:= select(isprime, [seq(6*n+1, n=1..floor((N-1)/6))]):
    A:= {1,3}:
    for p in P do
      A:= {seq(seq(a*p^k, k=0..floor(log[p](N/a))),a=A)}:
    od:
    sort(convert(A,list)); # Robert Israel, Nov 04 2015
  • Mathematica
    lst = {0}; maxLen = 331; Do[If[Reduce[m^2 + m*n + n^2 == k && m >= n >= 0 && GCD[m, n] == 1, {m, n}, Integers] === False, , AppendTo[lst, k]], {k, maxLen}]; lst (* Frank M Jackson, Jan 10 2013 *) (* simplified by T. D. Noe, Feb 05 2013 *)
  • PARI
    is(n)=my(f=factor(n));for(i=1,#f[,1],if(f[i,1]%3!=1 && (f[i,1]!=3 || f[i,2]>1), return(n==0))); 1 \\ Charles R Greathouse IV, Jan 10 2013
    
  • PARI
    list(lim)=if(lim<7, return(select(n->n<=lim, [0,1,3]))); my(v=List([0,1,3])); for(x=1,sqrtint(lim\=1), my(y,t); while(y++Charles R Greathouse IV, Jan 20 2022
    

Formula

Extensions

Extended by Ray Chandler, Jan 29 2009

A034022 Numbers that are primitively or imprimitively represented by x^2+xy+y^2, but not both.

Original entry on oeis.org

0, 1, 3, 4, 7, 9, 12, 13, 16, 19, 21, 25, 27, 28, 31, 36, 37, 39, 43, 48, 52, 57, 61, 63, 64, 67, 73, 75, 76, 79, 81, 84, 91, 93, 97, 100, 103, 108, 109, 111, 112, 117, 121, 124, 127, 129, 133, 139, 144, 148, 151, 156, 157, 163, 171, 172, 175, 181, 183, 189, 192, 193, 196, 199, 201, 208, 211, 217, 219, 223, 225
Offset: 1

Views

Author

Keywords

Comments

a(n) = A198772(n) for n <= 32. - Reinhard Zumkeller, Oct 30 2011

Crossrefs

Symmetric difference of A034017 and A034019.
After the initial 0, differs from A329963 next time at a(63) = 196, term which is not present in the latter.

Programs

  • PARI
    prim(f)=for(i=1, #f~, if(f[i, 1]%3!=1 && (f[i, 1]!=3 || f[i, 2]>1), return(factorback(f)==0))); 1
    imprim(f)=my(t); for(i=1, #f~, if(f[i, 1]%3<2 && f[i, 2]>1, t=1); if(f[i, 1]%3==2, if(f[i, 2]%2, return(0), t=1))); t
    is(n)=my(f=factor(n)); prim(f)+imprim(f)==1 \\ Charles R Greathouse IV, Nov 04 2015

Extensions

Extended by Ray Chandler, Jan 29 2009
Data section further extended up to a(71), to better differentiate from nearby sequences, by Antti Karttunen, Jul 04 2024

A045612 Number of different energy states between two stacked necklaces, each of n positive and n negative charges.

Original entry on oeis.org

1, 2, 5, 18, 71, 538, 2619, 33052, 233613, 3153926, 24581647
Offset: 0

Views

Author

Keywords

Comments

For this sequence, only the distances between the two necklaces are taken into account (that is, distances within a necklace are ignored). - Sean A. Irvine, Mar 18 2021

Crossrefs

Cf. A045611.

Extensions

a(6)-a(10) from Sean A. Irvine, Mar 18 2021

A103441 Triangle read by rows: T(n,k) = number of bracelets of n beads (necklaces that can be flipped over) with exactly two colors and k white beads for which the set of distances among the white beads are different.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 3, 3, 3, 1, 1, 3, 4, 4, 3, 1, 1, 4, 5, 7, 5, 4, 1, 1, 4, 7, 10, 10, 7, 4, 1, 1, 5, 8, 16, 13, 16, 8, 5, 1, 1, 5, 10, 20, 26, 26, 20, 10, 5, 1, 1, 6, 12, 28, 35, 35, 35, 28, 12, 6, 1, 1, 6, 14, 34, 57, 74, 74, 57, 34, 14, 6, 1, 1, 7, 16, 47, 73, 120, 85, 120, 73
Offset: 2

Views

Author

Wouter Meeussen, Feb 06 2005

Keywords

Comments

If two bracelets can be made to coincide by rotation or flipping over they necessarily have the same set of distances, but the reverse is obviously not true.
Offset is 2, since exactly two colors are required, ergo at least two beads.
T[2n,n] equals A045611. Row sums equal A103442.
Same as A052307, except for bracelets such as {0,0,0,1,1,0,1,1} and{0,0,1,0,0,1,1,1}, that both have the same set of distances between the "1" beads: 4 d[0]+ 4 d[1]+ 2 d[2]+ 4 d[3]+ 2 d[4], where d[k] represents the unidirectional distance between two beads k places apart.

Examples

			Table starts as
  1;
  1,1;
  1,2,1;
  1,2,2,1;
  ...
		

Crossrefs

Programs

  • Mathematica
    Needs[DiscreteMath`NewCombinatorica`]; f[bi_]:=DeleteCases[bi Range[Length[bi]], 0]; dist[li_, l_]:=Plus@@Flatten[Outer[d[Min[ #, l-# ]&@Mod[Abs[ #1-#2], l, 0]]&, li, li]]; Table[Length[Union[(dist[f[ #1], n]&)/@ListNecklaces[n, Join[1+0*Range[i], 0*Range[n-i]], Dihedral]]], {n, 2, 16}, {i, 1, n-1}]
Showing 1-4 of 4 results.