hdu 5755 Gambler Bo （高斯消元法解同余方程组）

http://acm.hdu.edu.cn/showproblem.php?pid=5755

题意：

n*m矩阵，每个格有数字0/1/2

每选择一个格子，这个格子+2，4方向相邻格子+1

如何选择格子，可以使每个格子的数最后 %3=0

#include<cstdio>
#include<cstring>
#include<algorithm>

using namespace std;

const int mod=;

int n,m,t;

int a[][];

int x[];

int turn(int i,int j)
{
    return (i-)*m+j-;
}

int getgcd(int a,int b) { return !b ? a : getgcd(b,a%b); }

int getlcm(int a,int b) { return a*b/getgcd(a,b); }

void gauss()
{
    int equ=n*m,var=n*m;
    int i,j,k;
    int max_r,col;
    int ta,tb;
    int lcm;
    int tmp;
    for(k=,col=;k<equ && col<var;++k,++col)
    {
        max_r=k;
        for(i=k+;i<equ;++i)
            if(abs(a[i][col])>abs(a[max_r][col])) max_r=i;
        if(!a[max_r][col]) { --k; continue; }
        if(max_r!=k) swap(a[k],a[max_r]);
        for(i=k+;i<equ;++i)
            if(a[i][col])
            {
                lcm=getlcm(abs(a[i][col]),abs(a[k][col]));
                ta=lcm/abs(a[i][col]);
                tb=lcm/abs(a[k][col]);
                if(a[i][col]*a[k][col]<) tb=-tb;
                for(j=col;j<var+;++j) a[i][j]=((a[i][j]*ta-a[k][j]*tb)%mod+mod)%mod;
            }
    }
    for(int i=var-;i>=;--i)
    {
        int tmp=a[i][var];
        for(int j=i+;j<var;++j)
            if(a[i][j])
            {
                tmp-=a[i][j]*x[j];
                tmp=(tmp%mod+mod)%mod;
            }
        x[i]=a[i][i]*tmp%mod;
    }
}

int main()
{
    int T;
    scanf("%d",&T);
    int xi,tmp;
    int ans;
    while(T--)
    {
        memset(a,,sizeof(a));
        scanf("%d%d",&n,&m);
        t=n*m;
        for(int i=;i<=n;++i)
            for(int j=;j<=m;++j)
            {
                scanf("%d",&xi);
                a[turn(i,j)][t]=-xi;
            }
        for(int i=;i<=n;++i)
            for(int j=;j<=m;++j)
            {
                tmp=turn(i,j);
                a[tmp][tmp]=;
                if(i>) a[turn(i-,j)][tmp]=;
                if(i<n) a[turn(i+,j)][tmp]=;
                if(j>) a[turn(i,j-)][tmp]=;
                if(j<m) a[turn(i,j+)][tmp]=;
            }
        gauss();
        ans=;
        for(int i=;i<t;++i) ans+=x[i];
        printf("%d\n",ans);
        for(int i=;i<=n;++i)
            for(int j=;j<=m;++j)
            {
                tmp=turn(i,j);
                if(x[tmp])
                    while(x[tmp]--) printf("%d %d\n",i,j);
            }
    }
}