（哈斯克尔）
a051451 n=a051451列表！！（n-1）
a051451_list=扫描1 lcm a000961_list
（PARI）do（lim）=my（v=素数（primepi（lin）），u=列表（[1]））；forprime（p=2，sqrtint（lim\1），for（e=2，log（lim+.5）\log（p），listput（u，p^e））；v=vecsort（concat（v，Vec（u）））；对于（i=2，#v，v[i]=lcm（v[i]，v[i-1]））；v（v）\\查尔斯·格里特豪斯四世2012年11月20日
（PARI）{lim=100；n=1；i=1；j=1；直到（n==lim，直到（a！=j，a=lcm（j，i+1）；i++；）；j=a；n++；打印（n“”a）；）；}\\迈克·温克勒2013年9月7日
（PARI）x=1；对于（i=1100，如果（ω（i）==1，x*=系数（i）[1,1]））\\弗洛里安·鲍尔2022年4月11日
（鼠尾草）
a=[]
L=[1]
对于（1..n）中的i：
a.附录（i）
如果（is_prime_power（i）==1）：
追加（lcm（a））
返回（L）
（Python）
从数学导入prod
从sympy导入primepi，integer_nthroot，integer日志，primerange
def f（x）：对于范围（1，x.bit_length（））中的k，返回int（n+x-1-sum（primepi（整数_nthroot（x，k）[0]））
m、 k=n，f（n）
而m！=克：
m、 k=k，f（k）
返回prod（素数范围（m+1）中p的p**integer_log（m，p）[0]）#柴华武2024年8月15日
