有N个点（N<=40）标记为0,1,2，...N-1,每个点i有个价值val[i],如果val[i]=-1那么这个点被定义为bad，否则如果val[i] >=0那么这个点为定义为good。现在给这N个点间连上N-1条边，使它们构成一个生成树，定义树中的点为great点当且仅当这个点本身是good点且与其相邻的点中至少有另一个good点。树的价值等于树中所有great点的价值和。定义限制价值树是指价值不大于maxVal的树，问对给定的val[]与maxVal，一共有多少种不同的限制价格树？由于答案太大，可取
modulo 1,000,000,007后的结果。
说明：两棵树是不同的，指两棵树的边集不同，注意这里的边都是无向边。
题解
首先我们发现我们只需要求一个数组f表示有i个点是(sweet)的方案数。
然后我们再去搜出所有组合情况来乘一下就好了，这个就是折半搜索+排序+双指针。
然后我们再去设(g[i])表示钦点i个点不能是(sweet)的方案数，这个用(matrix-tree)高斯消元就可以求出来。
于是：
[f_i=g_i-sum_{j=i+1}^{x}f_j*inom{x-i}{j-i}
]
代码
#include<bits/stdc++.h>
using namespace std;
#define N 43
typedef long long ll;
const int mod=1e9+7;
const int maxn=40;
ll a[N][N],c[N][N],sum[N],f[N],g[N],tong[N],cnt1,cnt2,mxval,x,val[N],num[N];
int n;
inline ll rd(){
ll x=0;char c=getchar();bool f=0;
while(!isdigit(c)){if(c=='-')f=1;c=getchar();}
while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
return f?-x:x;
}
inline void MOD(ll &x){x=x>=mod?x-mod:x;}
struct node{
int sum,cnt;
inline bool operator <(const node &b)const{
return sum<b.sum;
}
}a1[1<<20],a2[1<<20];
void dfs1(int now,int sum,int cnt){
if(now>x){
a1[++cnt1]=node{sum,cnt};
return;
}
dfs1(now+1,sum,cnt);
dfs1(now+1,sum+num[now],cnt+1);
}
void dfs2(int now,int sum,int cnt){
if(now>num[0]){
a2[++cnt2]=node{sum,cnt};
return;
}
dfs2(now+1,sum,cnt);
dfs2(now+1,sum+num[now],cnt+1);
}
inline ll power(ll x,ll y){
ll ans=1;
while(y){
if(y&1)ans=ans*x%mod;x=x*x%mod;y>>=1;
}
return ans;
}
inline ll ni(ll x){return power(x,mod-2);}
inline ll guass(int n){
ll ans=1;
for(int i=1;i<=n;++i)
for(int j=1;j<=n;++j)MOD(a[i][j]=a[i][j]+mod);
for(int i=1;i<=n;++i){
for(int j=i+1;j<=n;++j)if(i!=j){
ll x=a[j][i]*ni(a[i][i])%mod;
for(int k=1;k<=n;++k)a[j][k]=(a[j][k]-x*a[i][k]%mod+mod)%mod;
}
}
for(int i=1;i<=n;++i)ans=ans*a[i][i]%mod;
return ans;
}
inline void solve(){
num[0]=0;
memset(tong,0,sizeof(tong));
memset(f,0,sizeof(f));
memset(g,0,sizeof(g));
memset(sum,0,sizeof(sum));
n=rd();mxval=rd();
for(int i=1;i<=n;++i)val[i]=rd();
sort(val+1,val+n+1);
for(int i=1;i<=n;++i)if(val[i]!=-1){
num[++num[0]]=val[i];
}
for(int i=0;i<=num[0];++i){ //you i ge bu xiang
memset(a,0,sizeof(a));
for(int j=1;j<=i;++j)
for(int k=num[0]+1;k<=n;++k){
a[j][k]--;a[k][j]--;
a[j][j]++;a[k][k]++;
}
for(int j=i+1;j<=n;++j)
for(int k=j+1;k<=n;++k){
a[j][k]--;a[k][j]--;
a[j][j]++;a[k][k]++;
}
g[i]=guass(n-1);
}
x=num[0]/2;
cnt1=cnt2=0;
dfs1(1,0,0);
dfs2(x+1,0,0);
sort(a1+1,a1+cnt1+1);
sort(a2+1,a2+cnt2+1);
for(int i=num[0];i>=0;--i){
f[i]=g[i];
for(int j=i+1;j<=num[0];++j)MOD(f[i]=f[i]-f[j]*c[num[0]-i][j-i]%mod+mod);
}
int p=1;
for(int i=cnt2;i>=1;--i){
while(p<=cnt1&&a1[p].sum+a2[i].sum<=mxval){
tong[a1[p].cnt]++;
p++;
}
for(int j=0;j<=num[0]-a2[i].cnt;++j)MOD(sum[j+a2[i].cnt]+=tong[j]);
}
ll ans=0;
/* for(int i=0;i<=num[0];++i)cout<<sum[i]<<" ";puts("");
for(int i=0;i<=num[0];++i)cout<<g[i]<<" ";puts("");
for(int i=0;i<=num[0];++i)cout<<f[i]<<" ";puts("");*/
for(int i=0;i<=num[0];++i)MOD(ans+=f[num[0]-i]*sum[i]%mod);
printf("%lld
",ans);
}
int main(){
// freopen("1.in","r",stdin);
c[0][0]=1;
for(int i=1;i<=maxn;++i){
c[i][0]=1;
for(int j=1;j<=maxn;++j)MOD(c[i][j]=c[i-1][j]+c[i-1][j-1]);
}
int T=rd();
while(T--)solve();
return 0;
}