2021“MINIEYE杯”中国大学生算法设计超级联赛（5）Random Walk 2（推式子+矩阵逆+矩阵乘）

#include<bits/stdc++.h>
using namespace std;
using ll=long long;
const ll mod=998244353;
const int N=310;
int W[N][N];
ll P[N][N],P0[N][N];
ll I[N][N],A[N][N];
int n;
ll qmi(ll a,ll b)
{ll v=1;while(b){if(b&1) v=v*a%mod;b>>=1;a=a*a%mod;}return v;
}
void Inv()
{for(int c=1,r=1;c<=n;c++){int t=r;for(int i=c;i<=n;i++)if(abs(P[i][c])>abs(P[t][c])) t=i;if(P[t][c]==0) continue;for(int j=1;j<=n;j++) {swap(I[t][j],I[r][j]);swap(P[t][j],P[r][j]);}ll inv=qmi(P[r][c],mod-2);for(int j=1;j<=n;j++){I[r][j]=I[r][j]*inv%mod;P[r][j]=P[r][j]*inv%mod;}for(int i=1;i<=n;i++)if(i!=r&&P[i][c]!=0){ll v=P[i][c];for(int j=1;j<=n;j++){I[i][j]=(I[i][j]-v*I[r][j]%mod+mod)%mod;P[i][j]=(P[i][j]-v*P[r][j]%mod+mod)%mod;}}r++;}}
void Mul() // I-1 × P0
{for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) A[i][j]=0;for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)   for(int k=1;k<=n;k++) A[i][j]=(A[i][j]+I[i][k]*P0[k][j]%mod)%mod;
}
int main()
{ios::sync_with_stdio(false);cin.tie(nullptr);cout.tie(nullptr);int Tc=1;cin>>Tc;while(Tc--){cin>>n;for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) cin>>W[i][j];for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) P[i][j]=P0[i][j]=0,I[i][j]=(i==j);for(int i=1;i<=n;i++){ll tot=0;for(int j=1;j<=n;j++) tot+=W[i][j];ll inv=qmi(tot,mod-2);for(int j=1;j<=n;j++){if(i!=j) P[i][j]=1ll*W[i][j]*inv%mod;else P0[i][j]=1ll*W[i][j]*inv%mod;}for(int j=1;j<=n;j++){P[i][j]=mod-P[i][j];if(i==j) P[i][j]=1;}}Inv();Mul();for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) cout<<A[i][j]<<" \n"[j==n];}return 0;
}

2021“MINIEYE杯”中国大学生算法设计超级联赛（5）Random Walk 2（推式子+矩阵逆+矩阵乘）相关推荐

【2021杭电多校赛】2021“MINIEYE杯”中国大学生算法设计超级联赛（3）签到题3题
2021"MINIEYE杯"中国大学生算法设计超级联赛(3) Start Time : 2021-07-27 12:00:00 End Time : 2021-07-27 17:0 ...
【2021杭电多校赛】2021“MINIEYE杯”中国大学生算法设计超级联赛（1）签到题15869
2021"MINIEYE杯"中国大学生算法设计超级联赛(1) Start Time : 2021-07-20 12:10:00 End Time : 2021-07-20 17:1 ...
2021“MINIEYE杯”中国大学生算法设计超级联赛
2021"MINIEYE杯"中国大学生算法设计超级联赛 1006 Given a sequence of integers of length n, find the shorte ...
2021“MINIEYE杯”中国大学生算法设计超级联赛（2）
2021"MINIEYE杯"中国大学生算法设计超级联赛(2) 1008 I love exam (类背包DP) 1010 I love permutation (数学构造,剩余系) ...
7068 Dota2 Pro Circuit 杭电多校(2021“MINIEYE杯”中国大学生算法设计超级联赛9) [贪心+双指针]
题目 Dota2 Pro Circuit *Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Oth ...
【2021杭电多校赛】2021“MINIEYE杯”中国大学生算法设计超级联赛（5）签到题4题
Solved Pro.ID Title Ratio(Accepted / Submitted) 1001 Miserable Faith 33.33%(19/57) 1002 String Mod 2 ...
【2021杭电多校赛】2021“MINIEYE杯”中国大学生算法设计超级联赛（8）
1006 GCD Game (博弈论,nim游戏,质因子个数) 题意:有n个数a1,a2....an两个人玩游戏,Alice先动,Bob后动每次可以将一个数变为他的因子,直到不能动就输了. 题解:将问 ...
2021“MINIEYE杯”中国大学生算法设计超级联赛（1）个人解题报告
文章目录 HDU6950 Mod, Or and Everything HDU6954 Minimum spanning tree HDU6958 KD-Graph HDU6957 Maximal s ...
2021“MINIEYE杯”中国大学生算法设计超级联赛(8)（1002,1004,1006,1009）
前言依旧是白嫖账号,只打了一些题/kk 正题 1002 Buying Snacks 题目大意 nnn个物品,每个可以买一次也可以不买,如果买需要选择1/21/21/2块钱的,然后也可以相邻两个一起买 ...
2021“MINIEYE杯”中国大学生算法设计超级联赛(7)部分题解
前言找大佬嫖到个号来划水打比赛了,有的题没写或者不是我写的就不放了. 目前只有:1004,1005,1007,1008,1011 正题题目链接:https://acm.hdu.edu.cn/con ...

2021“MINIEYE杯”中国大学生算法设计超级联赛（5）Random Walk 2（推式子+矩阵逆+矩阵乘）

Random Walk 2

2021“MINIEYE杯”中国大学生算法设计超级联赛（5）Random Walk 2（推式子+矩阵逆+矩阵乘）相关推荐

最新文章

热门文章