AcWing 227. 小部件厂 （高斯消元）打卡

题目：https://www.acwing.com/problem/content/description/229/

题意：有很多个零件，每个零件的生产时间都在3-9天之间，现在只知道每个工人的生产部件有哪些，还有生产日期的星期几和完成日期的星期几，求每个部件的具体生产日期

思路：首先我们根据两个星期，我们实际上可以计算出具体时间，然后我们到每个工人生产零件的数目还有时间，相当于我们可以建立这么多个方程组，然后我们就可以消元得出答案

注意处理3-9天范围的细节即可

#include<bits/stdc++.h>
using namespace std;
#define N 305
char day[][]={"","MON","TUE","WED","THU","FRI","SAT","SUN"};
int n,m,k,mat[N][N],ans[N];

int gcd(int a,int b)
{
    if(b==)
        return a;
    else return gcd(b,a%b);
}
bool mul_solu;
bool gauss(int mat[N][N],int m,int n)
{
    mul_solu=false;
    int row,col;
    for(row=,col=;row<m && col<n;++row,++col)
    {
        int p=row;
        for(int j=row+;j<m;++j)
            if(abs(mat[j][col])>abs(mat[p][col]))
                p=j;
        if(p!=row)
        {
            for(int j=;j<=n;++j)
                swap(mat[row][j],mat[p][j]);
        }
        if(mat[row][col]==) // 最大的都为0，说明这一列以下全是0
        {
            row--;
            continue;
        }
        for(int j=row+;j<m;++j)
        {
            int gg=gcd(mat[row][col],mat[j][col]);
            int muli=mat[j][col]/gg%;
            int mulj=mat[row][col]/gg%;
            for(int k=col;k<=n;++k)
            {
                mat[j][k]=mat[j][k]*mulj-mat[row][k]*muli;
                mat[j][k]=(mat[j][k]%+)%;
            }
        }
    }
    for(int i=row;i<m;++i) // inconsistent必须先与multiple判断
        if(mat[i][n])
            return false;
    if(row<n) // free variable 少
    {
        mul_solu=true;
        return true;
    }

    for(int i=n-;i>=;--i)
    {
        bool flag=false;
        for(int j=;j<=;++j)
        {
            int rr=;
            for(int k=i+;k<n;++k)
                rr+=mat[k][n]*mat[i][k]%,rr%=;
            rr+=j*mat[i][i];
            if(rr%==mat[i][n])
            {
                mat[i][n]=j;
                flag=true;
                break;
            }
        }
        if(!flag)
            return false;
    }

    return true;
}

int get_time(char *t1,char *t2)
{
    int st,ed;
    for(int i=;i<=;++i)
    {
        if(strcmp(t1,day[i])==)
            st=i;
        if(strcmp(t2,day[i])==)
            ed=i;
    }
    if(ed-st>=)
        return ed-st+;
    else return ed+-st;
}
char t1[],t2[];
int tt;
int main ()
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        if(n== && m==) break;
        memset(mat,,sizeof(mat));
        for(int i=;i<m;++i)
        {
            scanf("%d%s%s",&k,t1,t2);
            mat[i][n]=get_time(t1,t2)%;
            for(int j=;j<=k;++j)
            {
                scanf("%d",&tt);
                mat[i][tt-]++;
                mat[i][tt-]%=;
            }
        }
        if(gauss(mat,m,n))
        {
            if(mul_solu)
            {
                printf("Multiple solutions.\n");
            }
            else
            {
                for(int i=;i<n-;++i)
                    printf("%d ",mat[i][n]);
                printf("%d\n",mat[n-][n]);
            }
        }
        else printf("Inconsistent data.\n");
    }
    return ;
}