高斯消元模板

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];
}