|
%I #14 Sep 12 2023 12:29:49
%S 0,1,1,2,1,2,1,4,2,2,1,4,1,2,3,8,1,6,1,4,2,2,1,8,2,2,5,4,1,6,1,6,2,2,
%T 3,12,1,2,3,8,1,6,1,4,2,2,1,16,3,10,3,4,1,18,3,6,2,2,1,8,1,2,6,18,3,6,
%U 1,4,3,4,1,14,1,2,9,4,5,6,1,16,2,2,1,12,2,2
%N Least k such that floor(n/k) + (n mod k) is a prime, or 0 if no such k exists.
%C a(n) = 1 only if n is a prime.
%C a(2m) <= m, because with k=m, floor(2m/m)+(2m mod m) = 2.
%C a(2m+1) <= 2m: floor((2m+1)/2m) + ((2m+1) mod 2m) = 1 + 1 = 2.
%H Antti Karttunen, <a href="/A234541/b234541.txt">Table of n, a(n) for n = 1..10000</a>
%o (Python)
%o primes = [2, 3]
%o primFlg = [0]*100000
%o primFlg[2] = primFlg[3] = 1
%o def appPrime(k):
%o for p in primes:
%o if k%p==0: return
%o if p*p > k: break
%o primes.append(k)
%o primFlg[k] = 1
%o for n in range(5, 100000, 6):
%o appPrime(n)
%o appPrime(n+2)
%o for n in range(1, 100000):
%o a = 0
%o for k in range(1, n):
%o c = n//k + n%k
%o if primFlg[c]: # if c in primes:
%o a = k
%o break
%o print(str(a), end=', ')
%o (Scheme)
%o ;; MIT/GNU Scheme, with Aubrey Jaffer's SLIB Scheme library and function A234575bi as defined in A234575
%o (require 'factor) ;; For predicate prime? from SLIB-library.
%o (define (A234541 n) (let loop ((k 1)) (cond ((prime? (A234575bi n k)) k) ((> k n) 0) (else (loop (+ 1 k))))))
%o ;; _Antti Karttunen_, Dec 29 2013
%Y Cf. A000040, A234575.
%K nonn,easy
%O 1,4
%A _Alex Ratushnyak_, Dec 27 2013
|
