poj解题报告——1914

        这就是一道简单的线性代数问题，三阶矩阵求值。

代码如下

#include<stdio.h>
#include<string.h>
__int64 a[3][3],b[3],c[3][3];
__int64 cal()
{
	return a[0][0]*(a[1][1]*a[2][2]-a[1][2]*a[2][1])
		-a[0][1]*(a[1][0]*a[2][2]-a[1][2]*a[2][0])
		+a[0][2]*(a[1][0]*a[2][1]-a[1][1]*a[2][0]);
}
double val(double a)
{
	if(a<0.0005&&-0.0005<a)
        return 0;
    return a;
}
void main()
{
	int i,j,t;
	__int64 A,ans[3];
	scanf("%d",&t);
	while(t--)
	{
		for(i=0;i<3;i++)
		{
			for(j=0;j<3;j++)
			{
				scanf("%I64d",&a[i][j]);
				c[i][j]=a[i][j];
			}
			scanf("%I64d",&b[i]);
		}
		A=cal();
		for(i=0;i<3;i++)
		{
			for(j=0;j<3;j++)
				a[j][i]=b[j];
			ans[i]=cal();
			for(j=0;j<3;j++)
				a[j][i]=c[j][i];
		}
		for(i=0;i<3;i++)
			printf("%I64d ",ans[i]);
		printf("%I64d\n",A);
		if(A==0)
			printf("No unique solution\n");
		else
		{
			printf("Unique solution:");
			for(i=0;i<3;i++)
				printf(" %.3f",val(ans[i]*1.0/A));
			printf("\n");
		}
		printf("\n");
	}
}