poj2947Widget Factory

对于同余方程的高斯消元啊。

其实也差不多吧。先同一位通分，然后减一下就好了。

主要是判无解和多解的麻烦，需要注意即使有自由元也可能先无解

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>

#define MOD(x) (x%7+7)%7
using namespace std;

int n,m;
char ss[];
int scs()
{
    scanf("%s",ss+);
         if(ss[]=='M')return ;
    else if(ss[]=='T'&&ss[]=='U')return ;
    else if(ss[]=='W')return ;
    else if(ss[]=='T')return ;
    else if(ss[]=='F')return ;
    else if(ss[]=='S'&&ss[]=='A')return ;
    else return ;
}

int gcd(int a,int b)
{
    if(a==)return b;
    return gcd(b%a,a);
}
int exgcd(int a,int b,int &x,int &y)
{
    if(a==)
    {
        x=,y=;
        return b;
    }
    else
    {
        int tx,ty;
        int d=exgcd(b%a,a,tx,ty);
        x=ty-b/a*tx;
        y=tx;
        return d;
    }
}

int qg[][],qc[];
int as[];
void gauss()
{
    int j,jj;
    for(j=,jj=;j<=n&&jj<=m;j++)
    {
        for(int i=jj;i<=m;i++)
            if(qg[i][j]!=)
            {
                for(int k=j;k<=n;k++)swap(qg[i][k],qg[jj][k]);
                swap(qc[i],qc[jj]);
                break;
            }
        if(qg[jj][j]==)continue;

        for(int i=;i<=m;i++)
        {
            if(i==jj)continue;

            int gg=gcd(qg[i][j],qg[jj][j]);
            int ml1=qg[jj][j]/gg,ml2=qg[i][j]/gg;

            for(int k=;k<=n;k++)
                qg[i][k]=MOD(qg[i][k]*ml1-qg[jj][k]*ml2);
            qc[i]=MOD(qc[i]*ml1-qc[jj]*ml2);
        }
        jj++;
    }
    for(int i=jj;i<=m;i++)
        if(qc[i]>){printf("Inconsistent data.\n");return ;}
    if(jj!=j||j<=n){printf("Multiple solutions.\n");return ;}

    for(int j=;j<=n;j++)
    {
        int A=qg[j][j],B=,K=qc[j],x,y;
        int d=exgcd(A,B,x,y);

        x=(x*(K/d)%(B/d)+(B/d))%(B/d);
        if(x<)x+=;
        as[j]=x;
    }
    for(int j=;j<n;j++)printf("%d ",as[j]);
    printf("%d\n",as[n]);
}

int v[];
int main()
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        if(n==&&m==)break;

        for(int i=;i<=m;i++)
        {
            int k,x;
            scanf("%d",&k);int s1=scs(),s2=scs();
            qc[i]=MOD(s2-s1+);

            memset(v,,sizeof(v));
            for(int j=;j<=k;j++)
                scanf("%d",&x), v[x]++;

            for(int j=;j<=n;j++)qg[i][j]=MOD(v[j]);
        }

        gauss();
    }
    return ;
}