本文共 1252 字，大约阅读时间需要 4 分钟。
可以设n*m个未知量，建立n*m个方程。位置i，j可以建立方程 (2*x[i*m+j]+x[(i-1)*m+j]+x[(i+1)*m+j]+x[i*m+j-1]+x[i*m+j+1])%3=3-b[i][j]; 套了个高斯消元的板子过了。
#pragma comment(linker, "/STACK:1024000000,1024000000")#include #include #include #include #include #include #include #include #include #include using namespace std;typedef long long LL;const double pi=acos(-1.0),eps=1e-8;void File(){ freopen("D:\\in.txt","r",stdin); freopen("D:\\out.txt","w",stdout);}inline int read(){ char c = getchar(); while(!isdigit(c)) c = getchar(); int x = 0; while(isdigit(c)) { x = x * 10 + c - '0'; c = getchar(); } return x;}const int mod = 3;int exgcd(int a,int b,int &x,int &y){ if(!b){x = 1; y = 0; return a;} else{ int r = exgcd(b,a%b,y,x); y -= x * (a/b); return r; }}int lcm(int a,int b){ int x = 0, y =0; return a / exgcd(a,b,x,y) * b;}const int MAXN=1000;int A[MAXN][MAXN],free_x[MAXN],x[MAXN];void Gauss(int n,int m){ int r,c; for(r=0,c=0;r abs(A[maxr][c])) maxr = i; if(maxr != r) for(int i=c;i<=m;i++) swap(A[r][i],A[maxr][i]); if(!A[r][c]) continue; for(int i=r+1;i =0;i--){ x[i] = A[i][m]; for(int j=i+1;j =0&&a =0&&b
转载于:https://www.cnblogs.com/zufezzt/p/5712819.html