Paysages Mathématiques · @paysmaths
135 followers · 619 posts · Server mathstodon.xyz

"Le raisonnement par récurrence [...] est un instrument qui permet de passer du fini à l’infini [...]" – Henri Poincaré (1854-1912)

#math #maths #recurrence #mathematiques #citation

Last updated 1 year ago

Paysages Mathématiques · @paysmaths
135 followers · 619 posts · Server mathstodon.xyz

"Le raisonnement par récurrence [...] est un instrument qui permet de passer du fini à l’infini [...]" – Henri Poincaré (1854-1912)

#math #maths #recurrence #mathematiques #citation

Last updated 1 year ago

Poetry News · @haikubot
600 followers · 4693 posts · Server mastodon.cloud
Paysages Mathématiques · @paysmaths
39 followers · 65 posts · Server mathstodon.xyz

"Le caractère essentiel du raisonnement par récurrence c’est qu’il contient, condensés pour ainsi dire en une formule unique, une infinité de syllogismes." – Henri Poincaré (1854-1912)

#recurrence #mathematiques #citation

Last updated 2 years ago

Matt McIrvin · @mattmcirvin
94 followers · 150 posts · Server mathstodon.xyz

So, the first four pseudoprimes generated by the sequence
a(n) = a(n-1) - a(n-2) - a(n-4), starting with (4, 1, -1, -2),
(that is, composite numbers passing the probable-prime test a(p) ≡ 1 mod p)
are, unless I have slipped up somewhere, 2706, 5205530, 42161779, 1146180898. That's a big gap before number four!

#recurrence

Last updated 2 years ago

Matt McIrvin · @mattmcirvin
91 followers · 136 posts · Server mathstodon.xyz

Here's a new one of these prime-detecting sequences for you:

a(n) = a(n-1) - a(n-2) - a(n-4)
Start with 4, 1, -1, -2.
For prime p, a(p) ≡ 1 (mod p).

My search in 2013 passed over this because the first pseudoprime is just 2706, but the second one is 5205530 and the third is > 40 million.

The other good fourth-order ones with coefficients in {+-1, 0} all seem to be basically the Perrin case again (characteristic polynomial is a multiple of the Perrin one).

#recurrence

Last updated 2 years ago

Matt McIrvin · @mattmcirvin
91 followers · 136 posts · Server mathstodon.xyz

@johncarlosbaez @ColinTheMathmo These methods could definitely be applied to the whole family of sequences I've been messing with. The fastest method I knew of was Joerg Arndt's approach of calculating large powers of matrices over Z_n, since we only need the moduli--PARI/GP is way better at this than sagemath--this is pretty good but does get harder as the numbers get truly immense.

#recurrence

Last updated 2 years ago

Matt McIrvin · @mattmcirvin
46 followers · 106 posts · Server mathstodon.xyz

But the analysis having to do with PV numbers doesn't apply to Lehmer's number, since it is a "Salem number": most of the other roots are ON the unit circle and will keep wreaking havoc with the moduli arbitrarily far out. So I understand that less. It is conjectured to be the smallest Salem number.

#recurrence

Last updated 2 years ago

Matt McIrvin · @mattmcirvin
46 followers · 105 posts · Server mathstodon.xyz

Anyway, what got me started on all this again was that a post here that I can no longer find mentioned "Lehmer's number", the leading root of x^10+x^9-x^7-x^6-x^5-x^4-x^3+x+1. Naturally I wondered if the recurrence corresponding to this polynomial was a good prime detector. It is: the only pseudoprime less than, I think, 4 million is 2014. But it turns out that 10th-order recursions that good are not all that hard to find, so it isn't super special in that respect.

#recurrence

Last updated 2 years ago

Matt McIrvin · @mattmcirvin
46 followers · 99 posts · Server mathstodon.xyz

So I guess all of this gives us a generalization of the Pisot-Vijayaraghavan numbers beyond the positive reals, and I'm sure someone has studied this already, but I haven't seen it.

#recurrence

Last updated 2 years ago

Matt McIrvin · @mattmcirvin
46 followers · 98 posts · Server mathstodon.xyz

I also found a lot where the leading roots are a conjugate pair of complex numbers, and there, if you add their powers, you get this sort of exponentially exploding sinusoidal oscillation--but (if the other roots have norm < 1) one where all the values get exponentially close to integers! That's real weird.

#recurrence

Last updated 2 years ago

Matt McIrvin · @mattmcirvin
46 followers · 97 posts · Server mathstodon.xyz

It is, by the way, profoundly strange to see this exact expression where we raise a bunch of generally irrational and/or complex numbers to large powers, add them together, and it turns out to just be *rounding* the leading term to an integer.

But that's a remarkable property of these PV numbers--they get close to integers when you raise them to large powers. But the PV numbers are just the positive real ones--for my fifth-order recurrence, the leading root is negative.

#recurrence

Last updated 2 years ago

Matt McIrvin · @mattmcirvin
46 followers · 96 posts · Server mathstodon.xyz

Anyway, if you look at that Binet-like formula, you can see that the sum should pretty quickly be dominated by the root(s) of largest norm.

But we're doing modulo arithmetic on the sums, so the smaller terms will still be important--but *if* they all have norm less than 1, pretty soon all they're doing is rounding the sum to the nearest integer. But early on, they can be doing more than that. I think that is what causes some of the small-number weirdness.

#recurrence

Last updated 2 years ago

Matt McIrvin · @mattmcirvin
46 followers · 95 posts · Server mathstodon.xyz

But as yet, I have no good feel for what makes a "good" sequence whose probable-prime test passes few composites, preferably with a few small-integer coefficients in the recurrence.

The polynomials associated with the Lucas and Perrin numbers are famous: their leading eigenvalues are the golden ratio and the "plastic number" respectively, small Pisot-Vijayaraghavan numbers as I was saying. But this doesn't lead to a simple pattern for "good" sequences.

#recurrence

Last updated 2 years ago

Matt McIrvin · @mattmcirvin
46 followers · 94 posts · Server mathstodon.xyz

Some other relevant facts: Sequences of this type all have a simple closed-form "Binet-like" exact expression for their elements:

\[ a_n = \sum_{j} c_j^n \]

where the c_j are the eigenvalues of the matrix M. Those are also the roots of a monic polynomial that comes from det (Ix - M) = 0. Basically the leading coefficient is always 1 and the others are minus the recurrence coefficients. So these recurrences seem to have a lot to do with the study of algebraic integers.

#recurrence

Last updated 2 years ago

· @xpil
606 followers · 4865 posts · Server fosstodon.org

You don’t understand “You won’t understand until you’re older” until you’re older.

#theshowerthoughts #recurrence

Last updated 3 years ago