行索引稀疏矩阵查找某一列的元素没那么方便,所以在做矩阵乘法时(这里以M乘N=Q为例),严书的做法是:在求Q的某一行的值是,用M的一行去遍历N的每一行,对结果中同列的值进行累加,最后稀疏化存入Q的当前行中,这个过程还原成正常矩阵比较容易理解:
求Q(2,2)的第一行时,肯定是M的第一行和N的第一列逐乘再累加,然后再M的第一行和N的第二列逐乘累加
M(2,3)* N(3,2):
直接算的话就是:
Q[1,1] = 1*1 + 2*0 + 3*1 = 1+0+3 = 4
Q[1,2] = 1*0 + 2*1 + 3*1 = 0+2+3 = 5
这回我们直接同时求两列的值,其实就是改变一下计算顺序:
这个就相当于严书代码中ctemp的作用,只不过ctemp是直接累加.
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MAX_E = 1e+4+5, MAX_R = 105;
//严版特色,下标从1开始
struct Triple{
int i,j,e;};
struct RLSMatrix{
Triple data[MAX_E];
int rpos[MAX_R];//行首索引, rpos[i]指向第i行中的首元素在data中的序号, 则rpos[i+1]-1指向第i行中最后一个非0元素在N.data中的序号.
int rows, cols, elems;
};
int ctemp[MAX_R];//Q的第i行元素结果累加器,算完一行就压缩存储进Q.data中
RLSMatrix mul(RLSMatrix M, RLSMatrix N){
RLSMatrix Q;
Q.rows = M.rows;Q.cols = N.cols;Q.elems = 0;
if(M.elems * N.elems != 0){
// Q是非零矩阵
for(int arow = 1; arow <= M.rows; arow++){
memset(ctemp, 0, sizeof(ctemp)); //到了下一行,清零
Q.rpos[arow] = Q.elems+1; //新一行的行首索引是当前data中元素个数+1,从该处继续向data中存三元组
int tp; //用tp指向data中该行行末的下一个位置,方便遍历
if(arow < M.rows) tp = M.rpos[arow+1]; //如果不是最后一行, tp指向下一行行首
else tp = M.elems+1; //最后一行, tp直接指向data中最后一个+1
for(int p = M.rpos[arow];p < tp; p++){
//对当前行找到的每一个非零元
int brow = M.data[p].j; //找到对应元在N中的行号(M中某行第三列的元素, 肯定和N中第三行某列的元素相乘)
int t;//和tp同样的套路
if(brow < N.rows) t = N.rpos[brow+1];
else t = N.elems+1;
for(int q = N.rpos[brow]; q < t; q++){
//这里不是直接算出M的某行乘N的某列的值
//而是用M的某行,去遍历N的所有行,然后对每行的应该在同列的值进行累加
int ccol = N.data[q].j;
ctemp[ccol] += M.data[p].e * N.data[q].e;//累加器
}
}//求得Q中第crow(=arow)行的非零元
for(int ccol = 1; ccol <= Q.cols; ++ccol){
//用M的一行遍历完了整个N矩阵
if(ctemp[ccol]){
Q.data[++Q.elems] = {
arow, ccol, ctemp[ccol]};//压缩存储
}
}
}
}
return Q;
}
RLSMatrix makeMat(){
/*
* 输入rows, cols, 三元组个数
* 输入三元组
*/
RLSMatrix A;
cin>>A.rows>>A.cols>>A.elems;
for(int i = 1;i <= A.elems;i++){
int x,y,e;
cin>>x>>y>>e;
A.data[i] = {
x,y,e};
}
/*
根据乘法中这段代码写的初始化rpos数组:
int tp; //用tp指向data中该行行末的下一个位置,方便遍历
if(arow < M.rows) tp = M.rpos[arow+1]; //如果不是最后一行, tp指向下一行行首
else tp = M.elems+1; //最后一行, tp直接指向data中最后一个+1
for(int p = M.rpos[arow];p < tp; p++) ...
可以看出如果某行没有元素,则让tp = M.rpos[arow]则可不进行遍历,也就是M.rpos[arow] = M.rpos[arow+1]
而当最后几行为空时,并不会执行else tp = M.elems+1;这时应该把rpos最后几行空的填成M.elems+1
*/
int row = 1;
for(int e = 1;e <= A.elems; e++){
int arow = A.data[e].i;
while(row <= arow){
A.rpos[row++] = e;
}
}
while(row <= A.rows){
//如果最后几行没有元素,让索引指到elems+1的位置
A.rpos[row++] = A.elems+1;
}
return A;
}
void printMat(RLSMatrix A){
cout<<"rows:"<<A.rows<<" cols:"<<A.cols<<" elems:"<<A.elems<<endl;
cout<<"-------------------------"<<endl;
cout<<"data:"<<endl;
for(int i = 1;i <= A.elems;i++) {
cout<<setw(4)<<A.data[i].i<<setw(4)<<A.data[i].j<<setw(4)<<A.data[i].e<<endl;
}
cout<<"-------------------------"<<endl;
cout<<"rpos: ";
for(int j = 1;j <= A.rows; j++){
cout<<A.rpos[j]<<' ';
}
cout<<endl<<endl;
}
int main(){
ios::sync_with_stdio(false);
RLSMatrix A = makeMat();
printMat(A);
RLSMatrix B = makeMat();
printMat(B);
RLSMatrix C = mul(A, B);
printMat(C);
return 0;
}
/*
严书上的例子
3 4 4
1 1 3
1 4 5
2 2 -1
3 1 2
4 2 4
1 2 2
2 1 1
3 1 -2
3 2 4
*/
文章浏览阅读2.6k次。解决Python词云库wordcloud不显示中文的问题2018-11-25背景:wordcloud是基于Python开发的词云生成库,功能强大使用简单。github地址:https://github.com/amueller/word_cloudwordcloud默认是不支持显示中文的,中文会被显示成方框。安装:安装命令:pip install wordcloud解决:经过测试发现不支持显示中文..._词云python代码无法输出文字
文章浏览阅读807次。扩展。_jmap 在线分析
文章浏览阅读1.1w次。随着炎热夏季的到来,当玩游戏正爽的时候,电脑突然死机了,自动关机了,是不是有想给主机一脚的冲动呢?这个很大的原因是因为CPU温度过高导致的。很多新手玩家可能都有一个疑虑,cpu温度多少以下正常?有些说是60,有些说是70,到底多高CPU温度不会死机呢?首先我们先看看如何查看CPU的温度。下载鲁大师并安装,运行鲁大师软件,即可进入软件界面,并点击温度管理,即可看到电脑各个硬件的温度。鲁大师一般情况下..._台式机玩游戏温度多少正常
文章浏览阅读243次。Day2-打印打印打印!我终于更新了!(哭腔)一、 最简单的打印最最简单的打印语句: print(“打印内容”)注意:python是全英的,符号记得是半角下面是我写的例子:然后进入power shell ,注意:你需要使用cd来进入你保存的例子的文件夹,保存时名字应该取为xxx.py我终于知道为什么文件夹取名都建议取英文了,因为进入的时候是真的很麻烦!如果你没有进入正确的文件夹..._puthon打印任务收获
文章浏览阅读1k次。centos8问题参考CentOS 8 EOL如何切换源? - 云服务器 ECS - 阿里云_"cenerrors during download metadata for repository \"appstream"
文章浏览阅读2.7k次,点赞3次,收藏11次。SpringBoot+Maven+MabatisPlusmaven在新建springboot项目引入RELEASE版本出错maven在新建springboot项目引入RELEASE版本出错maven详解maven就是通过pom.xml中的配置,就能够从仓库获取到想要的jar包。仓库分为:本地仓库、第三方仓库(私服)、中央仓库springframework.boot:spring-boot-starter-parent:2.2.1.RELEASE’ not found若出现jar包下载不了只有两_基于微服务的在线教育平台尚硅谷
文章浏览阅读316次。路由的概念路由器它称之为网关设备。路由器就是用于连接不同网络的设备路由器是位于OSI模型的第三层。路由器通过路由决定数据的转发。网关的背景:当时每家计算机厂商,用于交换数据的通信程序(协议)和数据描述格式各不相同。因此,就把用于相互转换这些协议和格式的计算机称为网关。路由器与三层交换器的对比路由协议对比路由器的作用:1.路由寻址2.实现不同网络之间相连的功能3.通过路由决定数据的转发,转发策略称为 路由选择。VLAN相关技术什么是VLAN?中文名称叫:虚拟局域网。虚_路由和vlan
文章浏览阅读2.8w次,点赞6次,收藏22次。设置div背景颜色透明度,内部元素不透明:.demo{ background-color:rgba(255,255,255,0.15) } 错误方式:.demo{ background-color:#5CACEE;opacity:0.75;} 这样会导致div里面的元素内容和背景颜色一起变透明只针对谷歌浏览器的测试_div设置透明度,里面的内容不透明
文章浏览阅读563次。1.[ u]文字:在文字的位置可以任意加入您需要的字符,显示为下划线效果。2.[ align=center]文字:在文字的位置可以任意加入您需要的字符,center位置center表示居中,left表示居左,right表示居右。5.[ color=red]文字:输入您的颜色代码,在标签的中间插入文字可以实现文字颜色改变。6.[ SIZE=数字]文字:输入您的字体大小,在标签的中间插入文..._discuzcode 大全
文章浏览阅读2.6k次。iOS中定时器有三种,分别是NSTimer、CADisplayLink、dispatch_source,下面就分别对这三种计时器进行说明。一、NSTimerNSTimer这种定时器用的比较多,但是特别需要注意释放问题,如果处理不好很容易引起循环引用问题,造成内存泄漏。1.1 NSTimer的创建NSTimer有两种创建方法。方法一:这种方法虽然创建了NSTimer,但是定时器却没有起作用。这种方式创建的NSTimer,需要加入到NSRunLoop中,有NSRunLoop的驱动才会让定时器跑起来。_ios nstimer
文章浏览阅读4.8k次,点赞17次,收藏51次。Linux的命令有几百个,对程序员来说,常用的并不多,考虑各位是初学者,先学习本章节前15个命令就可以了,其它的命令以后用到的时候再学习。1、开机 物理机服务器,按下电源开关,就像windows开机一样。 在VMware中点击“开启此虚拟机”。2、登录 启动完成后,输入用户名和密码,一般情况下,不要用root用户..._ls-lmore
文章浏览阅读4.1k次。1.登录MYSQL系统命令打开DOS命令框shengfen,以管理员的身份运行命令1:mysql -u usernae -p password命令2:mysql -u username -p password -h 需要连接的mysql主机名(localhost本地主机名)或是mysql的ip地址(默认为:127.0.0.1)-P 端口号(默认:3306端口)使用其中任意一个就OK,输入命令后DOS命令框得到mysql>就说明已经进入了mysql系统2. 查看mysql当中的._mysql -u user