fzu1704（高斯消元法解异或方程组+高精度输出）

题目链接：https://vjudge.net/problem/FZU-1704

题意：经典开关问题，求使得灯全0的方案数。

思路：题目保证至少存在一种方案，即方程组一定有解，那么套上高斯消元法的板子，求出自由变元的个数t，方案总数即2^t，t可能大于64，要用到高精度计算。

AC代码：

#include<iostream>
#include<string>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;

int compare(string str1,string str2){
    int len1=str1.length(),len2=str2.length();
    if(len1<len2)
        return -;
    else if(len1>len2)
        return ;
    else
        return str1.compare(str2);
}

//高精度加法，只能是两个正数相加
string add(string str1,string str2){
    string str;
    int len1=str1.length(),len2=str2.length();
    if(len1<len2)
        for(int i=;i<=len2-len1;i++)
            str1=""+str1;
    else
        for(int i=;i<=len1-len2;i++)
            str2=""+str2;
    int cf=,tmp;
    len1=str1.length();
    for(int i=len1-;i>=;i--){
        tmp=str1[i]-''+str2[i]-''+cf;
        cf=tmp/;
        tmp%=;
        str=char(tmp+'')+str;
    }
    if(cf)
        str=""+str;
    return str;
}

//高精度减法，只能是两个正数相减，而且要大减小
string sub(string str1,string str2){
    string str;
    int len1=str1.length(),len2=str2.length();
    for(int i=;i<=len1-len2;i++)
        str2=""+str2;
    int cf=;
    for(int i=len1-;i>=;i--){
        if(str1[i]-cf>=str2[i]){
            str=char(str1[i]-cf-str2[i]+'')+str;
            cf=;
        }
        else{
            str=char(str1[i]-cf-str2[i]++'')+str;
            cf=;
        }
    }
    str.erase(,str.find_first_not_of(''));
    if(str.empty())
        str="";
    return str;
}
//高精度乘法，只能是两个正数相乘
string mul(string str1,string str2){
    string str;
    int len1=str1.length(),len2=str2.length();
    for(int i=len2-;i>=;i--){
        string tmpstr;
        int tmp=str2[i]-'',cf=,t;
        if(tmp){
            for(int j=;j<=len2--i;j++)
                tmpstr+="";
            for(int j=len1-;j>=;j--){
                t=(tmp*(str1[j]-'')+cf)%;
                cf=(tmp*(str1[j]-'')+cf)/;
                tmpstr=char(t+'')+tmpstr;
            }
            if(cf)
                tmpstr=char(cf+'')+tmpstr;
        }
        str=add(str,tmpstr);
    }
    str.erase(,str.find_first_not_of(''));  //删除前面多余的0，因为如果是0×10，结果将会是00
    if(str.empty())
        str="";
    return str;
}
//高精度除法，只能是两个正数相除
void div(string str1,string str2,string& con,string &rem){
    if(str2=="")
        return;
    else if(str1==""){
        con="",rem="";
        return;
    }
    else if(compare(str1,str2)<){
        con="",rem=str1;
        return;
    }
    else if(compare(str1,str2)==){
        con="",rem="";
        return;
    }
    else{
        int len1=str1.length(),len2=str2.length();
        string tmpstr;
        for(int i=;i<len2-;i++)
            tmpstr=tmpstr+str1[i];
        for(int i=len2-;i<len1;i++){
            tmpstr=tmpstr+str1[i];
            tmpstr.erase(,tmpstr.find_first_not_of(''));
            if(tmpstr.empty())
                tmpstr="";
            for(char j='';j>='';j--){
                string str,tmp;
                str=str+j;
                tmp=mul(str,str2);
                if(compare(tmp,tmpstr)<=){
                    con=con+j;
                    tmpstr=sub(tmpstr,tmp);
                    break;
                }
            }
        }
        rem=tmpstr;
    }
    con.erase(,con.find_first_not_of(''));
    if(con.empty())
        con="";
} 

bool JudgeZero(string s1){
    for(int i=;i<s1.length();++i)
        if(s1[i]!='') return false;
    return true;
}

string gcd(string a,string b){
    if(!JudgeZero(b)){
        string s1,s2;
        div(a,b,s1,s2);
        return gcd(b,s2);
    }
    return a;
}

string lcm(string a,string b){
    string t=gcd(a,b),s1,s2;
    a=mul(a,b);
    div(a,t,s1,s2);
    return s1;
}
const int maxn=;
int T,equ,var,a[maxn][maxn],x[maxn];

int Gauss(){
    int k=,max_r;
    for(int col=;k<equ&&col<var;++k,++col){
        max_r=k;
        for(int i=k+;i<equ;++i){
            if(abs(a[i][col])>abs(a[max_r][col]))
                max_r=i;
        }
        if(max_r!=k){
            for(int i=col;i<var+;++i)
                swap(a[max_r][i],a[k][i]);
        }
        if(!a[k][col]){
            --k;
            continue;
        }
        for(int i=k+;i<equ;++i){
            if(!a[i][col]) continue;
            for(int j=col;j<var+;++j)
                a[i][j]^=a[k][j];
        }
    }
    for(int i=k;i<equ;++i)
        if(a[i][var])
            return -;
    return var-k;
}

int main(){
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&equ,&var);
        for(int i=;i<equ;++i)
            for(int j=;j<var+;++j)
                a[i][j]=;
        for(int i=;i<var+;++i)
            x[i]=;
        for(int i=;i<equ;++i)
            scanf("%d",&a[i][var]);
        for(int i=;i<var;++i){
            int t1,t2;
            scanf("%d",&t1);
            while(t1--){
                scanf("%d",&t2);
                --t2;
                a[t2][i]=;
            }
        }
        int t=Gauss();
        string s1,s2;
        s1=s1+"",s2=s2+"";
        while(t--)
            s1=mul(s1,s2);
        cout<<s1<<endl;
    }
    return ;
}