cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A272353 Numbers n such that the number of divisors of n+1 divides n and the number of divisors of n divides n+1.

Original entry on oeis.org

3, 15, 1023, 6399, 10815, 15375, 26895, 53823, 55695, 65535, 80655, 107583, 118335, 262143, 309135, 440895, 614655, 633615, 817215, 891135, 1236543, 1784895, 2676495, 2715903, 2849343, 2985983, 3182655, 3225615, 3268863, 4194303, 4326399, 4343055, 4596735, 5053503
Offset: 1

Views

Author

Waldemar Puszkarz, Apr 26 2016

Keywords

Comments

One may call such pairs {n, n+1} mutually (or amicably) one step refactorable numbers.
While most terms appear to be divisible by 3, some are not, the first being 2985983=7*11*13*19*157.
From Robert Israel, May 09 2016: (Start)
2^k-1 is a term if k+1 is an odd prime and 2^k-1 is squarefree.
If n is an odd term, then n+1 is a square, and n == 3 mod 4. (End)
There is no term such that last digit of it 1 or 7. Proof: If n is an odd term, then n+1 is a square. For any odd number k, last digit can be trivially 1, 3, 5, 7, 9, that is, the last digit of k+1 is 2, 4, 6, 8, 0 for corresponding odd k values. There cannot be square such that last digit of it 2 or 8. So in this sequence, there is no term such that the last digit of it 1 or 7. - Altug Alkan, May 11 2016

Examples

			15 is a term because the number of divisors of 16=15+1, which is 5, divides 15, and the number of divisors of 15, which is 4, divides 16.
		

Crossrefs

Cf. A000005 (number of divisors), A033950 (refactorable numbers), A268037, A269781 (related sequences).

Programs

  • Maple
    select(t -> (t+1) mod numtheory:-tau(t) = 0 and t mod numtheory:-tau(t+1) = 0, [$1..10^6]); # Robert Israel, May 09 2016
  • Mathematica
    lst={}; Do[ If[ Divisible[n, DivisorSigma[0, n+1]]&&Divisible[n+1, DivisorSigma[0, n]], AppendTo[lst, n]], {n, 7000000}]; lst
    Select[Range[7000000], Divisible[#, DivisorSigma[0, # + 1]] && Divisible[# + 1, DivisorSigma[0, #]] &]
  • PARI
    for(n=1, 7000000, (n%numdiv(n+1)==0) && ((n+1)%numdiv(n)==0)&& print1(n ", "))

A269781 a(n) is the smallest k different from n such that (n, k) is an amicably refactorable pair (see the comments).

Original entry on oeis.org

4, 3, 16, 4, 64, 24, 36, 16, 1024, 6, 4096, 64, 4, 5, 65536, 12, 262144, 6, 4, 1024, 4194304, 8, 81, 4096, 4, 6, 268435456, 16, 1073741824, 6, 4, 65536, 16, 9, 68719476736, 262144, 4, 8, 1099511627776, 32, 4398046511104, 6, 36, 4194304, 70368744177664, 10, 729, 48, 4, 6, 4503599627370496, 32, 16, 8
Offset: 3

Views

Author

Waldemar Puszkarz, May 01 2016

Keywords

Comments

Let m and k be distinct integers and numdiv(n) be the number of divisors of n (A000005(n)). We call m and k amicably refactorable if numdiv(m) divides k and numdiv(k) divides m.
For any n with no amicably refactorable partner, a(n) = 0.
Conjecture: the sequence contains no zeros.
1 does not have an amicable partner as all other numbers have more than one divisor and 2 does not have an amicable partner as all other numbers with two divisors are odd primes and cannot be divided by the number of divisors of 2, also 2. All other numbers may have an amicably refactorable partner, though for some, primes, semiprimes and squares of primes in particular, this number can be quite large.
For primes and semiprimes, a(n) = 2^(f(n) - 1), (see A061286), where f(n) is their largest prime factor. For squares of primes, a(n) = 3^(|sqrt(n)| - 1), except for n = 9 for which this formula yields 9; this forces us to choose the next best candidate: 36.

Examples

			For n=5, a(5)=16 as the number of divisors of n (2) divides a(n) while the number of divisors of a(n) (5) divides 5 and 16 is the smallest number for which this happens.
		

Crossrefs

Cf. A000005 (number of divisors), A033950 (refactorable numbers), A061286 (subsequence for odd prime indices and semiprime indices), A268037, A272353 (related sequences).

Programs

  • Mathematica
    A269781 = {}; Do[k = 1; If[PrimeQ[n] || PrimeNu[n] == 2 && PrimeOmega[n] == 2, AppendTo[A269781, 2^(First[Last[FactorInteger[n]]] - 1)], If[PrimeQ @ Sqrt @ n && (n > 9), AppendTo[A269781, 3^(Sqrt[n] - 1)],While[k != n && !(Divisible[n, DivisorSigma[0, k]] && Divisible[k, DivisorSigma[0, n]]), k++]; If[k == n, k = n + 1; While[!(Divisible[n, DivisorSigma[0, k]] && Divisible[k, DivisorSigma[0, n]]), k++]]; AppendTo[A269781, k]]], {n, 3, 56}]; A269781
Showing 1-2 of 2 results.