%I #27 Sep 22 2023 08:19:13
%S 2,3,9,10,11,12,13,14,15,21,27,33,39,45,51,52,53,54,55,56,57,58,59,60,
%T 61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,
%U 84,85,86,87,93,99,105,111,117,123,129,135,141,147,153,159,165
%N Unique monotonic sequence of positive integers satisfying a(a(n)) = k*(n-1) + 3, where k = 6.
%C Numbers m such that the base-6 representation of (5*m-3) starts with 11 or 12 or 13 or 14 or 15 or ends with 0.
%C First differences give a run of 6^i 1's followed by a run of 6^i 6's, for i >= 0.
%H Yifan Xie, <a href="/A353653/b353653.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/Aa#aan">Index entries for sequences of the a(a(n)) = 2n family</a>
%F For n in the range (2*6^i + 3)/5 <= n < (7*6^i + 3)/5, for i >= 0:
%F a(n) = n + 6^i.
%F a(n+1) = 1 + a(n).
%F Otherwise, for n in the range (7*6^i + 3)/5 <= n < (12*6^i + 3)/5, for i >= 0:
%F a(n) = 6*(n - 6^i) - 3.
%F a(n+1) = 6 + a(n).
%e a(5) = 11 because (2*6^1 + 3)/3 <= 5 < (7*6^1 + 3)/5, hence a(5) = 5 + 6^1 = 11;
%e a(10) = 21 because (7*6^1 + 3)/5 <= 10 < (12*6^1 + 3)/5, hence a(10) = 6*(10 - 6^1) - 3 = 21.
%t okQ[m_] := With[{id = IntegerDigits[5 m - 3, 6] }, MatchQ[id[[1 ;; 2]], {1, 1}|{1, 2}|{1, 3}|{1, 4}|{1, 5}] || id[[-1]] == 0];
%t Join[{2}, Select[Range[3, 1000], okQ]] (* _Jean-François Alcover_, Sep 22 2023 *)
%o (C++)
%o /* program used to generate the b-file */
%o #include<iostream>
%o using namespace std;
%o int main(){
%o int cnt1=1, flag=0, cnt2=1, a=2;
%o for(int n=1; n<=10000; n++) {
%o cout<<n<<" "<<a<<endl;
%o if(cnt2==cnt1) {
%o flag=1-flag;
%o cnt1=1;
%o if(flag) a+=1;
%o else {
%o a+=6;
%o cnt2*=6;
%o }
%o }
%o else {
%o cnt1++;
%o a+=(flag?6:1);
%o }
%o }
%o return 0;
%o }
%Y For other values of k: A080637 (k=2), A003605 (k=3), A353651 (k=4), A353652 (k=5), this sequence (k=6).
%K nonn,easy
%O 1,1
%A _Yifan Xie_, Jul 15 2022