杨氏矩阵与钩子公式
杨氏矩阵又叫杨氏图表，它是这样一个矩阵，满足条件：
(1)如果格子(i,j)没有元素，则它右边和上边的相邻格子也一定没有元素。
(2)如果格子(i,j)有元素a[i][j]，则它右边和上边的相邻格子要么没有元素，要么有元素且比a[i][j]大。
1 ~ n所组成杨氏矩阵的个数可以通过下面的递推式得到：
如图就是n=3时的杨氏矩阵。
下面介绍一个公式，那就是著名的钩子公式。
对于给定形状，不同的杨氏矩阵的个数为：n！除以每个格子的钩子长度加1的积。其中钩子长度定义为该格子
右边的格子数和它上边的格子数之和。
题目：http://poj.org/problem?id=1825
介绍完了钩子公式，那么我们可以来做一道基础题了。
题目：给四行，第一行放5个数字，第二行放三个数字，第三行放3个数字，第四行放1个数字，都是左对齐的排列，
现有1~12共12个数字，要求放到这四行中，从上到下，从左到右都是按小到大排列，问你共有几种排法？
() () () () ()
() () ()
() () ()
()
这个问题直接利用钩子公式解决即可。
杨氏矩阵既可以用来当堆，又可以当成平衡树。通常杨氏矩阵会涉及到两个问题：
(1)在杨氏矩阵中查找值为x的元素 (2)在杨氏矩阵中找第K大的元素
对于第一个问题，其实有两种方法，第一种方法就是二分查找法，这种方法的时间效率不是很好。第二种方法就是类
堆查找法。方法是这样的：从矩阵的右上角出发，对于元素a[i][j]，如果a[i][j]==x，则找到元素x，直接返
回； 如果a[i][j]> x，则向下移动，即继续比较a[i+1][j]与x；如果a[i][j] < x，则向左移动，即继续比
较a[i][j-1]与x。该算法的时间复杂度是O(m+n)。
bool Find(int a[][N],int n,int m,int x) {assert(a != NULL && n > 0 && m > 0);int row = 0;int col = m - 1;while(row <= n - 1 && col >= 0){if(a[row][col] == x) return true;else if(a[row][col] > x) col--;else row++;}return false; }
对于第二个问题，首先，二分枚举找到一个数x，它比杨氏矩阵中k个数大；然后，利用类堆查找法找到刚好小于x的
元素。该算法的时间复杂度为O((m+n)log(mn))，但不需要额外存储空间。
int get_order(int a[][N],int n,int m,int k) {int row = 0;int col = m - 1;int order = 0;while(row <= n - 1 && col >= 0){if(a[row][col] < k){order += col + 1;row++;}else col--;}return order; }int Find_Kth_Num(int a[][N],int n,int m,int k) {int low = a[0][0];int high = a[n-1][m-1];int order = 0;int mid = 0;do{mid = (low + high) >> 1;order = get_order(a,n,m,mid);if(order == k) break;else if(order > k) high = mid - 1;else low = mid + 1;}while(1);int row = 0;int col = m - 1;int ret = mid;while(row <= n - 1 && col >= 0){if(a[row][col] < mid){ret = max(ret,a[row][col]);row++;}else col--;}return ret; }
总结
- 上一篇： codeforces 数论分析题
- 下一篇： 杨氏矩阵的基本操作