(PARI)
up_to=65537;
rgs_transform(invec)={my(om=Map(),outvec=vector(length(invec)),u=1);对于(i=1,长度(invesc),如果(mapisdefined(om,invec[i]),my(pp=mapget(om,invec[i];
A003415号(n) =如果(n<=1,0,my(f=系数(n));n*和(i=1,#f~,f[i,2]/f[i,1]);
v353522=rgs_transform(向量(up_to,n,Aux353522(n)));
|