参考:http://blog.csdn.net/WeYoungg/article/details/77073140?locationNum=6&fps=1
这个地方把统计超集的代码改为了统计子集的代码。。
所以学的就是高维前缀和统计子集和超集的方法J
#include<cmath>
#include<algorithm>
#include<cstring>
#include<string>
#include<set>
#include<map>
#include<time.h>
#include<cstdio>
#include<vector>
#include<list>
#include<stack>
#include<queue>
#include<iostream>
#include<stdlib.h>
using namespace std;
#define LONG long long
const LONG INF=0x3f3f3f3f;
const LONG MOD=1e9;
const double PI=acos(-1.0);
#define clrI(x) memset(x,-1,sizeof(x))
#define clr0(x) memset(x,0,sizeof x)
#define clr1(x) memset(x,INF,sizeof x)
#define clr2(x) memset(x,-INF,sizeof x)
#define EPS 1e-10
#define lson l , mid , rt<< 1
#define rson mid + 1 ,r , (rt<<1)+1
#define root 1, n , 1
LONG dp[55][(1<<15)+1000] ;
LONG Sum[(1<<15)+1000] ;
LONG C[55] ;
int take[(1<<15) + 1000] ;
int main()
{
int T;
cin >> T;
while(T --)
{
int n ;
int m ;
cin>> n >> m ;
for(int i = 1; i<= n ; ++ i)
scanf("%lld",&C[i]) ;
clr0(take) ;
clr0(dp) ;
clr0(Sum) ;
for(int i = 0 ;i < (1<<m) ; ++ i) if(i%C[1] != 0)dp[1][i] = 1 ,Sum[i] = 1;
for(int i = 0; i < m ; ++i)
for(int j = 0 ; j < (1<<m) ; ++ j)
if(j&(1<<i))Sum[j] += Sum[j^(1<<i)] ;
for(int i = 2;i <= n ; ++ i)
{
for(int j = 0;j< (1<<m) ; ++ j)
{
if(j%C[i]!=0)
dp[i][j] = Sum[j^((1<<m)-1)]; //统计所有超集的补集
}
clr0(Sum) ;
for(int j = 0 ; j <(1<<m) ; ++ j)Sum[j] = dp[i][j] ;
for(int k =0; k< m ; ++ k)
{
for(int j = 0 ; j < (1<<m) ; ++ j)
{
if( (j&(1<<k)) ) Sum[j] += Sum[j^(1<<k)] ;
Sum[j] %= MOD ;
}
}
}
LONG ans = 0;
for(int i = 0 ; i< (1<<m) ; ++ i)
ans = (ans + dp[n][i]) %MOD ;
cout<<ans<<endl;
}
}