zoukankan      html  css  js  c++  java
  • zz

    #include<iostream>
    #include<cstdio>
    #include<algorithm>
    #include<vector> 
    using namespace std;
    #define pp 100000007
    #define N 35 
    vector<int> pri;
    #define num(i,j) ((((i)-1)*n)+(j))
    int n,m,c[N][N],Map[5000][N][N],sz;
    int mat[5065][N*N];
    long long ans;
    void prepare()
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
            {
                scanf("%d",&c[i][j]);
                int x=c[i][j];
                for(int k=2;k*k<=n;k++)
                if(x%k==0)
                {
                    pri.push_back(k);
                    while(x%k==0) x/=k;
                }
                if(x>1)
                pri.push_back(x);
            }
        sort(pri.begin(),pri.end());
        pri.erase(unique(pri.begin(),pri.end()),pri.end());
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
            {
                int x=c[i][j];
                sz=pri.size()-1;
                for(int k=0;k<=sz;k++)
                    while(x%pri[k]==0) Map[k][i][j]^=1,x/=pri[k];
            }
    }
    void makematrix()
    {
        for(int i=1;i<=n;i++)
        {
            mat[i][0]=1;
            for(int j=1;j<=n;j++) mat[i][num(i,j)]=1;
        }
        for(int i=1;i<=n;i++)
        {
            mat[n+i][0]=1;
            for(int j=1;j<=n;j++) mat[n+i][num(j,i)]=1;
        }
        m=n+n;
        for(int k=0;k<=sz;k++)
        {
            mat[++m][0]=0;
            for(int i=1;i<=n;i++)
                for(int j=1;j<=n;j++) mat[m][num(i,j)]=Map[k][i][j];
        }
    }
    void gauss()
    {
        int num=0;
        ans=1;
        int pos=1;
        for(int i=1;i<=n*n;i++)
        {
            bool q=0;
            for(int j=pos;j<=m;j++)
            {
                if(mat[i][j])
                {
                    swap(mat[pos],mat[j]),q=1;
                    break;
                }
                if(!q)num++;
                else pos++;
                for(int j=pos;j<=m;j++)
                {
                    if(mat[j][i])
                        for(int k=0;k<=n*n;k++)
                        mat[j][k]^=mat[pos-1][k]; 
                    }
                }
            }
            for(int i=1;i<=m;i++)
            {
                bool q=0;
                for(int j=1;j<=n*n;j++)
                    if(mat[i][j]) {q=1;break;}
                if(q) continue;
                if(mat[i][0]) ans=0;
            }
            for(int i=1;i<=num;i++)
                 ans=ans*2%pp;
        }
    }
    int main()
    {
        prepare();
        makematrix();
        gauss();
        printf("%lld
    ",ans);
     } 

    zz

  • 相关阅读:
    httpclient妙用一 httpclient作为客户端调用soap webservice(转)
    WebService学习总结(转)
    C++的override和final
    C++类const和static成员初始化
    C++数据存储方式
    C++类成员存储大小
    内联函数
    C++接口的概念
    C++深拷贝和浅拷贝
    C++构造函数以及何时被调用
  • 原文地址:https://www.cnblogs.com/dancer16/p/7571300.html
Copyright © 2011-2022 走看看