A141063 - OEIS

%I #16 Feb 19 2021 05:55:56

%S 1,2,3,4,1,2,7,8,4,10,3,2,3,0,4,0,3,4,6,4,1,6,1,7,9,0,3,4,9,0,1,2,0,8,

%T 7,1,11,8,11,7,7,2,10,5,11,8,3,6,5,11,3,10,4,6,13,1,6,8,11,5,9,6,3,4,

%U 2,0,4,3,13,6,5,4,9,9,18,6,17,2,4,2,11,5,3,4,5,9,2,8,1,12,6,12,17,10,7,0,13

%N a(n) = n mod (sum of digits of prime(n)).

%e Prime(1) = 2, so a(1) = 1 mod 2 = 1;

%e prime(5) = 11, so a(5) = 5 mod (1+1) = 5 mod 2 = 1.

%p sd:=proc(n) options operator, arrow: add(convert(n,base,10)[j], j=1..nops(convert(n,base,10))) end proc: seq(`mod`(n, sd(ithprime(n))),n=1..100); # _Emeric Deutsch_, Aug 16 2008

%o (PARI) a(n) = n % sumdigits(prime(n)); \\ _Michel Marcus_, Feb 18 2021

%Y Cf. A007605.

%K nonn,base,less

%O 1,2

%A _Juri-Stepan Gerasimov_, Aug 04 2008

%E Corrected and extended by _Emeric Deutsch_, Aug 16 2008