login
Search: keyword:new
     Sort: relevance | references | number | modified | created      Format: long | short | data
Number of n X n binary matrices with exactly 10 1's such that there are no odd number of 1's in a row and no odd number of 1's in a column.
+0
0
0, 0, 0, 96, 12840, 239040, 2163840, 12888960, 58042656, 213615360, 674570160, 1888587360, 4797335400, 11245730496, 24642268560, 50981199360, 100370578560, 189249592320, 343524515616, 602905753440, 1026789137640, 1702093525440, 2753544282720, 4356978950016, 6756347796000, 10285188576000
OFFSET
1,4
COMMENTS
Also weight A_10 of the weight distribution for a square n X n Single Parity Check Product (SPCP) code when the 1's represent bit errors, and 0's represent unchanged data bits. For a SPCP code, A_odd = 0, A_4 = C(n,2)^2 = A000537(n-1), A_6 = 3!*C(n,3) = A179058(n), A_8 = A396317(n). See the rectangular formula for other message data formats.
REFERENCES
E. R. Berlekamp, Algebraic Coding Theory, McGraw-Hill, New York, 1968, pp. 397-400.
R. E. Blahut, Theory and Practice of Error Control Codes, Addison-Wesley, 1983, p. 431.
C. Leung, Evaluation of the Error Probability of Single Parity-Check Codes, IEEE Transactions on Communications Vol. COM-31, NO. 2, Feb 1983, pp. 250-253.
F. J. MacWilliams, and N. J. A. Sloane The Theory of Error Correcting Codes, Elsevier/North Holland 1977, p. 21.
W. W. Peterson, and E. J. Weldon, Error Correcting Codes, MIT Press, 1994, pp. 64-70, pp. 131-132.
LINKS
G. Gordon Thompson, Probability of Undetected Error for Some Simple Product Codes, MS Thesis, UT Arlington TX, May 1984, pp. 35-37.
FORMULA
a(n) = 2040*C(n,5)^2 + 1440*C(n,5)*C(n,4) + 120*C(n,5)*C(n,3) + 96*C(n,4)^2.
a(n) = 1/120*(n - 3)*(n - 2)^2 * (n - 1)^2 * n^2 * (17*n^3 - 127*n^2 + 300*n - 236).
The rectangular formula for the n X k matrices is a(n,k) = 2040*C(n,5)*C(k,5) + 720*C(n,5)*C(k,4) + 720*C(n,4)*C(k,5) + 60*C(n,3)*C(k,5) + 60*C(n,5)*C(k,3) + 96*C(n,4)*C(k,4). This is shown in table form below, where each entry is multiplied by its column header and left row header, then all entries are summed for a(n,k). The square case n X n, is when k = n.
A_10 C(k,2) C(k,3) C(k,4) C(k,5)
C(n,2) 0 0 0 0
C(n,3) 0 0 0 60
C(n,4) 0 0 96 720
C(n,5) 0 60 720 2040
EXAMPLE
The first formula given above adds the counts of 4 subsets that are categorized by the same number of 1's in rows and columns (with transpositions in the same category).
For n=4, a(4) = 90*0 + 1440*0 + 120*0 + 96 = 96 by the formula. See also the figures below, showing an example for each subset:
.
Subset 1, 5 X 5 exactly 2 1's per row and 2 1's per column.
1 1 0 0 0
0 1 1 0 0 The number of patterns of 2 1's per row and column
0 0 1 1 0 is 2040, see A001499(5).
0 0 0 1 1
1 0 0 0 1
.
Subset 2, 5 X 4, one row with 4 1's. We count the number of patterns in nonzero rows.
1 1 1 1 0 The number of patterns in the first row is 5. The 0 in the first row
0 0 0 1 1 requires a pair of 1's in that column and there are C(3,2) = 3 choices of
0 0 1 0 1 those. The second 0 down from the top in the last right column requires two
1 1 0 0 0 1's in that row, there are C(4,2) = 6 patterns of those. That same row has
0 0 0 0 0 two 0's that have 2 remaining choices in respective columns to add 1's.
The row of 4 can be placed in 4 ways. Finally, a doubling of choices comes
from transposition. Total choices = 5*3*6*2*4*2 = 1440.
.
Subset 3, 5 X 3, two rows with 4 1's. We count the number of patterns in nonzero rows.
1 1 0 0 0 Number of patterns of 2 1's in the first row is C(5,2) = 10.
1 0 1 1 1 The next row of 4 requires a 0 in a shared column of 1 from first row.
0 1 1 1 1 The last row requires a 0 in the shared column of the other 1 from the
0 0 0 0 0 first row. The number of row selections is 3! Finally, a doubling of
0 0 0 0 0 choices comes from transposition. Total choices = 2*3!*C(5,2) = 120.
.
Subset 4, 4 X 4, one row of 4 1's and one column of 4 1's. We count the nonzero row
1 1 1 1 0 and nonzero columns. There are 4 choices for row of 4, and 4 choices for
1 1 0 0 0 column of 4. The remaining part has 6 patterns. Total choices = 4*4*6 = 96.
1 0 1 0 0
1 0 0 1 0
0 0 0 0 0
.
For n=5, a(5) = 2040*1 + 1440*5 + 120*10 + 96*25 = 12840.
CROSSREFS
KEYWORD
nonn,easy,new
AUTHOR
G. Gordon Thompson, May 21 2026
STATUS
approved
Numbers m such that rad(m)^omega(m) = m.
+0
0
1, 2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 36, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 100, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 196, 197, 199, 211, 223, 225, 227, 229, 233, 239
OFFSET
1,2
COMMENTS
A term m > 1 is in the sequence if and only if the exponent of every prime factor in its prime factorization is exactly equal to omega(m). Thus, m must be of the form (p_1 * p_2 * ... * p_k)^k for distinct primes p_i.
Terms are mostly primes: for any p, omega(p) = 1 and rad(p) = p, so rad(p)^omega(p) = p ^ 1.
Except for 1 and the prime numbers (A000040), all terms are perfect powers.
Composite terms are the k-th powers of squarefree numbers having exactly k prime factors. For k = 2, these are the squares of squarefree semiprimes.
LINKS
EXAMPLE
m = 36 is a term since rad(36)^omega(36) = 6 ^ 2 = 36.
m = 8 is not a term since rad(8)^omega(8) = 2 ^ 1 = 2 != 8.
MATHEMATICA
q[m_]:=(Times@@First/@FactorInteger[m])^PrimeNu[m]==m; Select[Range[239], q]
PROG
(Python)
from math import prod
from sympy import primefactors
def ok(m): return prod(f := primefactors(m)) ** len(f) == m
print([m for m in range(1, 250) if ok(m)])
(PARI) isok(m) = {my(p = factor(m)[, 1]); m == vecprod(p)^#p; } \\ Amiram Eldar, Jul 02 2026
CROSSREFS
KEYWORD
nonn,new
AUTHOR
STATUS
approved
a(n) = denominator of A069835(n)/2^n.
+0
0
1, 1, 2, 1, 8, 4, 16, 8, 128, 64, 256, 128, 1024, 512, 2048, 1024, 32768, 16384, 65536, 32768, 262144, 131072, 524288, 262144, 4194304, 2097152, 8388608, 4194304, 33554432, 16777216, 67108864, 33554432, 2147483648, 1073741824, 4294967296, 2147483648
OFFSET
0,3
LINKS
H. Bateman, Some problems in potential theory, Messenger Math., 52 (1922), 71-78. [Annotated scanned copy]
FORMULA
Conjecture: a(2*n) = 2^A005187(n), a(2*n+1) = max(1,a(2*n)/2).
PROG
(PARI) a(n)=denominator(sum(k=0, n, binomial(n, k)^2*3^k)/2^n)
CROSSREFS
Cf. A069835, A397547 (numerators).
KEYWORD
frac,nonn,new
AUTHOR
Zhuorui He, Jul 02 2026
STATUS
approved
a(n) = numerator of A069835(n)/2^n.
+0
0
1, 2, 11, 17, 443, 743, 10159, 17593, 984467, 1734443, 24591493, 43793863, 1252842991, 2247613027, 32350749719, 58349238473, 6750044324867, 12224356333619, 177406570353769, 322314014477851, 9382052284364629, 17089477261939129, 249306090098352089
OFFSET
0,2
LINKS
H. Bateman, Some problems in potential theory, Messenger Math., 52 (1922), 71-78. [Annotated scanned copy]
PROG
(PARI) a(n)=numerator(sum(k=0, n, binomial(n, k)^2*3^k)/2^n)
CROSSREFS
Cf. A069835, A397606 (denominators).
KEYWORD
frac,nonn,new
AUTHOR
Zhuorui He, Jun 30 2026
STATUS
approved
Primes in A161601 whose binary reversal is prime.
+0
1
11, 23, 37, 43, 47, 67, 71, 83, 131, 151, 163, 167, 173, 199, 223, 263, 269, 277, 283, 307, 331, 349, 359, 383, 431, 463, 479, 521, 571, 599, 601, 619, 631, 643, 653, 661, 683, 691, 701, 727, 739, 823, 827, 839, 911, 1033, 1049, 1087, 1091, 1093, 1109, 1117
OFFSET
1,1
LINKS
EXAMPLE
The pairs (A161601(k), A161603(k)) for k=1..5 are (11,13), (19,25), (23,29), (35,49), (37,41), so this sequence begins with 11,23,37.
For (37,41), we have (37 base 2) = 100101 and (41 base 2) = 10101, which is the reverse of 100101, as in A161601(5) and A161603(5).
MATHEMATICA
p[f_] := Module[{n = 1, c = 0}, Reap[While[c < f, With[{d = IntegerDigits[n, 2]}, With[{r = Reverse[d]}, If[d =!= r && n < FromDigits[r, 2],
Sow[{n, FromDigits[r, 2]}]; c++]]]; n += 2]][[2, 1]]];
u = Select[p[300], PrimeQ[#[[1]]] && PrimeQ[#[[2]]] &];
Map[First, u]
(* Peter J. C. Moses, Jun 24 2026 *)
KEYWORD
nonn,easy,base,new
AUTHOR
Clark Kimberling, Jun 24 2026
STATUS
approved
Composites in A161601 whose binary reversal is prime.
+0
1
55, 77, 91, 143, 203, 275, 295, 323, 355, 395, 407, 415, 515, 517, 533, 535, 553, 559, 575, 595, 637, 671, 703, 731, 755, 767, 847, 895, 943, 1037, 1043, 1067, 1079, 1081, 1111, 1147, 1157, 1159, 1183, 1189, 1243, 1267, 1271, 1331, 1333, 1357, 1363, 1391, 1411, 1435, 1463, 1507
OFFSET
1,1
LINKS
MATHEMATICA
p[f_] := Module[{n = 1, c = 0}, Reap[While[c < f, With[{d = IntegerDigits[n, 2]}, With[{r = Reverse[d]}, If[d =!= r && n < FromDigits[r, 2],
Sow[{n, FromDigits[r, 2]}]; c++]]]; n += 2]][[2, 1]]];
u = Select[p[300], !PrimeQ[#[[1]]] && PrimeQ[#[[2]]] &];
Map[First, u]
(* Peter J. C. Moses, Jun 24 2026 *)
KEYWORD
nonn,easy,base,new
AUTHOR
Clark Kimberling, Jun 24 2026
EXTENSIONS
a(40) corrected and more terms from Vincenzo Librandi, Jul 07 2026
STATUS
approved
Composites in A161601 whose binary reversal is composite.
+0
1
35, 39, 69, 75, 87, 95, 111, 133, 135, 141, 147, 155, 159, 171, 175, 183, 187, 207, 215, 259, 261, 265, 267, 279, 285, 287, 291, 299, 301, 303, 309, 315, 319, 327, 333, 335, 339, 343, 351, 363, 371, 375, 391, 399, 411, 423, 447, 519, 525, 527, 529, 531, 5373
OFFSET
1,1
MATHEMATICA
p[f_] := Module[{n = 1, c = 0}, Reap[While[c < f, With[{d = IntegerDigits[n, 2]}, With[{r = Reverse[d]}, If[d =!= r && n < FromDigits[r, 2],
Sow[{n, FromDigits[r, 2]}]; c++]]]; n += 2]][[2, 1]]];
u = Select[p[300], !PrimeQ[#[[1]]] && !PrimeQ[#[[2]]] &];
Map[First, u]
(* Peter J. C. Moses, Jun 24 2026 *)
KEYWORD
nonn,base,easy,new
AUTHOR
Clark Kimberling, Jun 26 2026
STATUS
approved
Number of ways to tile a 2 X n strip with one corner missing, using squares, dominoes, and T-shaped tetrominoes.
+0
0
0, 1, 3, 11, 38, 132, 458, 1589, 5513, 19127, 66360, 230232, 798776, 2771305, 9614875, 33358227, 115734350, 401533324, 1393095570, 4833260781, 16768705809, 58178010095, 201845085552, 700289310256, 2429611385968, 8429389682769, 29245257424435, 101464650942619
OFFSET
0,3
COMMENTS
Compare to A397582 which counts these tilings for the full 2 X n strip (without a missing corner).
FORMULA
G.f.: x/(1 - 3*x - 2*x^2 + x^3 + x^4).
a(n) = 3*a(n-1) + 2*a(n-2) - a(n-3) - a(n-4).
a(n) = a(n-1) + a(n-2) + A397582(n-1).
EXAMPLE
Here is one of the a(6)=458 ways to tile this 2 X 6 strip (with one missing corner) using squares, dominoes, and T-shaped tetrominoes:
_________
|_ _| |_|_
|_|_|_|_|___|.
MATHEMATICA
LinearRecurrence[{3, 2, -1, -1}, {0, 1, 3, 11}, 50]
CROSSREFS
Cf. A397582.
KEYWORD
nonn,easy,new
AUTHOR
Greg Dresden and Yutang (Daniel) Bai, Jul 01 2026
STATUS
approved
Expansion of ordinary generating function A(x) satisfying (A(x)^2-A(x^2))/2 = x/(1-x)^2.
+0
0
1, 1, 2, 1, 2, 1, 1, 2, 2, 0, 2, 2, 0, 3, 2, -3, 7, 2, -10, 15, 3, -24, 37, -5, -60, 106, -29, -159, 293, -115, -396, 824, -441, -983, 2375, -1569, -2455, 6849, -5281, -5928, 19665, -17604, -13629, 56730, -57761, -29324, 163281, -186758, -53979, 467872, -598545, -60594, 1334821, -1904160, 119258, 3784805
OFFSET
0,3
COMMENTS
This sequence was inspired by the analysis of A397174 and is a variant of A397525.
FORMULA
Sum_{k=0..n} a(k)*a(n-k) = 2*n + [n even] a(n/2).
MAPLE
a:= proc(n) option remember; `if`(n=0, 1, n+`if`(n::even, (t->
t*(1-t)/2)(a(n/2)), 0)-add(a(k)*a(n-k), k=1..(n-1)/2))
end:
seq(a(n), n=0..55); # Alois P. Heinz, Jun 30 2026
PROG
(PARI)
h(v) = (Ser(v)^2-subst(Ser(v), x, x^2))/2
ListA(max_n) = {my(v=[1]); while(#v<max_n, v=concat(v, #v-polcoeff(Vec(h(concat(v, y)))[#v], 0))); vector(#v, k, v[k])}
CROSSREFS
KEYWORD
sign,new
AUTHOR
Thomas Scheuerle, Jun 29 2026
STATUS
approved
Number of pairs (P, I) where P is a partially ordered set on n labeled elements and I is an order ideal (down-set) of P; equivalently, sum over labeled posets P on [n] of the number of antichains of P.
+0
0
1, 2, 10, 98, 1678, 46922, 2049550, 135499898, 13243258318, 1878894285002, 381362574101710, 109430505055180058, 43954986389427881998, 24508712037068391623402, 18835419731708228516982670, 19828793128398628087131589178, 28441643117705315333254490986318
OFFSET
0,2
COMMENTS
a(n) = Sum_{P} d(P) over the A001035(n) labeled posets P on [n], where d(P) is the number of antichains (equivalently order ideals / down-sets) of P; so a(n) counts the pairs (poset, down-set). This is the k=1 member of the antichain-count moment family k=1..4 used to compute A001035(19).
CROSSREFS
Cf. A001035 (k=0, number of labeled posets), A000798.
KEYWORD
nonn,new
AUTHOR
Rafael Ayala, Jun 30 2026
STATUS
approved

Search completed in 0.060 seconds