矩阵幂(矩阵相乘)

矩阵幂(矩阵相乘)题目描述给定一个 n n 的矩阵 求该矩阵的 k 次幂 即 P k 输入描述 第一行 两个整数 n 2 lt n lt 10 k 1 lt k lt 5 两个数字之间用一个空格隔开 含义如上所示 接下来有 n 行 每行 n 个正整数 其中 第 i 行第 j 个整数表示矩阵中第 i 行第 j 列的矩阵元素 Pij 且 0 lt Pij lt 10 另外 数据保证最后结果不会超过 10 8 输出描述 对于每组测试数据 输

题目描述

给定一个n*n的矩阵,求该矩阵的k次幂,即P^k。

输入描述:

 第一行:两个整数n(2<=n<=10)、k(1<=k<=5),两个数字之间用一个空格隔开,含义如上所示。 接下来有n行,每行n个正整数,其中,第i行第j个整数表示矩阵中第i行第j列的矩阵元素Pij且(0<=Pij<=10)。另外,数据保证最后结果不会超过10^8。

输出描述:

对于每组测试数据,输出其结果。格式为: n行n列个整数,每行数之间用空格隔开,注意,每行最后一个数后面不应该有多余的空格。
示例1

输入

2 2 9 8 9 3 

输出

153 96 108 81 

#include<iostream> #include<string.h> using namespace std; struct matrix//用结构体封装一个二维数组 { int a[100][100]; }origin,res;//res.a存储结果 ,origin.a是需要运算的矩阵 ,k是次幂数 ,n是矩阵的行列数 matrix multiply(matrix x,matrix y,int N)//矩阵相乘函数 { matrix temp; memset(temp.a,0,sizeof(temp.a)); for(int i=0;i<N;i++) { for(int j=0;j<N;j++) { for(int k=0;k<N;k++) { temp.a[i][j]+=x.a[i][k]*y.a[k][j]; } } } return temp; } int main() { int n,k,i,j; while(~scanf("%d%d",&n,&k)) { for(i=0;i<n;i++) for(j=0;j<n;j++) scanf("%d",&origin.a[i][j]); k--; res=origin; while(k--) { res=multiply(res,origin,n); } for(i=0;i<n;i++){ for(j=0;j<n-1;j++) printf("%d ",res.a[i][j]); printf("%d\n",res.a[i][n-1]); } } return 0; } 

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/232746.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • pycharm R语言_pycharm安装后无解释器

    pycharm R语言_pycharm安装后无解释器最近被网上的一些R语言做出的图迷住,遂突发奇想学习R语言。按照教程装好了R语言解释器与Rstudio后,才得知pycharm中有R语言的插件,安装后可以在pycharm中写R语言。pycharm安装R语言插件打开File|Settings|Plugins,在搜索框里搜索RLanguageSupport,安装就可以了。安装结束后重启pycharm,再进入File|Settin…

    2022年8月28日
    3
  • 第十六届全国大学生智能车竞赛华东赛区_第十五届全国大学生智能汽车竞赛

    第十六届全国大学生智能车竞赛华东赛区_第十五届全国大学生智能汽车竞赛第十七届全国大学生智能车竞赛华南赛区奖项

    2022年9月30日
    3
  • vue的路由懒加载_vue图片懒加载

    vue的路由懒加载_vue图片懒加载vue打包后的js文件越来越大,这会是影响加载时间的重要因数。当构建的项目比较大的时候,懒加载可以分割代码块,提高页面的初始加载效率。下面是几种常见vue路由懒加载的方法。

    2022年10月6日
    2
  • python ip池(python 连接池)

    ,都说标题是文章的灵魂,想了半天没想到什么比较有创意的标题,只好拿一个去年的爆款标题套一下。啊哈哈哈哈哈哈,朕真是太机智了这是一篇介绍如何使用python搭建IP池的文章,如果爱卿对此不感兴趣,那很抱歉,标题耽误了你宝贵的时间。事情的起因是这样,前段时间我写了一篇介绍如何爬取小说的blog【python那些事.No2】,在爬取的过程中,发现同一个IP连续只能获取前几页小说内容,原本是想搭建…

    2022年4月11日
    39
  • python求最大公约数和最小公倍数的两种方法

    python求最大公约数和最小公倍数的两种方法最大公约数和最小公倍数的求解可以归结为求最大公约数,最小公倍数为两数乘积除以最大公倍数这里介绍两种求解方法,一种数常规易于理解的,一种是用辗转相除法实现的#求最大公倍数和最小公约数a=int(input(“输入一个数字”))b=int(input(“再输入一个数字”))m=[]ifa>b:smaller=belse:smaller=afori…

    2022年5月17日
    46
  • delphi xe datasnap 服务器显示客户端,Delphi xe datasnap[通俗易懂]

    delphi xe datasnap 服务器显示客户端,Delphi xe datasnap[通俗易懂]我想从客户端向服务端提交多个OleVariant内容.最初我想这样这实现functionSaveData(aDataArr:arrayofOleVariant;aTableArr:arrayofstring;aKeyArr:arrayofstring;varaErrorStr:string):Boolean;这样经测试不行,DATASNAP参数不能为数组.现在我用TJSONObje…

    2022年7月18日
    12

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号