%I#50 2023年10月13日19:29:07

%S 1,2,4,7,10,15,19,26,31,37,43,54,60,73,81,90,98115124143153165，

%电话：177200210226240255268297309340355373391411424461481，

%电话：50251855957661963966068473174877980182885190492695797910091039

%N不同产品i*j的数量，0<=i，j<=N。

%C a（n）=A027420（n，0）=A07420（n，n）_Reinhard Zumkeller_，2014年5月2日

%H Seiichi Manyama，n表，n（n）表示n=0..10000（术语0..1000来自T.D.Noe）

%H R.P.Brent、C.Pomerance、D.Purdum和J.Webster，<a href=“https://arxiv.org/abs/1908.04251“>乘法表的算法</a>，arXiv:1908.04251[math.NT]，2019-2021。

%F对于素数p，a（p）=a（p-1）+p.-达维德·a·科内斯，2019年1月1日

%p A027384:=程序（n）

%p局部L，i，j；

%p L：={}；

%p代表i从0到n do

%从i到ndo的j的p

%p L：=L联合{i*j}；

%p端do：

%p端do：

%p nops（L）；

%p结束程序：#_R.J.Mathar_，2016年5月6日

%t u={}；表[u=并集[u，n*范围[0，n]]；长度[u]，{n，0，100}]（*_T.D.Noe_，2012年1月7日*）

%o（哈斯克尔）

%o导入数据。列表（nub）

%o a027384 n=长度$nub[i*j|i<-[0..n]，j<-[0.n]]

%o——Reinhard Zumkeller，2012年1月1日

%o（PARI）a（n）={my（s=Set（））；对于（i=0，n，s=setunion（s，Set（vector（n+1，k，i*（k-1）））；）；#s；}\\_Michel Marcus_，2019年1月1日

%o（Python）

%o def A027384（n）：返回len（{i*j代表范围（1，n+1）中的i，代表范围（1,i+1）}中的j）+1#_Chai Wah Wu_，2023年10月13日

%Y等于A027424+1，n>0。

%Y参见A027417、A027427和A027425。

%K nonn，简单

%0、2

%弗雷德·施瓦布（Fred Schwab（AT）nrao.edu）