OEIS哀悼西蒙斯感谢西蒙斯基金会支持包括OEIS在内的许多科学分支的研究。
登录
OEIS由OEIS基金会的许多慷慨捐赠者.

 

徽标
提示
(来自的问候整数序列在线百科全书!)
A362344飞机 系数为0.1的n次多项式的最大不同实根数。 3
1, 2, 2, 2, 2, 2, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5 (列表;图表;参考;;历史;文本;内部格式)
抵消
1,2
链接
恩里科·贝托拉齐,Sturm公司
例子
对于n=7,系数为0,1的7次多项式的最大独立实根数为3;例如,多项式x^7+x^4+x^2+x有3个不同的实根。
MAPLE公司
f: =proc(n)局部i,j,L,p,v,m;
m: =0:
对于i从2^n到2^(n+1)-1 do
五十: =换算(i,基数,2);
p: =相加(L[j]*x^(j-1),j=1..n+1);
v: =sturm(sturmseq(p,x),x,-无穷大,无穷大);
如果v>m,则m:=v fi;
od;
结束进程:
地图(f,[1..19])#罗伯特·伊斯雷尔,2023年5月7日
数学
对于[n=1,n<=8,n++,打印[Max[删除重复项时的长度@n解算[Total[x^#]+x^n==0,x,Reals]&/@子集[Range[0,n-1]]]]
黄体脂酮素
(MATLAB)
%使用Sturm工具箱;查看链接
对于i=2:13
显示([i-1最大根(i)]);
结束
函数max_len=maxroots(n)
max_len=0;
组合=dec2bin(1:2:2^n-1)-“0”;
对于i=1:2^(n-1)
c=组合(i,:);
如果总和(c,2)==1
继续;
结束
len=distinct_real_roots(c);
如果len>max_len
max_len=长度;
结束
结束
结束
函数sign_var=distinct_real_roots(a)
S=Sturm();
S.build(Poly(a));
符号变量=0;
last_sign=0;
last_sign_parity=0;
奇偶校验=1;
对于i=1:长度(S.m_strum)
v=S.m_sturm{i}.m_coeffs(结束);
如果v>0
如果last_sign<0
sign_var=签名变量-1;
结束
如果last_sign_parity==-奇偶校验
签名变量=签名变量+1;
结束
last_sign=1;
last_sign_parity=奇偶校验;
elseif v<0
如果last_sign>0
sign_var=签名变量-1;
结束
如果last_sign_parity==奇偶校验
sign_var=sign_var+1;
结束
last_sign=-1;
last_sign_parity=-奇偶校验;
结束
奇偶校验=-奇偶校验;
结束
结束
(Python)
来自itertools导入产品
从sympy.abc导入x
来自sympy import sturm,oo,sign,nan,LT
定义A362344飞机(n) :
c=0
对于产品([0,1],repeat=n)中的s:
q=sturm(x**n+总和(int(s[i])*x**i表示范围(n)中的i)
a=[1如果(k:=LT(p).subs(x,-oo))==q中p的nan其他符号(k)[:-1]]+[符号(q[-1])]
b=[1如果(k:=LT(p).subs(x,oo))==q中p的其他符号(k)[:-1]]+[符号(q[-1])]
c=如果a[i]=a[i+1])-sum(1表示范围(len(b)-1)中的i,如果b[i]=b【i+1】)
返回c#柴华武2024年2月15日
(PARI)a(n)=我的(nb=0,k);对于(i=2^n,2^(n+1)-1,k=#集(polrootsreal(Pol(binary(i))));如果(k>nb,nb=k);nb\\米歇尔·马库斯2024年2月15日
交叉参考
的可能副本A090735号A090736号.
关键词
非n,更多
作者
李可阳2023年5月5日
扩展
a(14)根据修正罗伯特·伊斯雷尔,2023年5月7日
a(20)-a(22)来自米歇尔·马库斯2024年2月15日
状态
经核准的

查找|欢迎光临|维基|注册|音乐|地块2|演示|索引|浏览|更多|网络摄像头
贡献新序列。或评论|格式|样式表|转换|超级搜索|最近
OEIS社区|维护人OEIS基金会。

许可协议、使用条款、隐私政策。.

上次修改时间:美国东部夏令时2024年6月3日20:36。包含373088个序列。(在oeis4上运行。)