A171898 - OEIS

%I #33 Oct 08 2023 18:06:52

%S 1,2,6,2,2,5,1,6,42,5,2,4,5,9,14,3,9,3,15,2,4,6,17,3,6,32,56,5,3,131,

%T 5,11,5,3,20,6,2,8,15,31,170,3,31,18,3,3,33,5,1,11,46,56,4,37,152,307,

%U 3,7,92,4,7,62,52,3,42,3,6,2,19,6,8,3,9,3,650,2,23,8,223,7,206,3,21,25,5,8

%N Forward van Eck transform of A181391.

%C Given a sequence a, the forward van Eck transform b is defined as follows: If a(n) also appears again in a later position, let a(m) be the next occurrence, and set b(n)=m-n; otherwise b(n)=0.

%C This is a permutation of the positive terms in A181391, where each term m > 0 from that sequence is shifted backwards m+1 positions. - _Jan Ritsema van Eck_, Aug 16 2019

%C The backwards van Eck transform searches backwards for a repeated value: if a(n) also has appeared in earlier positions, a(m)=a(n) with m<n, then b(n) is the minimum n-m. - _R. J. Mathar_, Jun 24 2021

%H Alois P. Heinz, <a href="/A171898/b171898.txt">Table of n, a(n) for n = 1..10000</a>

%F From _Jan Ritsema van Eck_, Aug 16 2019: (Start)

%F A181391(i+a(i)+1) = a(i) for any i, a(i)>0.

%F Conversely, a(j-A(181391(j)-1) = A181391(j) for any j, A181391(j)>0. (End)

%p ECKf:=proc(a) local b,i,m,n;

%p if whattype(a) <> list then RETURN([]); fi:

%p b:=[];

%p for n from 1 to nops(a)-1 do

%p # does a(n) appear again?

%p m:=0;

%p for i from n+1 to nops(a) do

%p if (a[i]=a[n]) then m:=i-n; break; fi

%p od:

%p b:=[op(b),m];

%p od:

%p b:=[op(b),0];

%p RETURN(b);

%p end:

%t terms = 100;

%t m = 14 terms; (* Increase m until no zero appears in the output *)

%t ClearAll[b, last]; b[_] = 0; last[_] = -1; last[0] = 2; nxt = 1;

%t Do[hist = last[nxt]; b[n] = nxt; last[nxt] = n; nxt = 0; If[hist > 0, nxt = n - hist], {n, 3, m}];

%t A181391 = Array[b, m];

%t ECKf[a_List] := Module[{b = {}, i, m, n}, For[n = 1, n <= Length[a]-1, n++, m = 0; For[i = n+1, i <= Length[a], i++, If[a[[i]] == a[[n]], m = i-n; Break[]]]; b = Append[b, m]]; b = Append[b, 0]; Return[b]];

%t ECKf[A181391][[;; terms]] (* _Jean-François Alcover_, Oct 30 2020, after Maple *)

%Y Cf. A181391 (van Eck's sequence), A171899, A171942.

%K nonn

%O 1,2

%A _N. J. A. Sloane_, Oct 22 2010