%I#17 2023年8月30日02:00:40
%S 1,5,2935831499819682426373122736784270089158028202000,
%电话2853152994118553866904840963028755855938535876921233326141376,
%电话:111916280261483009775739858744965571138908168272402828091266521771625215094103508198780449024
%N a(N)=(N-1)!*求和{d|n}d^n/(d-1)!。
%如果p是素数,a(p)=p^p+(p-1)!。
%例如:-求和{k>0}log(1-(k*x)^k)/k!。
%t a[n]:=(n-1)!*除数和[n,#^n/(#-1)!&];阵列[a,18](*_Amiram Eldar_,2023年8月30日*)
%o(PARI)a(n)=(n-1)*sumdiv(n,d,d^n/(d-1)!);
%o(PARI)我的(N=20,x='x+o('x^N));Vec(serlaplace(-sum(k=1,N,log(1-(k*x)^k)/k!))
%Y参考A087906、A354890、A356487。
%K nonn公司
%O 1,2号机组
%A _Seiichi Manyama,2022年8月9日
|