poj The Clocks 高斯消元

由于数据量不大，所以这题有很多解法。

我用的是高斯消元化为逆矩阵解决的……

代码如下：

#include<stdio.h>
#include<iostream>
using namespace std;
int an[][]={
{,,,,,,,,,},
{,,,,,,,,,},
{,,,,,,,,,},
{,,,,,,,,,},
{,,,,,,,,,},
{,,,,,,,,,},
{,,,,,,,,,},
{,,,,,,,,,},
{,,,,,,,,,}
},ans[];
int main(){
    int n,i,j,k,sum,ma,mb;
    for(i=;i<;i++){
        cin>>an[i][];
        an[i][]=(-an[i][])%;
    }
    for(i=;i<;i++){
        if(an[i][i]==){
            for(k=i+;k<;k++){
                if(an[k][i]){
                    for(j=;j<=;j++)
                        swap(an[i][j],an[k][j]);
                    break;
                }
            }
        }
        for(k=i+;k<;k++){
            if(an[k][i]){
                ma=an[i][i];
                mb=an[k][i];
                for(j=i;j<=;j++){
                    an[k][j]=ma*an[k][j]-mb*an[i][j];
                    an[k][j]=(an[k][j]%+)%;
                }
            }
        }
    }
    sum=;
    for(i=;i>=;i--){
        for(j=i+;j<;j++){
            an[i][]-=ans[j]*an[i][j];
            an[i][]=(an[i][]%+)%;
        }
        for(ans[i]=;ans[i]<=;ans[i]++)
            if((ans[i]*an[i][i]%+)%==an[i][])
                break;
        ans[i]%=;
        sum+=ans[i];
    }
    j=;k=;
    while(j<){
        if(ans[j]){
            cout<<j+;
            k++;
            if(k>=sum) cout<<endl;
            else cout<<' ';
            ans[j]--;
        }
        else j++;
    }
    return ;
}