[go: up one dir, main page]

login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Search: a365406 -id:a365406
     Sort: relevance | references | number | modified | created      Format: long | short | data
Numbers j for which the symmetric representation of sigma(j) has two parts, each of width one.
+10
21
3, 5, 7, 10, 11, 13, 14, 17, 19, 22, 23, 26, 29, 31, 34, 37, 38, 41, 43, 44, 46, 47, 52, 53, 58, 59, 61, 62, 67, 68, 71, 73, 74, 76, 79, 82, 83, 86, 89, 92, 94, 97, 101, 103, 106, 107, 109, 113, 116, 118, 122, 124, 127, 131, 134, 136, 137, 139, 142, 146, 148, 149, 151, 152, 157, 158, 163, 164, 166, 167, 172, 173, 178, 179, 181, 184, 188, 191, 193, 194, 197, 199
OFFSET
1,1
COMMENTS
The sequence is the intersection of A239929 (sigma(j) has two parts) and of A241008 (sigma(j) has an even number of parts of width one).
The numbers in the sequence are precisely those defined by the formula for the triangle, see the link. The symmetric representation of sigma(j) has two parts, each part having width one, precisely when j = 2^(k - 1) * p where 2^k <= row(j) < p, p is prime and row(j) = floor((sqrt(8*j + 1) - 1)/2). Therefore, the sequence can be written naturally as a triangle as shown in the Example section.
The symmetric representation of sigma(j) = 2*j - 2 consists of two regions of width 1 that meet on the diagonal precisely when j = 2^(2^m - 1)*(2^(2^m) + 1) where 2^(2^m) + 1 is a Fermat prime (see A019434). This subsequence of numbers j is 3, 10, 136, 32896, 2147516416, ...[?]... (A191363).
The k-th column of the triangle starts in the row whose initial entry is the first prime larger than 2^(k+1) (that sequence of primes is A014210, except for 2).
Observation: at least the first 82 terms coincide with the numbers j with no middle divisors whose largest divisor <= sqrt(j) is a power of 2, or in other words, coincide with the intersection of A071561 and A365406. - Omar E. Pol, Oct 11 2023
FORMULA
Formula for the triangle of numbers associated with the sequence:
P(n, k) = 2^k * prime(n) where n >= 2, 0 <= k <= floor(log_2(prime(n)) - 1).
EXAMPLE
We show portions of the first eight columns, 0 <= k <= 7, of the triangle.
0 1 2 3 4 5 6 7
3
5 10
7 14
11 22 44
13 26 52
17 34 68 136
19 38 76 152
23 46 92 184
29 58 116 232
31 62 124 248
37 74 148 296 592
41 82 164 328 656
43 86 172 344 688
47 94 188 376 752
53 106 212 424 848
59 118 236 472 944
61 122 244 488 976
67 134 268 536 1072 2144
71 142 284 568 1136 2272
. . . . . .
. . . . . .
127 254 508 1016 2032 4064
131 262 524 1048 2096 4192 8384
137 274 548 1096 2192 4384 8768
. . . . . . .
. . . . . . .
251 502 1004 2008 4016 8032 16064
257 514 1028 2056 4112 8224 16448 32896
263 526 1052 2104 4208 8416 16832 33664
Since 2^(2^4) + 1 = 65537 is the 6543rd prime, column k = 15 starts with 2^15*(2^(2^16) + 1) = 2147516416 in row 6542 with 65537 in column k = 0.
For an image of the symmetric representations of sigma(m) for all values m <= 137 in the triangle see the link.
The first column is the sequence of odd primes, see A065091.
The second column is the sequence of twice the primes starting with 10, see A001747.
The third column is the sequence of four times the primes starting with 44, see A001749.
For related references also see A033676 (largest divisor of n less than or equal to sqrt(n)).
MATHEMATICA
(* functions path[] and a237270[ ] are defined in A237270 *)
atmostOneDiagonalsQ[n_]:=SubsetQ[{0, 1}, Union[Flatten[Drop[Drop[path[n], 1], - 1] - path[n - 1], 1]]]
(* data *)
Select[Range[200], Length[a237270[#]]==2 && atmostOneDiagonalsQ[#]&]
(* function for computing triangle in the Example section through row 55 *)
TableForm[Table[2^k Prime[n], {n, 2, 56}, {k, 0, Floor[Log[2, Prime[n]] - 1]}], TableDepth->2]
KEYWORD
nonn
AUTHOR
Hartmut F. W. Hoft, Sep 08 2014
STATUS
approved
Numbers j whose largest divisor <= sqrt(j) is not a power of 2.
+10
2
9, 12, 15, 18, 21, 25, 27, 30, 33, 35, 36, 39, 40, 42, 45, 48, 49, 50, 51, 54, 55, 56, 57, 60, 63, 65, 66, 69, 70, 75, 77, 78, 81, 84, 85, 87, 90, 91, 93, 95, 98, 99, 100, 102, 105, 108, 110, 111, 114, 115, 117, 119, 120, 121, 123, 125, 126, 129, 130, 132, 133, 135, 138, 140, 141, 143, 144, 145, 147, 150, 153
OFFSET
1,1
COMMENTS
Also indices of terms in A033676 that are not a power of 2.
Differs from A342082 in not having {24, 72, 80, 96, 112, ...}. - Hugo Pfoertner, Oct 15 2023
MATHEMATICA
A365408Q[n_]:=With[{d=Divisors[n]}, !IntegerQ[Log2[d[[Ceiling[Length[d]/2]]]]]];
Select[Range[200], A365408Q] (* Paolo Xausa, Oct 19 2023 *)
PROG
(Python)
from itertools import count, islice
from sympy import divisors
def A365408_gen(startvalue=1): # generator of terms >= startvalue
return filter(lambda i:(a:=(d:=divisors(i))[len(d)-1>>1])!=1<<a.bit_length()-1, count(max(startvalue, 1)))
A365408_list = list(islice(A365408_gen(), 30)) # Chai Wah Wu, Oct 18 2023
CROSSREFS
Complement of A365406.
KEYWORD
nonn
AUTHOR
Omar E. Pol, Oct 14 2023
STATUS
approved
a(n) = 1 if the largest divisor of n <= sqrt(n) is a power of 2, otherwise 0.
+10
2
1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 0, 1, 0, 0, 1, 1, 0, 1, 1, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0
OFFSET
1
FORMULA
a(n) = A209229(A033676(n)).
a(n) >= A365605(n).
MATHEMATICA
Table[If[IntegerQ[Log2[Max[Select[Divisors[n], #<=Sqrt[n]&]]]], 1, 0], {n, 110}] (* Harvey P. Dale, May 26 2024 *)
PROG
(PARI)
A209229(n) = (n && !bitand(n, n-1));
A033676(n) = { my(d); if(n<2, 1, d=divisors(n); d[(length(d)+1)\2]) }; \\ From A033676.
CROSSREFS
Characteristic function of A365406.
KEYWORD
nonn
AUTHOR
Antti Karttunen, Oct 22 2023
STATUS
approved
Numbers k whose largest divisor <= sqrt(k) is a power of 2, listing only the first such number with any given prime signature.
+10
0
1, 2, 4, 6, 8, 16, 20, 24, 32, 64, 72, 80, 96, 128, 256, 288, 320, 336, 384, 512, 1024, 1056, 1152, 1280, 1344, 1536, 2048, 4096, 4224, 4608, 4800, 5120, 5376, 6144, 8192, 16384, 16896, 17280, 18432, 18816, 19200, 20480, 21504, 24576, 32768, 65536, 67584, 69120, 69888
OFFSET
1,2
COMMENTS
This sequence is a primitive sequence related to A365406 in the sense that it can be used to find the smallest term k in A365406 such that tau(k), omega(k) or bigomega(k) has some particular value.
Not every prime signature produces a term. For example no term has prime signature (3, 2, 1). Proof: any number with prime signature (3, 2, 1) has 24 divisors. Hence the 12th divisor must be a power of 2. But the largest power of 2 such number can have as a divisor is 8. 8 can never be the 12th divisor of a number. Therefore (3, 2, 1) can never be the prime signature of a term.
EXAMPLE
k = 20 = 2^2 * 5 is in the sequence as it has prime signature (2, 1) and its largest divisor <= sqrt(k) is 4, a power of 2. It is the smallest such number since smaller numbers with prime signature (2, 1), namely 12 and 18, do not have the relevant divisor being a power of 2.
PROG
(PARI)
upto(n) = {
my(res = List([1]), m = Map());
forstep(i = 2, n, 2,
if(isok(i),
s = sig(i);
sb = sigback(s);
if(!mapisdefined(m, sb),
listput(res, i);
mapput(m, sb, i)
)
)
);
res
}
sig(n) = {
vecsort(factor(n)[, 2], , 4)
}
sigback(v) = {
my(pr = primes(#v));
prod(i = 1, #v, pr[i]^v[i])
}
isok(n) = my(d = divisors(n)); hammingweight(d[(#d + 1)\2]) == 1
CROSSREFS
KEYWORD
nonn
AUTHOR
David A. Corneth, Oct 21 2023
EXTENSIONS
Edited by Peter Munn, Oct 26 2023
STATUS
approved
Square array read by upward antidiagonals: T(n,k) is the n-th number j with the property that the parts of the symmetric representation of sigma(j) are two s-gon of width 1, where s = 2^(k+1), n >= 1, k >= 1.
+10
0
3, 5, 10, 7, 14, 44, 11, 22, 52, 136, 13, 26, 68, 152, 592, 17, 34, 76, 184, 656
OFFSET
1,1
COMMENTS
For column k = 1, 2, 3, 4, 5, ... the number of sides of the mentioned s-gon are respectively 4, 8, 16, 32, 64, ...
Conjecture 1: column k gives the row numbers of the triangle A364639 where the rows are [1, A036563(k+1) zeros, -1, 1] or where the rows start with [1, A036563(k+1) zeros, -1, 1] and the remaining terms are zeros.
Conjecture 2: every column gives a subsequence of A246955.
Conjecture 3: the sequence is infinite.
Observation 1: at least the terms <= 199 in increasing order coincide with at least the first 82 terms of the intersection of A071561 and A365406.
Observation 2: in the Example section of A246955 there is an irregular triangle. It seems that the terms sorted of the triangle give the sequence A246955. At least the first r(k) terms in the column (k - 1) of the triangle coincide with the first r(k) terms of the column k of this square array, where r(k) are 19, 18, 16, 14, 7 for k = 1..5 respectively.
Observation 3: at least the first five terms of the row 1 coincide with the first five terms of A246956.
EXAMPLE
The corner of the square array is as shown below:
3, 10, 44, 136, 592, ...
5, 14, 52, 152, 656, ...
7, 22, 68, 184, 688, ...
11, 26, 76, 232, 752, ...
13, 34, 92, 248, 848, ...
17, 38, 116, 296, 944, ...
19, 46, 124, 328, 976, ...
...
KEYWORD
nonn,tabl,more
AUTHOR
Omar E. Pol, Sep 25 2023
STATUS
approved

Search completed in 0.008 seconds