51nod1446 Kirchhoff矩阵+Gauss消元+容斥+折半DFS

思路：

//By SiriusRen
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
const int mod=;
int cases,n,maxval,a[][],C[][],f[],val[],X,g[];
int top2,top,bgn,half,cnts2[],allnum[],Ans,T;
struct Node{int wei,num;Node(){}Node(int x,int y){wei=x,num=y;}}s[+(<<)],s2[+(<<)];
bool operator<(Node a,Node b){return a.wei<b.wei;}
int pow(ll x,int y){
    ll res=;
    while(y){
        if(y&)res=res*x%mod;
        x=x*x%mod,y>>=;
    }return (int)res;
}
int Gauss(int n){
    int f=;
    for(int i=;i<=n;i++){
        int j=i;while(!a[i][j]&&j<=n)j++;
        if(j==n+)return ;
        if(j!=i){
            for(int k=;k<=n;k++)swap(a[i][k],a[j][k]);
            f*=-;
        }
        int t=pow(a[i][i],mod-);
        for(int j=i+;j<=n;j++){
            int ww=1ll*a[j][i]*t%mod;
            for(int k=i;k<=n;k++)a[j][k]=(a[j][k]-1ll*ww*a[i][k]%mod+mod)%mod;
        }
    }
    ll ans=;
    for(int i=;i<=n;i++)ans=ans*a[i][i]%mod;
    if(f==-)ans=(mod-ans)%mod;
    return ans;
}
void Kirchhoff(int x){
    memset(a,,sizeof(a));
    for(int i=;i<=n;i++)
        if(i<=x)for(int j=max(i+,X+);j<=n;j++)a[i][j]--,a[j][i]--,a[i][i]++,a[j][j]++;
        else for(int j=i+;j<=n;j++)a[i][j]--,a[j][i]--,a[i][i]++,a[j][j]++;
}
bool cmp(int x,int y){return x>y;}
void dfs(int x,int wei,int deep){
    !T?s[top++]=Node(wei,deep):s2[top2++]=Node(wei,deep);
    for(int i=x;i>bgn;i--)if(wei+val[i]<=maxval)dfs(i-,wei+val[i],deep+);
}
signed main(){
    for(int i=;i<=;i++){
        C[i][]=C[i][i]=;
        for(int j=;j<i;j++)C[i][j]=(C[i-][j-]+C[i-][j])%mod;
    }
    scanf("%d",&cases);
    while(cases--){
        memset(allnum,,sizeof(allnum));
        memset(cnts2,,sizeof(cnts2));
        top=top2=X=bgn=Ans=T=;
        scanf("%d%d",&n,&maxval);
        for(int i=;i<=n;i++){
            scanf("%d",&val[i]);
            if(~val[i])X++;
        }half=(X+)/;
        for(int i=;i<=X;i++)Kirchhoff(i),f[i]=g[i]=Gauss(n-);
        for(int i=X;~i;i--)
            for(int j=i+;j<=X;j++)f[i]=((f[i]-1ll*C[X-i][j-i]*f[j])%mod+mod)%mod;
        sort(val+,val++n,cmp),random_shuffle(val+,val++X);
        dfs(half,,),sort(s,s+top);
        T=,bgn=half,dfs(X,,),sort(s2,s2+top2);
        for(int i=;i<top2;i++)cnts2[s2[i].num]++;
        for(int i=,j=top2-;i<top;i++){
            while(s[i].wei+s2[j].wei>maxval)cnts2[s2[j].num]--,j--;
            for(int k=;k<=X-half;k++)(allnum[k+s[i].num]+=cnts2[k])%=mod;
        }
        for(int i=;i<=X;i++)Ans=(Ans+1ll*f[X-i]*allnum[i])%mod;
        printf("%d\n",(Ans+mod)%mod);
    }
}