传送门
大意:有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;
}