A340349 a(n) is the smallest k such that A292849(k) = 2n-1.

1, 3, 13, 5, 57, 35, 21, 9, 241, 219, 49, 45, 169, 83, 73, 17, 993, 59, 941, 53, 3197, 51, 185, 93, 209, 81, 349, 85, 41, 89, 105, 33, 4033, 491, 4749, 247, 449, 227, 429, 363, 3249, 401, 193, 259, 233, 107, 117, 189, 697249, 1355, 173, 517, 473, 1091, 101, 231, 725, 305

Offset: 1

Author

Thomas Scheuerle, Jan 05 2021

Keywords

• nonn
• base

Comments

This implies that a(n)=k is a solution to A000120(2*n-1) = A000120((2*n-1)*k), where A000120 is the Hamming weight. If no such number exists we define a(n) = 2.

Does this sequence consist only of odd numbers? It appears so.

This will be the case if A292849 contains all odd numbers, because a(n) will then never become 2.

A292849 must contain all odd numbers if two conditions are satisfied:

First: For every odd number 2n-1 there must be an odd k > 1 that satisfies A000120(2n-1) = A000120((2n-1)*k). To prove that this condition is satisfied, it may be helpful that if k*m = 2^j+r, we know that A000120(k*m) = A063787(r) and for each Hamming weight there exists an r such that A063787(r) = A063787(r+1). This allows us to choose r such that the Hamming weight becomes A000120(m) = A063787(r). For a given r, k*m = 2^j+r may have no solutions if k or m are divisors of r, but there may still exist a solution for k*m = 2^j+r+1. Of course this is not a complete proof.

Second: For every n there needs to be a number k such that 2n-1 is the smallest solution to A000120(2n-1) = A000120((2n-1)*k). This is satisfied if no row in A340441 is a subset of a previous row.

No odd number can appear more than once in this sequence, but not all odd numbers will appear, so this sequence is not a permutation of the odd numbers.

This sequence can be constructed from A340441. Start with a(1) = 1, then a(n) is the least number in row n-1 of A340441 that has not already appeared in A340441.

Crossrefs

Cf. A000120, A292849, A340069, A263132, A077459, A295827, A340441, A340351 table of multiplications.

Programs

MATLAB

function a = A340349(maxA292849)
c = A340351(maxA292849,1);
n = 1; run = 1;
while run == 1
     i = find(c==(n*2)-1);
     if ~isempty(i);
         a(n) = i(1);
         n = n+1;
     else
         run = 0;
     end
end
end
function a = A340351(max_n,max_m)
    for n = 1:max_n
        m = 1; k = 1;
        while m < max_m+1
            c = length(find(bitget(k,1:32)== 1));
            if c == length(find(bitget(n*k,1:32)== 1))
                a(n,m) = k;
                m = m+1;
            end
            k = k +1;
        end
    end
end

PARI

f(n) = my(k=1); while ((hammingweight(k)) != hammingweight(k*n), k++); k; \\ A292849
a(n) = my(k=1); while(f(k) != 2*n-1, k++); k; \\ Michel Marcus, Jan 09 2021