bzoj4589 Hard Nim（博弈论+fwt+快速幂）

搞出生成函数，把它异或卷积n次即可，最后答案就是0的系数。
fwt之后直接快速幂即可。再ifwt回来。
复杂度$O(mlogm+mlogn)$

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 50010
#define mod 1000000007
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(T==S){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,prime[N>>2],tot=0,inv2=mod+1>>1;
bool notprime[N];
struct Icefox{
    int a[N<<1];
    inline void init(bool f){
        memset(a,0,sizeof(a));
        if(f) for(int i=0;i<n;++i) a[i]=1;
    }int& operator[](int x){return a[x];}
    friend Icefox operator*(Icefox a,Icefox b){
        Icefox res;
        for(int i=0;i<n;++i) res[i]=(ll)a[i]*b[i]%mod;return res;
    }
    friend Icefox operator^(Icefox a,int k){
        Icefox res;res.init(1);
        for(;k;k>>=1,a=a*a) if(k&1) res=res*a;return res;
    }
}a;
inline void init(){
    notprime[1]=1;
    for(int i=2;i<=50000;++i){
        if(!notprime[i]) prime[++tot]=i;
        for(int j=1;prime[j]*i<=50000;++j){
            notprime[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }
}
inline void fwt(int *a,int f){
    for(int i=1;i<n;i+=i)
        for(int j=0,p=i<<1;j<n;j+=p)
            for(int k=0;k<i;++k){
                int x=a[j+k],y=a[j+k+i];
                if(f==1){
                    a[j+k]=(x+y)%mod;
                    a[j+k+i]=(x-y+mod)%mod;
                }else{
                    a[j+k]=(ll)(x+y)%mod*inv2%mod;
                    a[j+k+i]=(ll)(x-y+mod)%mod*inv2%mod;
                }
            }
}
int main(){
//  freopen("a.in","r",stdin);
    init();
    while(~scanf("%d%d",&n,&m)){
        a.init(0);
        for(int i=1;i<=m;++i)
            if(!notprime[i]) a[i]=1;
        int k=n;for(n=1;n<=m;n+=n);fwt(a.a,1);
        a=a^k;fwt(a.a,-1);
        printf("%d\n",a[0]);
    }return 0;
}