A263095 -id:A263095 - cp's OEIS frontend

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

David W. Wilson

nonn

Complement of A236562. - Jaroslav Krizek, Feb 09 2014
Positions of zeros in A060990, leaf-nodes in the tree generated by edge-relation A049820(child) = parent. - Antti Karttunen, Oct 06 2015
Since A000005(x) <= 1 + x/2, k is in the sequence if there are no x <= 2*(k+1) with k = x - d(x). - Robert Israel, Oct 12 2015
This can be improved as: k is in the sequence if there are no x <= k + A002183(2+A261100(k)) with k = x - d(x). Cf. also A070319, A262686. - Antti Karttunen, Oct 12 2015
Luca (2005) proved that this seqeunce is infinite. - Amiram Eldar, Jul 26 2025

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.

Top row of A262898.
Cf. A000005, A002183, A049820, A060990, A070319, A236562, A236565, A261100, A262511, A262686, A262901, A262902, A262903, A262909, A263081.
Cf. A263091 (primes in this sequence), A263095 (squares).
Cf. A259934 (gives the infinite trunk of the same tree, conjectured to be unique).

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

lim = 10000; Take[Complement[Range@ lim, Sort@ DeleteDuplicates@ Table[n - DivisorSigma[0, n], {n, lim}]], 57] (* Michael De Vlieger, Oct 13 2015 *)

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

(define A045765 (ZERO-POS 1 1 A060990))
;; Using also IntSeq-library of Antti Karttunen, Oct 06 2015

A263093 Numbers whose squares are in A045765.

Original entry on oeis.org

5, 6, 7, 8, 10, 14, 16, 18, 20, 22, 26, 27, 28, 34, 35, 37, 46, 47, 50, 54, 56, 58, 59, 60, 62, 67, 73, 78, 82, 85, 89, 90, 94, 95, 98, 100, 103, 104, 106, 110, 114, 116, 118, 122, 124, 125, 126, 127, 128, 130, 135, 140, 141, 142, 148, 150, 155, 158, 161, 164, 170, 172, 174, 177, 178, 182, 184, 188, 190, 199, 202, 205, 207

Offset: 1

Antti Karttunen, Oct 11 2015

nonn

Numbers n such that there is no such k for which k - d(k) = n^2, where d(k) is the number of divisors of k (A000005).

Numbers n for which A060990(n^2) = A263087(n) = 0.

Antti Karttunen, Table of n, a(n) for n = 1..12443
A. Karttunen, Ratio a(n)/A263092(n) drawn with the help of OEIS Plot2-utility

Cf. A000005, A000196, A045765, A049820, A060990, A263098.
Complement: A263092.
Positions of zeros in A263087 and positions of ones in A263088.
Cf. A263095 (the squares of these numbers).

\\ Compute A263093 and A263095 at the same time:
A060990(n) = { my(k = n + 1440, s=0); while(k > n, if(((k-numdiv(k)) == n),s++); k--;); s}; \\ Hard limit 1440 is good for at least up to A002182(67) = 1102701600 as A002183(67) = 1440.
n = 1; k = 0; while((n^2)<1102701600, if((0 == A060990(n*n)), k++; write("b263093.txt", k, " ", n); write("b263095.txt", k, " ", (n*n)); ); n++; if(!(n%8192),print1(n,",k=", k, ", ")); );

;; With Antti Karttunen's IntSeq-library.
(define A263093 (MATCHING-POS 1 1 (lambda (n) (zero? (A060990 (* n n))))))
(define A263093 (ZERO-POS 1 0 A263087))

a(n) = A000196(A263095(n)).

A263094 Squares in A236562; numbers n^2 such that there is at least one such k for which k - d(k) = n^2, where d(k) is the number of divisors of k (A000005).

Original entry on oeis.org

0, 1, 4, 9, 16, 81, 121, 144, 169, 225, 289, 361, 441, 529, 576, 625, 841, 900, 961, 1024, 1089, 1296, 1444, 1521, 1600, 1681, 1764, 1849, 1936, 2025, 2304, 2401, 2601, 2704, 2809, 3025, 3249, 3721, 3969, 4096, 4225, 4356, 4624, 4761, 4900, 5041, 5184, 5476, 5625, 5776, 5929, 6241, 6400, 6561, 6889, 7056, 7396, 7569, 7744, 8281, 8464, 8649, 9216, 9409, 9801, 10201, 10404, 11025

Offset: 0

Antti Karttunen, Oct 11 2015

nonn

Starting offset is zero, because a(0)=0 is a special case in this sequence.

Antti Karttunen, Table of n, a(n) for n = 0..20763

Cf. A000005, A049820, A060990.
Intersection of A000290 and A236562.
Cf. A263092 (gives the square roots of these terms).
Cf. A263095 (complement among squares).
Cf. A262514 (a subsequence).
Cf. also A263090, A263098.

Take[Select[Sort@ DeleteDuplicates@ Table[n - DivisorSigma[0, n], {n, 20000}], IntegerQ@ Sqrt@ # &], 68] (* Michael De Vlieger, Oct 13 2015 *)

PARI
```
\\ See code in A263092.
```

(define (A263094 n) (A000290 (A263092 n)))

a(n) = A000290(A263092(n)).

cp's OEIS Frontend

A045765 k - d(k) never takes these values, where d(k) = A000005(k).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Scheme

A263093 Numbers whose squares are in A045765.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Scheme

Formula

A263094 Squares in A236562; numbers n^2 such that there is at least one such k for which k - d(k) = n^2, where d(k) is the number of divisors of k (A000005).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Scheme

Formula