%I#9 2019年8月22日23:11:14
%S 1,1,1,5,11,33,57157303683135829745932125602532852400106256,
%电话:217875441278899551822703370140174911731517825330691135,
%电话:620858461254356892534141432651547323103242763520825519314199956099846685251706477665
%N正整数求和到N的非周期数组的数目。
%C一维外壳为A000740。
%如果n×k矩阵的行序列和列序列的所有n×k旋转都是不同的,则该矩阵是非周期的。
%H Andrew Howroyd,<a href=“/A323867/b323867.txt”>n的表,a(n)表示n=0..200</a>
%e a(5)=33个阵列:
%电子邮箱5 14 23 32 41 113 122 131 212 221 311 1112 1121 1211 2111
%e、。
%e 1 2 3 4 11 12 21
%e 4 3 2 1 12 21 11 11
%e、。
%e 1 1 1 2 2 3
%e 1 2 3 1 2 1
%e 3 2 1 2 1 1
%e、。
%e 1 1 1 2
%e 1 1 2 1
%e 1 2 1 1
%e 2 1 1 1
%t素数MS[n_]:=如果[n==1,{},扁平[Cases[FactorInteger[n],{p_,k_}:>表[PrimePi[p],{k}]]];
%t facs[n_]:=如果[n<=1,{{}},连接@@表[Map[Prepend[#,d]&,选择[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}];
%t ptnmats[n_]:=并集@@置换/@Select[Union@@(Tuples[Permutations/@#]&/@Map[primeMS,facs[n],{2}]),SameQ@@Length/@#&];
%t apematQ[m_]:=不命名Q@@Join@@Table[RotateLeft[m,{i,j}],{i、Length[m]},{j、Length[First[m]]}];
%t表[Length[Union@@Table[Select[ptnmats[k],apromatQ],{k,Times@@Prime/@#&/@IntegerPartitions[n]}],{n,15}]
%o(GAP)列表([0..30],A323867);#代码见A323861_安德鲁·霍罗伊,2019年8月21日
%Y参见A000670、A000740、A027375、A101509。
%Y参见A323860、A323861、A32386、A323963、A3238.64、A323886、A323669。
%K nonn公司
%0、4
%A _Gus Wiseman_,2019年2月4日
%E 2019年8月21日_Andrew Howroyd_的第a(16)条及以后条款
|