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-3 of 3 results.

A349308 Numbers k such that A321167(k) = A321167(k+1) > 1.

Original entry on oeis.org

80, 135, 296, 343, 375, 624, 728, 1160, 1431, 1592, 1624, 2240, 2295, 2456, 2511, 2624, 2727, 2888, 3429, 3591, 3624, 3752, 3992, 4023, 4184, 4671, 4887, 4913, 5048, 5144, 5264, 5319, 5480, 5696, 6183, 6344, 6375, 6591, 6615, 6776, 6858, 6859, 7479, 7624, 7640
Offset: 1

Views

Author

Amiram Eldar, Nov 14 2021

Keywords

Comments

Without the restriction that A321167(k) > 1, all the terms of A340152 would be in this sequence.
In contrast to A001274, which has only one known pair of consecutive terms (5186 and 5187), this sequence seems to have many pairs of consecutive terms. The smaller members of these pairs are 6858, 13375, 22625, ...

Examples

			80 is a term since A321167(80) = A321167(81) = 3.
		

Crossrefs

Subsequence of A068140.
Similar sequences: A001274, A287055, A293184, A326403, A349307.

Programs

  • Mathematica
    f[p_, e_] := p^e - 1; uphi[1] = 1; uphi[n_] := Times @@ f @@@ FactorInteger[n]; fe[p_, e_] := uphi[e]; euphi[n_] := Times @@ fe @@@ FactorInteger[n]; Select[Range[8000], euphi[#] == euphi[# + 1] > 1 &]

A349309 Numbers k such that A254926(k) = A254926(k+1).

Original entry on oeis.org

7, 26, 124, 342, 1330, 2196, 12166, 24388, 29790, 79506, 103822, 148876, 205378, 226980, 300762, 357910, 493038, 571786, 1030300, 1092726, 1225042, 2248090, 2685618, 3307948, 3442950, 3869892, 4657462, 5177716, 5735338, 6967870, 7645372, 9393930, 11089566, 11697082
Offset: 1

Views

Author

Amiram Eldar, Nov 14 2021

Keywords

Examples

			7 is a term since A254926(7) = A254926(8) = 7.
		

Crossrefs

Cf. A254926.

Programs

  • Mathematica
    f[p_, e_] := If[e < 3, p^e, p^e - p^(e - 3)]; s[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[10^6], s[#] == s[# + 1] &]
  • Python
    from math import prod
    from itertools import count, islice
    from sympy import factorint
    def A349309_gen(startvalue=1): # generator of terms >= startvalue
        a = prod(p**e - (p**(e-3) if e >= 3 else 0) for p, e in factorint(max(startvalue,1)).items())
        for k in count(max(startvalue,1)):
            b = prod(p**e - (p**(e-3) if e >= 3 else 0) for p, e in factorint(k+1).items())
            if a == b:
                yield k
            a = b
    A349309_list = list(islice(A349309_gen(),10)) # Chai Wah Wu, Jan 24 2022

A385743 Numbers k such that A384247(k) = A384247(k+1).

Original entry on oeis.org

1, 20, 27, 35, 63, 64, 104, 143, 194, 208, 740, 836, 1220, 1299, 1419, 1803, 1892, 2625, 3255, 3705, 3716, 3843, 4096, 5184, 5186, 5635, 5695, 7868, 10659, 13365, 16904, 17948, 18507, 18914, 21007, 22935, 25388, 25545, 27675, 30380, 31599, 32304, 32864, 34595
Offset: 1

Views

Author

Amiram Eldar, Jul 08 2025

Keywords

Comments

63 is the only number k below 10^11 such that A384247(k) = A384247(k+1) = A384247(k+2). Are there any other such terms?

Examples

			1 is a term since A384247(1) = A384247(2) = 1.
20 is a term since A384247(20) = A384247(21) = 12.
		

Crossrefs

Cf. A384247.
Similar sequences: A001274, A287055, A293184, A301866, A326403, A349307.

Programs

  • Mathematica
    f[p_, e_] := p^e*(1 - 1/p^(2^(IntegerExponent[e, 2]))); iphi[1] = 1; iphi[n_] := iphi[n] = Times @@ f @@@ FactorInteger[n]; Select[Range[35000], iphi[#] == iphi[# + 1] &]
  • PARI
    iphi(n) = {my(f = factor(n)); n * prod(i = 1, #f~, (1 - 1/f[i, 1]^(1 << valuation(f[i, 2], 2)))); }
    list(lim) = {my(s1 = iphi(1), s2); for(k = 2, lim, s2 = iphi(k); if(s1 == s2, print1(k-1, ", ")); s1 = s2);}
Showing 1-3 of 3 results.