%I#10 2023年3月10日10:02:02

%电话：48498455870865656155574023957912905827326585804984385，

%电话：94444351001500510140585104654551055555451070569510818885，

%电话：109759651156501511696685119969851226725512777765127859551309654513408395134951851513677451383405

%N正好是7个不同奇素数的乘积。

%H Karl-Heinz Hofmann，n的表格，n=1..10000的a（n）</a>

%e a（1）=4849845=3*5*7*11*13*17*19

%e a（9663）=253808555=5*7*11*13*17*19*157

%e a（9961）=258573315=3*5*7*11*13*17*1013

%e a（10000）=259173915=3*5*7*11*13*41*421

%o（Python）

%o导入数字

%o从sympy导入nextprime、sieve、primepi

%o k_upto=14*10**6

%o数组=numpy.zeros（k_upto，dtype=“i4”）

%o sieve_max_number=素数（下一素数（k_upto//255255））

%o对于范围（2，sieve_max_number）内的s：

%o数组[筛[s]：k_upto][：：筛[s]]+=1

%o对于范围（2，sieve_max_number）内的s：

%o数组[筛[s]**2:k_upto][：：筛[s]**2]=0

%o打印（如果数组[k]==7]，则[k代表范围（1，k_upto，2）中的k）

%Y参考A065091，A046388（2个不同的奇素数）。

%Y参见A046389（3个不同的奇素数），A046390（4个不同的奇数素数）。

%Y参见A046391（5个不同的奇素数），A168352（6个不同的奇数素数）。

%K nonn公司

%O 1,1号机组

%A _Karl-Heinz Hofmann，2023年3月1日