POJ 2947 Widget Factory 解题报告（高斯消元法）

    题目大意：某工厂，每个零件的加工时间是3天到9天。现在只有一些记录，记录中有开始是星期几，结束是星期几，加工了哪些零件，问你能不能求出加工每个零件所需要的天数。

    解题报告：高斯消元法解线性方程组。这题可以直接用模板做的。直接贴代码~

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;

const int maxn=400;
int a[maxn][maxn];
int x[maxn];
bool free_x[maxn];

int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}

int lcm(int a,int b)
{
    return a/gcd(a,b)*b;
}

int Guass(int equ,int var)
{
    memset(x,0,sizeof(x));

    int k=0;
    int col=0;
    for(;k<equ && col<var;k++,col++)
    {
        int max_r=k;
        for(int i=k+1;i<equ;i++)
            if(abs(a[i][col])>abs(a[max_r][col]))
                max_r=i;
        if(max_r!=k)
            for(int j=k;j<=var;j++)
                swap(a[k][j],a[max_r][j]);
        if(a[k][col]==0)
        {
            k--;
            continue;
        }
        for(int i=k+1;i<equ;i++) if(a[i][col])
        {
            int LCM=lcm(abs(a[i][col]),abs(a[k][col]));
            int ta =LCM/abs(a[i][col]);
            int tb =LCM/abs(a[k][col]);

            if(a[i][col]*a[k][col]<0)
                tb=-tb;

            for(int j=col;j<=var;j++)
                a[i][j]=((a[i][j]*ta-a[k][j]*tb)%7+7)%7;
        }
    }

    for(int i=k;i<equ;i++) if(a[i][col])
        return -1;

    if(k<var)
        return var-k;

    for(int i=var-1;i>=0;i--)
    {
        int temp=a[i][var];
        for(int j=i+1;j<var;j++) if(a[i][j])
        {
            temp-=a[i][j]*x[j];
            temp=(temp%7+7)%7;
        }
        while(temp%a[i][i])
            temp+=7;
        x[i]=((temp/a[i][i])%7+7)%7;
    }
    return 0;
}

int trans(char* s)
{
    switch(s[0])
    {
    case 'M':   return 1;
    case 'T':   return s[1]=='U'?2:4;
    case 'W':   return 3;
    case 'F':   return 5;
    case 'S':   return s[1]=='A'?6:7;
    }
}

int main()
{
    int equ,var;
    while(~scanf("%d%d",&var,&equ) && equ|var)
    {
        memset(a,0,sizeof(a));

        for(int i=0;i<equ;i++)
        {
            int k;
            char s1[4],s2[4];
            scanf("%d%s%s",&k,s1,s2);
            a[i][var]=((trans(s2)-trans(s1)+1)%7+7)%7;
            while(k--)
            {
                int t;
                scanf("%d",&t);
                t--;
                a[i][t]++;
                a[i][t]%=7;
            }
        }
        int ans=Guass(equ,var);
        if(ans==0)
        {
            for(int i=0;i<var;i++) if(x[i]<=2)
                x[i]+=7;    // 3-9
            printf("%d",x[0]);
            for(int i=1;i<var;i++)
                printf(" %d",x[i]);
            puts("");
        }
        else if(ans<0)
            puts("Inconsistent data.");
        else
            puts("Multiple solutions.");
    }
}