login
"Flag numbers": number of dots that can be arranged in successive rows of K, K-1, K, K-1, K, ..., K-1, K (assuming there is a total of L > 1 rows of size K > 1).
11

%I #50 Oct 23 2024 00:44:00

%S 5,8,11,13,14,17,18,20,23,25,26,28,29,32,33,35,38,39,41,43,44,46,47,

%T 48,50,53,56,58,59,60,61,62,63,65,67,68,71,72,73,74,77,78,80,81,83,85,

%U 86,88,89,92,93,94,95,98,101,102,103,104,105,107,108,109,110,111,113,116

%N "Flag numbers": number of dots that can be arranged in successive rows of K, K-1, K, K-1, K, ..., K-1, K (assuming there is a total of L > 1 rows of size K > 1).

%C Numbers of the form F(K, L) = KL+(K-1)(L-1), K, L > 1, i.e. 2KL - (K+L) + 1, sorted and duplicates removed.

%C If K=1, L=1 were allowed, this would contain all positive integers.

%C Positive numbers > 1 but not of the form (odd primes plus one)/2. - Douglas Winston (douglas.winston(AT)srupc.com), Sep 11 2003

%C In other words, numbers n such that 2n-1, or equally, A064216(n) is a composite number. - _Antti Karttunen_, Apr 17 2015

%C Note: the following comment was originally applied in error to the numerically similar A246371. - _Allan C. Wechsler_, Aug 01 2022

%C From _Matthijs Coster_, Dec 22 2014: (Start)

%C Also area of (over 45 degree) rotated rectangles with sides > 1. The area of such rectangles is 2ab - a - b + 1 = 1/2((2a-1)(2b-1)+1).

%C Example: Here a = 3 and b = 5. The area = 23.

%C *

%C ***

%C *****

%C *****

%C *****

%C ***

%C *

%C (End)

%C The smallest integer > k/2 and coprime to k, where k is the n-th odd composite number. - _Mike Jones_, Jul 22 2024

%C Numbers k such that A193773(k-1) > 1. - _Allan C. Wechsler_, Oct 22 2024

%H Antti Karttunen, <a href="/A053726/b053726.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = A008508(n) + n + 1.

%F From _Antti Karttunen_, Apr 17 2015: (Start)

%F a(n) = n + A000720(A071904(n)). [The above formula reduces to this. A000720(k) gives number of primes <= k, and A071904 gives the n-th odd composite number.]

%F a(n) = A104275(n+1). (End)

%F a(n) = A116922(A071904(n)). - _Mike Jones_, Jul 22 2024

%F a(n) = A047845(n+1)+1. - _Amiram Eldar_, Jul 30 2024

%o (Scheme) ;; with Antti Karttunen's IntSeq-library.

%o (define A053726 (MATCHING-POS 1 1 (lambda (n) (and (> n 1) (not (prime? (+ n n -1)))))))

%o ;; _Antti Karttunen_, Apr 17 2015

%o (Scheme) ;; with Antti Karttunen's IntSeq-library.

%o (define (A053726 n) (+ n (A000720 (A071904 n))))

%o ;; _Antti Karttunen_, Apr 17 2015

%o (Python)

%o from sympy import isprime

%o def ok(n): return n > 1 and not isprime(2*n-1)

%o print(list(filter(ok, range(1, 117)))) # _Michael S. Branicky_, May 08 2021

%o (Python)

%o from sympy import primepi

%o def A053726(n):

%o if n == 1: return 5

%o m, k = n, (r:=primepi(n)) + n + (n>>1)

%o while m != k:

%o m, k = k, (r:=primepi(k)) + n + (k>>1)

%o return r+n # _Chai Wah Wu_, Aug 02 2024

%o (PARI) select( {is_A053726(n)=n>4 && !isprime(n*2-1)}, [1..115]) \\ _M. F. Hasler_, Aug 02 2022

%Y Essentially same as A104275, but without the initial one.

%Y A144650 sorted into ascending order, with duplicates removes.

%Y Cf. A006254 (complement, apart from 1, which is in neither sequence).

%Y Cf. also A000720, A008508, A064216, A071904, A199593, A250474.

%Y Differs from its subsequence A246371 for the first time at a(8) = 20, which is missing from A246371.

%K nonn,easy

%O 1,1

%A Dan Asimov, asimovd(AT)aol.com, Apr 09 2003

%E More terms from Douglas Winston (douglas.winston(AT)srupc.com), Sep 11 2003