bzoj 4589 Hard Nim

http://www.elijahqi.win/archives/3544
Description

Claris和NanoApe在玩石子游戏，他们有n堆石子，规则如下：
1. Claris和NanoApe两个人轮流拿石子，Claris先拿。
2. 每次只能从一堆中取若干个，可将一堆全取走，但不可不取，拿到最后1颗石子的人获胜。
不同的初始局面，决定了最终的获胜者，有些局面下先拿的Claris会赢，其余的局面Claris会负。
Claris很好奇，如果这n堆石子满足每堆石子的初始数量是不超过m的质数，而且他们都会按照最优策略玩游戏，那么NanoApe能获胜的局面有多少种。
由于答案可能很大，你只需要给出答案对10^9+7取模的值。

Input
输入文件包含多组数据，以EOF为结尾。
对于每组数据：
共一行两个正整数n和m。
每组数据有1<=n<=10^9, 2<=m<=50000。
不超过80组数据。

Output

Sample Input
3 7
4 13
Sample Output
6
120
HINT

Source
Topcoder SRM 518 Div1 Hard Nim By Tangjz

设生成函数 每个素数的位置都给1 那么最后的答案就相当于是给 这个生成函数在异或下做n次卷积 最后求0位置上的值就好 考虑直接做会T 其实我们可以求出正变换序列之后直接对这个序列进行快速幂即可 并不需要每次fwt回来再搞

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int mod=1e9+7;
inline int ksm(ll b,int t){static ll tmp;
    for (tmp=1;t;b=b*b%mod,t>>=1) if (t&1) tmp=tmp*b%mod;return tmp;
}
inline int inc(int x,int v){return x+v>mod?x+v-mod:x+v;}
inline int dec(int x,int v){return x-v<0?x-v+mod:x-v;}
const int N=50010;
int inv2,a[N<<2],prime[N],tot,n,m,len;bool not_prime[N];
inline void fwt(int *x,int f){
    for (int i=1;i<len;i<<=1){
        for (int j=0;j<len;j+=i<<1){
            for (int k=0;k<i;++k){
                int t1=x[j+k],t2=x[j+k+i];
                if (f==1){
                    x[j+k]=inc(t1,t2);x[j+k+i]=dec(t1,t2);
                }else{
                    x[j+k]=(ll)inc(t1,t2)*inv2%mod;x[j+k+i]=(ll)dec(t1,t2)*inv2%mod;
                }
            }
        }
    }
}
inline void init(){not_prime[1]=1;
    for (int i=2;i<=5e4;++i){
        if (!not_prime[i]) prime[++tot]=i;
        for (int j=1;prime[j]*i<=5e4;++j){
            not_prime[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }
}
int main(){
    freopen("bzoj4589.in","r",stdin);
    init();inv2=ksm(2,mod-2);
    while(~scanf("%d%d",&n,&m)){
        memset(a,0,sizeof(a));for (len=1;len<=m;len<<=1);
        for (int i=1;i<=m;++i) if (!not_prime[i]) a[i]=1;
        fwt(a,1);
        for (int i=0;i<len;++i) a[i]=ksm(a[i],n);
        fwt(a,-1);
        printf("%d\n",a[0]);
    }
    return 0;
}