计算n阶行列式的模板

最新推荐文章于 2022-09-04 09:00:20 发布

weixin_30556161

最新推荐文章于 2022-09-04 09:00:20 发布

阅读量297

点赞数

原文链接：http://www.cnblogs.com/mjy0724/p/4442373.html

版权

　　之前在学习计数问题的时候也在网上找了很多关于行列式的资料

　　但是发现很多地方都只介绍2\3阶的情况

　　一些论文介绍的方法又看不懂

　　然后就一直耽搁着

　　今天恰好出到这样的题目发现标算的代码简介明了还挺开心的

function det(a:arr):int64;
var i,j,k:longint;
    t,ans:int64;

begin
    ans:=1;
    for i:=1 to n do
    begin
        for j:=i+1 to n do
        while a[j,i]<>0 do
        begin
            t:=a[i,i] div a[j,i];
            for k:=i to n do a[i,k]:=(a[i,k]-t*a[j,k]) mod p;
            for k:=1 to n do swap(a[i,k],a[j,k]);
            ans:=-ans;
        end;
        if a[i,i] = 0 then exit(0);
        ans:=(ans * a[i,i]) mod p;
    end;
    while ans<0 do inc(ans,p);
    exit(ans);
end;