N: =7 ^6:
cv7:=proc(n)局部L;五十: =换算(n,基数,7);
添加(L[i]*10^(i-1),i=1..nops(L))结束进程:
V: =矩阵(N,7,数据类型=integer[8]):
计数:=0:结果:=空:
对于i从1到N do
五十: =换算(i,基数,7);
M: =向量[行]([seq(numbeccurse(d,L),d=0..6)],数据类型=整数[8]);
对于从i到N的r,通过i do V[r,..]:=V[r、..]+M od;
如果nops(转换(V[i,..],设置)=1,则
计数:=计数+1;
w: =cv7(i);
Res:=Res,w;
fi(菲涅耳)
日期: