W dniu 2013-05-13 22:03, Miroslaw Kwasniak pisze:
> bartekltg <b...@g...com> wrote:
>> Hehe, spryciarz;> A dla 39x39?
>
> No, tak 2^32 rozwiązań nie przejrzę :O
>
> ? matrank(matker(m))
> time = 84 ms.
> %16 = 32

To wspomniane pari-gp?

Ile wychodzi dla większych? powiedzmy 55 czy 64.


> Ale przypadkowe poniżej 1s
>
> ? m=Mod(M(39,39),2); l=m;l=l[,1]+l[,39]+l[,7]; s=matinverseimage(m,l);m*s==l
> time = 889 ms.
> %15 = 1
>
> Z czego większość zajmuje kiepska procedura budowy m:

Spróbuj może tak:
N = macierz pasmowa n x n o paśmie [1,0.5,1]
(na diagonali 0.5. Nad i pod diagonalą 1).
M = to samo, ale rozmiaru m x m

Teraz poszukiwana macierz to

A = M (*) Id[n] + Id[m] (*) N

gdzie (*) to iloczyn kroneckera, pewnie jest zaimplementowany.

Mathematice zajmuje to 1.2s, ale matlabowi (w bebechach pewnie
bliższe pari) wyszło 0.007236s:)

pzdr
bartekltg