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.

A248219 Indices where A248218(n) = 1.

This page as a plain text file.
%I A248219 #17 Mar 22 2017 03:39:25
%S A248219 1,3,7,19,21,31,37,43,57,93,111,129,133,157,217,259,301,307,313,399,
%T A248219 463,471,499,589,613,651,703,777,817,903,921,939,967,1099,1147,1333,
%U A248219 1389,1497,1591,1767,1839,2109,2149,2191,2451,2683,2901,2983,3241,3297,3441
%N A248219 Indices where A248218(n) = 1.
%C A248219 n is in the sequence if there exists a number k0 such that A003095(k) mod n = A003095(k0) mod n for all k >= k0.
%C A248219 Conjecture: a(n) ~ n^2.
%H A248219 Vaclav Kotesovec, <a href="/A248219/b248219.txt">Table of n, a(n) for n = 1..717</a>
%e A248219 7 is in the sequence since A003095 mod 7 is {1, 2, 5, 5, 5, 5, ...}, k0 = 3, A003095(k) mod 7 = 5 for all k >= k0. The period is 1.
%e A248219 10 is not in the sequence since A003095 mod 10 is {1, 2, 5, 6, 7, 0, 1, 2, 5, 6, 7, 0, 1, 2, ...} and the period is 6, not 1.
%t A248219 nmax=10000; periods=Table[m=Rest[NestList[Mod[#^2+1,n]&,0,nmax]]; period=0; j=1; While[j<=Length[m]&&period==0,If[m[[Length[m]-j]]==m[[Length[m]]],period=j]; j++]; period,{n,1,nmax}]; Select[Range[nmax],periods[[#]]==1&]
%Y A248219 Cf. A003095, A247981, A248218.
%K A248219 nonn
%O A248219 1,2
%A A248219 _Vaclav Kotesovec_, Oct 04 2014