|
|
A064446号 |
| a(n)=gcd(n!,n^n,lcm(1,2,…,n))或gcd(n^n、lcm(1,2,…、n))。 |
|
三
|
|
|
1, 2, 3, 4, 5, 12, 7, 8, 9, 40, 11, 72, 13, 56, 45, 16, 17, 144, 19, 80, 63, 176, 23, 144, 25, 208, 27, 112, 29, 10800, 31, 32, 297, 544, 175, 864, 37, 608, 351, 800, 41, 6048, 43, 352, 675, 736, 47, 864, 49, 800, 459, 416, 53, 864, 275, 1568, 513, 928, 59, 21600, 61
(列表;图表;参考文献;听;历史;文本;内部格式)
|
|
|
抵消
|
1,2
|
|
评论
|
gcd(n^n,lcm(1..n,当考虑gcd(n^n,lcm(1..n))时,我们只需要考虑r,可以忽略s。因为r是每个不同素数除数p的最大幂e_1的乘积,并且由于划分n^n的相应素数的幂e_2必须始终是e_2>=e_1,所以计算r就足以确定a(n)-迈克尔·德弗利格2015年10月26日
|
|
链接
|
|
|
配方奶粉
|
a(n)=产品{prime p|n}p^floor(log_p(n))-罗伯特·伊斯雷尔2015年10月26日
a(n)=e^(总和{k=1..n}(floor(n^k/k)-floor((n^k-1)/k))*Mangoldt(k)),其中Mangoldt是Mangold函数-安东尼布朗2016年6月16日
|
|
例子
|
n=6:a(6)=gcd(720,60,46656)=12。
由于只有1和5是相对于6的质数,a(6)=lcm(1,2,3,4,5,6)/lcm(1,5)=60/5=12。
|
|
MAPLE公司
|
A064446号:=n->ilcm(seq(i,i=1..n))/ilcm(op(select(k->igcd(n,k)=1,[$1..n]));
N: =1000:#将a(1)转换为a(N)
素数:=选择(isprime,[2,seq(2*i+1,i=1..floor((N-1)/2))]):
A: =矢量(N,1):
对于素数do中的p
对于从1到地板的d(log[p](N))do
对于从p^d到min(N,p^(d+1)-p)的j,通过p do
A[j]:=A[j]*p^d
od od od日期:
|
|
数学
|
表[GCD[n!,n^n,LCM@@Range[n]],{n,70}](*哈维·P·戴尔2011年6月25日*)
f[n_]:=块[{p=First/@FactorInteger@n},时间@@Power@@Transpose[{p,Floor@Log[#,n]&/@p}]];{1} ~连接~表[f@n,{n,2,10000}](*迈克尔·德弗利格2015年10月26日*)
|
|
黄体脂酮素
|
(PARI)L=1;对于(n=11000,L=lcm(L,n));写入(“b064446.txt”,n,“”,gcd(n^n,L))\\哈里·史密斯2009年9月14日
(PARI)a(n)=我的(f=系数(n));对于(i=1,#f~,f[i,2]=logint(n,f[i,1]));因子回收(f)\\查尔斯·格里特豪斯四世2015年11月19日
(PARI)a(n)=gcd(n^n,lcm(矢量(n,k,k))\\米歇尔·马库斯2018年3月18日
(GAP)列表([1.70],n->Gcd(因子(n),n^n,Lcm([1.n]))#穆尼鲁·A·阿西鲁2018年3月20日
|
|
交叉参考
|
|
|
关键字
|
非n
|
|
作者
|
|
|
状态
|
经核准的
|
|
|
|