POJ2947 Widget Factory（高斯消元模板题）

传送门
大意：有n个变量，m个等式，求在 模7 意义下的解。
其实就是模板题，消完元之后如果有自由元就是有多个解，如果消完元之后有等式两边不等的情况，就是无解，否则有解。
代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <assert.h>
using namespace std;
#define MAXN 305
const int MOD = 7;
int n, m, a[MAXN][MAXN], ans[MAXN];
int inv(int n) {
    if(n == 1) return 1;
    return (MOD-MOD/n) * inv(MOD%n) % MOD;
}
int gcd(int a, int b) { return b ? gcd(b, a%b) : a; }
int Guass() {
    int i, j, k, l;
    for(i = 0, j = 0; i < n && j < m; ++ i, ++ j) {
        for(k = i; k < n; ++ k) if(a[k][j]) break;
        if(k == n) { -- i; continue; }
        if(i != k) for(int l = 0; l <= m; ++ l) swap(a[k][l], a[i][l]);
        for(k = 0; k < n; ++ k) {
            if(!a[k][j] || i == k) continue;
            int t = a[k][j] * a[i][j] / gcd(a[k][j], a[i][j]);
            int t1 = t / a[k][j], t2 = t / a[i][j];
            for(l = 0; l <= m; ++ l)
                a[k][l] = ((a[k][l] * t1 - a[i][l] * t2) % MOD + MOD) % MOD;
        }
    }
    for(k = i; k < n; ++ k) if(a[k][m]) return puts("Inconsistent data."), 0;
    if(i < m) return puts("Multiple solutions."), 0;
    for(i = m-1; ~i; -- i) {
        ans[i] = a[i][m]*inv(a[i][i]%7)%MOD;
        if(ans[i] < 3) ans[i] += 7;
    }
    for(int i = 0; i < m; ++ i)
        printf("%d%c", ans[i], i == m-1 ? '\n' : ' ');
    return 0;
}
inline int idx(char *s) {
    if(s[0] == 'M') return 0;
    if(s[0] == 'T' && s[1] == 'U') return 1;
    if(s[0] == 'W') return 2;
    if(s[0] == 'T' && s[1] == 'H') return 3;
    if(s[0] == 'F') return 4;
    if(s[0] == 'S' && s[1] == 'A') return 5;
    return 6;
}
int main() {
    while(scanf("%d%d", &m, &n), n+m) {
        memset(a, 0, sizeof a); memset(ans, 0, sizeof ans);
        int k; char s1[10], s2[10];
        for(int i = 0; i < n; ++ i) {
            scanf("%d%s%s", &k, s1, s2);
            int p = (idx(s2) - idx(s1) + 1) % MOD;
            if(p < 0) p += MOD; a[i][m] = p;
            while(k --) { scanf("%d", &p); (++ a[i][p-1]) %= MOD; }
        } Guass();
    }
    return 0;
}