cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 12 results. Next

A048674 Fixed points of A048673 and A064216: Numbers n such that if n = product_{k >= 1} (p_k)^(c_k), then Product_{k >= 1} (p_{k+1})^(c_k) = (2*n)-1, where p_k indicates the k-th prime, A000040(k).

Original entry on oeis.org

1, 2, 3, 25, 26, 33, 93, 1034, 970225, 8550146, 325422273, 414690595, 1864797542, 2438037206
Offset: 1

Views

Author

Antti Karttunen, Jul 14 1999

Keywords

Comments

Equally: after 1, numbers n such that, if the prime factorization of 2n-1 = Product_{k >= 1} (p_k)^(c_k) then Product_{k >= 1} (p_{k-1})^(c_k) = n.
Factorization of the initial terms: 1, 2, 3, 5^2, 2*13, 3*11, 3*31, 2*11*47, 5^2*197^2, 2*11*47*8269, 3*11*797*12373, 5*11^2*433*1583, 2*23*59*101*6803, 2*11*53*1201*1741.
The only 3-cycle of permutation A048673 in range 1 .. 402653184 is (2821 3460 5639).
For 2-cycles, take setwise difference of A245449 and this sequence.
Numbers k for which A336853(k) = k-1. - Antti Karttunen, Nov 26 2021

Examples

			25 is present, as 2*25 - 1 = 49 = p_4^2, and p_3^2 = 5*5 = 25.
26 is present, as 2*26 - 1 = 51 = 3*17 = p_2 * p_8, and p_1 * p_7 = 2*13 = 26.
Alternatively, as 26 = 2*13 = p_1 * p_7, and ((p_2 * p_8)+1)/2 = ((3*17)+1)/2 = 26 also, thus 26 is present.
		

Crossrefs

Fixed points of permutation pair A048673/A064216.
Positions of zeros in A349573.
Subsequence of the following sequences: A245449, A269860, A319630, A349622, A378980 (see also A379216).
This sequence is also obtained as a setwise difference of the following pairs of sequences: A246281 \ A246351, A246352 \ A246282, A246361 \ A246371, A246372 \ A246362.
Cf. also A348514 (fixed points of map A108228, similar to A048673).

Programs

  • Maple
    A048673 := n -> (A003961(n)+1)/2;
    A048674list := proc(upto_n) local b,i; b := [ ]; for i from 1 to upto_n do if(A048673(i) = i) then b := [ op(b), i ]; fi; od: RETURN(b); end;
  • Mathematica
    Join[{1}, Reap[For[n = 1, n < 10^7, n++, ff = FactorInteger[n]; If[Times @@ Power @@@ (NextPrime[ff[[All, 1]]]^ff[[All, 2]]) == 2 n - 1, Print[n]; Sow[n]]]][[2, 1]]] (* Jean-François Alcover, Mar 04 2016 *)
  • PARI
    A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    isA048674(n) = ((n+n)==(1+A003961(n))); \\ Antti Karttunen, Nov 26 2021

Extensions

Entry revised and the names in Maple-code cleaned by Antti Karttunen, Aug 25 2014
Terms a(11) - a(14) added by Antti Karttunen, Sep 11-13 2014

A378980 Numbers k such that (A003961(k)-2*k) divides (A003961(k)-sigma(k)), where A003961 is fully multiplicative with a(prime(i)) = prime(i+1), and sigma is the sum of divisors function.

Original entry on oeis.org

1, 2, 3, 4, 6, 7, 10, 25, 26, 28, 33, 46, 55, 57, 69, 91, 93, 496, 1034, 1054, 1558, 2211, 2626, 4825, 8128, 11222, 12046, 12639, 28225, 32043, 68727, 89575, 970225, 1392386, 2245557, 8550146, 12371554, 16322559, 22799825, 33550336, 48980427, 51326726, 55037217, 60406599, 68258725, 142901438, 325422273, 342534446
Offset: 1

Views

Author

Antti Karttunen, Dec 12 2024

Keywords

Comments

Numbers k such that A252748(k) divides A286385(k).
Conjecture: Apart from a(5)=6, this is a subsequence of A319630, i.e., for all terms k<>6, gcd(k, A003961(k)) = 1. See also A372562, A372566.

Crossrefs

Positions of 0's in A378981.
Subsequence of A263837.
Subsequences: A000396, A048674, A348514, A326134, A349753 (odd terms of this sequence).
Cf. also A378983.

Programs

  • Mathematica
    f1[p_, e_] := (p^(e + 1) - 1)/(p - 1); f2[p_, e_] := NextPrime[p]^e; q[k_] := Module[{fct = FactorInteger[k], m, s}, s = Times @@ f1 @@@ fct; m = Times @@ f2 @@@ fct; Divisible[m - s, m - 2*k]]; q[1] = True; Select[Range[10^5], q] (* Amiram Eldar, Dec 19 2024 *)
  • PARI
    A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    A378981(n) = { my(u=A003961(n)); ((u-sigma(n))%((2*n)-u)); };
    isA378980(n) = !A378981(n);

A326134 Numbers k such that A326057(k) is equal to A252748(k) and A252748(k) is not 1.

Original entry on oeis.org

6, 28, 69, 91, 496, 2211, 4825, 8128, 12639, 22799825, 33550336, 60406599, 68258725, 569173299, 794579511, 984210266, 2830283326, 8589869056, 10759889913, 80295059913, 85871289682
Offset: 1

Views

Author

Antti Karttunen, Jun 11 2019

Keywords

Comments

No other terms below 3221225472.
Numbers k such that A252748(k) [= A003961(k) - 2*k] <> 1 (i.e., k is not in A348514), and A286385(k) [= A003961(k) - A000203(k)] = m*A252748(k) for some positive integer m. Note that this entails that k is nonabundant (A000203(k) <= 2*k) and primeshift-abundant (A252748(k) > 2), thus this is a subsequence of A341614. - revised Dec 13 2024
This is a subsequence of A378980, see further comments there. - Antti Karttunen, Dec 13 2024

Examples

			28 is a term as A252748(28) = 43 > 1 and A286385(28) = 43, which is a multiple of 43.
69 is a term as A252748(69) = 7 > 1 and A286385(69) = 49 is a multiple of 7.
91 is a term as A252748(91) = 5 > 1 and A286385(91) = 75 is a multiple of 5.
		

Crossrefs

Subsequence of the following sequences: A246282, A341614, A378980.
Odd terms form a subsequence of A349753.

Programs

  • Mathematica
    Select[Range[10^5], And[#3 - #1 != 1, GCD[#3 - #1, #3 - #2] == #3 - #1] & @@ {2 #, DivisorSigma[1, #], Times @@ Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}] - Boole[# == 1]} &] (* Michael De Vlieger, Feb 22 2021 *)
  • PARI
    A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961
    isA326134(n) = { my(s = A003961(n), t = (s-(2*n)), u = s-sigma(n)); ((1!=t)&&!(u%t)&&((u/t)>0)); };

Extensions

a(18) from Antti Karttunen, Dec 14 2024
a(19)..a(21) from Antti Karttunen (from the b-file of A378980 computed by Amiram Eldar), Dec 20 2024

A349753 Odd numbers k for which A003961(k)-2k divides A003961(k)-sigma(k), where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Original entry on oeis.org

1, 3, 7, 25, 33, 55, 57, 69, 91, 93, 2211, 4825, 12639, 28225, 32043, 68727, 89575, 970225, 2245557, 16322559, 22799825, 48980427, 55037217, 60406599, 68258725, 325422273, 414690595, 569173299, 794579511, 10056372275, 10475647197, 10759889913, 11154517557
Offset: 1

Views

Author

Antti Karttunen, Dec 01 2021

Keywords

Comments

Numbers k for which A326057(k) = gcd(A003961(k)-2k, A003961(k)-sigma(k)) is equal to abs(A252748(k)) = |A003961(k)-2k|.
The odd terms of A326134 form a subsequence of this sequence. Unlike in A326134, here we don't constrain the value of A252748(k) = A003961(k)-2k, thus allowing also values <= +1. Because of that, the odd terms of A048674 and A348514 are all included here, for example 57 and 68727 that occur in A348514, and 1, 3, 25, 33, 93, 970225, 325422273, 414690595 that occur in A048674.
Conjecture (1): This is a subsequence of A319630, in other words, for all terms k, gcd(k, A003961(k)) = 1.
Conjecture (2): Apart from 1, there are no common terms with A349169, which would imply that no odd perfect numbers exist.
None of the 36 initial terms is Zumkeller, in A083207, because all are deficient (in A005100). See also A337372. - Antti Karttunen, Dec 05 2024

Crossrefs

Subsequence of A378980 (its odd terms).

Programs

  • Mathematica
    f[p_, e_] := NextPrime[p]^e; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; q[n_] := Divisible[(sn = s[n]) - DivisorSigma[1, n], sn - 2*n]; Select[Range[1, 10^6, 2], q] (* Amiram Eldar, Dec 04 2021 *)
  • PARI
    A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    isA349753(n) = if(!(n%2), 0, my(s = A003961(n), t = (s-(2*n)), u = s-sigma(n)); !(u%t));

A387411 Numbers k such that the odd part of (1+k) divides (1+A003961(k)), where A003961 is fully multiplicative with a(p) = nextprime(p).

Original entry on oeis.org

1, 3, 4, 7, 10, 15, 18, 23, 27, 31, 47, 57, 63, 95, 119, 127, 255, 348, 383, 415, 447, 511, 575, 695, 767, 959, 1023, 1054, 1071, 1535, 1919, 2047, 2626, 3471, 3839, 4095, 4415, 6815, 8191, 8703, 13823, 16383, 31743, 32767, 39895, 42367, 48127, 64607, 65535, 68727, 74495, 81919, 92159, 98303, 113535, 124671, 131071
Offset: 1

Views

Author

Antti Karttunen, Sep 01 2025

Keywords

Comments

Like in many sequences of this type, the criterion seems to select for numbers with a long tail of trailing 1-bits. Terms that are not in A004767 are: 1, 4, 10, 18, 57, 348, 1054, 2626, 675348, 1869741, 12371554, 14070141, 1158654378, 1673018314, etc.

Crossrefs

Subsequences: A000225, A348514 (which is also a subsequence of A387414).
For similar sequences, see A336700, A387410, A387415, A387410, A387418, A387419.

Programs

  • Mathematica
    a3961[x_] := Apply[Times, Prime[PrimePi[#1] + 1]^#2 & @@@ FactorInteger[x]] - Boole[x == 1];
    a265[x_] := x/2^IntegerExponent[x, 2];
    Select[Range[2^17], Divisible[1 + a3961[#], a265[# + 1] ] &] (* Michael De Vlieger, Sep 01 2025 *)
  • PARI
    A000265(n) = (n>>valuation(n,2));
    A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    isA387411(n) = !((1+A003961(n))%A000265(1+n));

A378983 Numbers k such that (A003961(k)-2*k) divides (A003961(k)-(1+sigma(k))), where A003961 is fully multiplicative with a(prime(i)) = prime(i+1), and sigma is the sum of divisors function.

Original entry on oeis.org

1, 2, 3, 4, 5, 8, 10, 11, 15, 16, 17, 25, 26, 29, 32, 33, 35, 39, 41, 57, 59, 64, 71, 93, 101, 107, 125, 128, 137, 149, 161, 179, 191, 197, 227, 239, 256, 269, 281, 311, 347, 419, 431, 461, 512, 521, 569, 599, 617, 641, 659, 782, 809, 821, 827, 857, 881, 1019, 1024, 1030, 1031, 1034, 1049, 1054, 1061, 1091, 1151
Offset: 1

Views

Author

Antti Karttunen, Dec 13 2024

Keywords

Comments

Conjecture: A202274 gives all terms of A028982 that occur in this sequence.

Examples

			For k=16 we have A003961(16) = 81, A003961(k)-2*k = 49, and 49 divides (A003961(k)-(1+sigma(k))) = 81-32 = 49, therefore 16 is included in this sequence.
For k=25 we have A003961(25) = 49, A003961(k)-2*k = -1, and -1 divides (A003961(k)-(1+sigma(k))) regardless of what the latter is, therefore 25 is included.
		

Crossrefs

Positions of 0's in A378982.
Subsequences: A048674, A348514, A202274.
Cf. also A378980.

Programs

A379216 Difference 2*k - A003961(k) computed for k for which this difference divides difference (A003961(k)-sigma(k)), where A003961 is fully multiplicative with a(prime(i)) = prime(i+1).

Original entry on oeis.org

1, 1, 1, -1, -3, 3, -1, 1, 1, -43, 1, 5, 19, -1, -7, -5, 1, -2005, 1, -1, 149, -193, -1, -3, -79243, 1243, 1253, -7, 51, 581, -1, 3093, 1, 155491, 919, 1, -1, 15833, -877, -4295498497, 5129369, 31, 5779339, -69187, -29, 6745, 1, 181, 1, 69197, -397, -117433, -101, -1, 1, 2759, 1, -29479, 1, -5626288431709, 29669, -1, -132239, -1, -1, 14591, -2267959, -3187, 787250461
Offset: 1

Views

Author

Antti Karttunen, Dec 20 2024

Keywords

Comments

Among the initial 69 terms, there are eleven +1's and eleven -1's. The former correspond in A378980 with those of its terms that are in A048674 (1, 2, 3, 25, 26, 33, 93, 1034, ...), while the latter here correspond in A378980 with those of its terms that are in A348514 (4, 10, 57, 1054, 2626, ...).

Crossrefs

Programs

Formula

a(n) = -A252748(A378980(n)).

A379233 Numbers k such that A003961(k) = 2k +- 3, multiplied by the sign of difference A003961(k)-2k, where A003961 is fully multiplicative with a(prime(i)) = prime(i+1).

Original entry on oeis.org

-5, 6, -7, -161, 1045, -2525, 2795, 4825, 9725, -159115, 307993, -359315, -18377525, 25484825
Offset: 1

Views

Author

Antti Karttunen, Dec 23 2024

Keywords

Comments

6 is the only term that is in A104210, the absolute values of all other terms residing in its complement, A319630, thus 3 occurs only once in A379231. Proof: If k is not a multiple of 3 and k is in A104210, then there are primes p > 3 and q = nextprime(p) that both divide k and q also divides A003961(k). However, q does not divide 2k +- 3, therefore the equation 2k +- 3 = A003961(k) is unsolvable in these cases. So let's assume that k is a multiple of 3, which immediately entails that k must be also even, for A003961(k) to be a multiple of 3. Let x = k/6; then the equation can be rewritten as 2*6*x +- 3 = A003961(6)*A003961(x) <=> 12x +- 3 = 15*A003961(x) <=> 3*(4x +- 1) = 3*5*A003961(x). The only value of x that satisfies the equation is x=1 (as A003961(n)>n for all n>1), hence k=6.
If it exists, abs(a(15)) > 2^32.

Crossrefs

Formula

{sign(A252748(k)) * k, for k such that abs(A252748(k)) = 3}.

A379235 Numbers k such that A003961(k) = 2k +- 5, multiplied by the sign of difference A003961(k)-2k, where A003961 is fully multiplicative with a(prime(i)) = prime(i+1).

Original entry on oeis.org

14, 15, -22, -46, 91, -2782, -269434, -1056574, 14129726, -25652506, 26594126, 34233062, 147087493
Offset: 1

Views

Author

Antti Karttunen, Dec 23 2024

Keywords

Comments

15 is the only term that is in A104210, the absolute values of all other terms residing in its complement, A319630, thus 5 occurs only once in A379231. Proof: If k is not a multiple of 5 and k is in A104210, then there are primes p (either p=2 or p > 5 and q = nextprime(p) that both divide k and q also divides A003961(k). However, q does not divide 2k +- 5, therefore the equation 2k +- 5 = A003961(k) is unsolvable in these cases. So let's assume that k is a multiple of 5, which immediately entails that k must be also a multiple of 3, for A003961(k) to be a multiple of 5. Let x = k/15; then the equation can be rewritten as 2*15*x +- 5 = A003961(15)*A003961(x) <=> 30x +- 5 = 35*A003961(x) <=> 5*(6x +- 1) = 5*7*A003961(x). The only value of x that satisfies the equation is x=1 (as A003961(n)>n for all n>1), hence k=15.
If it exists, abs(a(14)) > 2^32.

Crossrefs

Formula

{sign(A252748(k)) * k, for k such that abs(A252748(k)) = 5}.

A379237 Numbers k such that A003961(k) = 2k +- 7, multiplied by the sign of difference A003961(k)-2k, where A003961 is fully multiplicative with a(prime(i)) = prime(i+1).

Original entry on oeis.org

9, 35, -38, 39, -51, 69, -374, -4521, 7869, 10426, 12639, -16094, -29354, 102579, -103881, 1295206, -3298514, 4267318, 478642449, -2120241621
Offset: 1

Views

Author

Antti Karttunen, Dec 23 2024

Keywords

Comments

35 is the only term that is in A104210, the absolute values of all other terms residing in its complement, A319630, thus 7 occurs only once in A379231. Proof: If k is not a multiple of 7 and k is in A104210, then there are primes p (either p=2, p=3 or p > 7 and q = nextprime(p) that both divide k and q also divides A003961(k). However, q does not divide 2k +- 7, therefore the equation 2k +- 7 = A003961(k) is unsolvable in these cases. So let's assume that k is a multiple of 7, which immediately entails that k must be also a multiple of 5, for A003961(k) to be a multiple of 7. Let x = k/35; then the equation can be rewritten as 2*35*x +- 7 = A003961(35)*A003961(x) <=> 70x +- 7 = 77*A003961(x) <=> 7*(10x +- 1) = 7*11*A003961(x). The only value of x that satisfies the equation is x=1 (as A003961(n)>n for all n>1), hence k=35.
If it exists, abs(a(21)) > 2^32.

Crossrefs

Formula

{sign(A252748(k)) * k, for k such that abs(A252748(k)) = 7}.
Showing 1-10 of 12 results. Next