SPOJTime Limit Exceeded 高维前缀和优化+ dp

参考：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;
    }

}