We prove by induction on k that
uk=0<i1<…<it<kij+1−ij⩾2∑ai1…ait(1)
Note that we have one trivial summand equal to 1 (which corresponds to t=0 and the empty sequence, whose product is 1).
For k=0,1 the sum on the right-hand side only contains the empty product, so (1) holds due to u0=u1=1. For k⩾1, assuming the result is true for 0,1,…,k, we have
uk+1=0<i1<…<it<k,ij+1−ij⩾2∑ai1…ait+0<i1<…<it<k−1,ij+1−ij⩾2∑ai1…ait⋅ak=0<i1<…<it<k+1,ij+1−ij⩾2,k∈/{i1,…,it}∑ai1…ait+0<i1<…<it<k+1,ij+1−ij⩾2,k∈{i1,…,it}∑ai1…ait=0<i1<…<it<k+1,ij+1−ij⩾2∑ai1…ait,
as required.
Applying (1) to the sequence b1,…,bn given by bk=an−k for 1⩽k⩽n, we get
vk=0<i1<…<it<kij+1−ij⩾2∑bi1…bit=n>i1>…>it>n−kij−ij+1⩾2∑ai1…ait(2)
For k=n the expressions (1) and (2) coincide, so indeed un=vn.
Solution 2:
Define recursively a sequence of multivariate polynomials by
P0=P1=1,Pk+1(x1,…,xk)=Pk(x1,…,xk−1)+xkPk−1(x1,…,xk−2),
so Pn is a polynomial in n−1 variables for each n⩾1. Two easy inductive arguments show that
un=Pn(a1,…,an−1),vn=Pn(an−1,…,a1)
so we need to prove Pn(x1,…,xn−1)=Pn(xn−1,…,x1) for every positive integer n. The cases n=1,2 are trivial, and the cases n=3,4 follow from P3(x,y)=1+x+y and P4(x,y,z)=1+x+y+z+xz.
Now we proceed by induction, assuming that n⩾5 and the claim hold for all smaller cases. Using F(a,b) as an abbreviation for P∣a−b∣+1(xa,…,xb) (where the indices a,…,b can be either in increasing or decreasing order),
F(n,1)=F(n,2)+x1F(n,3)=F(2,n)+x1F(3,n)=(F(2,n−1)+xnF(2,n−2))+x1(F(3,n−1)+xnF(3,n−2))=(F(n−1,2)+x1F(n−1,3))+xn(F(n−2,2)+x1F(n−2,3))=F(n−1,1)+xnF(n−2,1)=F(1,n−1)+xnF(1,n−2)=F(1,n)
as we wished to show.
Solution 3:
Using matrix notation, we can rewrite the recurrence relation as
(uk+1−ukuk+1)=(akuk−1uk+akuk−1)=(1+akak−ak−ak)(uk−uk−1uk)
for 1⩽k⩽n−1, and similarly
(vk+1;vk−vk+1)=(vk+an−kvk−1;−an−kvk−1)=(vk;vk−1−vk)(1+an−kan−k−an−k−an−k)
for 1⩽k⩽n−1. Hence, introducing the 2×2 matrices Ak=(1+akak−ak−ak) we have
(uk+1−ukuk+1)=Ak(uk−uk−1uk) and (vk+1;vk−vk+1)=(vk;vk−1−vk)An−k.
for 1⩽k⩽n−1. Since (u1−u0u1)=(01) and (v1;v0−v1)=(1;0), we get
(un−un−1un)=An−1An−2⋯A1⋅(01) and (vn;vn−1−vn)=(1;0)⋅An−1An−2⋯A1.
It follows that
(un)=(1;0)(un−un−1un)=(1;0)⋅An−1An−2⋯A1⋅(01)=(vn;vn−1−vn)(01)=(vn).