|
%I #10 Sep 15 2023 23:09:09
%S 0,0,0,1,-1,0,3,0,0,-3,3,3,12,6,-3,34,-1,0,16,8,9,103,17,0,75,6,0,-17,
%T -7,-9,24,12,9,36,21,9,12,60,36,135,99,18,207,36,-9,199,149,102,576,
%U 150,-3,448,11,0,22,54,48,217,29,24,289,50,27,279,425,309,808,212,51,1180,89,0,1152,318,225,3049,323,18
%N a(n) = A356867(n) - n.
%H Antti Karttunen, <a href="/A365462/b365462.txt">Table of n, a(n) for n = 1..19683</a>
%o (PARI)
%o up_to = 19683;
%o A356867list(up_to) = { my(v=vector(up_to),met=Map(),h=0,ak); for(i=1,#v,if(1==vecsum(digits(i,3)), v[i] = i; h = i, ak = v[i-h]; forprime(p=2,,if(3!=p && !mapisdefined(met,p*ak), v[i] = p*ak; break))); mapput(met,v[i],i)); (v); };
%o v356867 = A356867list(up_to);
%o A356867(n) = v356867[n];
%o A365462(n) = (A356867(n)-n);
%Y Cf. A356867, A364958 (positions of 0's), A365463.
%Y Cf. also A364499.
%K sign
%O 1,7
%A _Antti Karttunen_, Sep 15 2023
|
