高斯消元

double型 fabs

double  a[maxn][maxn] , x[maxn] ;  //a[i][j] 系数矩阵 ， a[i][n+1]  = y[i] ， x解
int     n ;   //n个方程
void    guass(){
int i , j , k ;
double  sum , rate ;
for(k = 1 ; k < n ;  k++){
for(i = j = k ; i <= n ; i++) j = fabs(a[i][k]) > fabs(a[j][k]) ? i : j ;
for(i = k ; i <= n+1 ; i++)  swap(a[j][i] , a[k][i]) ;
for(i = k+1 ; i <= n ; i++)
for(rate = a[i][k]/a[k][k] , j = k ; j <= n+1 ; j++)
a[i][j] -= a[k][j]*rate ;
}
for(i = n ; i >= 1 ; i--){
for(sum = 0 , j = i+1 ; j <= n ; sum += a[i][j]*x[j],j++) ;
x[i] = (a[i][n+1] - sum)/a[i][i] ;
}
}

double  Fix(double x){   // 将 -0.000 变成 0.00 ,结果输出Fix(x[])
if(fabs(x) < 1e-4) return 0 ;
return x ;
}

Xor 的高斯消元

int a[maxn][maxn] , x[maxn] ;
int     n ;
void    guass(){
int i , j , k ;
double  sum , rate ;
for(k = 1 ; k < n ;  k++){
for(i = k ; i <= n ; i++) if(a[i][k]){ j = i ; break ;}
for(i = k ; i <= n+1 ; i++)  swap(a[j][i] , a[k][i]) ;
for(i = k+1 ; i <= n ; i++){
if(a[i][k] == 0) continue ;
for(j = k ; j <= n+1 ; j++)
a[i][j] ^= a[k][j] ;
}
}
for(i = n ; i >= 1 ; i--){
x[i] = a[i][n+1] ;
for(j = i+1 ; j <= n ; x[i] ^= (a[i][j]&&x[j]),j++) ;
}
}