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

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

矩陣快速冪 非詳解

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


發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: h网站在线观看 | 日韩精品中文字幕在线播放 | 日本网站在线看 | 视频一区二区三区免费观看 | 欧美综合在线观看视频 | 欧美一级性 | 日日草夜夜操 | 日本aaaa片毛片免费观看视频 | 日韩视频一二三 | 91精品国产乱码久久桃 | 亚洲国产高清自拍 | 日本高清com | 日韩精品中文字幕一区二区三区 | 可以看逼的视频 | 日美黄色片 | 毛片网站视频 | a视频在线免费观看 | 免费黄色在线观看网站 | 羞羞的网址 | 日韩黄色片免费看 | 98国内自拍在线视频 | 亚洲精品永久视频 | 蜜桃视频在线免费播放 | 国产福利视频在线观看 | 狠狠操夜夜爱 | 3344永久免费| 国产精品刺激对白麻豆99 | 爽爽淫人网 | 国产精品久久久久久久娇妻 | 欧美a级在线免费观看 | 国产精品视频免费网站 | 天天曰夜夜操 | 草莓福利社区在线 | 成人aaaaa片毛片按摩 | 黄网站色成年大片免费高 | 久久久久久久一区二区三区 | 九九热精彩视频 | 嫩呦国产一区二区三区av | 黄色大片网 | a网站在线| 国产男女爽爽爽爽爽免费视频 |