bzoj 4589 Hard Nim——FWT

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=4589

一开始异或和为0的话先手必败。有 n 堆，每堆可以填那些数，求最后异或和为0的方案数，就是一个快速幂的异或FWT。

注意快速幂的过程中对那些数组直接乘就行，不用总是FWT！！！

为什么比Zinn慢了1008ms？

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int N=5e4+,mod=1e9+;
int n,m,pri[],cnt,len,inv;
int f[N<<],a[N<<],b[N<<];
bool vis[N];
int pw(int x,int k)
{int ret=;while(k){if(k&)ret=(ll)ret*x%mod;x=(ll)x*x%mod;k>>=;}return ret;}
void fwt(int *a,bool fx)
{
  for(int R=;R<=len;R<<=)
    for(int i=,m=R>>;i<len;i+=R)
      for(int j=;j<m;j++)
    {
      int x=a[i+j]+a[i+m+j],y=a[i+j]+mod-a[i+m+j];
      a[i+j]=(ll)x*(fx?inv:)%mod;
      a[i+m+j]=(ll)y*(fx?inv:)%mod;
    }
}
void init()
{
  for(int i=;i<=;i++)
    {
      if(!vis[i])pri[++cnt]=i;
      for(int j=;j<=cnt&&i*pri[j]<=;j++)
    {
      vis[i*pri[j]]=;
      if(i%pri[j]==)break;
    }
    }
  for(int i=;i<=cnt;i++)f[pri[i]]=;
  inv=pw(,mod-);
}
int main()
{
  init();
  while(scanf("%d%d",&n,&m)==)
    {
      for(len=;len<=m;len<<=);
      a[]=;for(int i=;i<len;i++)a[i]=;
      for(int i=;i<=m;i++)b[i]=f[i];
      for(int i=m+;i<len;i++)b[i]=;
      fwt(a,); fwt(b,);
      while(n)
    {
      if(n&)
        for(int i=;i<len;i++)a[i]=(ll)a[i]*b[i]%mod;
      for(int i=;i<len;i++)b[i]=(ll)b[i]*b[i]%mod;
      n>>=;
    }
      fwt(a,);
      printf("%d\n",a[]);
    }
  return ;
}