高斯消元模板

高斯消元模板

void gauss() {    fo(i,1,n-1) {        fo(j,i+1,n) if (abs(matrix[j].a[i])>abs(matrix[i].a[i]))         swap(matrix[i],matrix[j]);        fo(j,i+1,n) {            db f=matrix[j].a[i]/matrix[i].a[i];            fo(k,i,n) matrix[j].a[k]=matrix[i].a[k]*f-matrix[j].a[k];            matrix[j].b=matrix[i].b*f-matrix[j].b;        }    }    fd(i,n,2) {        fo(j,1,i-1) {            db f=matrix[j].a[i]/matrix[i].a[i];            fo(k,j,n) matrix[j].a[k]=matrix[i].a[k]*f-matrix[j].a[k];            matrix[j].b=matrix[i].b*f-matrix[j].b;        }    }    fo(i,1,n) an[i]=matrix[i].b/matrix[i].a[i];}