[51Nod1446] 限制价值树 (容斥+MT定理+折半搜索)

Description

有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后的结果。

说明：两棵树是不同的，指两棵树的边集不同，注意这里的边都是无向边。

Input

多组测试数据，第一行一个整数T，表示测试数据数量，1<=T<=5

每组测试数据有相同的结构构成：

每组数据第一行两个整数N与maxVal,满足1<=N<=40，0<=maxVal<=1,000,000,000。

第二行有N个整数，即val[0] ~ val[N-1]的数值，满足-1<=val[i]<=25,000,000。

Output

每组数据一行输出，即限制价格树的个数。

Sample Input

3

4 3

1 2 -1 3

4 5

1 2 -1 3

5 6

-1 -1 2 5 5

Sample Output

3

7

20

Solution

f[i]表示至少有i个good点不great的生成树个数（这i个good点只选一种情况）

g[i]表示恰有i个good点不great的生成树个数（这i个good点只选一种情况）

以上两个数组可以只用一个数组qwq

h[i]表示这tot个good点中使i个为great并且权值不超过maxval的方案数

最终答案就是$\sum_{i=1}^{tot}g[i]*h[tot-i]$

f[i]直接由矩阵树定理求出，g[i]显然可以由f[i]递推得到

h[i]直接折半搜索

PS:注意不要爆空间

Code

//By Menteur_Hxy

#include <cmath>

#include <cstdio>

#include <cstring>

#include <cstdlib>

#include <iostream>

#include <algorithm>

#define F(i,a,b) for(register int i=(a);i<=(b);i++)

#define R(i,a,b) for(register int i=(b);i>=(a);i--)

using namespace std;

typedef long long LL;

typedef pair<int,int> PII;

inline LL read() {

	LL x=0,f=1; char c=getchar();

	while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}

	while(isdigit(c)) x=(x<<1)+(x<<3)+c-48,c=getchar();

	return x*f;

}

const int N=41,MOD=1e9+7,MAX=1048580;

int n,m,tot,cnt,ans;

int val[N],g[N],h[N],mat[N][N],fac[N],inv[N],Cnt[MAX][21];//注意空间大小

PII da[MAX];

LL qpow(LL a,LL b) {

	LL t=1;

	while(b) {

		if(b&1) t=t*a%MOD;

		a=a*a%MOD; b>>=1;

	}

	return t;

}

void init() {

	fac[0]=1; F(i,1,40) fac[i]=(LL)fac[i-1]*i%MOD;

	inv[40]=qpow(fac[40],MOD-2);

	R(i,0,40-1) inv[i]=(LL)inv[i+1]*(i+1)%MOD;

}

LL C(int m,int n) {return (LL)fac[m]*inv[m-n]%MOD*inv[n]%MOD;}

void print() {

	F(i,1,n) {

		F(j,1,n) cout<<mat[i][j]<<" ";

		cout<<endl;

	}cout<<endl;

}

LL getf(int x) {

	LL res=1,fla=1;

	memset(mat,0,sizeof(mat));

	F(i,1,x) F(j,tot+1,n) mat[i][i]++,mat[j][j]++,mat[i][j]--,mat[j][i]--;

	F(i,x+1,n) F(j,i+1,n) mat[i][i]++,mat[j][j]++,mat[i][j]--,mat[j][i]--;//注意连边

	// print();

	// cout<<x<<endl;

	F(i,1,n-1) {//注意去掉一行一列!!

		int mx=i; F(j,i,n-1) if(mat[j][i]!=0) {mx=j;break;}

		if(mx!=i) {F(j,i,n-1) swap(mat[i][j],mat[mx][j]);fla=-fla;}

		if(!mat[mx][i]) return 0;

		res=(LL)res*mat[i][i]%MOD;

		LL Inv=qpow(mat[i][i],MOD-2);

		F(j,i,n-1) mat[i][j]=(LL)mat[i][j]*Inv%MOD;

		F(j,i+1,n-1) {

			LL tmp=mat[j][i]; mat[j][i]=0;

			F(k,i+1,n-1) mat[j][k]=(mat[j][k]-(LL)tmp*mat[i][k]%MOD+MOD)%MOD;

		}

		// print();

	}

	return res*fla;

}

void dfs1(int l,int r,int num,int sum) {

	if(l>r) {da[++cnt]=PII(sum,num);return ;}

	dfs1(l+1,r,num,sum);

	if(sum+val[l]<=m) dfs1(l+1,r,num+1,sum+val[l]);

}

void dfs2(int l,int r,int num,int sum) {

	if(l>r) {

		int p=upper_bound(da+1,da+1+cnt,PII(m-sum,n+1))-da-1;//!!!

		F(i,0,tot/2) h[i+num]=(h[i+num]+Cnt[p][i])%MOD;

		return ;

	}

	dfs2(l+1,r,num,sum);

	if(sum+val[l]<=m) dfs2(l+1,r,num+1,sum+val[l]);

}

bool cmp(int x,int y) {return x>y;}

void solve() {

	n=read(),m=read();tot=cnt=ans=0;

	F(i,1,n) val[i]=read(),tot+=(val[i]!=-1);

	sort(val+1,val+1+n,cmp);

	F(i,0,tot) g[i]=getf(i);

	// F(i,0,tot) cout<<g[i]<<" ";cout<<endl;

	R(i,0,tot) F(j,i+1,tot) g[i]=(g[i]-(LL)C(tot-i,j-i)*g[j]%MOD+MOD)%MOD;

	dfs1(1,tot/2,0,0);

	sort(da+1,da+1+cnt);

	F(i,1,cnt) {

		F(j,0,tot/2) Cnt[i][j]=Cnt[i-1][j];

		Cnt[i][da[i].second]++;

	}

	memset(h,0,sizeof(h));//初始化！！

	dfs2(tot/2+1,tot,0,0);

	F(i,0,tot) ans=(ans+(LL)g[i]*h[tot-i]%MOD)%MOD;

	// F(i,1,tot) cout<<val[i]<<" ";cout<<endl;cout<<endl;

	// F(i,1,cnt) cout<<da[i].first<<" "<<da[i].second<<endl;cout<<endl;

	// F(i,1,cnt) F(j,1,10) if(Cnt[i][j]) cout<<

	// F(i,0,tot) cout<<g[i]<<" "<<h[i]<<endl;

	printf("%lld\n",ans);

}

int main() {

	init();

	int T=read();

	while(T--) solve();

	return 0;

}