sequ.cpp

这是学校ACM程序设计大赛的题目,拿出来大家算是当练习

#include <stdio.h>
#include <stdlib.h>

void solv_e(double **a, int total)
{
/*	for(int i=0; i<total-1; i++)
	{
		if(a[i][i]==0)
			continue;
		else
		{
			for(int j=i+1; j<total; j++)
			{
				double temp=-a[j][i]/a[i][i];
				for(int k=i; k<=total; k++)
				{
					a[j][k]=a[i][k]*temp+a[j][k];
				}
			}
		}
	}*/
	for(int i=1; i<total; i++)
	{
		double temp=-a[i][0]/a[0][0];
		for(int k=0; k<=total; k++)
			a[i][k]=a[0][k]*temp+a[i][k];
	}
/*	for(i=2; i<total; i++)
	{
		double temp=-a[i][1]/a[1][1];
		for(int k=1; k<=total; k++)
			a[i][k]=a[1][k]*temp+a[i][k];
	}*/
}
void main()
{
	//读文件操作
	FILE *fp;
	if((fp=fopen("dao1.in","r"))==NULL)
	{
		printf("can't open file\n");
		exit(0);
	}
	int circle;
	fscanf(fp,"%d",&circle);
	for(int i=0; i<circle; i++)
	{

		int total;
		fscanf(fp,"%d",&total);
	//	printf("%d\n",total);
		double **array;
		array=(double **)malloc((total+5)*sizeof(double *));
		for(int j=0; j<total; j++)
			array[j]=(double *)malloc((total+5)*sizeof(double));
		for(int k=0; k<total; k++)
		{
			for(int s=0; s<=total; s++)
				fscanf(fp,"%lf",*(array+k)+s);
		}

		//进行行列式化简
		solv_e(array,total);
		//打印行列式
		for(k=0; k<10; k++)
		{
			for(int s=0; s<=total; s++)
				printf("%.0lf ",array[k][s]);
			printf("\n");
		}
		printf("***********\n");
	}
}