pku2947 Widget Factory

题目链接：http://162.105.81.212/JudgeOnline/problem?id=2947

题意简述：题目长，好难看。大致就是n个工人，从某个星期的x-day到某个星期的y-day，做了哪几个零件，求每个零件需要的时间。

解题思路：高斯消元。看到x和y ms 不好连接起来，但是它需要的天数一定是y-x+7*k，所有系数对7取模，因为这不会影响结果，这个可以理解一下。

代码：

#include<stdio.h>#include<string.h>#include<stdlib.h>using namespace std;#define swap(a,b) (a^=b^=a^=b)const int M=7;const int Len=320;const char S[8][4]={"","MON","TUE","WED","THU","FRI","SAT","SUN"};int map[Len][Len];int x[Len];int ee,vv;FILE *tp,*fp;inline int abs(int a){ if(a<0) return (-a); return a;}inline int gcd(int a,int b){ int temp; while(b!=0) { temp=b; b=a%b; a=temp; } return a;}inline int lcm(int a,int b){ return a*b/gcd(a,b);}void Gauss(){ int max_r; int col; int ta,tb; int LCM; int temp; col=1; int k; for(k=1;k<=ee&&col<=vv;++k,++col) { max_r=k; for(int i=k+1;i<=ee;++i) if(abs(map[i][col])>abs(map[max_r][col])) max_r=i; if(max_r!=k) for(int i=k;i<=vv+1;++i) swap(map[k][i],map[max_r][i]); if(map[k][col]==0) { k--; continue; } for(int i=k+1;i<=ee;++i) if(map[i][col]!=0) { LCM=lcm(abs(map[i][col]),abs(map[k][col])); ta=LCM/abs(map[i][col]),tb=LCM/abs(map[k][col]); if(map[i][col]*map[k][col]<0) tb=-tb; for(int j=col;j<=vv+1;++j) map[i][j]=((map[i][j]*ta)%7-(map[k][j]*tb)%7+7)%7; } } for(int i=k;i<=ee;++i) if(map[i][vv+1]%7!=0) { printf("Inconsistent data./n"); return; } if(k-1<vv) { printf("Multiple solutions./n"); return; } for(int i=vv;i>=1;--i) { temp=map[i][vv+1]; for(int j=i+1;j<=vv;++j) if(map[i][j]!=0) temp=temp-map[i][j]*x[j]; while(temp%map[i][i]) temp+=7; x[i]=temp/map[i][i]; if(x[i]<3) while(x[i]<3) x[i]+=7; if(x[i]>9) while(x[i]>9) x[i]-=7; } for(int i=1;i<=vv;++i) if(x[i]>9||x[i]<3) { printf("Inconsistent data./n"); return; } for(int i=1;i<=vv;++i) { printf("%d",x[i]); if(i<vv) printf(" "); } printf("/n");}int main(){ int n,m; while(scanf("%d%d",&n,&m)&&(n||m)) { memset(map,0,sizeof(map)); memset(x,0,sizeof(x)); int k; char s1[10],s2[10]; ee=m,vv=n; for(int i=1;i<=m;++i) { scanf("%d %s %s",&k,s1,s2); s1[3]=s2[3]='/0'; int num1,num2; for(int j=1;j<=7;++j) { if(strcmp(S[j],s1)==0) num1=j; if(strcmp(S[j],s2)==0) num2=j; } int value=(num2-num1+8)%7; map[i][vv+1]=value; for(int j=1;j<=k;++j) { scanf("%d",&num1); ++map[i][num1]; map[i][num1]%=7; } } Gauss(); } return 0;}