|
|
|
|
#46通过乔格·阿恩特2022年9月13日星期二11:14:38 EDT |
|
|
|
#45通过米歇尔·马库斯2022年9月13日星期二上午10:58:41 |
|
|
|
#44通过柴华武2022年9月13日星期二10:56:52 EDT |
|
|
|
#43通过柴华武2022年9月13日星期二美国东部夏令时10:49:04 |
| 黄体脂酮素
|
(Python)
从itertools导入组合
从sympy.utilities.iterables导入multiset_permutations
从sympy导入isprime
定义A075053号(n) :如果b[0]!='0'和isprime(int(''.join(b)))#柴华武2022年9月13日
|
| 状态
|
经核准的
编辑
|
|
|
|
#42通过N.J.A.斯隆2019年10月22日星期二21:31:17 EDT |
|
|
|
#41通过M.F.哈斯勒2019年10月18日星期五18:06:40 EDT |
|
|
讨论
|
10月19日星期六
| 06:07
| 大卫·A·科内斯:感谢您删除a()。现在这里的所有PARI代码都可以正常工作了。
|
|
|
|
#40通过M.F.哈斯勒2019年10月18日星期五18:05:15 EDT |
| 链接
|
M.F.Hasler,<a href=“https://oeis.org/wiki/Primeval_numbers网站“>原始数字</a>,OEIS wiki,2014年,2019年更新。
|
|
|
|
#39通过M.F.哈斯勒美国东部时间2019年10月18日星期五18:02:39 |
| 公式
|
a(n)=S(n)}中的和{kA039999号(k) ,如果n不超过一个数字0,其中S(n)是其 非零数字是的子序列n个'秒那些 非零属于 数字n个,跟着和 哪一个 通过包含数字0如果n 有 一做,具有多重性(如果 n个 有 重复 非零 数字).如果 n个 有 更多 比 一 数字 0,这 做 不 计数 具有 对的 多重性对于 素数数字使用一些但不是全部 数字 那个 是 重复的 属于在里面 他们_n个). - _M.F.Hasler,2019年10月15日
|
|
|
|
#38通过M.F.哈斯勒2019年10月18日星期五17:53:31 EDT |
| 黄体脂酮素
|
(PARI)a(n,D=vecsort(数字(n)),z=#D-hammingweight(D))={如果(z>1,A075053号(,D),总和(b=2,2^(#D-z),A039999号(,向量提取(D,b<<z-1))}\\M.F.哈斯勒,2019年10月17日
|
|
|
讨论
|
10月18日星期五
| 17:56
| M.F.哈斯勒:嗯,速度并没有显著加快(a()vs A075053()),但使用A0399999会缩短。但我会发布那个项目。代码没有更改,只有注释/公式。
|
|
|
|
#37通过M.F.哈斯勒2019年10月18日星期五17:27:16 EDT |
| 黄体脂酮素
|
(PARI)A075053号(n,D=vecsort(digits(n)),S=0)=对于(b=1,2^#D-1,对于perm(vecextract(D,b),p,p[1]&&isprime(from digities(Vec(p)))&&S++&&打印1(来自数字(Vec公司(第页))", ")));++));S公司\\ _\\第二 可选择的 参数 允许 到 给 这个 数字 如果 他们 是 已经 已知_M.F.Hasler,2019年10月14日,取代2014年的早期代码。
|
|
|
讨论
|
10月18日星期五
| 17:35
| M.F.哈斯勒:(实际上,函数a()的作用是*(我得到a(100)=0,a(101)=2,…)但当n的位数大于1时,使用另一个A075053()。我致力于实现一种更好的方法,即使用组合因子,使用n的一些但不是所有数字来正确计算素数的多重性。
|
|
|
|