|
（PARI）
\\广义Wythoff阵列W（h）的第i行，
\\其中h是介于1和2之间的无理数，
\\m是向量a和b中的项数。
行（h，i，m）={
如果（h<=1|h>=2，打印（“h的值无效”）；返回）；
我的(
a=矢量（m，n，地板（n*h）），
b=矢量（m，n，地板（n*h/（h-1））），
w=[a[a[i]]，b[a[i]]，
j=3
);
而（1，
如果（j%2==1，
如果（w[j-1]<=#a，w=concat（w，a[w[j-1]]），返回（w））
,
如果（w[j-2]<=#b，w=concat（w，b[w[j-2-]），返回（w））
);
j个++
)
}
分配（10^9）
行（exp（1）^Euler，1，10^7）
| 