看嚴蔚敏的數據結構關于稀疏矩陣相乘的一節，覺得有點復雜，自己按照個人理解，寫了一個相乘的程序，把相乘后的矩陣先轉化成了一維數組，代碼如下，大家指點指點（注：本人非計算機專業，若果有低級錯誤，見笑）
1 #include <stdio.h>
2
3 #include <stdlib.h>
4
5
6
7 #define max 12
8
9
10
11 typedef struct{
12
13 int i,j;
14
15 int e;
16
17 }triple;
18
19
20
21 typedef struct{
22
23 triple data[max+1];
24
25 int mu,tu,nu;
26
27 }tsmatrix;
28
29
30
31 void main()
32
33 {
34
35 tsmatrix M;
36
37 tsmatrix N;
38
39 int Q[7]={0};//數組下標從1開始
40
41 int p=1,q=1;//下面很長的一部分，是直接輸入三元組，
42
43 int s=1;
44
45
46
47 M.mu=3;//行數
48
49 M.nu=5;//列數
50
51 M.tu=6;//非零元的個數
52
53
54
55 M.data[1].i=1;
56
57 M.data[1].j=1;
58
59 M.data[1].e=3;
60
61
62
63 M.data[2].i=1;
64
65 M.data[2].j=4;
66
67 M.data[2].e=5;
68
69
70
71 M.data[3].i=1;
72
73 M.data[3].j=5;
74
75 M.data[3].e=-2;
76
77
78
79 M.data[4].i=2;
80
81 M.data[4].j=2;
82
83 M.data[4].e=-1;
84
85
86
87 M.data[5].i=3;
88
89 M.data[5].j=1;
90
91 M.data[5].e=2;
92
93
94
95 M.data[6].i=3;
96
97 M.data[6].j=5;
98
99 M.data[6].e=4;
100
101
102
103 N.mu=5;//行數
104
105 N.nu=2;//列數
106
107 N.tu=6;//非零元的個數
108
109
110
111 N.data[1].i=1;
112
113 N.data[1].j=2;
114
115 N.data[1].e=2;
116
117
118
119 N.data[2].i=2;
120
121 N.data[2].j=1;
122
123 N.data[2].e=1;
124
125
126
127 N.data[3].i=3;
128
129 N.data[3].j=1;
130
131 N.data[3].e=-2;
132
133
134
135 N.data[4].i=3;
136
137 N.data[4].j=2;
138
139 N.data[4].e=4;
140
141
142
143 N.data[5].i=5;
144
145 N.data[5].j=1;
146
147 N.data[5].e=3;
148
149
150
151 N.data[6].i=5;
152
153 N.data[6].j=2;
154
155 N.data[6].e=6;
156
157
158
159 for(p=1;p<=M.tu;p++) //開始進入正題
160
161 {
162
163 for(q=1;q<=N.tu;q++)
164
165 {
166
167 if(M.data[p].j==N.data[q].i)
168
169 {
170
171 printf("p=%d,q=%d \n",p,q);//打印測試
172
173
174
175 printf("%d,%d,%d \n",M.data[p].e,N.data[q].e,(M.data[p].i-1)*N.nu+N.data[q].j);//打印測試
176
177
178
179 Q[(M.data[p].i-1)*N.nu+N.data[q].j]+=M.data[p].e*N.data[q].e;//將相乘的結果先存入一維數組
180
181 }
182
183 else continue;
184
185 }
186
187 }
188
189
190
191 for(s=1;s<=6;s++)
192
193 {
194
195 printf("%d ",Q[s]);//打印一維數組
196
197 }
198
199 }
200
posted on 2010-08-04 13:06
ypp 閱讀(482)
評論(0) 編輯 收藏 引用 所屬分類:
c++