A213635 -id:A213635 - cp's OEIS frontend

A213636 Remainder when n is divided by its least nondivisor.

1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 1, 4, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 1, 3, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 1, 4, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 1, 4, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 2, 1

Offset: 1

Clark Kimberling, Jun 16 2012

nonn

Experimentation suggests that every positive integer occurs in this sequence and that
2 occurs only in even numbered positions,
3 occurs in only in positions that are multiples of 12,
4 occurs only in positions that are multiples of 12,
5 occurs only in positions that are multiples of 60,
6 occurs only in positions that are multiples of 60,
7 occurs only in positions that are multiples of 2520, etc.
See A213637 for positions of 1 and A213638 for positions of 2.
From Robert Israel, Jul 28 2017: (Start)
Given any positive number m, let q be a prime > m and r = A003418(q-1).  Then a(n) = m if n == m (mod q) and n == 0 (mod r).  By the Chinese Remainder Theorem, such n exists.
On the other hand, if a(n) = m, we must have A007978(n) > m, and then n must be divisible by A003418(q-1) where q = A007978(n) is a member of A000961 greater than m.  Moreover, if q=p^j with j>1, n is divisible by p^(j-1) so m must be divisible by p^(j-1).  Thus:
For m=2, A003418(2)=2.
For m=3, A007978(n) can't be 4 because m is odd, so A007978(n)>= 5 and n must be divisible by A003418(4)=12.
For m=4, A003418(4)=12.
For m=5 or 6, A003418(6)=60.
For m=7, A007978(n) can't be 8 because m is odd, and can't be 9 because m is not divisible by 3, so n must be divisible by A003418(10)=2520. (End)

			a(10) = 10-3*[10/3] = 1.

Antti Karttunen, Table of n, a(n) for n = 1..10000

Cf. A000961, A003418, A007978, A213633, A213635, A213637, A213638.

f:= proc(n) local k;
  for k from 2 do if n mod k <> 0 then return n mod k fi od
end proc:
map(f, [$1..100]); # Robert Israel, Jul 27 2017

y=120; z=2000;
t = Table[k := 1; While[Mod[n, k] == 0, k++];
   k, {n, 1, z}]  (*A007978*)
Table[Floor[n/t[[n]]], {n, 1, y}] (*A213633*)
Table[n - Floor[n/t[[n]]], {n, 1, y}] (*A213634*)
Table[t[[n]]*Floor[n/t[[n]]], {n, 1, y}] (*A213635*)
t1 = Table[n - t[[n]]*Floor[n/t[[n]]],
   {n, 1, z}] (* A213636 *)
Flatten[Position[t1, 1]] (* A213637 *)
Flatten[Position[t1, 2]] (* A213638 *)
rem[n_]:=Module[{lnd=First[Complement[Range[n],Divisors[n]]]},Mod[n,lnd]]; Join[{1,2},Array[rem,100,3]] (* Harvey P. Dale, Mar 26 2013 *)
Table[Mod[n, SelectFirst[Range[n + 1], ! Divisible[n, #] &]], {n, 105}] (* Michael De Vlieger, Jul 29 2017 *)

def a(n):
    k=2
    while n%k==0: k+=1
    return n%k
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jul 28 2017

def A213636(n): return next(filter(None, (n%d for d in range(2,n)))) if n>2 else n # Chai Wah Wu, Feb 22 2023

(define (A213636 n) (modulo n (A007978 n))) ;; Antti Karttunen, Jul 27 2017

a(n) = n - A213635(n).

a(n) = n - m(n)*floor(n/m(n)), where m(n) = A007978(n).

A213633 [A000027/A007978], where [ ] = floor.

Original entry on oeis.org

0, 0, 1, 1, 2, 1, 3, 2, 4, 3, 5, 2, 6, 4, 7, 5, 8, 4, 9, 6, 10, 7, 11, 4, 12, 8, 13, 9, 14, 7, 15, 10, 16, 11, 17, 7, 18, 12, 19, 13, 20, 10, 21, 14, 22, 15, 23, 9, 24, 16, 25, 17, 26, 13, 27, 18, 28, 19, 29, 8, 30, 20, 31, 21, 32, 16, 33, 22, 34, 23, 35, 14, 36, 24, 37

Offset: 1

Clark Kimberling, Jun 16 2012

a(n) = floor(n/m), where m is the least positive nondivisor of n, as in A007978.

			a(10) = floor(10/3) = 3.

Peter Kagey, Table of n, a(n) for n = 1..10000

Cf. A007978, A213634, A213635, A213636, A213637, A213638.

y=120; z=2000;
t = Table[k := 1; While[Mod[n, k] == 0, k++]; k, {n, 1, z}]  (* A007978 *)
Table[Floor[n/t[[n]]], {n, 1, y}] (* A213633 *)
Table[n - Floor[n/t[[n]]], {n, 1, y}] (* A213634 *)
Table[t[[n]]*Floor[n/t[[n]]], {n, 1, y}] (* A213635 *)
t1 = Table[n - t[[n]]*Floor[n/t[[n]]], {n, 1, z}] (* A213636 *)
Flatten[Position[t1, 1]] (* A213637 *)
Flatten[Position[t1, 2]] (* A213638 *)

cp's OEIS Frontend

A213636 Remainder when n is divided by its least nondivisor.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Python