%I#13 2018年6月7日22:01:23
%S 0,0,0,1,0,2,0,00,02,0,2,10,2,0,10,10,2,2,0,4,0,2,0,10,100,2,2,10,0,84,0,
%温度0,2,2,2122,0,2,44,0,84,0,10,10,2,0184,0,10,2,10,0,44,2,2,0,
%U 1590,0,2,10,0,2,84,0,10,2,84,01156,0,1,10,10,2
%N具有最小公共倍数N的相对素数正整数的有限集的数目。
%C自2018年6月6日起(开始)
%Ca(n)仅取决于n的素数签名。
%C如果n在A000961中,a(n)=0。
%C如果n在A006881中,a(n)=2。(完)
%C如果n=p^k*q,其中p和q是不同的素数,k>=1,那么a(n)=3*4^(k-1)-2^(k-1)_罗伯特·伊斯雷尔,2018年6月7日
%H Robert Israel,n的表,n=1..10000的a(n)</a>
%e a(12)=10组:
%e{3,4},
%e{2,3,4},{2,3,12},}3,4,6},
%e{2,3,4,6},{2,3,12},}2,3,6,12},
%e{2,3,4,6,12}。
%p f:=进程(n)g(排序(映射(t->t[2],ifactors(n)[2])))结束进程:
%p f(1):=0:
%p g:=proc(L)选项记忆;
%p本地nL、Cands、nC、Cons、i;
%p nL:=nops(L);
%p坎迪斯:=[[]];
%p代表i从1到nL do
%p Cands:=[seq(seq([op(s),t],t=0..L[i]),s=Cands)];
%日期:
%p Cands:=移除(t->max(t)=0,Cands);
%p nC:=nops(Cands);
%p Cons:=[seq(选择(t->Cands[t][i]=0,{$1..nC}),i=1..nL),
%p序列(选择(t->Cands[t][i]=L[i],{$1..nC}),i=1..nL)];
%p h(Cons,{$1..nC})
%p结束过程:
%p h:=进程(Cons,Cands)
%p本地t,i,Consi,Candsi;
%p如果Cons=[],则返回2^nops(Cands)fi;
%p t:=0;
%p代表i从1到nops(Cons[1])do
%p组成:=映射(proc(t)if成员(Cons[1][i],t),则为NULL,否则t减去Cons[1][1..i-1]fi结束进程,Cons[2..-1]);
%p如果成员({},Consi),则为下一个fi;
%p Candsi:=Cands减去Cons[1][1..i];
%p t:=t+程序名(Consi,Candsi)
%p od;
%价格
%p结束过程:
%p映射(f,[1..100]美元);#_罗伯特·伊斯雷尔,2018年6月7日
%t表[Length[Select[Subsets[Rest[Divisors[n]],And[GCD@@#==1,LCM@@#=n]&]],{n,100}]
%Y参考A000961、A006881、A076078、A181819、A281116、A285572、A290103、A304818、A305563、A305564、A305565、A305567。
%K nonn公司
%O 1,6型
%A _Gus Wiseman_,2018年6月5日
|