%I#62 2023年1月29日17:04:41

%S 0,3,2,1,12,15,14,13,8,11,10,9,4,7,6,5,48,51,50,49,60,63,62,61,56,59，

%电话58,57,52,55,54,53,32,35,34,33,44,47,46,45,40,43,42,41,36,39,38,37,16，

%U 19,18,17,28,31,30,29,24,27,26,25,20,23,221,2195194193204207206年

%N以4为基数写入N，然后将每个数字“1”替换为“3”，反之亦然，并转换回十进制。

%C[1..4^k]上的a（n）图类似于分形维数为1的平面分形。

%C自逆被认为是整数的置换。

%C序列的前4^n项构成了0..4^n-1的置换s（n），n>=1；s（n）的反演数为A115490（n）_Gheorghe Coserea，2018年4月23日

%H Reinhard Zumkeller，n表，n=0..16383的a（n）</a>

%H J.W.Layman，<a href=“http://intranet.math.vt.edu/people/layman/sequences/sequences.htm“>查看分形图形</a>

%H<a href=“/index/Per#IntegerPermutation”>自然数排列序列的索引项</a>

%F a（n）=如果n=0，则0其他4*a（楼层（n/4））+如果m=0，那么0其他4-m，其中m=n mod 4.-_Reinhard Zumkeller，2013年4月8日

%F G.F.G（x）满足：G（x）=4*（1+x+x^2+x^3）*G（x^4）+（3*x+2*x^2+x^3_罗伯特·伊斯雷尔，2014年11月3日

%e a（15）=5，因为15=33_4->11_4=5。

%p f：=进程（n）

%p选项记忆；

%p局部m，r；

%pm:=n模块4；

%p r：=4*进程名（（n-m）/4）；

%p，如果m＝0，则r为r+4m-fi；

%p结束过程：

%p f（0）：=0:

%p序列（f（n），n=0..100）；#_罗伯特·伊斯雷尔，2014年11月3日

%p#第二个Maple程序：

%p a:=proc（n）选项记忆`如果`（n=0，0，

%p a（iquo（n，4，'r'））*4+[0，3，2，1][r+1]）

%p端：

%p序列（a（n），n=0..70）；#_阿洛伊斯·海因茨，2022年1月25日

%t表[FromDigits[If[#==0,0,4-#]和/@IntegerDigits（n，4]，4]），{n，0,70}]（*哈维·P·戴尔，2012年7月23日*）

%o（哈斯克尔）

%o a048647 0=0

%o a048647 n=4*a048647n’+如果m==0，则0其他4-m

%o其中（n'，m）=divMod n 4

%o——Reinhard Zumkeller，2013年4月8日

%o（PARI）a（n）=来自数字（apply（d->if（d，4-d），digits（n，4）），4）\\查尔斯·格里特豪斯IV_，2017年6月23日

%o（Python）

%o来自sympy.theory.factor _导入数字

%o定义a（n）：

%o如果d！则返回int（“”.join（str（4-d））=以数字（n，4）[1:]），4）表示d的其他0“0”

%o打印（[a（n）代表范围（101）内的n）]#_Indranil Ghosh，2017年6月26日

%o（Python）

%o定义A048647（n）：返回n^（（n&（（1<<（m:=n.bit_length（））+（m&1））-1）//3）<<1）#_Chai Wah Wu_，2023年1月29日

%o（C）uint32_ta（uint32-tn）{return n^（（n&0x5555555）<<1）；}//_Falk Hüffner_，2022年1月22日

%Y参考A065256、A007090。

%A248813的Y列k=4。

%K nonn，简单，漂亮，基础，看

%O 0,2

%A _John W.Layman，1999年7月5日