Codeforces Round #448 (Div. 2) C. Square Subsets（状压）

题目链接：http://codeforces.com/contest/895/problem/C

70以内就19个素数，直接状压，0表示出现偶数次，1表示出现奇数次，然后对于每个数枚举选取奇数次还是偶数次，T了一发随便改改就过了，其实好像可以处理出2的幂次来搞啊？貌似那样就是稳定7e7的复杂度哦？绝望。。。

代码：

#include<bits/stdc++.h>using namespace std;typedef long long ll;const int MOD=1e9+7;const int MAXN=(1<<20)+5;int prime[100],pool[100];//保存素数 bool vis[100];//初始化 int cnt=0;ll qpow(ll a,ll b){ll ans=1;a%=MOD;    for(ll i=b;i;i>>=1,a=a*a%MOD)        if(i&1)ans=ans*a%MOD;    return ans;}void Prime(int n=MAXN){memset(vis,0,sizeof(vis));for(int i=2;i<n;i++){if(!vis[i])prime[cnt++]=i;for(int j=0;j<cnt&&i*prime[j]<n;j++){vis[i*prime[j]]=1;if(i%prime[j]==0)//关键 break;}}for(int i=1;i<=70;i++){int x=i;for(int j=0;j<cnt&&prime[j]<=x;j++){while(x%prime[j]==0){pool[i]^=(1<<j);x/=prime[j];}}}}int num[100];ll dp[2][MAXN];int main(){//freopen("in.txt","r",stdin);//freopen("out.txt","w",stdout);Prime(70);int n;scanf("%d",&n);for(int i=1;i<=n;i++){int x;scanf("%d",&x);num[x]++;}int now=0,nxt=1;dp[now][0]=1;for(int i=1;i<=70;i++){if(num[i]==0)continue;memset(dp[nxt],0,sizeof(dp[nxt]));for(int j=0;j<MAXN;j++){if(dp[now][j]){dp[nxt][j]=(dp[nxt][j]+qpow(2,num[i]-1)*dp[now][j])%MOD;dp[nxt][j^pool[i]]=(dp[nxt][j^pool[i]]+qpow(2,num[i]-1)*dp[now][j])%MOD;}}now^=1;nxt^=1;}printf("%lld\n",dp[now][0]-1);return 0;}