[poj2947][高斯消元]Widget Factory

传送门：
http://poj.org/problem?id=2947

这个代码是错的，我也不是很清楚为啥错了，高消的部分被我拿去交seti那题过了。。可能系数求错了？如果有人看我这个博文发现出了偏差请帮忙指正一下谢谢了。

sol:

题目就是给出了一系列同余方程组，解方程即可。

#include<iostream>
#include<algorithm>
#include<string>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<cmath>
using namespace std;
typedef long long ll;
typedef double s64;
int n,m;
inline int read()
{
    char c;
    int res,flag=0;
    while((c=getchar())>'9'||c<'0') if(c=='-')flag=1;
    res=c-'0';
    while((c=getchar())>='0'&&c<='9') res=(res<<3)+(res<<1)+c-'0';
    return flag?-res:res;
}
const int N=310;
char sr[N];
int c[N][N];
const int pyz=7;
inline int ksm(int s,int t)
{
    int res=1;
    while(t)
    {
        if(t&1) res=(ll)res*s%pyz;
        s=(ll)s*s%pyz;
        t>>=1;
    }
    return res;
}
int ans;
inline void guass()
{
    for(int i=1;i<=n;++i)
    for(int j=1;j<=m;++j)
    c[i][j]%=pyz;
    for(int i=1;i<m;++i)
    {
        int row;
        for(row=i;row<=n;++row) if(c[row][i]) break;
        if(row>n)
        {
            ans=1;
            continue;
        }
        if(row!=i)
            for(int j=1;j<=m;++j) swap(c[row][j],c[i][j]);
        int inv=ksm(c[i][i],5);
        for(int j=1;j<=m;++j) c[i][j]=c[i][j]*inv%pyz;
        for(int j=1;j<=n;++j)
        if(c[j][i]&&i!=j)
        {
            int tmp=c[j][i];
            for(int k=1;k<=m;++k)
            c[j][k]=(c[j][k]-c[i][k]*tmp%pyz+pyz)%pyz;
        }
    }
/*  for(int i=1;i<=n;++i)
    {
        for(int j=1;j<=m;++j)
        printf("%d ",c[i][j]);
        printf("\n");
    }*/
    for(int i=1;i<m;++i)
    if(!c[i][i]&&c[i][m])
    {
        ans=-1;
        return;
    }
}
inline int get(char *s)
{
    if(s[0]=='M') return 1;
    if(s[0]=='T'&&s[1]=='U') return 2;
    if(s[0]=='W') return 3;
    if(s[0]=='T'&&s[1]=='H') return 4;
    if(s[0]=='F') return 5;
    if(s[0]=='S'&&s[1]=='A') return 6;
    if(s[0]=='S'&&s[1]=='U') return 0;
}
int main()
{
//  freopen("2947.in","r",stdin);
//  freopen(".out","w",stdout);
    while(scanf("%d%d",&m,&n)!=EOF)
    {
        if(!n&&!m) break;
        ++m;
        ans=0;
        memset(c,0,sizeof(c));
        for(int i=1;i<=n;++i)
        {
            int k,a,b;
            k=read();
            scanf("%s",sr);
            a=get(sr);
            scanf("%s",sr);
            b=get(sr);
            c[i][m]=(b-a+1+pyz)%pyz;
            for(int j=1;j<=k;++j)
            c[i][read()]++;
        }
        guass();
        if(ans==-1) printf("Inconsistent data.\n");
        else if(ans==1) printf("Multiple solutions.\n");
        else
        {
            for(int i=1;i<m;++i)
            if(c[i][i]) printf("%d ",c[i][m]<3?c[i][m]+7:c[i][m]);
            else printf("7 ");
            printf("\n");
        }
    }
}