bzoj千題計劃308：bzoj4589: Hard Nim（倍增FWT+生成函數）

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

 

n*m*m 作法

dp[i][j] 前i堆石子，異或和爲j的方案數

 

 第一重循環能夠矩陣快速冪優化

後面求出序列的生成函數能夠FWT優化

 

作log次FWT也很慢（logn*m*logm）

兩個合併就是倍增FWT，即先對生成函數的序列作一次正變換，對正變換獲得的每一個結果快速冪，最後逆變換回去

時間複雜度O（logn*m+m*logm）

 

生成函數：是質數則係數爲1，不然爲0

 

#include<cstdio>
#include<cstring>

using namespace std;

#define N 50001

const int mod=1e9+7;

const int M=1<<16;

int inv;

int a[N];

int b[M+1];

void FWT(int *a,int n)
{
    int x,y;
    for(int d=1;d<n;d<<=1)
        for(int m=d<<1,i=0;i<n;i+=m)
            for(int j=0;j<d;++j)
            {
                x=a[i+j]; y=a[i+j+d];
                a[i+j]=x+y; a[i+j+d]=x-y;
                a[i+j]-=a[i+j]>=mod ? mod : 0;
                a[i+j+d]+=a[i+j+d]<0 ? mod : 0;
            }
}

void IFWT(int *a,int n)
{
    int x,y;
    for(int d=1;d<n;d<<=1)
        for(int m=d<<1,i=0;i<n;i+=m)
            for(int j=0;j<d;++j)
            {
                x=a[i+j]; y=a[i+j+d];
                a[i+j]=1LL*(x+y)*inv%mod;
                a[i+j+d]=1LL*(x-y+mod)%mod*inv%mod;
            }
}

int Pow(int a,int b)
{
    int res=1;
    for(;b;a=1LL*a*a%mod,b>>=1)
        if(b&1) res=1LL*res*a%mod;
    return res;
}

int main()
{
    for(int i=2;i<N;++i) a[i]=1;
    for(int i=2;i<N;++i)
        if(a[i])
            for(int j=2;i*j<N;++j) 
                a[i*j]=0;
    inv=Pow(2,mod-2);
    int n,m;
    int len;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        memset(b,0,sizeof(b));
        for(int i=1;i<=m;++i) b[i]=a[i];
        len=1;
        while(len<=m) len<<=1;
        FWT(b,len);
        for(int i=0;i<len;++i) b[i]=Pow(b[i],n);
        IFWT(b,len);
        printf("%d\n",b[0]);
    }
}