수다닷컴

질문 제목 :
쓰레기값 나오는 이유가 뭘까요 그리고 무엇을 고쳐야할까요

희소행렬의곱질문 내용 : #include stdio.h
#include stdlib.h
#define max_col 50 // 최대 열의 수 + 1
#define max_terms 101 //
#define max_size 100
typedef struct{
int row;
int col;
int value;
}term;

void transpose(term a[ ], term b[ ] )
// a는 입력 행렬, b는 출력 행렬. b = at
{
int count, i, j, currentb;
count = a[0].value; // a에서 0이 아닌 원소의 수
b[0].row = a[0].col;// b의 행의 수 = a의 열의 수
b[0].col = a[0].row; // b의 열의 수 = a의 행의 수
b[0].value = count; // 0이 아닌 원소 수는 a와 동일
if (count 0) { // a가 empty matrix가 아니라면…
currentb = 1;// 새로운 원소가 b에 저장될 위치.
for ( i = 0; i a[0].col; i++ ) {
// a의 열 순서로 전치 연산 실행 (i: 현재 열)
for (j = 1; j = count; j++) {
// a에서 현재 열을 찾자. (a는 행 순서로 저장)
if ( a[ j ].col == i ) {
// 현재 열의 원소 발견. b에 추가하자.
b[currentb].row = a[ j ].col;
b[currentb].col = a[ j ].row;
b[currentb].value = a[ j ].value;
currentb++;// b의 저장될 위치를 1 증가
}
}}}
}
void storesum(term d[], int *totald, int row, int column, int *sum)
{
/* sum이 0이 아니면, d 배열의 *totald+1 위치에 row, column 값과 함께 저장 */
if (*sum)
if (*totald max_terms) {
d[++*totald].row = row;
d[*totald].col = column;
d[*totald].value = *sum;
*sum = 0;
}
else {
fprintf(stderr, numbers of terms in product exceeds %d\n, max_terms);
exit(1);
}
}
int compare(int a,int b){
if(a==b){
return 0;}
else if(ab){
return 1;}
else if(ab){
return -1;}

}
void mmult(term a[], term b[], term d[])
{// a, b: 입력, d: 출력, d = a * b
int i, j, column, totald = 0;
int rows_a = a[0].row, cols_a = a[0].col, totala = a[0].value;
int cols_b = b[0].col, totalb = b[0].value;
int row_begin = 1, row = a[1].row, sum = 0; // row: a의 현재 행
term new_b[max_terms];
if (cols_a != b[0].row) {// a의 열의 수와 b의 행의 수는 동일
fprintf(stderr, incompatible matricses\n); exit(1);
}
transpose(b, new_b);
// 경계 조같麗?조건을 설정
a[totala+1].row = rows_a;
new_b[totalb+1].row = cols_b;
new_b[totalb+1].col = 0;
for (i = 1; i = totala; ) {
column = new_b[1].row;// b의 현재 열
for (j = 1; j = totalb + 1; ) {
// a의 현재 행과 b의 현재 열에 대해 곱셈 수행
if (a[i].row != row) { // a의 현재 행을 벗어남.
storesum(d, &totald, row, column, &sum);
i = row_begin;// b는 다음 열로. a는 원 위치로.
for (; new_b[j].row == column; j++) ;
column = new_b[j].row;
}
else if (new_b[j].row != column) { // b의 현재 열을 벗어남.
storesum(d, &totald, row, column, &sum);
i = row_begin;// a는 원 위치
column = new_b[j].row; // b는 다음 열로.
}
else switch (compare(a[i].col, new_b[j].col)){
case -1: // a[i].col new_b[j].col. a 증가
i++;break;
case 0: // 계산 후, a와 b를 모두 진행
sum += (a[i++].value * new_b[j++].value); break;
case 1: // a[i].col new_b[j].col. b 증가
j++;
}
}
for (; a[i].row == row; i++) ;// b의 모든 원소를 처리한 후,
row_begin = i;row = a[i].row;// a의 현재 행을 다음 행으로.
} /* end of for i = totala */
d[0].row = rows_a;
d[0].col = cols_b; d[0].value = totald;
}

void print(term d[]){
int i;
for(i=0; id[0].value+1; i++){
printf(%d %d %d\n,d[i].row,d[i].col,d[i].value);
}
}

void main(){
term a[max_terms]={{6,6,8},{0,0,15},{0,3,22},{0,5,-15},{1,1,11},{1,2,3},{2,3,-6},{4,0,91},{5,2,28}};
term b[max_terms]={{6,2,5},{0,0,1},{2,0,3},{4,0,2},{1,1,1},{3,1,1}};
term d[max_terms];
mmult(a,b,d);
print(d);
}

신청하기