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.
%I A295607 #48 Mar 31 2024 00:43:03 %S A295607 85,49,133,81,363,705,881,1023,417,653,773,1229,1985,273,275,585,1365, %T A295607 2505,3861,129,289,719,2069,2393,3113,4609,5549,5555,5789,6299,7517, %U A295607 7649,321,2321,2337,3567,6617,6993,9377,12957,13737,14505,15033,15225,15237,385,2177,2565,7097,8273,8897 %N A295607 a(n) = A001567(n) - 2^floor(log_2(A001567(n))). %C A295607 This sequence contains the distances from pseudoprime numbers (A001567) to the next smaller number of the form 2^n. Conjecture: It seems that these distances do not take all possible values. So, if we know that a certain distance does not appear with pseudoprime numbers, we are able to calculate these numbers using Fermat's little theorem and we know for sure that they are primes. %H A295607 Amiram Eldar, <a href="/A295607/b295607.txt">Table of n, a(n) for n = 1..10000</a> %H A295607 Jonas Kaiser, <a href="https://arxiv.org/abs/1608.00862">On the relationship between the Collatz conjecture and Mersenne prime numbers</a>, arXiv:1608.00862 [math.GM], 2016. %e A295607 There are no pseudoprimes detected by Fermat's little theorem for 2^k+m with m = {3,5,7,...,47} up to k = 10000 (checked using the Pari function ispseudoprime(k)). When this sequence is ordered for the first 10^5 pseudoprimes, the following first terms (up to 1000) result: 1, 49, 81, 85, 129, 133, 273, 275, 289, 321, 363, 385, 417, 585, 653, 705, 719, 773, 881. %t A295607 Map[# - 2^Floor@ Log2@ # &, Select[Range[3, 10^5, 2], And[! PrimeQ[#], PowerMod[2, (# - 1), #] == 1] &]] (* _Michael De Vlieger_, Nov 26 2017 *) %o A295607 (PARI) a(A001567)=A001567-2^(floor(log(A001567)/log(2))) \\ %Y A295607 Cf. A001567. %K A295607 nonn %O A295607 1,1 %A A295607 _Jonas Kaiser_, Nov 24 2017