国产探花免费观看_亚洲丰满少妇自慰呻吟_97日韩有码在线_资源在线日韩欧美_一区二区精品毛片,辰东完美世界有声小说,欢乐颂第一季,yy玄幻小说排行榜完本

首頁 > 學院 > 開發設計 > 正文

矩陣快速冪 非詳解

2019-11-14 08:59:18
字體：大 中 小
來源：轉載
供稿：網友
#include <cstdio>#include <cstring>int n, k;const int mod = 9973;struct matrix{ int tr[10][10]; matrix Operator * (const matrix &a) const{//重載運算符 matrix tmp; memset(tmp.tr, 0, sizeof(tmp.tr)); for(int i = 0; i < n; i++) for(int j = 0; j < n; j++){ for(int k = 0; k < n; k++) tmp.tr[i][j] += tr[i][k] * a.tr[k][j]; tmp.tr[i][j] %= mod;//這里要取模，不然可能溢出 } return tmp; }}ans, ori;matrix pow_mod(int k){ for(int i = 0; i < n; i++) ans.tr[i][i] = 1;//化為單位矩陣 while(k){ if(k&1) ans = ans*ori;//不能寫成ans *= ori，因為沒有重載*=運算符 k >>= 1; ori = ori*ori; }//核心代碼，下面有解釋}int main(){ int t; scanf("%d", &t); while(t--){ scanf("%d%d", &n, &k); memset(ans.tr, 0, sizeof(ans.tr)); memset(ori.tr, 0, sizeof(ori.tr)); for(int i = 0; i < n; i++) for(int j = 0; j < n; j++) scanf("%d", &ori.tr[i][j]); pow_mod(k); long long res = 0; for(int i = 0; i < n; i++){ res += ans.tr[i][i] % mod; } 核心代碼：

while(k){ if(k&1) ans = ans*ori; k >>= 1; ori = ori*ori; }

假設 k = 89，其二進制為 1011001， 顯然 a^k = ( a^1 ) * ( a^8 ) * ( a^16 ) * ( a^64 ); 也就是，k的二進制位為0時，可以跳過（右移） 。

每次判斷k的最后一位二進制位，若為1，則 ans = ans*ori ， 然后k右移一位，ori 乘以本身。

發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板： 磴口县| 东乌珠穆沁旗| 米林县| 池州市| 镇巴县| 平远县| 尚志市| 应用必备| 湟源县| 邵东县| 浦江县| 邛崃市| 磐石市| 虹口区| 宁武县| 鸡东县| 长汀县| 诏安县| 永嘉县| 柳州市| 太湖县| 金塔县| 常德市| 阿尔山市| 墨玉县| 龙泉市| 平武县| 大渡口区| 广宗县| 建宁县| 项城市| 保山市| 三台县| 邳州市| 准格尔旗| 朝阳区| 宜宾市| 石林| 元朗区| 五原县| 崇文区|