2) A101880
Вот вам функция, которая очень загадочно, но при этом довольно быстро вычисляет

для
A101880Код:
a(n) =
{
my(k1 = n+1, kx = k1, a = 1, t = vector(k1, i, vector(n)));
t[k1][1] = 1;
for (i = 1, n-1,
my(ky = a);
t[a] = vector(n, j, if (!(ky--), ky = k1);
if (j == 1, t[ky][1]
, t[kx][j-1] + t[ky][j]));
kx = a; if (a++ > k1, a = 1));
a=if(a == 1, k1, a-1);
return(sum(i=1,n,t[a][i]*i!));
}
a(1000) у меня вычисляется за одну секунду.
Но поля этого форума слишком узкие для результата (

)