A317756 Number of distinct primes obtained by cyclically shifting the decimal digits of the n-th prime.
1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1, 2, 2, 2, 2, 1, 3, 2, 3, 1, 1, 2, 1, 2, 2, 1, 2, 1, 2, 2, 1, 3, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 3, 2, 2, 2, 3, 1, 1, 1, 2, 2, 3, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2
Offset: 1
Links
- Robert Israel, Table of n, a(n) for n = 1..10000
Programs
-
Mathematica
f[n_] := Block[{len = IntegerLength@n, s = {n}}, Do[AppendTo[s, FromDigits@RotateRight@IntegerDigits@s[[k - 1]]], {k, 2, len}]; DeleteDuplicates@Select[s, PrimeQ]] (* after Michael De Vlieger in A262988 *); Array[Length@f@Prime@# &, 105] (* Robert G. Wilson v, Aug 06 2018 *) Table[Count[Union[FromDigits/@Table[RotateRight[IntegerDigits[p],n],{n,IntegerLength[p]}]],?PrimeQ],{p,Prime[Range[120]]}] (* _Harvey P. Dale, Jan 18 2025 *)
-
PARI
eva(n) = subst(Pol(n), x, 10) rot(n) = if(#Str(n)==1, v=vector(1), v=vector(#n-1)); for(i=2, #n, v[i-1]=n[i]); u=vector(#n); for(i=1, #n, u[i]=n[i]); v=concat(v, u[1]); v count_primes(n) = my(d=digits(n), i=0); while(1, if(ispseudoprime(eva(d)), i++); d=rot(d); if(d==digits(n), return(i))) a(n) = my(p=prime(n)); count_primes(p) \\ Felix Fröhlich, Aug 06 2018
Comments