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

时空隧道

---

题意：
有n种工作，m个人，每个人在某段时间中完成k个任务，求出n个任务的单个完成时间…（每个任务完成时间在[3,7]区间中）…
如果无解输出Inconsistent data.
如果有无穷多个解输出Inconsistent data.

---

分析：
很容易就可以得出方程….然后消元就好了…

---

代码如下：

#include<algorithm>#include<iostream>#include<cstring>#include<cstdio>#include<map>//by NeighThornusing namespace std;const int maxn=300+5,MOD=7;int n,m,coe[maxn][maxn];map<string,int> mp;string str[2];inline int gauss(void){    int i,j;    for(i=1,j=1;i<=m&&j<=n;j++){//i枚举第几个方程，j枚举第几个主元         int k=i;        while(!coe[k][j]&&k<=m)            k++;        if(coe[k][j]){//消元             for(int s=1;s<=n+1;s++)                swap(coe[k][s],coe[i][s]);            for(int l=i+1;l<=m;l++)                if(coe[l][j]){                    int x=coe[i][j],y=coe[l][j];                    for(int s=j;s<=n+1;s++)                        coe[l][s]=((coe[l][s]*x-coe[i][s]*y)%MOD+MOD)%MOD;                }            i++;        }    }    for(int x=i;x<=m;x++)        if(coe[x][n+1])            return -1;    if(i<=n)        return n-i+1;    for(int x=n;x>=1;x--){        for(int y=x+1;y<=n;y++)            if(coe[x][y])                coe[x][n+1]=((coe[x][n+1]-(coe[y][n+1]*coe[x][y]%MOD))%MOD+MOD)%MOD;        while(coe[x][n+1]%coe[x][x]!=0)            coe[x][n+1]+=MOD;        coe[x][n+1]=(coe[x][n+1]/coe[x][x])%MOD;    }    return 0;}signed main(void){    mp["MON"]=1,mp["TUE"]=2,mp["WED"]=3,mp["THU"]=4,mp["FRI"]=5,mp["SAT"]=6,mp["SUN"]=7;    while(scanf("%d%d",&n,&m)&&!(!n&&!m)){        memset(coe,0,sizeof(coe));        for(int i=1,num;i<=m;i++){            cin>>num>>str[0]>>str[1];            coe[i][n+1]=(mp[str[1]]-mp[str[0]]+1+7)%7;            for(int j=1,x;j<=num;j++)                scanf("%d",&x),coe[i][x]++;            for(int j=1;j<=n;j++)                coe[i][j]%=7;        }        int lala=gauss();        if(lala==-1)            puts("Inconsistent data.");//没有解         else if(lala)            puts("Multiple solutions.");//无穷解         else{            for(int i=1;i<=n;i++)                if(coe[i][n+1]<3)                    coe[i][n+1]+=7;            for(int i=1;i<n;i++)                printf("%d ",coe[i][n+1]);            printf("%d\n",coe[n][n+1]);        }    }    return 0;}

---

by >_< NeighThorn