[atARC084F]XorShift

如果异或变为加法和减法，那么根据扩欧，$k$合法当且仅当$k|\gcd_{i=1}^{n}a_{i}$

换一种方式定义约数：$x$是$y$的约数当且仅当存在$p_{i}\in \{0,1\}$使得$\sum_{i=0}^{\infty}2^{i}x=y$，那么类似的，再把加法改为异或，我们就得到了本题中关于约数的定义

如何求$d=\gcd(x,y)$：假设$x$的最高位为$2^{p}$，$y$的最高位为$2^{q}$（二进制下，且不妨假设$p\ge q$），那么有$d|x$和$d|2^{p-q}y$，又因为$d|x$，所以$d|(x\oplus 2^{p-q}y)$，即$\gcd(x,y)=\gcd(y,x-2^{p-q}y)$

对于求gcd的过程，每一次必然会使得最高位-1，可以通过bitset优化到$o(\frac{nL^{2}}{64})$

令$d=\gcd_{i=1}^{n}a_{i}$，考虑$k\le C$等价于$k\oplus C$的最高位上的1是$C$的1或$k\oplus C$为0，因此枚举$k\oplus C$第一个非0的位置（通过$p_{i}$来控制），最后再判断所有$p_{i}$都确定了（即$k\oplus C$在$d$最高位即以上都为0时）能否即可

#include<bits/stdc++.h>
using namespace std;
#define N 4005
#define mod 998244353
struct ji{
    int l;
    bitset<N>a;
}m,a[11];
int n,ans,mi[N];
char s[N];
void read(ji &a){
    scanf("%s",s);
    a.l=strlen(s);
    for(int i=0;i<a.l;i++)a.a[a.l-i]=s[i]-'0';
}
ji gcd(ji x,ji y){
    if (!y.l)return x;
    x.a^=(y.a<<x.l-y.l);
    while ((x.l)&&(!x.a[x.l]))x.l--;
    if (x.l<y.l)swap(x,y);
    return gcd(x,y);
}
void write(ji a){
    for(int i=a.l;i;i--){
        int p=a.a[i];
        printf("%d",p);
    }
    printf("\n");
}
int main(){
    mi[0]=1;
    for(int i=1;i<N-4;i++)mi[i]=mi[i-1]*2%mod;
    scanf("%d",&n);
    read(m);
    read(a[1]);
    for(int i=2;i<=n;i++){
        read(a[i]);
        if (a[1].l<a[i].l)swap(a[1],a[i]);
        a[1]=gcd(a[1],a[i]);
    }
    if (m.l<a[1].l){
        printf("1");
        return 0;
    }
    ji s=m;
    for(int i=m.l;i>=a[1].l;i--){
        if (m.a[i])ans=(ans+mi[i-a[1].l])%mod;
        if (s.a[i])s.a^=(a[1].a<<i-a[1].l);
    }
    ans=(ans+1)%mod;
    for(int i=a[1].l-1;i;i--)
        if (s.a[i]){
            if (!m.a[i])ans=(ans+mod-1)%mod;
            break;
        }
    printf("%d",ans);
}