10.12.2012 Views

Prime Numbers

Prime Numbers

Prime Numbers

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

9.8 Research problems 539<br />

9.82. A certain prime suggested in [Craig-Wood 1998], namely<br />

p =2 64 − 2 32 +1,<br />

has advantageous properties in regard to CRT-based convolution. Investigate<br />

some of these advantages, for example by stating the possible signal lengths<br />

for number-theoretical transforms modulo p, exhibiting a small-magnitude<br />

element of order 64 (such elements might figure well into certain FFT<br />

structures), and so on.<br />

9.83. Here is a surprising result: Length-8 cyclic convolution modulo a<br />

Mersenne prime can be done via only eleven multiplies. It is surprising because<br />

the Winograd bound would be 2 · 8 − 4 = 12 multiplies, as in Exercise 9.39.<br />

Of course, the resolution of this paradox is that the Mersenne mod changes<br />

the problem slightly.<br />

To reveal the phenomenon, first establish the existence of an 8-th<br />

root of unity in F p 2,withp being a Mersenne prime and the root being<br />

symbolically simple enough that DGTs can be performed without explicit<br />

integer multiplications. Then consider the length-8 DGT, used to cyclically<br />

convolve two integer signals x, y. Next, argue that the transforms X, Y have<br />

sufficient symmetry that the dyadic product X ∗Y requires two real multiplies<br />

and three complex multiplies. This is the requisite count of 11 muls.<br />

An open question is: Are there similar “violations” of the Winograd bound<br />

for lengths greater than eight?<br />

9.84. Study the interesting observations of [Yagle 1995], who notes that<br />

matrix multiplication involving n×n matrices can be effected via a convolution<br />

of length n 3 . This is not especially surprising, since we cannot do an<br />

arbitrary length-n convolution faster than O(n ln n) operations. However,<br />

Yagle saw that the indicated convolution is sparse, and this leads to interesting<br />

developments, touching, even, on number-theoretical transforms.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!