OFFSET
1,1
MAPLE
A076382 := proc(n) local i, j, t, t1, sod, sopd; t := NULL; for i from 2 to n do t1 := i; sod := 0; while t1 <> 0 do sod := sod + (t1 mod 4); t1 := floor(t1/4); od; sopd := 0; j := 1; while ithprime(j) <= i do if i mod ithprime(j) = 0 then sopd := sopd+ithprime(j); fi; j := j+1; od; if sopd mod sod = 0 then t := t, i; fi; od; t; end;
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Floor van Lamoen, Oct 08 2002
STATUS
approved