HDU 5833

2016CCPC网赛第二题。现场没出，好可惜。。

当时找到书上的板子了，不敢抄怕查重，然后就手写了一个高斯消元，数据量也挺小没用模版，直接模拟的，结果最后一直wa，找了一晚上原因竟然是数组开小了。。数组开小了，一直wa。。。不明嚼栗。。

蓝书161页原题。。

附代码：

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
using namespace std;
string s[3005];
int p[3005];
int nn = 0;
int a[3005];
int n;
int m;
const int MOD = 1000000007;
void init()
{
memset(p, 0, sizeof(p));
for(int i = 2; i <= 2000; ++i)
{

if(!p[i])
{
for(int j = i + i; j <= 2000; j += i)
{
p[j] = 1;
}
}
}
for(int i = 2; i <= 2000; ++i)
{
if(!p[i])
{
a[nn++] = i;
}
}
for(int i = 0; i <= 2001; ++i)
for(int j = 0; j <= 305; ++j)
s[i].push_back('0');

}
void solve(int j, long long x)
{
int i = 0;
while(i < nn)
{
if(x % a[i] == 0)
{
x /= a[i];
int tmp = s[i][j] - '0';
tmp = !tmp;
s[i][j] = '0' + tmp;
}
else
{
++i;
if(x == 1)
{
m = max(m, i);
return ;
}
}
}
}
bool cmp(const string &t1, const string &t2)
{
return t1 > t2;
}
void Xor(int ii, int res, int tar)
{
for(int i = ii; i < m; ++i)
{
s[tar][i] = '0' + ((s[tar][i] - '0') ^ (s[res][i] - '0'));
}
}
int Rank(int m, int n)
{
int i = 0, j = 0, k, r, u;
while(i < m && j < n)
{
r = i;
for(k = i; k < m; ++k)
if(s[k][j] == '1')
{
r = k;
break;
}
if(s[r][j] == '1')
{
if(r != i)
for(k = 0; k <= n; ++k)
{
swap(s[r][k], s[i][k]);
}
for(u = i + 1; u < m; ++u)
if(s[u][j] == '1')
{
for(k = i; k <= n; ++k)
{
s[u][k] = '0' + ((s[u][k] - '0') ^ (s[i][k] - '0'));
}
}
++i;
}
++j;
}
return i;
}
int main()
{
int T;
init();
scanf("%d", &T);
int kase = 0;
while(T--)
{
m = 0;
for(int i = 0; i <= 2001; ++i)
for(int j = 0; j <= 305; ++j)
s[i][j] = '0';
scanf("%d", &n);
for(int i = 0; i < n; ++i)
{
long long x;
scanf("%I64d", &x);
solve(i, x);
}
swap(n, m);
sort(s, s + n, cmp);
int ii = 0;
int ans;
// cout << "----------------" << endl;
while(true)
{
int flag = 0;
int jj;
for(int j = 0; j < m; ++j)
if(s[ii][j] != '0')
{
flag = 1;
jj = j;
break;
}
if(!flag || ii == n)
{
ans = ii;
break;
}
for(int i = ii + 1; i < n; ++i)
if(s[i][jj] != '0')
{
Xor(jj, ii, i);
}
++ii;
sort(s + ii, s + n, cmp);
// for(int i = 0; i < n; ++i)
// {
// for(int j = 0; j < m; ++j)
// cout << s[i][j] << " ";
// cout << endl;
// }
// cout << "----------------------------" << endl;
}
long long aa = 1;
for(int i = 1; i <= (m - ans); ++i)
{
aa *= 2;
aa %= MOD;
}
aa = (aa - 1 + MOD) % MOD;
printf("Case #%d:\n", ++kase);
printf("%I64d\n", aa);
}
}