UVA 12849 Mother’s Jam Puzzle( 高斯消元 )

题目： http://uva.onlinejudge.org/external/128/12849.pdf

#include <bits/stdc++.h>
using namespace std;
typedef long long LL ;
typedef unsigned long long ULL ;
typedef pair<int,int> pii;
#define X first
#define Y second
const int N =  ;
const int mod = 1e9+ ;
const double K = ;
double A[N][N] ;

void solve( int n ) {
    int i , k , j , r ;
    for( i =  ; i < n ; ++i ) {
        r = i ;
        for( j = i +  ; j < n ; ++j )
            if( fabs(A[j][i]) > fabs( A[r][i] ) ) r = j ;

        if( r != i ) for( j =  ; j <= n ; ++j ) swap( A[r][j] , A[i][j] ) ;

        for( j = n ; j >= i ; j-- ) {
            for( k = i +  ; k < n ; ++k ) {
                A[k][j] -= A[k][i] / A[i][i] * A[i][j] ;
            }
        }
    }
    for( i = n- ; i >= ; i-- ) {
        for( j = i +  ; j < n ; ++j )
            A[i][n] -= A[j][n] * A[i][j] ;
        A[i][n] /= A[i][i];
    }
}

int Run() {
    int _ , cas =  ; scanf("%d",&_);
    while( _-- ) {
        printf("Case %d: ",cas++);
        for( int i =  ; i <  ; ++i ) {
            for( int j =  ; j <  ; ++j ) {
                scanf("%lf",&A[i][j]);
                A[i][j] *= K ;
            }
        }
        solve() ;
        printf("%.2lf %.2lf %.2lf\n",A[][],A[][],A[][]);
    }
    return  ;
}

int main()
{
    #ifdef LOCAL
        freopen("in","r",stdin);
    #endif
    ios::sync_with_stdio();
    return Run();
}