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

A359402 Numbers whose binary expansion and reversed binary expansion have the same sum of positions of 1's, where positions in a sequence are read starting with 1 from the left.

Original entry on oeis.org

0, 1, 3, 5, 7, 9, 15, 17, 21, 27, 31, 33, 45, 51, 63, 65, 70, 73, 78, 85, 93, 99, 107, 119, 127, 129, 150, 153, 165, 189, 195, 219, 231, 255, 257, 266, 273, 282, 294, 297, 310, 313, 325, 334, 341, 350, 355, 365, 371, 381, 387, 397, 403, 413, 427, 443, 455, 471
Offset: 1

Views

Author

Gus Wiseman, Jan 05 2023

Keywords

Comments

Also numbers whose binary expansion and reversed binary expansion have the same sum of partial sums.
Also numbers whose average position of a 1 in their binary expansion is (c+1)/2, where c is the number of digits.
Conjecture: Also numbers whose binary expansion has as least squares fit a line of zero slope, counted by A222955.

Examples

			The binary expansion of 70 is (1,0,0,0,1,1,0), with positions of 1's {1,5,6}, while the reverse positions are {2,3,7}. Both sum to 12, so 70 is in the sequence.
		

Crossrefs

Binary words of this type appear to be counted by A222955.
For greater instead of equal sums we have A359401.
These are the indices of 0's in A359495.
A030190 gives binary expansion, reverse A030308.
A048793 lists partial sums of reversed standard compositions, sums A029931.
A070939 counts binary digits, 1's A000120.
A326669 lists numbers with integer mean position of a 1 in binary expansion.

Programs

  • Mathematica
    Select[Range[0,100],#==0||Mean[Join@@Position[IntegerDigits[#,2],1]]==(IntegerLength[#,2]+1)/2&]
  • Python
    from functools import reduce
    from itertools import count, islice
    def A359402_gen(startvalue=0): # generator of terms
        return filter(lambda n:(r:=reduce(lambda c, d:(c[0]+d[0]*(e:=int(d[1])),c[1]+e),enumerate(bin(n)[2:],start=1),(0,0)))[0]<<1==(n.bit_length()+1)*r[1],count(max(startvalue,0)))
    A359402_list = list(islice(A359402_gen(),30)) # Chai Wah Wu, Jan 08 2023

Formula

A230877(a(n)) = A029931(a(n)).

A359676 Least positive integer whose weakly increasing prime indices have zero-based weighted sum n (A359674).

Original entry on oeis.org

1, 4, 6, 8, 14, 12, 16, 20, 30, 24, 32, 36, 40, 52, 48, 56, 100, 72, 80, 92, 96, 104, 112, 124, 136, 148, 176, 152, 214, 172, 184, 188, 262, 212, 272, 236, 248, 244, 304, 268, 346, 284, 328, 292, 386, 316, 398, 332, 376, 356, 458, 388, 478, 404, 472, 412, 526
Offset: 1

Views

Author

Gus Wiseman, Jan 14 2023

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.
The zero-based weighted sum of a sequence (y_1,...,y_k) is Sum_{i=1..k} (i-1)*y_i.

Examples

			The terms together with their prime indices begin:
    1: {}
    4: {1,1}
    6: {1,2}
    8: {1,1,1}
   14: {1,4}
   12: {1,1,2}
   16: {1,1,1,1}
   20: {1,1,3}
   30: {1,2,3}
   24: {1,1,1,2}
   32: {1,1,1,1,1}
   36: {1,1,2,2}
   40: {1,1,1,3}
   52: {1,1,6}
   48: {1,1,1,1,2}
		

Crossrefs

First position of n in A359674, reverse A359677.
The sorted version is A359675, reverse A359680.
The reverse one-based version is A359679, sorted A359754.
The reverse version is A359681.
The one-based version is A359682, sorted A359755.
The version for standard compositions is A359756, one-based A089633.
A053632 counts compositions by zero-based weighted sum.
A112798 lists prime indices, length A001222, sum A056239.
A124757 gives zero-based weighted sum of standard compositions, rev A231204.
A304818 gives weighted sums of prime indices, reverse A318283.
A320387 counts multisets by weighted sum, zero-based A359678.
A358136 lists partial sums of prime indices, ranked by A358137, rev A359361.

Programs

  • Mathematica
    nn=20;
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    wts[y_]:=Sum[(i-1)*y[[i]],{i,Length[y]}];
    seq=Table[wts[primeMS[n]],{n,1,Prime[nn]^2}];
    Table[Position[seq,k][[1,1]],{k,0,nn}]

A359677 Zero-based weighted sum of the reversed (weakly decreasing) prime indices of n.

Original entry on oeis.org

0, 0, 0, 1, 0, 1, 0, 3, 2, 1, 0, 3, 0, 1, 2, 6, 0, 4, 0, 3, 2, 1, 0, 6, 3, 1, 6, 3, 0, 4, 0, 10, 2, 1, 3, 7, 0, 1, 2, 6, 0, 4, 0, 3, 6, 1, 0, 10, 4, 5, 2, 3, 0, 9, 3, 6, 2, 1, 0, 7, 0, 1, 6, 15, 3, 4, 0, 3, 2, 5, 0, 11, 0, 1, 7, 3, 4, 4, 0, 10, 12, 1, 0, 7, 3
Offset: 1

Views

Author

Gus Wiseman, Jan 13 2023

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.
The zero-based weighted sum of a sequence (y_1,...,y_k) is Sum_{i=1..k} (i-1)*y_i.

Examples

			The reversed prime indices of 12 are (2,1,1), so a(12) = 0*2 + 1*1 + 2*1 = 3.
		

Crossrefs

Positions of 0's are A008578.
Positions of 1's are A100484.
The version for standard compositions is A231204, reverse of A124757.
The one-based version is A318283, unreversed A304818.
The one-based version for standard compositions is A359042, rev of A029931.
This is the reverse version of A359674.
First position of n is A359679(n), reverse of A359675.
Positions of first appearances are A359680, reverse of A359676.
A053632 counts compositions by weighted sum.
A112798 lists prime indices, length A001222, sum A056239.
A358136 lists partial sums of prime indices, ranked by A358137, rev A359361.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    wts[y_]:=Sum[(i-1)*y[[i]],{i,Length[y]}];
    Table[wts[Reverse[primeMS[n]]],{n,100}]

A359681 Least positive integer whose reversed (weakly decreasing) prime indices have zero-based weighted sum (A359677) equal to n.

Original entry on oeis.org

1, 4, 9, 8, 18, 50, 16, 36, 100, 54, 32, 72, 81, 108, 300, 64, 144, 400, 216, 600, 243, 128, 288, 800, 432, 486, 1350, 648, 256, 576, 729, 864, 2400, 3375, 1296, 3600, 512, 1152, 1944, 1728, 4800, 9000, 2187, 2916, 8100, 1024, 2304, 6400, 3456, 4374, 12150
Offset: 0

Views

Author

Gus Wiseman, Jan 15 2023

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.
The zero-based weighted sum of a sequence (y_1,...,y_k) is Sum_{i=1..k} (i-1)*y_i.

Examples

			The terms together with their prime indices begin:
    1: {}
    4: {1,1}
    9: {2,2}
    8: {1,1,1}
   18: {1,2,2}
   50: {1,3,3}
   16: {1,1,1,1}
   36: {1,1,2,2}
  100: {1,1,3,3}
   54: {1,2,2,2}
   32: {1,1,1,1,1}
   72: {1,1,1,2,2}
   81: {2,2,2,2}
  108: {1,1,2,2,2}
  300: {1,1,2,3,3}
		

Crossrefs

The unreversed version is A359676.
First position of n in A359677, reverse A359674.
The one-based version is A359679, sorted A359754.
The sorted version is A359680, reverse A359675.
The unreversed one-based version is A359682, sorted A359755.
A053632 counts compositions by zero-based weighted sum.
A112798 lists prime indices, length A001222, sum A056239.
A124757 gives zero-based weighted sum of standard compositions, rev A231204.
A304818 gives weighted sum of prime indices, reverse A318283.
A320387 counts multisets by weighted sum, zero-based A359678.

Programs

  • Mathematica
    nn=20;
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    wts[y_]:=Sum[(i-1)*y[[i]],{i,Length[y]}];
    seq=Table[wts[Reverse[primeMS[n]]],{n,1,Prime[nn]^2}];
    Table[Position[seq,k][[1,1]],{k,0,nn}]

A066571 Number of sets of positive integers with arithmetic mean n.

Original entry on oeis.org

1, 3, 9, 31, 117, 479, 2061, 9183, 42021, 196239, 931457, 4480531, 21793257, 107004891, 529656765, 2640160039, 13241371629, 66771501151, 338333343825, 1721768732423, 8796192611917, 45096680384635, 231945566136129, 1196461977291959, 6188390166782849
Offset: 1

Views

Author

Amarnath Murthy, Dec 19 2001

Keywords

Comments

From Franklin T. Adams-Watters, Sep 13 2011: (Start)
If we use nonnegative integers instead of positive integers, we get this sequence shifted left (i.e., with offset 0).
The largest number that can be included in set of positive integers with mean n is the triangular number n*(n+1)/2 = A000217(n).
All values are odd. Sets including n are paired with the same set with n removed, with exception of {n}, as the empty set has no average.
(End)

Examples

			a(2) = 3 as there are three sets viz. {2}, {1,3}, {1,2,3}, each of which has the arithmetic mean 2.
a(3) = 9: the nine sets are {3}, {1, 5}, {2, 4}, {1, 2, 6}, {1, 3, 5}, {2, 3, 4}, {1, 2, 3, 6}, {1, 2, 4, 5}, {1, 2, 3, 4, 5}.
		

Crossrefs

Programs

  • Haskell
    a066571 n = f [1..] 1 n 0 where
       f (k:ks) l nl x
         | y > nl  = 0
         | y < nl  = f ks (l + 1) (nl + n) y + f ks l nl x
         | otherwise = if y `mod` l == 0 then 1 else 0
         where y = x + k
    -- Reinhard Zumkeller, Feb 13 2013
  • Maple
    g := k->expand(mul(1+t*x^i,i=1..k)); A066571 := proc(n) local k; add(coeff(coeff(g(n*k),t,k),x,n*k),k=1..2*n-1); end;
  • Mathematica
    g[k_] := Expand[Product[1 + t*x^i, {i, 1, k}]]; a[n_] := Sum[Coefficient[ Coefficient[g[n*k], t, k], x, n*k], {k, 1, 2*n - 1}]; Table[an = a[n]; Print["a(", n, ") = ", an]; an, {n, 1, 10}] (* Jean-François Alcover, Feb 10 2018, translated from Maple *)

Formula

Sum of coefficient of t^k x^(n*k) in Product_{i=1..n*k} (1+t*x^i) for k <= 2*n-1. - N. J. A. Sloane
From Martin Fuller, Sep 14 2023: (Start)
Constant term in formal Laurent series (Product_{i=1-n..n*(n-1)/2} (1+x^i)) - 1.
a(n) = (Sum_{i=0..n*(n-1)/2} A053632(n-1,i)*A000009(i))*2-1. (End)

Extensions

Corrected and extended by N. J. A. Sloane, Dec 19 2001
More terms from Naohiro Nomoto, Jun 19 2002
More terms from David Wasserman, Sep 10 2002
More terms from Martin Fuller, Sep 14 2023

A126024 Number of subsets of {1,2,3,...,n} whose sum is a square integer (including the empty subset).

Original entry on oeis.org

1, 2, 2, 3, 5, 7, 12, 20, 34, 60, 106, 190, 346, 639, 1183, 2204, 4129, 7758, 14642, 27728, 52648, 100236, 191294, 365827, 700975, 1345561, 2587057, 4981567, 9605777, 18546389, 35851756, 69382558, 134414736, 260658770, 505941852, 982896850
Offset: 0

Views

Author

John W. Layman, Feb 27 2007

Keywords

Examples

			The subsets of {1,2,3,4,5} that sum to a square are {}, {1}, {1,3}, {4}, {2,3,4}, {1,3,5} and {4,5}. Thus a(5)=7.
		

Crossrefs

Cf. A181522. - Reinhard Zumkeller, Oct 27 2010
Row sums of A281871.

Programs

  • Haskell
    import Data.List (subsequences)
    a126024 = length . filter ((== 1) . a010052 . sum) .
                              subsequences . enumFromTo 1
    -- Reinhard Zumkeller, Feb 22 2012, Oct 27 2010
  • Maple
    b:= proc(n, i) option remember; (m->
          `if`(n=0 or n=m, 1, `if`(n<0 or n>m, 0, b(n, i-1)+
          `if`(i>n, 0, b(n-i, i-1)))))(i*(i+1)/2)
        end:
    a:= proc(n) option remember; `if`(n<0, 0, a(n-1)+
          add(b(j^2-n, n-1), j=isqrt(n)..isqrt(n*(n+1)/2)))
        end:
    seq(a(n), n=0..50);  # Alois P. Heinz, Feb 02 2017
  • Mathematica
    g[n_] := Block[{p = Product[1 + z^i, {i, n}]},Sum[Boole[IntegerQ[Sqrt[k]]]*Coefficient[p, z, k], {k, 0, n*(n + 1)/2}]];Array[g, 35] (* Ray Chandler, Mar 05 2007 *)

Extensions

Extended by Ray Chandler, Mar 05 2007
a(0)=1 prepended by Alois P. Heinz, Jan 30 2017

A359495 Sum of positions of 1's in binary expansion minus sum of positions of 1's in reversed binary expansion, where positions in a sequence are read starting with 1 from the left.

Original entry on oeis.org

0, 0, -1, 0, -2, 0, -2, 0, -3, 0, -2, 1, -4, -1, -3, 0, -4, 0, -2, 2, -4, 0, -2, 2, -6, -2, -4, 0, -6, -2, -4, 0, -5, 0, -2, 3, -4, 1, -1, 4, -6, -1, -3, 2, -5, 0, -2, 3, -8, -3, -5, 0, -7, -2, -4, 1, -9, -4, -6, -1, -8, -3, -5, 0, -6, 0, -2, 4, -4, 2, 0, 6
Offset: 0

Views

Author

Gus Wiseman, Jan 05 2023

Keywords

Comments

Also the sum of partial sums of reversed binary expansion minus sum of partial sums of binary expansion.

Examples

			The binary expansion of 158 is (1,0,0,1,1,1,1,0), with positions of 1's {1,4,5,6,7} with sum 23, reversed {2,3,4,5,8} with sum 22, so a(158) = 1.
		

Crossrefs

Indices of positive terms are A359401.
Indices of 0's are A359402.
A030190 gives binary expansion, reverse A030308.
A070939 counts binary digits.
A230877 adds up positions of 1's in binary expansion, reverse A029931.

Programs

  • Maple
    a:= n-> (l-> add(i*(l[-i]-l[i]), i=1..nops(l)))(Bits[Split](n)):
    seq(a(n), n=0..127);  # Alois P. Heinz, Jan 09 2023
  • Mathematica
    sap[q_]:=Sum[q[[i]]*(2i-Length[q]-1),{i,Length[q]}];
    Table[sap[IntegerDigits[n,2]],{n,0,100}]
  • Python
    def A359495(n):
        k = n.bit_length()-1
        return sum((i<<1)-k for i, j in enumerate(bin(n)[2:]) if j=='1') # Chai Wah Wu, Jan 09 2023

Formula

a(n) = A029931(n) - A230877(n).
If n = Sum_{i=1..k} q_i * 2^(i-1), then a(n) = Sum_{i=1..k} q_i * (2i-k-1).

A059529 For 1 < x, each c(i) is "multiply" (*) or "divide" (/); a(n) is number of choices for c(0),...,c(n-1) so that 1 c(0) x^1 c(1) x^2,.., c(n-1) x^n is an integer.

Original entry on oeis.org

1, 1, 2, 5, 9, 16, 32, 68, 135, 256, 512, 1059, 2110, 4096, 8192, 16745, 33425, 65536, 131072, 266254, 531924, 1048576, 2097152, 4244214, 8482454, 16777216, 33554432, 67741466, 135417620, 268435456, 536870912, 1082015434, 2163280087, 4294967296, 8589934592
Offset: 0

Views

Author

Naohiro Nomoto, Feb 16 2001

Keywords

Comments

From Gus Wiseman, Jul 04 2019: (Start)
Also the number of subsets of {1..n} whose sum is less than or equal to the sum of their complement. For example, the a(0) = 1 through a(5) = 16 subsets are:
{} {} {} {} {} {}
{1} {1} {1} {1}
{2} {2} {2}
{3} {3} {3}
{1,2} {4} {4}
{1,2} {5}
{1,3} {1,2}
{1,4} {1,3}
{2,3} {1,4}
{1,5}
{2,3}
{2,4}
{2,5}
{3,4}
{1,2,3}
{1,2,4}
(End)

Examples

			x = 3: for n = 2 there are 2 possibilities: 1*3*9=27 and 1/3*9=3. For n = 4 there are 9 possibilities: 1*3*9*27*81 1/3*9*27*81 1*3/9*27*81 1/3/9*27*81 1*3*9/27*81 1*3*9*27/81 1/3*9/27*81 1/3*9*27/81 1*3/9/27*81
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]],Plus@@Complement[Range[n],#]>=Plus@@#&]],{n,0,10}] (* Gus Wiseman, Jul 04 2019 *)

Formula

a(0)=1; for 0A058377(n)+2^(n-1).

Extensions

More terms from Alois P. Heinz, Jun 13 2019

A222955 Number of nX1 0..1 arrays with every row and column least squares fitting to a zero slope straight line, with a single point array taken as having zero slope.

Original entry on oeis.org

2, 2, 4, 4, 8, 8, 20, 18, 52, 48, 152, 138, 472, 428, 1520, 1392, 5044, 4652, 17112, 15884, 59008, 55124, 206260, 193724, 729096, 688008, 2601640, 2465134, 9358944, 8899700, 33904324, 32342236, 123580884, 118215780, 452902072, 434314138, 1667837680
Offset: 1

Views

Author

R. H. Hardin, Mar 10 2013

Keywords

Comments

Column 1 of A222959
Conjecture: A binary word is counted iff it has the same sum of positions of 1's as its reverse, or, equivalently, the same sum of partial sums as its reverse. - Gus Wiseman, Jan 07 2023

Examples

			All solutions for n=4
..0....1....1....0
..0....1....0....1
..0....1....0....1
..0....1....1....0
From _Gus Wiseman_, Jan 07 2023: (Start)
The a(1) = 2 through a(7) = 20 binary words with least squares fit a line of zero slope are:
  (0)  (00)  (000)  (0000)  (00000)  (000000)  (0000000)
  (1)  (11)  (010)  (0110)  (00100)  (001100)  (0001000)
             (101)  (1001)  (01010)  (010010)  (0010100)
             (111)  (1111)  (01110)  (011110)  (0011100)
                            (10001)  (100001)  (0100010)
                            (10101)  (101101)  (0101010)
                            (11011)  (110011)  (0110001)
                            (11111)  (111111)  (0110110)
                                               (0111001)
                                               (0111110)
                                               (1000001)
                                               (1000110)
                                               (1001001)
                                               (1001110)
                                               (1010101)
                                               (1011101)
                                               (1100011)
                                               (1101011)
                                               (1110111)
                                               (1111111)
(End)
		

Crossrefs

These words appear to be ranked by A359402.
A011782 counts compositions.
A359042 adds up partial sums of standard compositions, reversed A029931.

A359679 Least number with weighted sum of reversed (weakly decreasing) prime indices (A318283) equal to n.

Original entry on oeis.org

1, 2, 3, 4, 6, 10, 8, 12, 19, 18, 16, 24, 27, 36, 43, 32, 48, 59, 61, 67, 71, 64, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269
Offset: 0

Views

Author

Gus Wiseman, Jan 14 2023

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.
The weighted sum of a sequence (y_1,...,y_k) is Sum_{i=1..k} i*y_i.

Examples

			12 has reversed prime indices (2,1,1), with weighted sum 7, and no number < 12 has the same weighted sum of reversed prime indices, so a(7) = 12.
		

Crossrefs

The version for standard compositions is A089633, zero-based A359756.
First position of n in A318283, unreversed A304818.
The unreversed zero-based version is A359676.
The sorted zero-based version is A359680, unreversed A359675.
The zero-based version is A359681.
The unreversed version is A359682.
The greatest instead of least is A359683, unreversed A359497.
The sorted version is A359754, unreversed A359755.
A112798 lists prime indices, length A001222, sum A056239.
A320387 counts multisets by weighted sum, zero-based A359678.
A358136 lists partial sums of prime indices, ranked by A358137, rev A359361.

Programs

  • Mathematica
    nn=20;
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    ots[y_]:=Sum[i*y[[i]],{i,Length[y]}];
    seq=Table[ots[Reverse[primeMS[n]]],{n,1,Prime[nn]^2}];
    Table[Position[seq,k][[1,1]],{k,0,nn}]
Previous Showing 21-30 of 80 results. Next