关闭

杨氏矩阵与钩子公式

876人阅读 评论(0) 收藏 举报
分类:

杨氏矩阵又叫杨氏图表,它是这样一个矩阵,满足条件:

 

(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;
}


版权声明:本文为博主原创文章,未经博主允许不得转载。

0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:566171次
    • 积分:9669
    • 等级:
    • 排名:第1948名
    • 原创:400篇
    • 转载:76篇
    • 译文:6篇
    • 评论:111条
    博客专栏
    相关链接
    最新评论