（矩阵行列式求模）Find The Determinant III

最新推荐文章于 2019-03-14 12:31:21 发布

原创

最新推荐文章于 2019-03-14 12:31:21 发布 · 2.8k 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#acm

该博客探讨了如何计算矩阵行列式的模，特别是在模不是质数时的处理方法。通过分类讨论，当模为质数时应用费马小定理，否则采用特定公式进行计算。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

https://www.spoj.com/problems/DETER3/en/

题意即为输出矩阵的行列式求模的答案
分类讨论模为不为质数的情况即可，若为质数，可以使用费马小定理；否则利用公式a/b%m=(a%(b*m))/b

#include <iostream>
#include <algorithm>
#include <cstdlib>
#include <cstring>
#include <vector>
#include <string>
#include <cmath>
using namespace std;
typedef long long ll;

const int maxn = 200 + 3;
ll a[maxn][maxn], p;
int isprime;
ll qpow(ll a, ll b) {
    ll ret = 1;
    while(b) {
        if(b & 1) ret = ret * a % p;
        b >>= 1;
        a = a * a % p;
    }
    return ret;
}
ll Dot_is_prime(int n) {
    ll ret = 1;
    int flg = 0;
    for(int i = 1; i <= n; ++i) {
        for(int j = i + 1; j <= n; ++j) {
            while(a[j][i]) {
                ll t = a[i][i] * qpow(a[j][i], p - 2) % p;
                for(int k = 1; k <= n; ++k) 
                    a[i][k] = (a[i][k