码迷,mamicode.com
首页 > 其他好文 > 详细

计算n阶行列式的模板

时间:2015-04-20 20:29:39      阅读:142      评论:0      收藏:0      [点我收藏+]

标签:

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

  但是发现很多地方都只介绍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;

 

计算n阶行列式的模板

标签:

原文地址:http://www.cnblogs.com/mjy0724/p/4442373.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!