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
,否则它将在该术语之前溢出。
更新(2019年7月12日):已更新,性能更高。现在可以在30秒内在我的计算机上计算第13个术语,而不是5分钟。
更新(2019年7月17日):修复了该g
函数的for循环边界以及该函数底部的数组长度边界的错误f
。这些错误最终应该导致了问题,但还不够早,仅通过检查输出就可以发现它们。无论哪种情况,由于游戏中存在5天的错误可能使某些人困惑不解,所以我完全可以将“安全”截止日期延长至7月24日。
更新(2019年7月18日):经过一些测试,我确认溢出在序列的第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
为0。该项目工程通过计算第n泰勒膨胀系数凡原函数是多项式为代表的商u
以及v
这是我从拿到这里,除了链接文档中的分母不出来倍增,他们无处不说您必须计算泰勒级数,我偶然发现了泰勒级数,然后通过另一个来源进行了确认。
通过对导数重复应用商规则来完成计算。,整个数组
的不正确的第一项以及诸如带有任何参数的函数之类的其他一些东西都被扔进了人们的视线。v
w
f