Solution:
It is possible, as long as the sum S desired is a multiple of 48, with N=S/6, which in this case is 3360336, and the ak repeats the 8-term pattern −1,1,1,−1,1,−1,−1,1.
Use the observation that if f(x) is a degree-k polynomial, then for any constant h, the difference f(x+h)−f(x) will be a degree-(k−1) polynomial. If we iterate this process three times, we can find a way to manipulate consecutive cubes to always get a constant.
More precisely, let
c0,c1,c2,c3,c4,c5,c6,c7,…
be consecutive cubes. In other words, cm=(m+u)3 for some fixed starting integer u. Then the differences
c1−c0,c3−c2,c5−c4,c7−c6,…
will be quadratic functions; i.e., if we define am:=cm+1−cm, then am is a quadratic function of m (depending on the parameter u, as well), and the differences are a0,a2,a4,a6,… Continuing, we see that the differences
a2−a0,a6−a4,…
will be a linear sequence; i.e., bm:=am+2−am is a linear function of m (with parameter u), and our differences are b0,b4,… Finally, the sequence
b4−b0,b12−b8,…
is constant, no matter what the parameter u equals! We have
b4−b0=a6−a4−(a2−a0)=a6−a4−a2+a0=c7−c6−c5+c4−c3+c2+c1−c0,
and since this is constant, we can compute it using any value of u. Taking u=−3, the constant must equal
43−33−22+13−03+(−1)3+(−2)3−(−3)3=48.
In other words, if we define
su:=−u3+(u+1)3+(u+2)3−(u+3)3+(u+4)3−(u+5)3−(u+6)3+(u+7)3,
then su=48 for all values of u.
Since 20162016=420042⋅48, we can easily write 20162016 as a sum of 420042 8-element sum/differences of 420042⋅8=3360336 consecutive cubes:
20162016=k=0∑420041s8k+1= −13+23+33−43+53−63−73+83−93+103+113−123+133−143−143+153⋯−33603293+33603303+⋯−33603353+33603363.