坑爹的钩子公式！POJ 2279 （《算法竞赛进阶指南》）

猛戳我 传送至POJ 2279！
参考网站：http://blog.sina.com.cn/s/blog_67bccccb0100vwfa.html

题意：有n个学生，编号1~n，站成k排，给出每排的人数，然后每排身高从左往右递增，每列身高从前往后递增。问有多少种排列方式。

思路：题目中的排列方式称为杨氏矩阵，它有一个公式来计算排列的数目，该数目设为F(n)函数的话，那么有递推式：
F(1) = 1
F(2) = 2
F(n) = F(n-1) + (n-1) * F(n-2) ，( n > 2 )

如果是规定了形状的杨氏矩阵，那么就有一个钩子公式来计算排列的数目
对于n个元素，公式为：
ans = n! * （每个格子**【下方】**和右边格子的数量+1 的乘积）
（PS：网上别的钩子公式是说上方和右方的格子数量，但实际上代码写是下方和右方的格子数量，所以这里擅自更改，有错误还请提出）

举例：如k=3 ，各行人数为 row[1] = 3 ，row[2] = 2 ，row[3] = 1，每个格子的值如图所示：

AC代码如下：

#include<iostream>
#include<stdio.h>
#include<cmath>
#include<queue>
#include<algorithm>
#include<cstring>
#include<map>
using namespace std;
int main()
{
    int n,i,j,k,cnt,tmp,sum[2005],row[15];
    long long ans,a,b;
    while(scanf("%d",&n),n)
    {
        memset(sum,0,sizeof(sum));
        cnt=0;
        for(i=1;i<=n;i++) scanf("%d",&row[i]);
        for(i=1;i<=n;i++) // i 表示第i行
        {
            for(j=1;j<=row[i];j++) // j 表示第i行的第j个
            {
                cnt++; // cnt 表示所有元素中第 cnt 个元素
                for(k=i+1;k<=n;k++) // k 表示第i行的下面的行
                {
                    if(row[k]>=j) //说明第cnt元素的下面存在元素
                        sum[cnt]++;
                    else break;
                }
                sum[cnt]+=row[i]-j+1; //加上右边剩下的元素 并且 +1
            }
        }
        a=1; //阶乘，公式的分子
        b=1; //公式的分母
        // 这里不是直接套公式，因为30!比较大会爆long long，
        // 所以一边乘，一边看能不能约分
        for(i=1;i<=cnt;i++)
        {
            a*=i;
            b*=sum[i];
            tmp=__gcd(a,b);
            if(tmp!=1) //如果乘到这里可以约分
            {
                a/=tmp;
                b/=tmp;
            }
        }
        ans=a/b;
        printf("%lld\n",ans);
    }
    return 0;
}