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.

A333590 a(n) = a(n-1) if half of the previous term pairs are inverted. a(n) = a(n-1) + 1 if more than half of the previous term pairs are inverted. a(n) = a(n-1) - 1 if fewer than half of the previous term pairs are inverted. a(1) = 0.

This page as a plain text file.
%I A333590 #25 Dec 20 2024 23:37:37
%S A333590 0,0,-1,0,-1,-1,0,-1,-2,-1,0,-1,-2,-1,-2,-1,0,-1,-2,-3,-2,-1,0,-1,-2,
%T A333590 -3,-2,-1,0,-1,-2,-3,-2,-1,0,-1,-2,-3,-4,-3,-2,-1,0,1,0,-1,-2,-3,-4,
%U A333590 -5,-4,-3,-2,-1,0,1,0,-1,-2,-3,-4,-3,-2,-1,0,-1,-2,-3,-4,-3,-2,-1
%N A333590 a(n) = a(n-1) if half of the previous term pairs are inverted. a(n) = a(n-1) + 1 if more than half of the previous term pairs are inverted. a(n) = a(n-1) - 1 if fewer than half of the previous term pairs are inverted. a(1) = 0.
%C A333590 An inversion occurs when some term is greater than some later term.
%C A333590 A sequence of length n can have at most n * (n - 1) / 2 inversions.
%H A333590 Samuel B. Reid, <a href="/A333590/b333590.txt">Table of n, a(n) for n = 1..10000</a>
%H A333590 Samuel B. Reid, <a href="/A333590/a333590.png">Plot of one billion terms</a>
%H A333590 Samuel B. Reid, <a href="/A333590/a333590_1.c.txt">C program for A333590</a>
%H A333590 Rémy Sigrist, <a href="/A333590/a333590_1.png">Scatterplot of the ordinal transform of the first 10000000 terms</a>
%e A333590 There are zero inversions in the first two terms of this sequence. The maximum possible number of inversions in a sequence of length 2 is 1. 0 is less than half of 1, so the third term is 0 - 1 or -1.
%e A333590 There are two inversions in the first 3 terms of this sequence. The maximum possible number of inversions in a sequence of length 3 is 3. 2 is more than half of 3, so the fourth term is -1 + 1 or 0.
%o A333590 (C) // See Links section.
%o A333590 (PARI) { my(v=0, f=vector(2*M=100), s=0, inv=0); for (n=1, 72, f[M+v]++; inv+=s; print1 (v", "); if (2*inv<t=n*(n-1)/2, s+=f[M+v]; v--, 2*inv>t, v++; s-=f[M+v])) } \\ _Rémy Sigrist_, Mar 28 2020
%Y A333590 Cf. A323186.
%K A333590 sign
%O A333590 1,9
%A A333590 _Samuel B. Reid_, Mar 27 2020