login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A227502 Numbers n such that Sum_{i=1..n} i^(i') == 0 (mod n), where i' is the arithmetic derivative of i. 8
1, 3, 7, 19, 32, 57, 99, 103, 439, 540, 2656, 18156, 179171, 235056 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

a(14) > 200000. - Giovanni Resta, Jul 15 2013

a(15) > 1000000. - Bert Dobbelaere, Dec 24 2018

LINKS

Table of n, a(n) for n=1..14.

EXAMPLE

1^1' + 2^2' + 3^3' = 1^0 + 2^1 + 3^1 = 6 and 6 == 0 (mod 3).

MAPLE

with(numtheory); ListA227502:=proc(q)  local a, n, p;  a:=0;

for n from 1 to q do a:=a+n^(n*add(op(2, p)/op(1, p), p=ifactors(n)[2]));

if a mod n=0 then print(n);  fi; od; end: ListA227502(10^6);

MATHEMATICA

d[n_] := d[n] = n*Total[Apply[#2/#1 &, FactorInteger[n], {1}]]; Reap[For[n = 1, n <= 2*10^5, n++, If[Mod[Sum[k^d[k], {k, 1, n}], n] == 0, Print[n]; Sow[n]]]][[2, 1]] (* Jean-Fran├žois Alcover, Feb 21 2014 *)

CROSSREFS

Cf. A003415, A227848.

Sequence in context: A145039 A112633 A113916 * A268065 A049490 A112391

Adjacent sequences:  A227499 A227500 A227501 * A227503 A227504 A227505

KEYWORD

nonn,more

AUTHOR

Paolo P. Lava, Jul 13 2013

EXTENSIONS

a(13) from Giovanni Resta, Jul 15 2013

a(14) from Bert Dobbelaere, Dec 24 2018

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 3 22:03 EDT 2021. Contains 346441 sequences. (Running on oeis4.)