void guess()//n行m列的矩阵a,且都为整数的
{
rep(i,n)
{
int k=-1;
repf(j,i,n-1)
if(a[j][i]!=0)
{
k=j;break;
}
if(k==-1) continue;//这一列都为空的
if(k!=i)
rep(j,m) swap(a[i][j],a[k][j]);
rep(j,n)
{
if(j!=i)
{
int x=a[i][i],y=a[j][i];
rep(k,m)
a[j][k]=a[j][k]*x-a[i][k]*y;
}
}
}
}
高斯消元
最新推荐文章于 2023-07-24 18:50:40 发布