title
给定 \(n\) 和 \(n\) 个数 \(c[i]\) ,要求构造出一串序 列 \(a[i]\) ,要求:
\(a[i] \in [0,2^m)\) ;
\(a[i] \%c[i]!=0\) ;
\(a[i]\&a[i+1]=0\) ;
求方案数(模数 \(10^9\) )。
analysis
对 \(3\) 的处理:用 \(f[i]\) 表示以 \(i\) 结尾的方案数。
要满足 \(a[i]\&a[i+1]=0\) ,那么 \(f[i]\) 是从上一次结果中 \(i\&j=0\) 的结果转移过来的。
而 \(i\&j=0\) 又可等价于 \(i\&(\sim j)=i\) ,所以每次要对上一次的结果的下标进行取反。
所以,求 \(f[i]\) ,就相当于求以 \(i\) 为子集的上一次求出的 \(f\) 的高维前缀和。
对于 \(2\) 的处理:就把下标为 \(c[i]\) 的倍数的 \(f[i](i\%c[i]=0)\) 的取零即可。
code
#include<bits/stdc++.h>
const int maxn=55,maxm=15,mod=1e9;
namespace IO
{
char buf[1<<15],*fs,*ft;
inline char getc() { return (ft==fs&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),ft==fs))?0:*fs++; }
template<typename T>inline void read(T &x)
{
x=0;
T f=1, ch=getchar();
while (!isdigit(ch) && ch^'-') ch=getchar();
if (ch=='-') f=-1, ch=getchar();
while (isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48), ch=getchar();
x*=f;
}
char Out[1<<24],*fe=Out;
inline void flush() { fwrite(Out,1,fe-Out,stdout); fe=Out; }
template<typename T>inline void write(T x,char str)
{
if (!x) *fe++=48;
if (x<0) *fe++='-', x=-x;
T num=0, ch[20];
while (x) ch[++num]=x%10+48, x/=10;
while (num) *fe++=ch[num--];
*fe++=str;
}
}
using IO::read;
using IO::write;
int c[maxn],f[1<<maxm];
int main()
{
int T;read(T);
while (T--)
{
int n,m;read(n);read(m);
for (int i=1; i<=n; ++i) read(c[i]);
for (int i=0; i< (1<<m); ++i) f[i]=(i==0);//赋初值,相当于 f[0]=1, f[1~(1<<m)]=0
for (int i=1; i<=n; ++i)
{
for (int j=0; j< (1<<m); j+=2) std::swap(f[j],f[j^((1<<m)-1)]);//把下标互补的值互换,下标之异或和等于 (1<<m)-1
for (int j=0; j< m; ++j)
for (int k=0; k< (1<<m); ++k)
if (~k&(1<<j)) f[k]=(f[k]+f[k|(1<<j)])%mod;//高维前缀和
for (int j=0; j< (1<<m); j+=c[i]) f[j]=0;//把下标为 c[i] 的倍数的清零
}
int ans=0;
for (int i=0; i< (1<<m); ++i) ans=(ans+f[i])%mod;//方案数叠加取模
write(ans,'\n');
}
IO::flush();
return 0;
}