%2018年10月24日02:33:05
%S 2,2,3,1,0,3,2,1,2,0,0,0,1,2,3,4,3,3,,3,4,1,3,0,2,2,2,0,1,0,1,
%T 0,0,0,4,0,0,0,1,0,2,3,1,4,0,3,4,4,3,1,3,2,4,4,3,0,3,1,1,2,0,0,3,3,1,2,0,3,
%U 0,2,3,4,3,0,2,2,3,1,4,3,1,2,0,1,0,3,4
%N 3^(1/3)的五元展开式中的系数。
%H Seiichi Manyama,<a href=“/A290563/b290563.txt”>n表,n=0..10000的a(n)</a>
%o(红宝石)
%o需要“OpenSSL”
%o定义fa(ary,a)
%o(0..ary.size-1).inject(0){|s,i|s+ary[i]*a**i}
%o结束
%o定义df(ary)
%o(1..ary.size-1).map{i|i*ary[i]}
%o端
%o定义A(c_ary,k,m,n)
%o x=OpenSSL::BN.new((-f_a(df(c_ary),k)).to_s).mod_inverse(m).to_i%m
%o f元=c元映射{i|x*i}
%o f元[1]+=1
%o _ ary=[]
%o元=[0]
%o a,mod=k,m
%o(n+1).乘以{|i|
%o b=a%mod
%o d元<<(b元[-1])/m**i
%o ary<<b
%o a=fα(f元,b)
%o模块*=m
%o}(o)
%o _ ary(_A)
%o端
%o定义A290563(n)
%o A([-3,0,0,1],2,5,n)
%o端
%o第A290563页(100)
%o(PARI)Vecrev(数字(截断((3+o(5^100))^(1/3)),5))\\_Joerg Arndt_,2017年8月6日
%Y参见A290566和A290568。
%K nonn公司
%0、1
%A _Seiichi Manyama,2017年8月6日
|