LG4516/LOJ2546 「JSOI2018」潜入行动 树上背包

问题描述

LG4516

LOJ2546

---

题解

好一个毒瘤题。

hkk：JSOI的签到题

设(opt[i][j][0/1][0/1])代表结点(i)的子树，放置(j)个，(i)放不放，(i)是否覆盖的方案数。

DP方程太长，无力打出（真·原因：我要睡觉！）。

---

(mathrm{Code})

#include<bits/stdc++.h>
using namespace std;

template <typename Tp>
void read(Tp &x){
	x=0;char ch=1;int fh;
	while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
	if(ch=='-'){
		fh=-1;ch=getchar();
	}
	else fh=1;
	while(ch>='0'&&ch<='9'){
		x=(x<<1)+(x<<3)+ch-'0';
		ch=getchar();
	}
	x*=fh;
}

const int maxn=100003;
const int mod=1000000007;

int sum(long long x,long long y){
	x=x%mod,y=y%mod;
	return (int)((x+y)%mod);
}

int opt[maxn][103][2][2];
int tmp[103][2][2],size[maxn];

int n,k;

int Head[maxn],to[maxn*2],Next[maxn*2],tot;

void add(int x,int y){
	to[++tot]=y,Next[tot]=Head[x],Head[x]=tot;
}

void dp(int x,int fa){
	size[x]=opt[x][0][0][0]=opt[x][1][1][0]=1;
	for(int i=Head[x];i;i=Next[i]){
		int y=to[i];
		if(y==fa) continue;
		dp(y,x);
		for(int j=0;j<=min(size[x],k);j++){
			tmp[j][0][0]=opt[x][j][0][0],opt[x][j][0][0]=0;
			tmp[j][0][1]=opt[x][j][0][1],opt[x][j][0][1]=0;
			tmp[j][1][0]=opt[x][j][1][0],opt[x][j][1][0]=0;
			tmp[j][1][1]=opt[x][j][1][1],opt[x][j][1][1]=0;
		}
		for(int j=0;j<=min(size[x],k);j++){
			for(int p=0;p<=min(size[y],k-j);p++){
				opt[x][j+p][0][0]=sum((long long)opt[x][j+p][0][0],(long long)tmp[j][0][0]*(long long)opt[y][p][0][1]);
			    opt[x][j+p][0][1]=sum((long long)opt[x][j+p][0][1],(long long)tmp[j][0][0]*(long long)opt[y][p][1][1]+(long long)tmp[j][0][1]*((long long)opt[y][p][1][1]+(long long)opt[y][p][0][1]));
                opt[x][j+p][1][0]=sum((long long)opt[x][j+p][1][0],(long long)tmp[j][1][0]*((long long)opt[y][p][0][0]+(long long)opt[y][p][0][1]));
                opt[x][j+p][1][1]=sum((long long)opt[x][j+p][1][1],(long long)tmp[j][1][0]*((long long)opt[y][p][1][0]+(long long)opt[y][p][1][1])+(long long)tmp[j][1][1]*((long long)opt[y][p][0][0]+(long long)opt[y][p][0][1]+(long long)opt[y][p][1][0]+(long long)opt[y][p][1][1]));
			}
		}
		size[x]+=size[y];
	}
}

int main(){
	read(n);read(k);
	for(int i=1,x,y;i<n;i++){
		read(x);read(y);
		add(x,y);add(y,x);
	}
	dp(1,0);
	printf("%d
",(int)(opt[1][k][0][1]+opt[1][k][1][1])%mod);
	return 0;
}