保镖排队
(p3.pas/cpp/in/out)
【问题背景】
教主LHX作为知名人物,时刻会有恐怖分子威胁他的生命。于是教主雇佣了一些保镖来保障他的人生安全。
【题目描述】
教主一共雇佣了N个保镖,编号为1~N。每个保镖虽然身手敏捷武功高强,但是他在其余N-1个保镖里,都会有一个“上司”,他会对他的上司言听计从。但一号保镖例外,他武功盖世,不惧怕其余任何保镖,所以他没有上司。
教主LHX会对这N个保镖进行定期视察。每次视察的时候,首先会让所有保镖排队。
对于每个保镖,在他心目中会对他的所有下属的武功实力排个队。
现在教主要求排出来的队伍满足:①互为上司-下属的两个保镖,上司在前,下属在后 ②对于一个保镖的所有下属,武功实力较强的在前,较弱的在后。
教主想知道,总的排队方法数除以10007的余数是多少。
【输入格式】
输入的第一行为一个正整数T,表示了数据组数。
对于每组数据:
第一行为一个正整数N。
接下来N行,每行描述一个保镖。
第i+1行,会有一个整数K,代表第i个保镖的下属个数,接下来K个数,代表第i个保镖的下属按照武功实力从高到低的编号。
【输出格式】
输出包括C行,每行对于每组数据输出方案数mod 10007后的结果。
【样例输入】
2
5
2 2 3
2 4 5
0
0
0
7
2 2 3
2 4 5
2 6 7
0
0
0
0
【样例输出】
3
10
【样例说明】
对于第1组数据,有以下3种排列是合法的:
1 2 4 3 5
1 2 3 4 5
1 2 4 5 3
同时满足了1在2与3之前且2在3之前,2在4与5之前且4在5之前
【数据规模】
对于20%的数据,有N ≤ 9;
对于40%的数据,有对于所有K,有K ≤ 2;
对于60%的数据,有N ≤ 100;
对于100%的数据,有T ≤ 10,N ≤ 1000,K ≤ N。
这是一道比较水的树形DP。
但是又爆零了,只是因为统计树的度的时候统计错了!!!就是DFS错了!!!我想杀人呀。
因为到了空节点没有跳出,还有统计度忘了加上根的1。。。
我是先将树多叉转二叉,因为根据题意,父亲,和较大的兄弟必须排在自己之前,所以我用这种方法统一表示——只要是父亲就必须在自己之前。
多叉转二叉很容易实现。right[s] = left[i];left[i] = s;
而这道题的难点在数学上,即两个序列,要把它们合并成一个序列的方法总数。。我是一个数学弱菜,所以没有用排列组合公式推出来,于是我写的朴素找的规律。
F(i,j)表示两个序列的长度分别为i和j时,合并的方法总数,F(i,j) = F(i-1,j)+F(i,j-1)。所以f[u] = f[left[u]] * f[right[u]] * F[sum[left[u]]][sum[right[u]]]
其实这个很好推的,而且一开始我想到了,不知为何把它给否定了。(i + j)! / (i! * j!),也就是先随便排,再分别除以i和j的全排列,因为两个序列是毫不相干的。转换一下就是C(i+j,j)。
#include <cstdio>
#include <string>
#include <cstring>
long func[1100][1100];
long left[1110];
long right[1110];
long sum[1110];
long f[1110];
void dp(long u)
{
if (!left[u]&&!right[u])
{
f[u] = 1;
return;
}
else if(left[u]&&right[u])
{
dp(left[u]);
dp(right[u]);
f[u] = ((f[left[u]]*f[right[u]])%10007*func[sum[left[u]]][sum[right[u]]])%10007;
return;
}
else if (left[u])
{
dp(left[u]);
f[u] = f[left[u]];
return;
}
else if (right[u])
{
dp(right[u]);
f[u] = f[right[u]];
return;
}
}
void dfs(long u)
{
if (!u) return;
if (!left[u] && !right[u])
{
sum[u] = 1;
return;
}
dfs(left[u]);
dfs(right[u]);
sum[u] = sum[left[u]] + sum[right[u]] + 1;
}
long getint()
{
long rs=0;bool sgn=1;char tmp;
do tmp = getchar();
while (!isdigit(tmp)&&tmp-'-');
if (tmp == '-'){tmp=getchar();sgn=0;}
do rs=(rs<<3)+(rs<<1)+tmp-'0';
while (isdigit(tmp=getchar()));
return sgn?rs:-rs;
}
long tmp[1100];
int main()
{
freopen("p3.in","r",stdin);
freopen("p3.out","w",stdout);
long T = getint();
for (long i=1;i<1001;i++)
func[0][i] = 1;
for (long i=1;i<1001;i++)
{
func[i][0] = 1;
for (long j=1;j<1001;j++)
{
func[i][j] = (func[i-1][j]+func[i][j-1])%10007;
}
}
while (T--)
{
memset(left,0,sizeof left);
memset(right,0,sizeof right);
memset(sum,0,sizeof sum);
memset(f,0,sizeof f);
long n = getint();
for (long i=1;i<n+1;i++)
{
long c = getint();
for (long j=c;j>0;j--)
{
tmp[j] = getint();
}
for (long j=1;j<c+1;j++)
{
long s = tmp[j];
right[s] = left[i];
left[i] = s;
}
}
dfs(1);
dp(1);
printf("%ld\n",f[1]);
}
return 0;
}