POJ 2947 Widget Factory <高斯消元同余线性方程>

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

题意：一个部件工厂要生产一批不同种类的零件，每种零件生产时间（3<=day<=9）,给你n种零件和m工人生产数据，求每个零件生产时间？

分析：高斯消元同余线性方程。

有n种装饰物，m个已知条件,每个已知条件的描述如下：
p start end
a1,a2......ap(1<=ai<=n)
第一行表示从星期start到星期end一共生产了p件装饰物(工作的天数为end-start+1+7*x,加7*x是因为它可能生产很多周)，第二行表示这p件装饰物的种类(可能出现相同的种类，即ai=aj)。规定每件装饰物至少生产3天，最多生产9天。问每种装饰物需要生产的天数。如果没有解，则输出“Inconsistentdata.”，如果有多解，则输出“Multiple solutions.”，如果只有唯一解，则输出每种装饰物需要生产的天数。

 

Day = [(strat – end+ 1)+7]%7

思路：依旧是高斯消元。这里的m个已知条件就是m个方程组，n为n个未知数。设生产第i个装饰物需要的天数为x(i)，生产第i个装饰物的件数是A(i)件，那么对于每一个方程有x(1)*A(1)+x(2)*A(2)  + ... x(n)*A(n)= end-start+1+7*x。

代码下次附上：

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

using namespace std;
const int maxn = 305;
#define INF 0x3f3f3f3f
int matrix[maxn][maxn];
int X[maxn];
int n,m;

int gcd(int a,int b)
{
    return b?gcd(b,a%b):a;
}
int change(char s[7])
{
    if(strcmp(s,"MON")==0)
        return 1;
    else if(strcmp(s,"TUE")==0)
        return 2;
    else if(strcmp(s,"WED")==0)
        return 3;
    else if(strcmp(s,"THU")==0)
        return 4;
    else if(strcmp(s,"FRI")==0)
        return 5;
    else if(strcmp(s,"SAT")==0)
        return 6;
    else return 7;
}
int LCM(int a,int b)
{
    return a*b/gcd(a,b);
}
int Guass()
{
    int i,j,k,col;
    memset(X,0,sizeof(X));
    for (k = 0,col = 0; k < m && col < n; ++k, ++col)
    {
        int max_r = k;
        for (i = k + 1; i < m; ++i)
        {
            if (abs(matrix[i][col]) > abs(matrix[max_r][col])) max_r = i;
        }
        if (max_r != k) //½»»»
        {
            for (i = col; i < n + 1; ++i) swap(matrix[k][i],matrix[max_r][i]);
        }

        if (matrix[k][col] == 0)
        {
            k--;
            continue;
        }
        for (i = k + 1; i < m; ++i)
        {
            if (matrix[i][col] != 0)
            {
                int lcm = LCM(matrix[k][col],matrix[i][col]);
                int ta = lcm/abs(matrix[i][col]);
                int tb = lcm/abs(matrix[k][col]);
                if (matrix[i][col]*matrix[k][col] < 0) tb = -tb;
                for (j = col; j < n + 1; ++j)
                {
                    matrix[i][j] =((ta*matrix[i][j] - tb*matrix[k][j])%7+7)%7;
                }
            }
        }
    }
    //printf("k = %d col = %d matrix = %d\n",k,col,matrix[k][col]);
    for (i = k; i < m; ++i)
    {
        if (matrix[i][col] != 0) return -1;
    }
    if (k < n)
    {
        return n - k;
    }
    for (i = k - 1; i >= 0; --i)
    {
        int tmp = matrix[i][n];
        for (j = i + 1; j < n; ++j)
        {
            tmp =((tmp- matrix[i][j]*X[j])%7+7)%7;
        }
        while(tmp%matrix[i][i]!=0) tmp+=7;
        X[i] = ((tmp/matrix[i][i])%7+7)%7;
        if(X[i]<3)X[i]+=7;
    }
    return 0;
}

int main()
{
    char s1[7],s2[7];
    int d,num;
    while(scanf("%d %d",&n,&m)&&(n&&m))
    {
        memset(matrix,0,sizeof(matrix));
        for(int i=0;i<m;i++)
        {
            scanf("%d %s %s",&d,s1,s2);
            matrix[i][n] = change(s2) - change(s1) + 1;
            matrix[i][n] = (matrix[i][n]%7 + 7)%7;
            for(int j=0;j<d;j++)
            {
                scanf("%d",&num);
                matrix[i][num-1] ++;
            }
            for(int j =0;j<n;j++)
                matrix[i][j] = (matrix[i][j]%7+7)%7;
        }
        int flag = Guass();
         if(flag==-1)
        {
            puts("Inconsistent data.");
        }
        else if(flag!=0)puts("Multiple solutions.");
        else
        {
            for(int i=0;i<n;i++)
            {
                if(i!=0)printf(" ");
                printf("%d",X[i]);
            }
           printf("\n");
        }

    }
    return 0;
}