麻豆小视频在线观看_中文黄色一级片_久久久成人精品_成片免费观看视频大全_午夜精品久久久久久久99热浪潮_成人一区二区三区四区

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

矩陣快速冪 非詳解

2019-11-14 09:34:20
字體:
來源:轉載
供稿:網友
#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 乘以本身。


發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 精品一区二区久久久久久按摩 | 精品国产乱码久久久久久久久 | 免费在线性爱视频 | 黄色片网页 | 久久久久久久.comav | 久久国产精品久久精品国产演员表 | 欧洲色阁中文字幕 | 日韩精品久久久久久久电影99爱 | 爱看久久| 国产一区日韩精品 | 中文字幕在线观看视频一区 | 精品黑人一区二区三区国语馆 | 九色中文 | 免费毛片观看 | 久久国产精品小视频 | 91九色精品 | h久久| 国产91丝袜在线播放 | 欧美中文字幕在线视频 | 私库av在线免费观看 | 国产乱淫av片免费 | 日美av在线| 牛牛热这里只有精品 | 一本视频在线观看 | 久久精品一区二区三区四区五区 | 成年免费视频黄网站在线观看 | 欧美一级不卡视频 | 成人在线免费观看视频 | 欧美另类在线视频 | 亚洲成人免费视频在线 | 欧美日韩在线看片 | 欧美爱爱视频免费看 | 九九精品免费 | 日韩视频在线观看免费视频 | av影片在线观看 | 亚洲性在线视频 | 久久精品国产99久久6动漫亮点 | lutube成人福利在线观看 | 天天夜夜操操 | 天天草夜夜| 18一20岁一级毛片 |