A045765 k - d(k) never takes these values, where d(k) = A000005(k).
7, 8, 13, 19, 20, 24, 25, 28, 33, 36, 37, 40, 43, 49, 50, 52, 55, 56, 63, 64, 66, 67, 68, 74, 75, 79, 85, 88, 98, 100, 103, 108, 109, 112, 113, 116, 117, 123, 124, 126, 131, 132, 133, 134, 136, 140, 143, 145, 150, 153, 156, 159, 160, 163, 164, 167, 168
Offset: 1
Keywords
Links
- Antti Karttunen, Table of n, a(n) for n = 1..10000
- Florian Luca, On numbers not of the form n-omega(n), Acta Mathematica Hungarica, Vol. 106, No. 1 (2005), pp. 117-135.
Crossrefs
Programs
-
Maple
N:= 1000: # to get all terms <= N sort(convert({$1..N} minus {seq(x - numtheory:-tau(x), x=1..2*(1+N))},list)); # Robert Israel, Oct 12 2015
-
Mathematica
lim = 10000; Take[Complement[Range@ lim, Sort@ DeleteDuplicates@ Table[n - DivisorSigma[0, n], {n, lim}]], 57] (* Michael De Vlieger, Oct 13 2015 *)
-
PARI
allocatemem((2^31)+(2^30)); uplim = 36756720 + 640; \\ = A002182(53) + A002183(53). v060990 = vector(uplim); for(n=3, uplim, v060990[n-numdiv(n)]++); A060990 = n -> if(!n,2,v060990[n]); uplim2 = 36756720; n=0; k=1; while(n <= uplim2, if(0==A060990(n), write("b045765_big.txt", k, " ", n); k++); n++;); \\ Antti Karttunen, Oct 09 2015
-
Scheme
(define A045765 (ZERO-POS 1 1 A060990)) ;; Using also IntSeq-library of Antti Karttunen, Oct 06 2015
Comments