[hdu6316]Odd shops

记$m=10$，即商品的种类

记$g(x)=1+sum_{i=1}^{m}a_{i}x_{i}$，问题即求$f_{n}(x)=g^{n}(x)$非0项数（模2意义下）

注意到$f^{2}(x)equiv f(x^{2})(mod 2)$，这是因为如果所选的项在两边不同，那么交换后即会抵消

令$F(n,A)$为$A(x)f_{n}(x)$非0项数，即有$egin{cases}F(n+1,A)=F(n,Acdot g)\F(2n,A)=F(n,E)+F(n,O)end{cases}$

（其中$E(x)$和$O(x)$满足$A(x)=E(x^{2})+xO(x^{2})$，即将$A(x)$按照奇偶次划分）

第一个式子显然成立，对于第二个式子，考虑有
$$
A(x)f_{2n}(x)=A(x)f_{n}^{2}(x)equiv E(x^{2})f_{n}(x^{2})+xO(x^{2})f_{n}(x^{2})(mod 2)
$$
注意到$E(x^{2})f_{n}(x^{2})$和$xO(x^{2})f_{n}(x^{2})$两者没有重复项，因此直接将两者的非0项数相加即可，同时显然两者的非0项数等于$E(x)f_{n}(x)$和$O(x)f_{n}(x)$的非0项数，递归即可

对于$F(n,A)$，注意到$n$有$o(log n)$种，$A$最高次数不超过10（归纳即可），在模2意义下只有$2^{m+1}$种，总状态数即为$o(2^{m+1}log n)$，直接递归计算即可

另外，关于$A$的存储可以使用int来表示，那么两数（多项式）相乘复杂度为$o(m)$

总复杂度为$o(m2^{m+1}log n)$，可以通过

#include<bits/stdc++.h>
using namespace std;
#define mod 998244353
int n,x,g,sum[(1<<11)],f[40][(1<<11)];
int mul(int x,int y){
    int ans=0;
    for(int i=0;i<=10;i++)
        if (x&(1<<i))ans^=(y<<i);
    return ans;
}
int dfs(int n,int a,int s){
    if (f[s][a]>=0)return f[s][a];
    if (!n)return sum[a];
    int aa=a;
    if (n&1)aa=mul(a,g);
    int E=0,O=0;
    for(int i=0;i<=20;i++)
        if (aa&(1<<i)){
            if (i&1)O|=(1<<(i>>1));
            else E|=(1<<(i>>1));
        }
    return f[s][a]=(dfs((n>>1),E,s+1)+dfs((n>>1),O,s+1))%mod;
}
int main(){
    for(int i=0;i<(1<<11);i++)sum[i]=sum[(i>>1)]+(i&1);
    while (scanf("%d",&n)!=EOF){
        g=1;
        for(int i=1;i<=10;i++){
            scanf("%d",&x);
            if (x&1)g|=(1<<i);
        }
        memset(f,-1,sizeof(f));
        printf("%d
",dfs(n,1,0));
    }
}