%I#15 2019年4月9日09:57:14

%S 9,1,3,5,2,4,4,3,7,9,10,11,5,19,22,26,8,17,16,19,9,8,13,7,17,4,17,3，

%电话11,18,13,5,23,17,18,7,17,15,9,18,16,17,9,7,12,28,6,23,9,24,23,13,18，

%U 11,7,14,4,18,14,13,19,11,25,17,6,8,14,27,11,26,8

%N最小的数字k>0，使得k^k的十进制展开式包含N。

%H Giovanni Resta，n的表格，n=0..10000的a（n）</a>

%e a（0）=9，因为9^9=387420489的子串为“0”；

%e a（1）=1，因为1^1=1有“1”作为子字符串；

%e a（2）=3，因为3^3=27有“2”作为子串；

%e a（3）=5，因为5^5=3125有“3”作为子串；

%e a（4）=2，因为2^2=4有“4”作为子字符串。

%t a[n_]：=（k=1；当[！MatchQ[IntegerDigits[k^k]，{___，Sequence@@IntegerDigits[n]，___}]，k++；k）；表[a[n]，{n，0，80}]（*Jean-Francois Alcover的程序适用于此序列-见A030001*）

%t snk[n_]：=模块[{k=1}，而[SequenceCount[IntegerDigits[k^k]，IntegerBigits[n]]<1，k++]；k] ；数组[snk，80,0]（*需要Mathematica版本10或更高版本*）（*Harvey P.Dale_，2019年4月9日*）

%o（PARI）重叠（长，短）=我的（D=10^#位数（短））；while（长>=短，如果（长%D==短，返回（1））；长=10）；0

%o a（n）=我的（k）；while（！重叠（k++^k，n））；2014年3月11日，科勒斯·R·格里特豪斯四世

%Y参考A030001。

%K nonn，基础

%0、1

%2014年3月10日，拉格瑙市