%我#47 2019年8月5日10:47:21
%S 0,1,1,2,1,1,1,1,2,3,4,2,1,2,2,3,5,6,1,1,1,1,1,1,2,1,2,3,1,2,2,
%温度3,4,5,6,7,8,9,10,2,1,2,1,1,1,2,3,5,6,7,8,10,11,2,5,6,12,3,7,14,
%U 13、4、11、2、1、8
%N按行读取的不规则三角形:T(N,k)=A109395(N)*k-A076512(N)*A038566(N,k)。
%C第n行的长度为A000010(n)。
%C行n>1的总和=n*A076512(n)/2。
%C第(n)行的第一个值=A076511(n)。
%C第(n)行上的最后一个值=A076512(n),对于n>1。
%C对于n>1,A109395(n)=最大(行)+最小(行)。
%C对于行n>1的位置a和b处的x和y值:
%C x+y=A109395(n),其中a=A000010(n)-(b-1)。
%C对于n>2,第A002110(n)行上的倒数第二个值由下式给出
%C A038110(n)*A000040(n)-A060753(n)。
%C来自Charlie Neder_,2019年6月5日:(开始)
%如果p是素数除以n,那么p*n行由n行的p个副本组成。
%C猜想:如果n是奇数,那么可以通过交换第一和第二部分从第n行得到第2n行。(结束)
%H Jamie Morken,n表,a(n)表示n=1..13413
%e序列为不规则三角形:
%不适用于1、2、3、4。。。
%e 1:0
%电子2:1
%e 3:1、2
%e 4:1、1
%e 5:1、2、3、4
%e 6:2、1
%e 7:1、2、3、4、5、6
%e 8:1、1、1和1
%e 9:1、2、1、2和1、2
%e 10:3、4、1、2
%e 11:1、2、3、4、5、6、7、8、9、10
%e 12:2、1、2、1
%e 13:1、2、3、4、5、6、7、8、9、10、11、12
%e 14:4、5、6、1、2、3
%e 15:7、14、13、4、11、2、1、8
%e。。。
%e行和:0、1、3、2、10、3、21、4、9、10、55、6、78、21、60。
%e T(14,5)=A109395(14)*5-A076512(14)*A038566(14.5)=7*5-3*11=2。
%e T(210,2)=A109395(210)*2-A076512(210)*A038566(210.2)=35*2-8*11=-18。
%t压扁@Table[与[{a=n/GCD[n,#],b=Numerator[#/n]},映射索引的[a First@2-b#1&,压扁@Position[GCD[Table[Mod[k,n],{k,n-1}],n],1]/。{}->{1}]]&@EulerPhi@n,{n,15}](*Michael De Vlieger_,2019年6月6日*)
%o(PARI)vtot(n)=选择(x->(gcd(n,x)==1),向量(n,k,k));
%o行(n)=my(q=eulerphi(n)/n,v=vtot(n));向量(#v,k,分母(q)*k-分子(q)*v[k]);\\_米歇尔·马库斯,2019年5月14日
%Y参见A000010、A109395、A038566、A076511、A076512。
%K符号,tabf
%O 1,4型
%A _杰米·莫肯(Jamie Morken),2019年5月13日
|