P4783-[模板]矩阵求逆

正题

题目链接:https://www.luogu.com.cn/problem/P4783


题目大意

给出一个矩阵,求它的逆矩阵。

1≤n≤4001\leq n\leq 4001n400


解题思路

记给出矩阵PPP,记单位矩阵EEE
P×P−1=E⇒P×(E×P−1)=EP\times P^{-1}=E\Rightarrow P\times (E\times P^{-1})=EP×P1=EP×(E×P1)=E

虽然看上去上面那个式子是废话,但是这是一个提示。

因为PPP进行初等变化变为EEE的过程中相当于乘上了一个P−1P^{-1}P1,而P−1×E=P−1P^{-1}\times E=P^{-1}P1×E=P1。所以如果我们拿一个EEEPPP做一样的初等变化就变为了P−1P^{-1}P1

写个高斯消元就好了,但是需要注意因为一般的消元会自动省略已经消掉的部分,但是因这里要处理P−1P^{-1}P1矩阵所以不能这么做。

时间复杂度O(n3)O(n^3)O(n3)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=410,P=1e9+7;
ll n,a[N][N],b[N][N];
ll power(ll x,ll b){ll ans=1;while(b){if(b&1)ans=ans*x%P;x=x*x%P;b>>=1;}return ans;
}
bool work(){for(ll i=1;i<=n;i++){for(ll j=i;j<=n;j++)if(a[j][i]){if(i!=j)swap(a[j],a[i]),swap(b[j],b[i]);break;}if(!a[i][i])return 1;ll inv=power(a[i][i],P-2);for(ll j=1;j<=n;j++)a[i][j]=a[i][j]*inv%P,b[i][j]=b[i][j]*inv%P;for(ll j=i+1;j<=n;j++){ll rate=P-a[j][i];for(ll k=1;k<=n;k++){(a[j][k]+=rate*a[i][k]%P)%=P;(b[j][k]+=rate*b[i][k]%P)%=P;}}}for(int i=n;i>=1;i--)for(int j=1;j<i;j++){for(int k=1;k<=n;k++)(b[j][k]+=P-a[j][i]*b[i][k]%P)%=P;a[j][i]=0;}return 0;
}
signed main()
{scanf("%lld",&n);for(ll i=1;i<=n;i++){for(ll j=1;j<=n;j++)scanf("%lld",&a[i][j]);b[i][i]=1;}if(work())return 0&puts("No Solution");for(ll i=1;i<=n;i++,putchar('\n'))for(ll j=1;j<=n;j++)printf("%lld ",b[i][j]);return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318435.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

.net core i上 K8S(六).netcore程序的service网络代理模式

正文上一章我们讲了pod的hostip模式&#xff0c;但在生产环境中&#xff0c;我们都是通过service来访问k8s集群的&#xff0c;service有两种模式来暴漏端口&#xff0c;今天我们来分享一下1.clusterIP模式我们在创建service的时候&#xff0c;默认创建的时clusterIP模式&#x…

CF183D-T-shirtx【dp,贪心】

正题 题目链接:https://www.luogu.com.cn/problem/CF183D 题目大意 nnn个人&#xff0c;mmm种衣服&#xff0c;给出每个人喜欢某件衣服的概率&#xff0c;你可以选择nnn件衣服带过去&#xff08;可以重复款式&#xff09;。求最大化能拿到喜欢衣服人的期望数量。 1≤n≤3000,1…

Jamie and Tree[CF916E]

Jamie and Tree[CF916E] 题意&#xff1a; 有一棵n个点的树&#xff0c;每个节点上有一个权值wi&#xff0c;最开始根为1号点&#xff0e;现在有3种 类型的操作&#xff1a; • 1 root, 表示将根设为root. • 2 u v x, 设u, v的最近公共祖先为p, 将p的子树中的所有点的权值加…

领域驱动设计,让程序员心中有码(三)

“正如西方古典哲学在现代社会逐渐式微&#xff0c;成为少数内心丰满者们填充自己精神世界的宝贵食物&#xff0c;UML也这样&#xff1b;互联网技术飞速发展的今天&#xff0c;各类软件设计思想层出不穷&#xff0c;正是站在UML和其他各种软件基础理论巨人的肩膀上&#xff0c;…

P4127 [AHOI2009]同类分布(数位dp)

洛谷传送门 文章目录题目描述解析代码题目描述 给出两个数a,b求出[a,b]中各位数字之和能整除原数的数的个数。 1<a<b<1018 解析 容易想到数位dp 但本题的难点是如果只记录数位和sum与取模的结果res&#xff0c;因为取模的除数发生改变&#xff0c;难以转移 如何解决…

线性代数(矩阵、高斯、线性基……)

矩阵 矩阵加法&#xff1a; 相同位置相加。 矩阵乘法&#xff1a; 满足分配率、结合律&#xff0c;不满足交换律(矩阵与逆矩阵之间除外) 。 矩阵转置&#xff1a; 记矩阵为 \(A\) &#xff0c;则 \(A\) 的转置记为 \(A^T\) 。 性质&#xff1a; \[{(A^T)}^TA \]\[{(AB)}^TA^TB^…

Tavan

Tavan–简单的k进制题 呵呵真的简单 也就是考试没做出来罢了&#xff0c;不慌不慌&#xff0c;奶我一口&#xff0c;还能活 【题目摘要】 题目描述 小 Zeljko 一直在阁楼里读他奶奶的旧信&#xff0c;并且发现了一个长度为 N 的单词。 不幸的是&#xff0c;由于溢出的墨水&…

.NET Core 2.x中使用Named Options处理多个强类型配置实例

来源&#xff1a; Using multiple instances of strongly-typed settings with named options in .NET Core 2.x作者&#xff1a; Andrew Lock译者&#xff1a; Lamond Lu.NET Core从1.0版本开始&#xff0c;就已经开始使用Options模式绑定强类型配置对象。从那时起到现在&…

CF235D-Graph Game【LCA,数学期望】

正题 题目链接:https://www.luogu.com.cn/problem/CF235D 题目大意 给出一棵基环树&#xff0c;每次随机选择一个点让权值加上这个点的连通块大小然后删掉这个点。 求删光所有点时期望权值。 1≤n≤30001\leq n\leq 30001≤n≤3000 解题思路 先找到环&#xff0c;然后考虑暴…

YBTOJ:魔法数字(数位dp)

文章目录题目描述解析题目描述 解析 迷惑。。。 首先&#xff0c;比较容易想到用二进制状态压缩记录1-9是否在十进制中出现过 然后就是整除的问题 如果记录余数&#xff0c;它的模数又有9个 开九维余数直接爆炸。。。 怎么办嘞&#xff1f; 有一个结论&#xff1a; 若&#xf…

斯特林数

第一类斯特林数 咕咕咕 第二类斯特林数 定义&#xff1a;把 \(n\) 个不同的球放入 \(r\) 个相同的盒子的方案数(盒子不能为空&#xff0c;记为&#xff1a;\(S(n,r)\) 或 \(\begin{Bmatrix}n\\r\end{Bmatrix}\) 。 递推式&#xff1a; \[\begin{Bmatrix}n\\r\end{Bmatrix}\begi…

Prosjecni(构造)

Prosjecni 【题目摘要】 描述 Slavko很无聊&#xff0c;所以他把正整数填到N*N的方阵中。 如果他填出来的方阵满足以下条件&#xff0c;他会特别高兴&#xff1a; ●每行中的数字的平均值是一个位于同一行的整数。 ●每列中的数字的平均值是一个位于同一列的整数。 ●表中的所…

撒花!中文翻译仓库链接已加入 ML.NET 官方示例网站首页

从2018年12月02日决定开始做ML.NET 示例中文版https://github.com/feiyun0112/machinelearning-samples.zh-cn&#xff0c;然后以每天一篇的速度进行翻译&#xff0c;总共耗时15天&#xff0c;将现有的官方实例全部翻译成了中文&#xff0c;并提交了添加中文链接PR&#xff0c;…

P3803 【模板】多项式乘法(FFT)

P3803 【模板】多项式乘法&#xff08;FFT&#xff09; 题目描述 给定一个 n 次多项式 F(x)&#xff0c;和一个 m 次多项式 G(x)。 请求出 F(x)和 G(x)的卷积。 从低到高输出F(x)*G(x)的系数 另一种问法&#xff1a; 如果有两个无限序列a和b&#xff0c;那么它们卷积的结果是…

CF666E-Forensic Examination【广义SAM,线段树合并】

正题 题目链接:https://www.luogu.com.cn/problem/CF666E 解题思路 给出一个串SSS和nnn个串TiT_iTi​。mmm次询问Sa∼bS_{a\sim b}Sa∼b​在Tl∼rT_{l\sim r}Tl∼r​中出现的最多次数并且输出这个串的编号。 1≤∣s∣≤5105,∑Ti≤5104,1≤m≤51051\leq |s|\leq 5\times 10^5…

简简单单组合数学

简简单单组合数学 P3158 [CQOI2011]放棋子 \(\uparrow\) 假组合数学&#xff0c;真 \(\text{DP}\) 。 \(f[i][j][k]\) : 用了 \(i\) 行 \(j\) 列&#xff0c;涂了前 \(k\) 种颜色的方案数。 \(g[i][j][k]\) : 用了 \(i\) 行 \(j\) 列&#xff0c;涂了第 \(k\) 种颜色的方案数(用…

YBTOJ:单词频率(AC自动机)

解析 我对力量一无所知 通过这题&#xff0c;可以看出我对AC自动机还是完全没有理解 qwq 首先容易想到&#xff1a; 建一课trie树&#xff0c;然后建树时记录每个串s的终点&#xff0c;这个点后面每被经过一次&#xff0c;就相当于出现一次该单词s 但是&#xff0c;这种“出现”…

一元一次方程

一元一次方程–逆波兰栈 【题目摘要】 题目描述 SLON是一个调皮的学生&#xff0c;为了让他静下心来&#xff0c;老师给他出了一道数学题&#xff1a; 给定表达式A&#xff0c;A中含有变量x和,-,*,(,)这些符号&#xff0c;括号成对出现&#xff0c;一个算术运算符均对应两个操…

.net core上 K8S(七).netcore程序的服务发现

正文上一章我们分享了k8s的网络代理模式&#xff0c;今天我们来分享一下k8s中的服务发现。1.环境变量模式的服务发现k8s默认为我们提供了通过环境变量来实现服务发现的功能&#xff0c;前提是1.需要service在pod之前创建2.适用于同一命名空间1.1创建servicekubectl create -f n…

质数与合数

质数与合数 题意&#xff1a; FFF和GGG正在玩一个质数与合数的游戏 一开始有N个石头 FFF和GGG轮流对这堆石头进行操作&#xff0c;FFF每次选择1到K之间的一个数x&#xff0c;并拿走x个石头&#xff0c;拿走之后剩下的石头数量必须是质数 接着GGG进行同样的操作&#xff0c;但…