A324107 -id:A324107 - cp's OEIS frontend

A324106 Multiplicative with a(p^e) = A005940(p^e).

1, 2, 3, 4, 5, 6, 9, 8, 7, 10, 15, 12, 25, 18, 15, 16, 11, 14, 21, 20, 27, 30, 45, 24, 49, 50, 75, 36, 125, 30, 81, 32, 45, 22, 45, 28, 55, 42, 75, 40, 77, 54, 105, 60, 35, 90, 135, 48, 121, 98, 33, 100, 245, 150, 75, 72, 63, 250, 375, 60, 625, 162, 63, 64, 125, 90, 39, 44, 135, 90, 99, 56, 91, 110, 147, 84, 135, 150, 189, 80, 143, 154, 231, 108, 55

Offset: 1

Antti Karttunen, Feb 15 2019

Question: are there any other numbers n besides 1 and those in A070776, for which a(n) = A005940(n)? At least not below 2^25. This is probably easy to prove.

			For n = 85 = 5*17, a(85) = A005940(5) * A005940(17) = 5*11 = 55. Note that A005940(5) is obtained from the binary expansion of 5-1 = 4, which is "100", and A005940(17) is obtained from the binary expansion of 17-1 = 16, which is "1000".

Antti Karttunen, Table of n, a(n) for n = 1..16384
Michael De Vlieger, Fan style binary tree of a(n), n = 1..2^12, color coded to show the smallest values in the range r = (2^r - 1)..2^(r+1) in blue and highlighting the largest with red.
Index entries for sequences related to binary expansion of n

Cf. A005940, A070776, A324107 (fixed points), A324108, A324109.

nn = 128; Array[Set[a[#], #] &, 2]; Do[If[EvenQ[n], Set[a[n], 2 a[n/2]], Set[a[n], Times @@ Power @@@ Map[{Prime[PrimePi[#1] + 1], #2} & @@ # &, FactorInteger[a[(n + 1)/2]]]]], {n, 3, nn}]; Array[Times @@ Map[a, Power @@@ FactorInteger[#]] &, nn] (* Michael De Vlieger, Sep 18 2022 *)

A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); t }; \\ From A005940
A324106(n) = { my(f=factor(n)); prod(i=1, #f~, A005940(f[i,1]^f[i,2])); };

A324109 Numbers n such that A324108(n) = A324054(n-1).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 16, 17, 18, 19, 20, 22, 23, 24, 25, 26, 27, 28, 29, 31, 32, 34, 36, 37, 38, 40, 41, 43, 44, 46, 47, 48, 49, 50, 52, 53, 54, 56, 58, 59, 61, 62, 64, 67, 68, 71, 72, 73, 74, 76, 79, 80, 81, 82, 83, 86, 87, 88, 89, 92, 94, 96, 97, 98, 100, 101, 103, 104, 106, 107, 108, 109, 112, 113, 116, 118, 121

Offset: 1

Antti Karttunen and David A. Corneth, Feb 15 2019

nonn

Numbers n such that A324054(n-1) is equal to A324108(n), which is a multiplicative function with A324108(p^e) = A324054((p^e)-1).
Prime powers (A000961) is a subsequence by definition.
Also A070776 is a subsequence. This follows because for every n of the form 2^i * p^j (where p is an odd prime, and i >= 0, j >= 0), we have A324108(2^i * p^j) = A324054(2^i - 1)*A324054(p^j - 1) = sigma(A005940(2^i)) * sigma(A005940(p^j)). Because A005940(1) = 1, and A005940(2n) = 2*A005940(n), the powers of two are among the fixed points of A005940 (cf. A029747), thus the left half of product is sigma(2^i), while on the other hand, we know that A005940(p^j) is odd (because A005940 also preserves parity), and thus the whole product is equal to sigma(2^i * A005940(p^j)) = sigma(A005940(2^i * p^j)) = A324054((2^i * p^j)-1).
See subsequence A324111 for less regular solutions.

Antti Karttunen, Table of n, a(n) for n = 1..10001

Union of A070776 and A324111.

Cf. A000961 (a subsequence), A029747, A324054, A324107, A324108, A324110 (complement).

A324054(n) = { my(p=2,mp=p*p,m=1); while(n, if(!(n%2), p=nextprime(1+p); mp = p*p, if(3==(n%4),mp *= p,m *= (mp-1)/(p-1))); n>>=1); (m); };
A324108(n) = { my(f=factor(n)); prod(i=1, #f~, A324054((f[i,1]^f[i,2])-1)); };
isA324109(n) = (A324054(n-1)==A324108(n));
for(n=1,121,if(isA324109(n), print1(n,", ")));

cp's OEIS Frontend

A324106 Multiplicative with a(p^e) = A005940(p^e).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI