poj2279排队——杨氏矩阵与钩子公式(DP爆内存)

题目：http://poj.org/problem?id=2279

书上的DP做法会爆内存，尝试写了一个，过了样例。

转载：

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#define ll unsigned long long
using namespace std;
ll k,n[],f[][][][][];
void cl(ll a1,ll a2,ll a3,ll a4,ll a5)
{
    if(a1)f[a1][a2][a3][a4][a5]+=f[a1-][a2][a3][a4][a5];
    if(a2)f[a1][a2][a3][a4][a5]+=f[a1][a2-][a3][a4][a5];
    if(a3)f[a1][a2][a3][a4][a5]+=f[a1][a2][a3-][a4][a5];
    if(a4)f[a1][a2][a3][a4][a5]+=f[a1][a2][a3][a4-][a5];
    if(a5)f[a1][a2][a3][a4][a5]+=f[a1][a2][a3][a4][a5-];
}
int main()
{
    while(scanf("%lld",&k)==)
    {
        memset(f,,sizeof f);
        memset(n,,sizeof n);//!!!
        if(!k)return ;
        for(ll i=;i<=k;i++)
            scanf("%lld",&n[i]);
        f[][][][][]=;
        for(ll a1=;a1<=n[];a1++)
        {
            if(k>=)
                for(ll a2=;a2<=n[]&&a2<=a1;a2++)
                {
                    if(k>=)
                        for(ll a3=;a3<=n[]&&a3<=a2;a3++)
                        {
                            if(k>=)
                                for(ll a4=;a4<=n[]&&a4<=a3;a4++)
                                {
                                    if(k>=)
                                        for(ll a5=;a5<=n[]&&a5<=a4;a5++)
                                            cl(a1,a2,a3,a4,a5);
                                    else cl(a1,a2,a3,a4,);
                                }
                            else cl(a1,a2,a3,,);
                        }
                    else cl(a1,a2,,,);
                }
            else cl(a1,,,,);
        }
        printf("%lld\n",f[n[]][n[]][n[]][n[]][n[]]);
    }
    return ;
}

然后看到了“杨氏矩阵”与“钩子公式”。

转载：

---

杨氏矩阵定义（需满足的条件/特征）：（1）若格子（i，j），则该格子的右边和上边一定没有元素；

　　　　　　　　　　　　　　　　     （2）若格子（i，j）有元素data[i][j]，则该格子右边和上边相邻的格子

　　　　　　　　　　　　　　　　　　　　要么没有元素，要么有比data[i][j]大的元素。

显然有同已写元素组成的杨氏矩阵不唯一，1~n组成杨氏矩阵的个数可以写出：F[1]=1，F[2]=2，

　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　  F[n]=F[n-1]+(n-1)*F[n-2] (n>2)。

 

钩子长度的定义：该格子右边的格子数和它上边的格子数之和;

钩子公式：对于给定形状，不同的杨氏矩阵的个数为(n！/(每个格子的钩子长度加1的积))。

---

于是直接用了，注意要步步约分。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#define ll unsigned long long
using namespace std;
int k,n[];
ll h[][],s1,s2;
ll gcd(ll x,ll y){
return x%y==?y:gcd(y,x%y);}
int main()
{
    while(scanf("%d",&k)==)
    {
        if(!k)return ;
        memset(h,,sizeof h);
        memset(n,,sizeof n);//!
        for(int i=;i<=k;i++)
            scanf("%d",&n[i]);
        for(int i=;i<=k;i++)
            for(int j=;j<=n[i];j++)
            {
                h[i][j]=n[i]-j+;
                int k=i+;
                while(n[k]>=j)h[i][j]++,k++;
            }
        int t=;s1=;s2=;
        for(int i=;i<=k;i++)
            for(int j=;j<=n[i];j++)
            {
                t++;
                s1*=t;s2*=h[i][j];
                ll tmp=gcd(s1,s2);
                s1/=tmp;s2/=tmp;
            }
        printf("%lld\n",s1/s2);
    }
    return ;
}