Problem 2 树(tree.cpp/c/pas)

【题目描述】

L发明了一种与树有关的游戏（友情提醒：树是一个没有环的连通图）：他从树中删除任意数量（可以为0）的边，计算删除后所有连通块大小的乘积，L将得到这么多的分数。你的任务就是对于一颗给定的树，求出L能得到的最大分数。

【输入格式】

第一行一个整数n，表示树的节点个数。
　　接下来n-1行，每行两个整数a[i],b[i](1<=a[i],b[i]<=n)，表示a[i]与b[i]之间连边。
　　保证输入的图是一棵树。

【输出格式】

输出一个整数，表示L能得到的最大分数。

 【样例输入】

样例1：
5
1 2
2 3
3 4
4 5
样例2：
8
1 2
1 3
2 4
2 5
3 6
3 7
6 8
样例3：
3
1 2
1 3

【样例输出】

样例1：
6
样例2：
18
样例3：
3

【数据范围】

    对于10%的数据，1<=n<=5；
　　对于30%的数据，1<=n<=100；
　　另有30%的数据，保证数据是一条链。
　　对于100%的数据，1<=n<=700；

树上背包

f[i][j]表示i的父亲的连通块在子树i中有j个的最大的最大值。

于是这就是树形Dp+背包合并了、

背包合并2个先合并，再与第三个……

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cstdlib>
#include<cmath>
using namespace std;
#define MAXN (700+10)
#define ll long long
#define F (100000000)
int n,edge[MAXN*2],pre[MAXN],next[MAXN*2],size=0,son[MAXN];
struct bign
{
	ll a[40];
	bign(){memset(a,0,sizeof(a));a[0]=1;}
	bign(int x){memset(a,0,sizeof(a)); a[0]=1;a[1]=x;	}
	ll& operator[](const int i){return a[i];	}
	friend bign operator*(bign a,bign b)
	{
		bign c;
		for (int i=1;i<=a[0];i++)
			for (int j=1;j<=b[0];j++)
			{
				c[i+j-1]+=a[i]*b[j];
				if (c[i+j-1]>F)
				{
					c[i+j]+=c[i+j-1]/F;
					c[i+j-1]%=F;
				}
			}
		c[0]=min(a[0]+b[0],(long long)39);while (!c[c[0]]&&c[0]>1) c[0]--;
		return c;
	}
	friend bool operator>(bign a,bign b)
	{
		if (a[0]!=b[0]) return a[0]>b[0];
		for (int i=a[0],j=b[0];i>0;i--,j--) if (a[i]!=b[j]) return a[i]>b[j];
		return false;
	}
	void print()
	{
		printf("%I64d",a[a[0]]);
		for (int i=a[0]-1;i;i--)
		{
			printf("%.8I64d",a[i]);
		}
	}
}f[MAXN][MAXN];
bign max(bign a,bign b)
{
	if (a>b) return a;
	return b;
}
void addedge(int u,int v)
{
	edge[++size]=v;
	next[size]=pre[u];
	pre[u]=size;
}
void dfs(int x,int father)
{
	son[x]=1;//f[x][1]=1;
	f[x][1]=1;
	for (int p=pre[x];p;p=next[p])
	{
		int &v=edge[p];
		if (v!=father)
		{
			dfs(v,x);
			/*
			for (int i=son[x]+son[v];i>0;i--)
			{
				if (i<son[x]) f[x][i]=f[x][i]*son[v];
				for (int k=son[v];k>=0;k--)
					if (i-k-1>=0) f[x][i]=max(f[x][i],f[x][i-k-1]*f[v][k]);

			}
			son[x]+=son[v];
			bign maxv=son[v];
			for (int k=0;k<=son[v]-1;k++) maxv=max(maxv,f[v][k]*(k+1));
			f[x][0]=f[x][0]*maxv;
			*/
			for (int i=son[x];i;i--)
				for (int j=son[v];j>=0;j--)
					f[x][i+j]=max(f[x][i+j],f[x][i]*f[v][j]);
			son[x]+=son[v];
		}
	}
	f[x][0]=bign(son[x]);
	for (int i=1;i<=son[x];i++) f[x][0]=max(f[x][0],f[x][i]*bign(i));




	return;
}
int main()
{
//	freopen("tree.in","r",stdin);
//	freopen("tree.out","w",stdout);
	scanf("%d",&n);
	memset(pre,0,sizeof(pre));
	memset(next,0,sizeof(next));
	for (int i=1;i<n;i++)
	{
		int u,v;
		scanf("%d%d",&u,&v);
		addedge(u,v);addedge(v,u);
	}
	addedge(n+1,1);
	dfs(n+1,0);	//n+1 is ans
	/*
	for (int i=1;i<=n+1;i++)
	{
		for (int j=0;j<=son[i]-1;j++)
		{
			f[i][j].print();printf(" ");
		}
		printf("n");
	}
	*/
//	f[n+1][1]=bign(123456789)*bign(234567899);

	f[1][0].print();
	printf("n");

	return 0;
}