西南交大峨眉校区INT杯2020-A.矩阵【行列式计算 + 快速幂】

星空皓月

于 2020-12-15 01:02:14 发布

阅读量305

点赞数

分类专栏： # 数学比赛

本文链接：https://blog.csdn.net/qq_45249273/article/details/111189572

版权

比赛同时被 2 个专栏收录

34 篇文章 3 订阅

订阅专栏

数学

4 篇文章 0 订阅

订阅专栏

这篇博客探讨了矩阵行列式的计算规律，通过手动计算1阶到3阶行列式，发现了一个公式：D=(a1*b1+a2*b2+…+an*bn)*xn-1+xn。然后利用快速幂算法实现高效计算，代码中展示了C++实现的快速幂求解行列式的方法。

摘要由CSDN通过智能技术生成

在这里插入图片描述

思路：

其实这个矩阵的行列式是有规律的，将1阶，2阶，3阶手动计算出来过后。
我们会发现一个规律，D = (a1 * b1 + a2 * b2+ … + an * bn) * x^n-1 + xⁿ.

AC代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MOD = 1e9 + 7;
const int MAXN = 1e5 + 5;
ll a[MAXN], b[MAXN];
ll fast_power(ll a, ll b) {
	ll ans=1;
	while (b > 0) {
		if(b & 1) {
           ans = ans * a % MOD;
		}
		a = a * a % MOD;
		b >>= 1;
	}
	return ans;
}
void solve() {
    ll n, x;
    while(~scanf("%lld%lld", &n, &x)) {
        for (int i = 0; i < n; ++i) {
            scanf("%lld", &a[i]);
        }
        for (int i = 0; i < n; ++i) {
            scanf("%lld", &b[i]);
        }
        ll ans = 0;
        for (int i = 0; i < n; ++i) {
            ans = (ans + a[i] * b[i] % MOD) % MOD;
        }
        ans = (ans * fast_power(x, n - 1) % MOD + fast_power(x, n) % MOD) % MOD;
        printf("%lld\n", ans);
    }
}
int main() {
    solve();
    return 0;
}