对角矩阵.cpp

矩阵的一些算法应用

#include "iostream.h"

void print(int a[],int n,int w)
{int i,j,k,m;
for(i=1;i<=n;i++)
{
 for(j=1;j<=n;j++)
 { m=i-j; 
   if(m<0) m=-m;
	if(m<(w-1))  
	{k=(w-1)*i+j-w;
	cout<<a[k]<<"  ";}
	else cout<<"   ";
}
 cout<<endl;
}
}
void duijiaoadd(int a[],int b[],int c[],int n,int w)
{
int i,j,k,m;
for(i=1;i<=n;i++)
 for(j=1;j<=n;j++)
 { m=i-j; 
   if(m<0) m=-m;
	if(m<(w-1))  
	{k=(w-1)*i+j-w;
	c[k]=a[k]+b[k];
	}

}

}
void duijiaosub(int a[],int b[],int c[],int n,int w)
{
int i,j,k,m;
for(i=1;i<=n;i++)
 for(j=1;j<=n;j++)
 { m=i-j; 
   if(m<0) m=-m;
	if(m<(w-1))  
	{k=(w-1)*i+j-w;
	c[k]=a[k]-b[k];
	}

}

}

void main()
{
	int b[]={1,2,3,4,5,6,7,8,9,10,11,12,13};
	int a[]={14,15,16,17,18,19,20,21,22,23,24,25,26};
    int c[25];
	print(b,5,3);
	print(a,5,3);
	duijiaoadd(a,b,c,5,3);
	print(c,5,3);
	duijiaosub(a,b,c,5,3);
	print(c,5,3);
}