Java 8+ , 1044 байта, последовательность A008008 (Safe)
class c{long[]u={1,4,11,21,35,52,74,102,136,172,212,257,306,354,400,445,488,529,563,587,595,592,584,575,558,530,482,421,354,292,232,164,85,0,-85,-164,-232,-292,-354,-421,-482,-530,-558,-575,-584,-592,-595,-587,-563,-529,-488,-445,-400,-354,-306,-257,-212,-172,-136,-102,-74,-52,-35,-21,-11,-4,-1},v={0,0,0,0,0,0,0,0,0,0,0,0,-1,0,-1,0,-1,0,0,0,0,0,0,0,-1,0,1,0,1,0,1,0,0,0,0,0,1,0,1,0,1,0,-1,0,0,0,0,0,0,0,-1,0,-1,0,-1,0,0,0,0,0,0,0,0,0,0,0,1},w={1,0,0,-1,5};long d=1,e=1;void f(long a,long b){long[]U=u,V=v,W,X;while(a-->0){U=g(U);w=h(v,w);}W=h(v,U);while(b-->0){V=g(V);v=h(v,v);}X=h(V,u);if(w[0]!=v[0]){int i,j,k=0;u=new long[i=(i=W.length)>(j=X.length)?i:j];for(;k<i;k++)u[k]=(k<i?W[k]:0)-(k<j?X[k]:0);d*=e++;}}long[]g(long[]y){int s=y.length,i=1;long[]Y=new long[s-1];for(;i<s;){Y[i-1]=y[i]*i++;}return Y;}long[]h(long[]x,long[]y){int q=x.length,r=y.length,i=0,j;long[]z=new long[q+r-1];for(;i<q;i++)if(x[i]!=0)for(j=0;j<r;)z[i+j]+=x[i]*y[j++];return z;}c(){f(3,0);System.out.println(u[0]/d);}public static void main(String[]args){new c();}}
Попробуйте онлайн!
Может быть решена с помощью скрытой строки размером 12 . Можно определенно больше играть в гольф, но это никоим образом не выигрывает. Я просто хотел внести свой вклад из уважения к числу 8008.
Примечание: прежде чем кто-либо будет жаловаться, что последовательность жестко запрограммирована, я проверил это до первого члена, который отличается от жесткого кода (13-й член = 307), и он получает его правильно, хотя и медленно. Это также почему он используетlong вместо int, в противном случае он переполняется до этого срока.
Обновление (12 июля 2019 г.) : обновлено, чтобы быть немного более производительным. Теперь вычисляет 13-й срок за 30 секунд на моем компьютере вместо 5 минут.
Обновление (17 июля 2019 г.) : исправлены ошибки для границ цикла для gфункции и границ длины массива в нижней частиf функции. Эти ошибки в конечном итоге должны были вызвать проблемы, но не настолько рано, чтобы их можно было обнаружить, просто проверив выходные данные. В любом случае, поскольку наличие этих ошибок в течение 5 дней в игре могло бы запутать некоторых людей настолько, что они не смогли решить эту загадку, я полностью согласен с продлением «безопасного» срока до 24 июля для этого представления.
Обновление (18 июля 2019 г.) : После некоторого тестирования я подтвердил, что переполнения начинаются после 4-го члена последовательности и начинают влиять на достоверность выходных данных после 19-го члена. Кроме того, в программе, как написано здесь, каждый последующий член занимает примерно в 5 раз больше времени, чем предыдущий. 15 семестр занимает около 14 минут на моем компьютере. Таким образом, на самом деле вычисление 19-го семестра с использованием программы, как написано, заняло бы более 6 дней.
Кроме того, вот код с нормальным интервалом / отступом, поэтому его немного легче читать, если у людей нет IDE с автоматическим форматированием под рукой.
class c {
long[] u = {1, 4, 11, 21, 35, 52, 74, 102, 136, 172, 212, 257, 306, 354, 400, 445, 488, 529, 563, 587, 595, 592, 584,
575, 558, 530, 482, 421, 354, 292, 232, 164, 85, 0, -85, -164, -232, -292, -354, -421, -482, -530, -558, -575,
-584, -592, -595, -587, -563, -529, -488, -445, -400, -354, -306, -257, -212, -172, -136, -102, -74, -52, -35,
-21, -11, -4, -1},
v = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, -1, 0, -1, 0, 0, 0, 0, 0, 0, 0, -1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0,
0, 1, 0, 1, 0, 1, 0, -1, 0, 0, 0, 0, 0, 0, 0, -1, 0, -1, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1},
w = {1, 0, 0, -1, 5};
long d = 1, e = 1;
void f(long a, long b) {
long[] U = u, V = v, W, X;
while (a-- > 0) {
U = g(U);
w = h(v, w);
}
W = h(v, U);
while (b-- > 0) {
V = g(V);
v = h(v, v);
}
X = h(V, u);
if (w[0] != v[0]) {
int i, j, k = 0;
u = new long[i = (i = W.length) > (j = X.length) ? i : j];
for (; k < i; k++)
u[k] = (k < i ? W[k] : 0) - (k < j ? X[k] : 0);
d *= e++;
}
}
long[] g(long[] y) {
int s = y.length, i = 1;
long[] Y = new long[s - 1];
for (; i < s;) {
Y[i - 1] = y[i] * i++;
}
return Y;
}
long[] h(long[] x, long[] y) {
int q = x.length, r = y.length, i = 0, j;
long[] z = new long[q + r - 1];
for (; i < q; i++)
if (x[i] != 0)
for (j = 0; j < r;)
z[i + j] += x[i] * y[j++];
return z;
}
c() {
f(3, 0);
System.out.println(u[0] / d);
}
public static void main(String[] args) {
new c();
}
}
Решение
f(1,v[0]=1);прямо перед тем, как System.out.println
программа работает путем вычисления n-го коэффициента расширения Тейлора в 0. Где исходная функция является частным от многочленов, представленных uи vкоторые я получил отсюда , за исключением того, что в связанном документе знаменатель не умножается, и нигде не говорят, что вы должны вычислить ряд Тейлора, я случайно наткнулся на это, а затем подтвердил из другого источника.
Расчет производится путем многократного применения коэффициента для производных.
Неправильный первый член v, весь массив wи некоторые другие вещи, такие как функция, fимеющая какие-либо аргументы, добавляются, чтобы связываться с людьми.