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.

Previous Showing 11-14 of 14 results.

A335905 Number of distinct integers encountered on all possible paths from n to any first encountered powers of 2 (that are excluded from the count), when using the transitions x -> x - (x/p) and x -> x + (x/p) in any order, where p is the largest prime dividing x.

Original entry on oeis.org

0, 0, 1, 0, 2, 1, 2, 0, 3, 2, 3, 1, 3, 2, 4, 0, 4, 3, 5, 2, 4, 3, 4, 1, 6, 3, 6, 2, 6, 4, 5, 0, 5, 4, 6, 3, 6, 5, 5, 2, 6, 4, 7, 3, 7, 4, 5, 1, 6, 6, 7, 3, 9, 6, 7, 2, 8, 6, 7, 4, 6, 5, 7, 0, 7, 5, 7, 4, 6, 6, 7, 3, 7, 6, 9, 5, 8, 5, 7, 2, 10, 6, 7, 4, 9, 7, 9, 3, 10, 7, 7, 4, 8, 5, 11, 1, 7, 6, 8, 6, 11, 7, 10, 3, 9
Offset: 1

Views

Author

Antti Karttunen, Jun 30 2020

Keywords

Examples

			From 9 one can reach with the transitions x -> A171462(x) (leftward arrow) and x -> A335876(x) (rightward arrow) the following three numbers, when one doesn't expand any power of 2 (in this case, 4, 8 and 16, that are not included in the count) further:
       9
      / \
     6   12
    / \ / \
  (4) (8) (16)
thus a(9) = 3.
From 10 one can reach with the transitions x -> A171462(x) and x -> A335876(x) the following two numbers (10 & 12), when one doesn't expand any powers of 2 (8 and 16 in this case, not counted) further:
  10
   |\
   | \
   | 12
   | /\
   |/  \
  (8)  (16)
thus a(10) = 2.
For n = 9, the numbers encountered are 6, 9, 12, thus a(9) = 3.
For n = 67, the numbers encountered are 48, 60, 66, 67, 68, 72, 96, thus a(67) = 7.
For n = 105, the numbers encountered are 48, 72, 90, 96, 105, 108, 120, 144, 192, thus a(105) = 9.
		

Crossrefs

Programs

  • PARI
    A171462(n) = if(1==n,0,(n-(n/vecmax(factor(n)[, 1]))));
    A335876(n) = if(1==n,2,(n+(n/vecmax(factor(n)[, 1]))));
    A209229(n) = (n && !bitand(n,n-1));
    A335905(n) = if(A209229(n),0,my(xs=Set([n]),allxs=xs,newxs,a,b,u); for(k=1,oo, newxs=Set([]); if(!#xs, return(#allxs)); allxs = setunion(allxs,xs); for(i=1,#xs,u = xs[i]; a = A171462(u); if(!A209229(a), newxs = setunion([a],newxs)); b = A335876(u); if(!A209229(b), newxs = setunion([b],newxs))); xs = newxs));

A336116 Primes of the form q*2^h - 1, where q is a Fermat prime.

Original entry on oeis.org

2, 5, 11, 19, 23, 47, 67, 79, 191, 271, 383, 1087, 1279, 4111, 5119, 6143, 16447, 20479, 81919, 262147, 263167, 786431, 1114111, 1310719, 16842751, 17825791, 1073758207, 4295032831, 4311744511, 17180131327, 21474836479, 51539607551, 824633720831, 1168231104511
Offset: 1

Views

Author

Antti Karttunen, Jul 09 2020

Keywords

Crossrefs

Programs

  • Mathematica
    NestList[NestWhile[NextPrime, #, ! (PrimeQ[#2] && With[{p = NestWhile[BitShiftRight, #2 + 1, EvenQ] - 1}, BitAnd[p, p - 1] == 0 && With[{b = BitLength[p]}, BitAnd[b - 1, b - 2] == 0]]) &, 2] &, 2, 25] (* Jan Mangaldan, Jul 14 2020 *)
  • PARI
    A000265(n) = (n>>valuation(n,2));
    isA019434(n) = ((n>2)&&isprime(n)&&!bitand(n-2,n-1));
    isA336116(n) = (isprime(n)&&isA019434(A000265(1+n)));

Formula

For all n >= 1, A335885(a(n)) <= 2.

Extensions

More terms from Jinyuan Wang, Jul 11 2020

A336117 Primes of the form q*2^h + 1, where q is a Mersenne prime.

Original entry on oeis.org

7, 13, 29, 97, 113, 193, 449, 509, 769, 7937, 12289, 114689, 520193, 786433, 7340033, 8388593, 33292289, 33550337, 469762049, 2130706433, 3221225473, 8588886017, 137438691329, 206158430209, 2199023254529, 6597069766657, 562948879679489, 7881299347898369, 9007182074871809
Offset: 1

Views

Author

Antti Karttunen, Jul 09 2020

Keywords

Crossrefs

Programs

Formula

For all n >= 1, A335885(a(n)) <= 2.

Extensions

More terms from Jinyuan Wang, Jul 11 2020

A336115 Semiprimes that are product of a Fermat prime and a Mersenne prime.

Original entry on oeis.org

9, 15, 21, 35, 51, 93, 119, 155, 381, 527, 635, 771, 1799, 2159, 7967, 24573, 32639, 40955, 139247, 196611, 393213, 458759, 655355, 1572861, 2031647, 2105087, 2228207, 2621435, 8323199, 8912879, 33685247, 134741759, 536813567, 6442450941, 8590000127, 10737418235
Offset: 1

Views

Author

Antti Karttunen, Jul 09 2020

Keywords

Comments

As 3 is both a Fermat prime and a Mersenne prime, A019434(1) * A000668(1) = 9 is also a term. It is the only square in this sequence.

Crossrefs

Programs

Formula

A335885(a(n)) = 2.

Extensions

Missing terms and more terms added by Jinyuan Wang, Jul 11 2020
Previous Showing 11-14 of 14 results.