传送门: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;
}