[seqfan] Re: A058055

[D. S. McNeil]

On Fri, Aug 10, 2012 at 6:31 PM,  <franktaw at netscape.net> wrote:
> Looking at A058056, it appears A058055(n) might be the least m such that
> m^2+n is the *first* prime greater than m^2.  Someone should check this.

I think you're right:

sage: old = [1, 3, 8, 5, 12, 11, 18, 51, 82, 49, 234, 23, 42, 75, 86,
231, 174, 107, 288, 63, 80, 69, 102, 325, 166, 765, 128, 143, 822,
727, 276, 597, 226, 835, 702, 461, 254, 693, 592, 797, 1284, 349, 370,
2337, 596, 645, 3012, 1033, 590, 4083, 1490, 757, 882, 833]
sage: new = [next(m for m in PositiveIntegers() if (m^2+n) ==
next_prime(m^2)) for n in [1..len(h)]]
sage: old == new
True